Ray-tracing for complex astrophysical high-opacity structures

Abstract

We present a ray-tracing technique for radiative transfer modeling of complex three-dimensional (3D) structures which include dense regions of high optical depth like in dense molecular clouds, circumstellar disks, envelopes of evolved stars, and dust tori around active galactic nuclei. The corresponding continuum radiative transfer problem is described and the numerical requirements for inverse 3D density and temperature modeling are defined. We introduce a relative intensity and transform the radiative transfer equation along the rays to solve machine precision problems and to relax strong gradients in the source term. For the optically thick regions where common ray-tracers are forced to perform small trace steps, we give two criteria for making use of a simple approximative solver crossing the optically thick region quickly. Using an example of a density structure with optical depth changes of 6 orders of magnitude and sharp temperature variations, we demonstrate the accuracy of the proposed scheme using a common 5th-order Runge-Kutta ray-tracer with adaptive step size control. In our test case, the gain in computational speed is about a factor of 870. The method is applied to calculate the temperature distribution within a massive molecular cloud core for different boundary conditions for the radiation field.

21 pages, 5 figures to appear in Astrophysical Journal

References in corpus (6)

• C^2-Ray: A new method for photon-conserving transport of ionizing radiation
• The 2D Continuum Radiative Transfer Problem: Benchmark Results for Disk Configurations
• 3D Continuum radiative transfer in complex dust configurations around young stellar objects and active nuclei II. 3D Structure of the dense molecular cloud core Rho Oph D
• Discovery of a massive SCUBA core with both inflow and outflow motions
• Obscuration of Active Galactic Nuclei by Circumnuclear Starbursts
• Three-dimensional Continuum Radiative Transfer Images of a Molecular Cloud Core Evolution

Cited by in corpus (12)

• Comparing the statistics of interstellar turbulence in simulations and observations: Solenoidal versus compressive turbulence forcing
• HYPERION: An open-source parallelized three-dimensional dust continuum radiative transfer code
• Modelling the spectral energy distribution of galaxies. V. The dust and PAH emission SEDs of disk galaxies
• Three-Dimensional Dust Radiative Transfer
• Radiative transfer with POLARIS: I. Analysis of magnetic fields through synthetic dust continuum polarization measurements
• Modeling the Dust Properties of z ~ 6 Quasars with ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree
• ART^2 : Coupling Lyman-alpha Line and Multi-wavelength Continuum Radiative Transfer
• Hierarchical octree and k-d tree grids for 3D radiative transfer simulations
• Using 3D Voronoi grids in radiative transfer simulations
• The structured environments of embedded star-forming cores. PACS and SPIRE mapping of the enigmatic outflow source UYSO 1
• Effects of scattering and dust grain size on the temperature structure of protoplanetary discs: A three-layer approach
• Fast ray-tracing algorithm for circumstellar structures (FRACS) I. Algorithm description and parameter-space study for mid-IR interferometry of B[e] stars