Composite biasing in Monte Carlo radiative transfer

Abstract

Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.

12 pages, accepted for publication in A&A

References in corpus (22)

• Interpreting Spectral Energy Distributions from Young Stellar Objects. I. A grid of 200,000 YSO model SEDs
• Three-dimensional radiative transfer modeling of AGN dusty tori as a clumpy two-phase medium
• Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation
• Efficient 3D NLTE dust radiative transfer with SKIRT
• SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture
• Sunrise: Polychromatic Dust Radiative Transfer in Arbitrary Geometries
• Radiative transfer in very optically thick circumstellar disks
• Three-Dimensional Dust Radiative Transfer
• Radiative transfer in disc galaxies -- III. The observed kinematics of dusty disc galaxies
• The DIRTY Model. I. Monte Carlo Radiative Transfer Through Dust
• Three-Dimensional Radiation Transfer in Young Stellar Objects
• Efficient Monte Carlo methods for continuum radiative transfer
• Dust Extinction and Emission in a Clumpy Galactic Disk. An Application of the Radiative Transfer Code TRADING
• The distribution of interstellar dust in CALIFA edge-on galaxies via oligochromatic radiative transfer fitting
• The DIRTY Model II: Self-Consistent Treatment of Dust Heating and Emission in a 3-D Radiative Transfer Code
• Benchmarking the Calculation of Stochastic Heating and Emissivity of Dust Grains in the Context of Radiative Transfer Simulations
• Using hierarchical octrees in Monte Carlo radiative transfer simulations
• Self-consistent 2-phase AGN torus models: SED library for observers
• The 3-Dimensional Distribution of Dust in NGC 891
• On the Modified Random Walk for Monte-Carlo Radiation Transfer
• Radiative transfer on hierarchial grids
• NGC4370: a case study for testing our ability to infer dust distribution and mass in nearby galaxies

Cited by in corpus (25)

• The Origin of Massive Stars: The Inertial-Inflow Model
• Monte Carlo Radiative Transfer
• Radiative Transfer Model of Dust Attenuation Curves in Clumpy, Galactic Environments
• 3D radiative-transfer for exoplanet atmospheres. gCMCRT: a GPU accelerated MCRT code
• The Failure of Monte Carlo Radiative Transfer at Medium to High Optical Depths
• Simplified 3D GCM modelling of the irradiated brown dwarf WD0137-349B
• TRUST I: A 3D externally illuminated slab benchmark for dust radiative transfer
• SOC program for dust continuum radiative transfer
• X-ray radiative transfer in full 3D with SKIRT
• Ly-alpha Radiative Transfer: Monte-Carlo Simulation of the Wouthuysen-Field Effect
• UV to submillimetre luminosity functions of TNG50 galaxies
• High-resolution, 3D radiative transfer modelling. II. The early-type spiral galaxy M81
• Modeling H$\alpha$ and He 10830 transmission spectrum of WASP-52b
• Self-consistent dust and non-LTE line radiative transfer with SKIRT
• Unbiased Monte Carlo continuum radiative transfer in optically thick regions
• Dynamic mineral clouds on HD 189733b II. Monte Carlo radiative transfer for 3D cloudy exoplanet atmospheres: combining scattering and emission spectra
• The cosmic spectral energy distribution in the EAGLE simulation
• Polarised emission from aligned dust grains in nearby galaxies: predictions from the Auriga simulations
• Three-dimensional continuum radiative transfer of polarized radiation in exoplanetary atmospheres
• The scattering order problem in Monte Carlo radiative transfer
• Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission
• Optical depth in polarised Monte Carlo radiative transfer
• Radiative transfer in disc galaxies $-$ V. The accuracy of the KB approximation
• Tetrahedral grids in Monte Carlo radiative transfer
• Improving Monte Carlo radiative transfer simulations: A shift of framework