A new versatile code for gamma-ray Monte-Carlo radiative transfer

Abstract

Ongoing MeV telescopes such as INTEGRAL/SPI and Fermi/GBM, and proposed telescopes including the recently accepted COSI and the e-ASTROGAM and AMEGO missions, provide another window in understanding transients. Their signals contain information about the stellar explosion mechanisms and their corresponding nucleosynthesis of short-lived radioactive isotopes. This raises the need of a radiative transfer code which may efficiently explore different types of astrophysical γ-ray sources and their dependence on model parameters and input physics. In view of this, we present our new Monte-Carlo Radiative Transfer code in Python. The code synthesizes the γ-ray spectra and light curves suitable for modeling supernova ejecta, including C+O novae, O+Ne novae, Type Ia and core-collapse supernovae. We test the code extensively for reproducing results consistent with analytic models. We also compare our results with similar models in the literature and discuss how our code depends on selected input physics and setting.