We study the PAH emission from protoplanetary disks. First, we discuss the dependence of the PAH band ratios on the hardness of the absorbed photons and the temperature of the stars. We show that the photon energy together with a varying degree of th
e PAH hydrogenation accounts for most of the observed PAH band ratios without the need to change the ionization degree of the molecules. We present an accurate treatment of stochastic heated grains in a vectorized three dimensional Monte Carlo dust radiative transfer code. The program is verified against results using ray tracing techniques. Disk models are presented for T Tauri and Herbig Ae stars. Particular attention is given to the photo-dissociation of the molecules. We consider beside PAH destruction also the survival of the molecules by vertical mixing within the disk. By applying typical X-ray luminosities the model accounts for the low PAH detection probability observed in T Tauri and the high PAH detection statistics found in Herbig Ae disks. Spherical halos above the disks are considered. We show that halos reduce the observed PAH band-to-continuum ratios when observed at high inclination. Finally, mid-IR images of disks around Herbig Ae disks are presented. We show that they are easier to resolve when PAH emission dominate.
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 p
ower 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.
We study the structure of passively heated disks around T Tauri and Herbig Ae stars, and present a vectorized Monte Carlo dust radiative transfer model of protoplanetary disks. The vectorization provides a speed up factor of 100 when compared to a sc
alar version of the code. Disks are composed of either fluffy carbon and silicate grains of various sizes or dust of the diffuse ISM. The IR emission and the midplane temperature derived by the MC method differ from models where the radiative transfer is solved in slab geometry of small ring segments. In the MC treatment, dusty halos above the disks are considered. Halos lead to an enhanced IR emission and warmer midplane temperature than do pure disks. Under the assumption of hydrostatic equilibrium we find that the disk in the inner rim puffs up, followed by a shadowed region. The shadow reduces the temperature of the midplane and decreases the height of the extinction layer of the disk. It can be seen as a gap in the disk unless the surface is again exposed to direct stellar radiation. There the disk puffs up a second time, a third time and so forth. Therefore several gaps and ring-like structures are present in the disk surface and appear in emission images. They result from shadows in the disks and are present without the need to postulate the existence of any companion or planet. As compared to Herbig Ae stars, such gaps and ring-like structures are more pronounced in regions of terrestrial planets around T Tauri stars.
Optical/UV photons and even harder radiation components in galaxies are absorbed and scattered by dust and re-emitted at infrared wavelengths. For a better understanding of the obscured regions of the galaxies detailed models of the interaction of ph
otons with dust grains and the propagation of light are required. A problem which can only be solved by means of numerical solution of the radiative transfer equation. As a prologue we present high angular mid IR observations of galactic nuclei in the spirit of future ELT instrumentation. Dust models are discussed, which are suited to fit the extinction curves and relevant to compute the emission of external galaxies. Self-consistent radiative transfer models have been presented in spherical symmetry for starburst nuclei, in two dimensions for disk galaxies (spirals) and, more recently, in three dimensional configuration of the dust density distribution. For the latter, a highlighting example is the clumpy dust tori around AGN. Modern advances in the field are reviewed which are either based on a more detailed physical picture or progress in computational sciences.
