ﻻ يوجد ملخص باللغة العربية
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.
The theory and numerical modelling of radiation processes and radiative transfer play a key role in astrophysics: they provide the link between the physical properties of an object and the radiation it emits. In the modern era of increasingly high-qu
We present a new algorithm for radiative transfer, based on a statistical Monte-Carlo approach, that does not suffer from teleportation effects on the one hand, and yields smooth results on the other hand. Implicit-Monte-Carlo (IMC) techniques for mo
Context: The Monte Carlo method is the most widely used method to solve radiative transfer problems in astronomy, especially in a fully general 3D geometry. A crucial concept in any Monte Carlo radiative transfer code is the random generation of the
The theory of radiative transfer provides the link between the physical conditions in an astrophysical object and the observable radiation which it emits. Thus accurately modelling radiative transfer is often a necessary part of testing theoretical m
A crucial aspect of 3D Monte Carlo radiative transfer is the choice of the spatial grid used to partition the dusty medium. We critically investigate the use of octree grids in Monte Carlo dust radiative transfer, with two different octree constructi