Parallel computing of radiative transfer in relativistic jets using Monte Carlo method

We present numerical attempts of radiative transfer in a relativistic scattering flow that can produce gamma rays using a three-dimensional Monte Carlo code. We prepared an initial background flowfield obtained from hydrodynamical simulation of a relativistic jet in which Thomson scattering dominates compared to absorption, and solved the radiative transfer equation for the background evolved by a simple expansion model. Since a large number of sample particles is required for an accurate computation, we have parallelized the Monte Carlo code in order to obtain solutions in a practical computational time even for a long-term simulation coupled with a time-dependent flowfield. Using this code, higher parallel efficiency is achieved with larger number of particles. The obtained light curve from the simple model shows a signal of the transition from the opaque post-shock flow to the transparent regime as the flow expands, and the high-energy photons are generated by not only the Doppler boosting but also the inverse Compton scattering.