We present a Monte Carlo (MC) radiative transfer code for complex three dimensional dust distributions and include transiently heated PAH. The correctness of the code is confirmed by comparison with benchmark results. The method makes use of the parallelization capabilities of modern vectorized computing units like graphic cards. The computational speed grows linearly with the number of graphical processing units (GPU). On a conventional desktop PC, our code is up to a factor 100 faster when compared to other MC algorithms. As an example, we compute the dust emission of proto-planetary disks. We simulate how a mid-IR instrument mounted at a future 42m ELT will detect such disks. Two cases are distinguished: a homogeneous disk and a disk with an outward migrating planet, producing a gap and a spiral density wave. We find that the resulting mid-IR spectra of both disks are almost identical. However, they can be distinguished at those wavelengths by coronographic, dual-band imaging. Finally, the emission of PAHs exposed to different radiation fields is computed. We demonstrate that PAH emission depends not only on the strength but also strongly on the hardness of the radiation, a fact which has often been neglected in previous models. We find that hard photons (>20eV) easily dissociate all PAHs in the disks of T Tauri stars. To explain the low, but not negligible detection rate (<10%) of PAHs in T Tau disks, we suggest that turbulent motions act as a possible path for PAH survival.

In Spitzer observations of Tauri stars and their disks, PAH features are detected in less than 10% of the objects, although the stellar photosphere is sufficiently hot to excite PAHs. To explain the deficiency, we discuss PAH destruction by photons assuming that the star has beside its photospheric emission also a FUV, an EUV and an X-ray component with fractional luminosity of 1%, 0.1% and 0.025%, respectively. As PAH destruction process we consider unimolecular dissociation and present a simplified scheme to estimate the location from the star where the molecules become photo-stable. We find that soft photons with energies below ~20eV dissociate PAHs only up to short distances from the star (r < 1AU); whereas dissociation by hard photons (EUV and X-ray) is so efficient that it would destroy all PAHs (from regions in the disk where they could be excited). As a possible path for PAH survival we suggest turbulent motions in the disk. They can replenish PAHs or remove them from the reach of hard photons. For standard disk models, where the surface density changes like 1/r and the mid plane temperature like 1/r^{0.5}, the critical vertical velocity for PAH survival is proportional to r^{-3/4} and equals ~5m/s at 10AU which is in the range of expected velocities in the surface layer. The uncertainty in the parameters is large enough to explain both detection and non-detection of PAHs. Our approximate treatment also takes into account the presence of gas which, at the top of the disk, is ionized and at lower levels neutral.

Dust enshrouded activity can ideally be studied by mid-infrared (MIR) observations. In order to explore the AGN versus star forming origin of the nuclear MIR emission of galaxies, observations of high spatial resolution are required. Here we report on 11.3mic. observations with VISIR at the VLT, reaching 0.35'' spatial resolution (FWHM). During the scientific verification of VISIR we have observed a sample of 36 nearby galaxies having a variety of optically classified nuclear activity: 17 black hole driven active galactic nuclei (AGN), 10 starbursts (SBs) and 9 quiet spirals. 16/17 AGN are detected and unresolved, 5/10 SBs are detected and resolved with structured emission up to a few arcsec, while for 5/10 SB and all 9 quiet nuclei low upper limits are provided. The morphology of the resolved SB nuclei follows that seen at radio frequencies. The compactness of AGN and the extent of the SB nuclei is consistent with predictions from radiative transfer models and with MIR spectra of lower spatial resolution. We explore the nuclear MIR surface brightness as a quantitative measure. While AGN and SB cannot be distinguished with MIR data from 4m class telescopes, our data provide evidence that, up to a distance of 100 Mpc, AGN and SB can well be separated by means of MIR surface brightness when using 8m class telescopes.

We determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33. We separate the thermal and nonthermal components of the radio continuum emission without the assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 $\mu$m and a standard dust model, we de-redden the H$\alpha$ emission. The extinction-corrected H$\alpha$ emission serves as a template for the thermal free-free radio emission. Subtracting from the observed 3.6 and 20cm emission (Effelsberg and the VLA) this free-free emission, we obtain maps of the nonthermal intensity and spectral index.

We determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33 at a linear resolution of 360 pc. We separate the thermal and nonthermal components of the radio continuum emission without the assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 $\mu$m and a standard dust model, we deredden the H$\alpha$ emission. The extinction corrected H$\alpha$ emission serves as a template for the thermal free-free radio emission. Subtracting from the observed 3.6 cm and 20 cm emission (Effelsberg and the VLA) this free-free emission, we obtain the nonthermal maps. A constant electron temperature used to obtain the thermal radio intensity seems appropriate for M~33 which, unlike the Milky Way, has a shallow metallicity gradient. For the first time, we derive the distribution of the nonthermal spectral index across a galaxy, M33. We detect strong nonthermal emission from the spiral arms and star-forming regions. Wavelet analysis shows that at 3.6 cm the nonthermal emission is dominated by contributions from star-forming regions, while it is smoothly distributed at 20 cm. For the whole galaxy, we obtain thermal fractions of 51% and 18% at 3.6 cm and 20 cm, respectively. The thermal emission is slightly stronger in the southern than in the northern half of the galaxy. We find a clear radial gradient of mean extinction in the galactic plane. The nonthermal spectral index map indicates that the relativistic electrons suffer energy-loss when diffusing from their origin in star-forming regions towards interarm regions and the outer parts of the galaxy. We also conclude that the radio emission is mostly nonthermal at R $>$ 5 kpc in M33.

We investigate the presence and properties of PAHs on the surface of circumstellar disks around Herbig Ae/Be stars by comparing the predictions of disk models with observations. We present results of a radiation transfer code for disks heated by the central star, in hydrostatic equilibrium in the vertical direction (flared disks). The dust is a mixture of large grains in thermal equilibrium, transiently heated small grains and PAHs. Special attention is given to the influence of the stellar, disk and PAHs properties on the strength of the PAH emission lines and their spatial distribution. The models predict an infrared SED showing PAH features at 3.3, 6.2, 7.7, and 11.3 \um clearly visible above the continuum, and with some of them very strong. The PAH emission, spatially extended, comes mostly from the outer disk region ($R\sim$100 AU) while the continuum emission at similar wavelengths, mostly due to warm large grains, is confined to the innermost disk regions ($R\sim$few AU). We compare the model results to infrared observations from ISO and ground-based telescopes of some thirty Herbig Ae/Be stars. Most of the observed PAH features in objects with spectral type later than B9 are well described by our disk models and we discuss in some detail the PAH characteristics one can derive from the existing data. Objects with strong radiation field (generally earlier than about B9) have the 3.3 $\mu$m feature (often the only one observed) much weaker than predicted, and we discuss possible explanations (dissipation of the disk, photoevaporation or modification of the PAH properties).