The High Energy Radiation Transport Code

Brief Code description
Examples of run products

We have devised a time-dependent blazar emission code that has the capability to follow the time evolution of injected relativistic particles and calculates the emitted radiation.

Considers a homogeneous spherical emission region of (jet-frame) size R' upon injection in a constant speed conical jet;
Emission region interfused by tangled magnetic field of strength B' upon injection;
Emission/loss processes: Photomeson production [using the SOPHIA code: Muecke et al 2000], Bethe-Heitler pair production, photon-photon pair production, inverse Compton scattering (incl. Klein-Nishina treatment), synchrotron radiation (incl. quantum regime treatment [Brainerd & Petrosian 1987] from protons, electrons/positrons, charged pions & muons, decay of unstable particles, adiabatic losses/charged particle and photon escape; 'on-off switch' for all processes.
Several injection scenarios (instantaneous, continuous) for relativistic particle population of various composition;
Photon/particle transport by matrix multiplication method [Protheroe & Stanev 1993]

The modular design of this code ensures high flexibility.

Examples of final state particle & photon distributions

Overall parameter setting:

As an example we simulate here ...

electrons / positrons

photons

protons

neutrons

neutrinos

For details please contact: A. Reimer anita.reimer@uibk.ac.at.

Resulting light curves depending on e/p ratio

Overall parameter setting:

As an example we simulate here the decay of light curves as a result of instantaneous injection of power law particle spectra (with index -2) into a emission region of size R' using no external target photon fields. We assume equipartition between particles and magnetic field (one cold proton per electron in the 'leptonic' case). Further fixed parameters: Doppler factor D=20, adiabatic loss time scale =10R'/c, electron injection between 0.05.... 40 GeV, proton injection between 2....10^10 GeV.

Leptonic Model