GPU-based Monte Carlo dust radiative transfer scheme applied to AGN
arXiv:1203.2941 · doi:10.1088/0004-637X/751/1/27
Abstract
A three dimensional parallel Monte Carlo (MC) dust radiative transfer code is presented. To overcome the huge computing time requirements of MC treatments, the computational power of vectorized hardware is used, utilizing either multi-core computer power or graphics processing units. The approach is a self-consistent way to solve the radiative transfer equation in arbitrary dust configurations. The code calculates the equilibrium temperatures of two populations of large grains and stochastic heated polycyclic aromatic hydrocarbons (PAH). Anisotropic scattering is treated applying the Heney-Greenstein phase function. The spectral energy distribution (SED) of the object is derived at low spatial resolution by a photon counting procedure and at high spatial resolution by a vectorized ray-tracer. The latter allows computation of high signal-to-noise images of the objects at any frequencies and arbitrary viewing angles. We test the robustness of our approach against other radiative transfer codes. The SED and dust temperatures of one and two dimensional benchmarks are reproduced at high precision. We utilize the Lucy-algorithm for the optical thin case where the Poisson noise is high, the iteration free Bjorkman & Wood method to reduce the calculation time, and the Fleck & Canfield diffusion approximation for extreme optical thick cells. The code is applied to model the appearance of active galactic nuclei (AGN) at optical and infrared wavelengths. The AGN torus is clumpy and includes fluffy composite grains of various sizes made-up of silicates and carbon. The dependence of the SED on the number of clumps in the torus and the viewing angle is studied. The appearance of the 10 micron silicate features in absorption or emission is discussed. The SED of the radio loud quasar 3C 249.1 is fit by the AGN model and a cirrus component to account for the far infrared emission.
12 pages, 6 figures, Accepted for publication in ApJ
References in corpus (13)
- AGN Dusty Tori: I. Handling of Clumpy Media
- Three-dimensional radiative transfer modeling of AGN dusty tori as a clumpy two-phase medium
- Efficient 3D NLTE dust radiative transfer with SKIRT
- Sunrise: Polychromatic Dust Radiative Transfer in Arbitrary Geometries
- Resolving the complex structure of the dust torus in the active nucleus of the Circinus galaxy
- Radiative transfer in very optically thick circumstellar disks
- Radiative Transfer Modeling of Three-Dimensional Clumpy AGN Tori and its Application to NGC 1068
- Three-dimensional radiative transfer models of clumpy tori in Seyfert galaxies
- Dust Extinction and Emission in a Clumpy Galactic Disk. An Application of the Radiative Transfer Code TRADING
- The Origin of the Silicate Emission Features in the Seyfert 2 Galaxy, NGC 2110
- Near- and mid-infrared photometry of high-redshift 3CR sources
- Shadows, gaps, and ring-like structures in protoplanetary disks
- Smart detectors for Monte Carlo radiative transfer
Cited by in corpus (30)
- Dust in the polar region as a major contributor to the IR emission of AGN
- Three-Dimensional Dust Radiative Transfer
- Dust in the diffuse interstellar medium: Extinction, emission, linear and circular polarisation
- Uncovering obscured luminous AGN with WISE
- Modeling optical and UV polarization of AGNs II. Polarization imaging and complex reprocessing
- Using hierarchical octrees in Monte Carlo radiative transfer simulations
- Hierarchical octree and k-d tree grids for 3D radiative transfer simulations
- Modeling optical and UV polarization of AGNs III. From uniform-density to clumpy regions
- Using 3D Voronoi grids in radiative transfer simulations
- 3-D Radiative Transfer Calculations of Radiation Feedback from Massive Black Holes: Outflow of Mass from the Dusty "Torus"
- Geometrically thick obscuration by radiation-driven outflow from magnetized tori of active galactic nuclei
- Herschel observations and a model for IRAS 08572+3915: a candidate for the most luminous infrared galaxy in the local (z < 0.2) Universe
- Synthetic observations of star formation and the interstellar medium
- Active galactic nucleus torus models and the puzzling infrared spectrum of IRAS F10214+4724
- SOC program for dust continuum radiative transfer
- A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics
- Near-Infrared Polarimetric Adaptive Optics Observations of NGC 1068: A torus created by a hydromagnetic outflow wind
- SKIRT: hybrid parallelization of radiative transfer simulations
- LOC program for line radiative transfer
- Extinction and dust properties in a clumpy medium
- PAHs in protoplanetary disks: emission and X-ray destruction
- Mid-infrared interferometry of Seyfert galaxies: Challenging the Standard Model
- Hypercubes of AGN Tori (HYPERCAT) -- II. Resolving the Torus with Extremely Large Telescopes
- High Spatial Resolution of the Mid-Infrared Emission of Compton-Thick Seyfert 2 Galaxy Mrk3
- Active Galactic Nuclei and their Role in Galaxy Evolution: The Infrared Perspective
- The XMM-Newton Bright Survey sample of absorbed quasars: X-ray and accretion properties
- Optical depth in polarised Monte Carlo radiative transfer
- SMART: spectral energy distributions Markov chain analysis with radiative transfer models
- Appearance of Dusty Filaments at Different Viewing Angles
- Paradigmatic examples for testing models of optical light polarization by spheroidal dust