Unbiased Monte Carlo continuum radiative transfer in optically thick regions

Abstract

Radiative transfer describes the propagation of electromagnetic radiation through an interacting medium. This process is often simulated by the use of the Monte Carlo method, which involves the probabilistic determination and tracking of simulated photon packages. In the regime of high optical depths, this approach encounters difficulties since a proper representation of the various physical processes can only be achieved by considering high numbers of simulated photon packages. As a consequence, the demand for computation time rises accordingly and thus practically puts a limit on the optical depth of models that can be simulated. Here we present a method that aims to solve the problem of high optical depths in dusty media, which relies solely on the use of unbiased Monte Carlo radiative transfer. For that end, we identified and precalculated repeatedly occuring and simulated processes, stored their outcome in a multidimensional cumulative distribution function, and immediately replaced the basic Monte Carlo transfer during a simulation by that outcome. During the precalculation, we generated emission spectra as well as deposited energy distributions of photon packages traveling from the center of a sphere to its rim. We carried out a performance test of the method to confirm its validity and gain a boost in computation speed by up to three orders of magnitude. We then applied the method to a simple model of a viscously heated circumstellar disk, and we discuss the necessity of finding a solution for the optical depth problem with regard to a proper temperature calculation. We find that the impact of an incorrect treatment of photon packages in highly optically thick regions extents even to optically thin regions, thus, changing the overall observational appearance of the disk.

arXiv comment17 pages, 11 figures, 1 table

References in corpus (13)

• Dust Grain Size Distributions and Extinction in the Milky Way, LMC, and SMC
• Protoplanetary Disk Structures in Ophiuchus
• Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer
• PIONIER: a 4-telescope visitor instrument at VLTI
• Radiative transfer in very optically thick circumstellar disks
• Radiative transfer with POLARIS: I. Analysis of magnetic fields through synthetic dust continuum polarization measurements
• Efficient radiative transfer in dust grain mixtures
• The Failure of Monte Carlo Radiative Transfer at Medium to High Optical Depths
• Composite biasing in Monte Carlo radiative transfer
• TRUST I: A 3D externally illuminated slab benchmark for dust radiative transfer
• Tracing planet-induced structures in circumstellar disks using molecular lines
• On the Modified Random Walk for Monte-Carlo Radiation Transfer
• Evolution of the T Tauri star population in the Lupus association

Cited by in corpus (16)

• A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
• Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
• Effects of Spatial Discretization in Lyman-alpha Line Radiation Transfer Simulations
• No self-shadowing instability in 2D radiation-hydrodynamical models of irradiated protoplanetary disks
• Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating
• Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions
• Feasibility of detecting and characterizing embedded low-mass giant planets in gaps in the VIS/NIR wavelength range
• Three-temperature radiation hydrodynamics with PLUTO: Thermal and kinematic signatures of accreting protoplanets
• The scattering order problem in Monte Carlo radiative transfer
• Can gap-edge illumination excite spirals in protoplanetary disks? Three-temperature radiation hydrodynamics and NIR image modelling
• Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks
• Monte Carlo post-processing for radiation hydro simulations of accreting planets in protoplanetary disks
• Improving Monte Carlo radiative transfer simulations: A shift of framework
• Starlight-driven flared-staircase geometry in radiation hydrodynamic models of protoplanetary disks
• Time-dependent Monte Carlo continuum radiative transfer
• Feasibility of interferometric observations and characterization of planet-induced structures at sub au to au scales in protoplanetary disks