A Monte-Carlo based relativistic radiation hydrodynamics code with a higher-order scheme

Understand

We develop a new relativistic radiation hydrodynamics code based on the Monte-Carlo algorithm.

• In this code, we implement a new scheme to achieve the second-order accuracy in time in the limit of a large packet number for solving the interaction between matter and radiation.
• This higher-order time integration scheme is implemented in the manner to guarantee the energy-momentum conservation to the precision of the geodesic integrator.
• The spatial dependence of radiative processes, such as the packet propagation, emission, absorption, and scattering, are also taken into account up to the second-order accuracy.