A New Fast Monte Carlo Code for Solving Radiative Transfer Equations based on Neumann Solution

Abstract

In this paper, we proposed a new Monte Carlo radiative transport (MCRT) scheme, which is based completely on the Neumann series solution of Fredholm integral equation. This scheme indicates that the essence of MCRT is the calculation of infinite terms of multiple integrals in Neumann solution simultaneously. Under this perspective we redescribed MCRT procedure systematically, in which the main work amounts to choose an associated probability distribution function (PDF) for a set of random variables and the corresponding unbiased estimation functions. We can select a relatively optimal estimation procedure that has a lower variance from an infinite possible choices, such as the term by term estimation. In this scheme, MCRT can be regarded as a pure problem of integral evaluation, rather than as the tracing of random walking photons. Keeping this in mind, one can avert some subtle intuitive mistakes. In addition the $\delta$-functions in these integrals can be eliminated in advance by integrating them out directly. This fact together with the optimal chosen random variables can remarkably improve the Monte Carlo (MC) computational efficiency and accuracy, especially in systems with axial or spherical symmetry. An MCRT code, Lemon (Linear Integral Equations' Monte Carlo Solver Based on the Neumann solution), has been developed completely based on this scheme. Finally, we intend to verify the validation of Lemon, a suite of test problems mainly restricted to flat spacetime have been reproduced and the corresponding results are illustrated in detail.

37 pages, 28 figures. The code can be download from: https://github.com/yangxiaolinyn/Lemon (or https://bitbucket.org/yangxiaolinsc/lemonsourcecode/src/main/) and https://doi.org/10.5281/zenodo.4686355. Comments are welcome

References in corpus (11)

• The Submillimeter Bump in Sgr A* from Relativistic MHD Simulations
• Thermal disc emission from a rotating black hole: X-ray polarization signatures
• Light Curves from an MHD Simulation of a Black Hole Accretion Disk
• Monte Carlo Radiative Transfer
• Constraining the size of the corona with fully relativistic calculations of spectra of extended corona. I - the Monte Carlo radiative transfer code
• Formal Solutions for Polarized Radiative Transfer. I. The DELO Family
• General relativistic polarized radiative transfer with inverse Compton scatterings
• Algorithms And Programs For Strong Gravitational Lensing In Kerr Space-time Including Polarization
• Formal Solutions for Polarized Radiative Transfer. II. High-order Methods
• Covariant polarized radiative transfer on cosmological scales for investigating large-scale magnetic field structures
• Discontinuities in numerical radiative transfer

Cited by in corpus (3)

• Modeling H$\alpha$ and He 10830 transmission spectrum of WASP-52b
• Polarization-sensitive Compton scattering by accelerated electrons
• X-ray Polarization of Inverse Compton Scattering by Thermal and Nonthermal Electrons