A New Fast Monte Carlo Code for Solving Radiative Transfer Equations based on Neumann Solution
arXiv:2104.07007 · doi:10.3847/1538-4365/abec73
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 -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.
arXiv comment37 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 (47)
- Cosmological parameters from CMB and other data: a Monte-Carlo approach
- The Two-Phase Pair Corona Model for Active Galactic Nuclei and X-ray Binaries: How to Obtain Exact Solutions
- Broad emission lines for negatively spinning black holes
- An extended scheme for fitting X-ray data with accretion disk spectra in the strong gravity regime
- Cosmic Microwave Background Polarization
- The Submillimeter Bump in Sgr A* from Relativistic MHD Simulations
- A Fast New Public Code for Computing Photon Orbits in a Kerr Spacetime
- A public code for general relativistic, polarised radiative transfer around spinning black holes
- X-ray Polarization from Accreting Black Holes: Coronal Emission
- GYOTO: a new general relativistic ray-tracing code
- Modeling optical and UV polarization of AGNs I. Imprints of individual scattering regions
- grmonty: a Monte Carlo Code for Relativistic Radiative Transport
- A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. I. Quasi-Kerr Black Holes
- General relativistic radiative transfer: formulation and emission from structured tori around black holes
- Thermal disc emission from a rotating black hole: X-ray polarization signatures
- Inferring the Inclination of a Black Hole Accretion Disk from Observations of its Polarized Continuum Radiation
- Light Curves from an MHD Simulation of a Black Hole Accretion Disk
- Geometric modeling of M87* as a Kerr black hole or a non-Kerr compact object
- Tests of General Relativity in the Strong Gravity Regime Based on X-Ray Spectropolarimetric Observations of Black Holes in X-Ray Binaries
- Monte Carlo Radiative Transfer
- RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes
- Covariant Magnetoionic Theory I: Ray Propagation
- Light bending scenario for accreting black holes in X-ray polarimetry
- Covariant Magnetoionic Theory II: Radiative Transfer
- Constraining the size of the corona with fully relativistic calculations of spectra of extended corona. I - the Monte Carlo radiative transfer code
- A parameter survey of Sgr A* radiative models from GRMHD simulations with self-consistent electron heating
- A Monte Carlo Code for Relativistic Radiation Transport Around Kerr Black Holes
- YNOGK: A new public code for calculating null geodesics in the Kerr spacetime
- General relativistic polarized radiative transfer: building a dynamics-observations interface
- The impact of Faraday effects on polarized black hole images of Sagittarius A*
- Faraday Conversion and Rotation In Uniformly Magnetized Relativistic Plasmas
- General purpose ray-tracing and polarized radiative transfer in General Relativity
- Images of the radiatively inefficient accretion flow surrounding a Kerr black hole: application in Sgr A*
- RAPTOR II: Polarized radiative transfer in curved spacetime
- General relativistic radiative transfer in hot astrophysical plasmas: a characteristic approach
- Formal Solutions for Polarized Radiative Transfer. I. The DELO Family
- Monte Carlo Simulation of Comptonization in Inhomogeneous Media
- Polarization Imaging of M87 Jets by General Relativistic Radiative Transfer Calculation based on GRMHD Simulations
- MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics
- General relativistic polarized radiative transfer with inverse Compton scatterings
- Algorithms And Programs For Strong Gravitational Lensing In Kerr Space-time Including Polarization
- Monte-Carlo simulations of thermal/nonthermal radiation from a neutron-star magnetospheric accretion shell
- Formal Solutions for Polarized Radiative Transfer. II. High-order Methods
- Covariant polarized radiative transfer on cosmological scales for investigating large-scale magnetic field structures
- Simulations of recoiling black holes: adaptive mesh refinement and radiative transfer
- Discontinuities in numerical radiative transfer
- Covariant Compton Scattering Kernel in General Relativistic Radiative Transfer