Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission

Abstract

Context: The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated in the Monte Carlo framework. Net stimulated emission, or absorption with a negative cross section, does not fit this method, however. Aims: We explore alterations to the standard photon packet life cycle in Monte Carlo radiative transfer that allow the treatment of net stimulated emission without loss of generality or efficiency. Methods: We present the explicit absorption technique that allows net stimulated emission to be handled efficiently. It uses the scattering rather than the extinction optical depth along a photon packet's path to randomly select the next interaction location, and offers a separate, deterministic treatment of absorption. We implemented the technique in a special-purpose Monte Carlo code for a two-stream 1D radiative transfer problem and in the fully featured 3D code SKIRT, and we studied its overall performance using quantitative statistical tests. Results: Our special-purpose code is capable of recovering the analytical solutions to the two-stream problem in all regimes, including the one of strong net stimulated emission. The implementation in SKIRT is straightforward, as the explicit absorption technique easily combines with the variance reduction and acceleration techniques already incorporated. In general, explicit absorption tends to improve the efficiency of the Monte Carlo routine in the regime of net absorption. Conclusions: Explicit absorption allows the treatment of net stimulated emission in Monte Carlo radiative transfer, it interfaces smoothly with... (abridged)

arXiv comment14 pages, accepted for publication in A&A

References in corpus (42)

• A computer program for fast non-LTE analysis of interstellar line spectra
• Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer
• Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation
• Efficient 3D NLTE dust radiative transfer with SKIRT
• Lyman Alpha Radiation From Collapsing Protogalaxies I: Characteristics of the Emergent Spectrum
• SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture
• LIME - a flexible, non-LTE line excitation and radiation transfer method for millimeter and far-infrared wavelengths
• 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
• The 2D Continuum Radiative Transfer Problem: Benchmark Results for Disk Configurations
• Ly-alpha Radiative Transfer in Cosmological Simulations and Application to a z~8 Emitter
• Far-infrared and dust properties of present-day galaxies in the EAGLE simulations
• Monte Carlo Radiative Transfer
• The dust-continuum size of TNG50 galaxies at $z=1-5$: a comparison with the distribution of stellar light, stars, dust and H${}_2$
• 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
• Hierarchical octree and k-d tree grids for 3D radiative transfer simulations
• Using 3D Voronoi grids in radiative transfer simulations
• Galaxy Evolution in all Five CANDELS Fields and IllustrisTNG: Morphological, Structural, and the Major Merger Evolution to $z\sim 3$
• Dust evolution in zoom-in cosmological simulations of galaxy formation
• Radiative equilibrium in Monte Carlo radiative transfer using frequency distribution adjustment
• High redshift JWST predictions from IllustrisTNG: III. Infrared luminosity functions, obscured star formation and dust temperature of high-redshift galaxies
• The Failure of Monte Carlo Radiative Transfer at Medium to High Optical Depths
• Infrared emission of $z\sim 6$ galaxies: AGN imprints
• Composite biasing in Monte Carlo radiative transfer
• TRUST I: A 3D externally illuminated slab benchmark for dust radiative transfer
• Large and small-scale structures and the dust energy balance problem in spiral galaxies
• Polarization in Monte Carlo radiative transfer and dust scattering polarization signatures of spiral galaxies
• High-resolution synthetic UV-submm images for simulated Milky Way-type galaxies from the Auriga project
• SKIRT: hybrid parallelization of radiative transfer simulations
• Geometry effects on dust attenuation curves with different grain sources at high redshift
• Effects of Spatial Discretization in Lyman-alpha Line Radiation Transfer Simulations
• The structure of the Milky Way based on unWISE 3.4$\mu$m integrated photometry
• Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes
• Relic galaxy analogues in Illustris TNG50: the formation pathways of surviving red nuggets in a cosmological simulation
• The dust attenuation law in $z\sim 6$ quasars
• Polarised emission from aligned dust grains in nearby galaxies: predictions from the Auriga simulations
• The complex infrared dust continuum emission of NGC1068: ground-based N- and Q-band spectroscopy and new radiative transfer models
• CosTuuM: polarized thermal dust emission by magnetically oriented spheroidal grains
• Optical depth in polarised Monte Carlo radiative transfer
• Evidence for the First Extragalactic Hydrogen Recombination Line Maser in NGC 253

Cited by in corpus (4)

• X-ray radiative transfer in full 3D with SKIRT
• Self-consistent dust and non-LTE line radiative transfer with SKIRT
• Tetrahedral grids in Monte Carlo radiative transfer
• Impact of a binary companion in AGB outflows on CO spectral lines