We assume that dust near active galactic nuclei (AGN) is distributed in a torus-like geometry, which may be described by a clumpy medium or a homogeneous disk or as a combination of the two (i.e. a 2-phase medium). The dust particles considered are fluffy and have higher submillimeter emissivities than grains in the diffuse ISM. The dust-photon interaction is treated in a fully self-consistent three dimensional radiative transfer code. We provide an AGN library of spectral energy distributions (SEDs). Its purpose is to quickly obtain estimates of the basic parameters of the AGN, such as the intrinsic luminosity of the central source, the viewing angle, the inner radius, the volume filling factor and optical depth of the clouds, and the optical depth of the disk midplane, and to predict the flux at yet unobserved wavelengths. The procedure is simple and consists of finding an element in the library that matches the observations. We discuss the general properties of the models and in particular the 10mic. silicate band. The AGN library accounts well for the observed scatter of the feature strengths and wavelengths of the peak emission. AGN extinction curves are discussed and we find that there is no direct one-to-one link between the observed extinction and the wavelength dependence of the dust cross sections. We show that objects of the library cover the observed range of mid IR colors of known AGN. The validity of the approach is demonstrated by matching the SEDs of a number of representative objects: Four Seyferts and two quasars for which we present new Herschel photometry, two radio galaxies, and one hyperluminous infrared galaxy. Strikingly, for the five luminous objects we find pure AGN models fit the SED without a need to postulate starburst activity.

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 the 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 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.

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 scalar 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.

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.

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 photons 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.

Using the Spitzer Space Telescope, we have obtained rest frame 9-16mu spectra of 11 quasars and 9 radio galaxies from the 3CRR catalog at redshifts 1.0<z<1.4. This complete flux-limited 178MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). The mean radio galaxy spectrum shows a clear silicate absorption feature (tau_9.7mu = 1.1) whereas the mean quasar spectrum shows silicates in emission. The mean radio galaxy spectrum matches a dust-absorbed mean quasar spectrum in both shape and overall flux level. The data for individual objects conform to these results. The trend of the silicate depth to increase with decreasing core fraction of the radio source further supports that for this sample, orientation is the main driver for the difference between radio galaxies and quasars, as predicted by AGN unification. However, comparing our high-z sample with lower redshift 3CRR objects reveals that the absorption of the high-z radio galaxy MIR continuum is lower than expected from a scaled up version of lower luminosity sources, and we discuss some effects that may explain these trends.

In the paradigm of hierarchical galaxy formation, luminous radio galaxies mark mass assembly peaks that should contain clusters of galaxies. Observations of the z=1.53 quasar 3C270.1 with the Spitzer Space Telescope at 3.6-24 micron and with the 6.5-m MMT in the z'- and Y-bands allow detection of potential cluster members via photometric redshifts. Compared with nearby control fields, there is an excess of 11 extremely red objects (EROs) at 1.33 < z_phot < 1.73, consistent with a proto-cluster around the quasar. The spectral energy distributions (SEDs) of 3/4 of the EROs are better fitted with passive elliptical galaxies than withdust-reddened starbursts, and of four sources well-detected on an archival HST snapshot image, all have undisturbed morphologies. However, one ERO, not covered by the HST image, is a double source with 0.8" separation on the z' image and a marginal (2sigma) 24 micron detection indicating a dust-enshrouded starburst. The EROs are more luminous than L* (H = -23.6 AB mag at z=1.5).

Using the Spitzer Space Telescope, we have obtained 3.6--24 micron photometry of 38 radio galaxies and 24 quasars from the 3CR catalog at redshift 1<z<2.5. This 178 MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). Quasar and radio galaxy subsamples matched in isotropic radio luminosity are compared. The quasars all have similar spectral energy distributions (SEDs), nearly constant in nu F_nu through the rest 1.6-10 micron range, consistent with a centrally heated dust distribution which outshines the host galaxy contribution. The radio galaxy SEDs show larger dispersion, but the mean radio galaxy SED declines from rest 1.6 to 3 micron and then rises from 3 to 8 micron. The radio galaxies are on average a factor 3-10 less luminous in this spectral range than the quasars. These characteristics are consistent with composite emission from a heavily reddened AGN plus starlight from the host galaxy. The mid-infrared colors and radio to mid-infrared spectral slopes of individual galaxies are also consistent with this picture. Individual galaxies show different amounts of extinction and host galaxy starlight, consistent with the orientation-dependent unified scheme.

To understand low- to intermediate-mass star-formation in the nearby R CrA molecular cloud, we try to identify the stellar content that is accessible with near-infrared observations. We obtained a JHK band mosaic of 10 x 60 arcmin square covering the entire RCrA molecular cloud with unprecedented sensitivity. We present a catalogue of about 3500 near-infrared sources fainter than the saturation limit K = 10 mag, reaching K = 18mag. We analysed the extended sources by inspecting their morphology and point sources by means of colour-colour and colour-magnitude diagrams. Additionally, we compared the extinction inferred from the NIR data with the line-of-sight dust emission at 1.2 mm. Sources towards high dust emission but relatively low H-K show a projected mm-exces; these sources are either immediately surrounded by cold circumstellar material or, if too red to be a true foreground object, they are embedded in the front layer of the 1.2 mm emitting dust cloud. In both cases they are most likely associated with the cloud.