Radiative equilibrium in Monte Carlo radiative transfer using frequency distribution adjustment

Abstract

The Monte Carlo method is a powerful tool for performing radiative equilibrium calculations, even in complex geometries. The main drawback of the standard Monte Carlo radiative equilibrium methods is that they require iteration, which makes them numerically very demanding. Bjorkman & Wood recently proposed a frequency distribution adjustment scheme, which allows radiative equilibrium Monte Carlo calculations to be performed without iteration, by choosing the frequency of each re-emitted photon such that it corrects for the incorrect spectrum of the previously re-emitted photons. Although the method appears to yield correct results, we argue that its theoretical basis is not completely transparent, and that it is not completely clear whether this technique is an exact rigorous method, or whether it is just a good and convenient approximation. We critically study the general problem of how an already sampled distribution can be adjusted to a new distribution by adding data points sampled from an adjustment distribution. We show that this adjustment is not always possible, and that it depends on the shape of the original and desired distributions, as well as on the relative number of data points that can be added. Applying this theorem to radiative equilibrium Monte Carlo calculations, we provide a firm theoretical basis for the frequency distribution adjustment method of Bjorkman & Wood, and we demonstrate that this method provides the correct frequency distribution through the additional requirement of radiative equilibrium. We discuss the advantages and limitations of this approach, and show that it can easily be combined with the presence of additional heating sources and the concept of photon weighting. However, the method may fail if small dust grains are included... (abridged)

17 pages, 2 figures, accepted for publication in New Astronomy

References in corpus (12)

• Self-similarity and scaling behavior of IR emission from radiatively heated dust: I. Theory
• Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer
• Mocassin: A fully three-dimensional Monte Carlo photoionization code
• The Spectral Energy Distribution of HH30 IRS: Constraining The Circumstellar Dust Size Distribution
• Radiative transfer in disc galaxies -- III. The observed kinematics of dusty disc galaxies
• Monte Carlo simulations of dusty spiral galaxies: extinction and polarization properties
• The DIRTY Model II: Self-Consistent Treatment of Dust Heating and Emission in a 3-D Radiative Transfer Code
• Modelling the Interstellar Medium of Low Surface Brightness Galaxies: Constraining Internal Extinction, Disk Color Gradients, and Intrinsic Rotation Curve Shapes
• Monte Carlo Radiative Transfer in Embedded Prestellar Cores
• Effects of Clumping on the Observed Properties of Dusty Galaxies
• Three-dimensional photoionization modelling of the planetary nebula NGC 3918
• Kinematics of elliptical galaxies with a diffuse dust component - III. A Monte Carlo approach to include the effects of scattering

Cited by in corpus (28)

• Monte Carlo radiative transfer in protoplanetary disks
• Three-dimensional radiative transfer modeling of AGN dusty tori as a clumpy two-phase medium
• Fitting the integrated Spectral Energy Distributions of Galaxies
• Efficient 3D NLTE dust radiative transfer with SKIRT
• Sunrise: Polychromatic Dust Radiative Transfer in Arbitrary Geometries
• The warm gas atmosphere of the HD 100546 disk seen by Herschel (Evidence of a gas-rich, carbon-poor atmosphere?)
• Modelling the spectral energy distribution of galaxies. V. The dust and PAH emission SEDs of disk galaxies
• Three-Dimensional Dust Radiative Transfer
• Dust Extinction and Emission in a Clumpy Galactic Disk. An Application of the Radiative Transfer Code TRADING
• SCUBA and Spitzer observations of the Taurus molecular cloud - pulling the bull's tail
• Modeling the Dust Properties of z ~ 6 Quasars with ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree
• FitSKIRT: genetic algorithms to automatically fit dusty galaxies with a Monte Carlo radiative transfer code
• The dust temperatures of the prestellar cores in the rho Oph main cloud and in other star forming regions: consequences for the core mass function
• Radiative transfer in disc galaxies IV - The effects of dust attenuation on bulge and disc structural parameters
• Panchromatic radiative transfer modeling of stars and dust in the Sombrero galaxy
• TRUST I: A 3D externally illuminated slab benchmark for dust radiative transfer
• GPU-based Monte Carlo dust radiative transfer scheme applied to AGN
• Panchromatic Spectral Energy Distributions of Dusty Galaxies with RADISHE. I. Predictions for Herschel: Correlating Colors with Galactic Energy Sources
• Monte Carlo radiative transfer in SPH density fields
• SKIRT: hybrid parallelization of radiative transfer simulations
• The initial conditions of isolated star formation -- IX. Akari mapping of an externally heated pre-stellar core
• Gravitational microlensing of AGN dusty tori
• Smart detectors for Monte Carlo radiative transfer
• Modelling Herschel observations of infrared-dark clouds in the Hi-GAL survey
• Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission
• Appearance of Dusty Filaments at Different Viewing Angles
• Diversity of dwarf galaxy IR-submm emission patterns: CLUES from hydrodynamical simulations
• Modelling the SEDs of spiral galaxies