The precise characteristics of clouds and nature of dust in the diffuse ISM can only be extracted by inspecting the rare cases of single-cloud sightlines, which may be identified by the lack of significant Doppler splitting of interstellar lines observed at resolving power$\lambda /\Delta \lambda > 75,000$. We have searched for such sightlines using high-resolution spectroscopy towards reddened OB stars for which far UV extinction curves are known. We have compiled a sample of 186 spectra with 100 obtained specifically for this project with UVES. From our sample we identified 66 single-cloud sightlines, about half of which were previously unknown. We have used the CH/CH$^+$ line ratio of our targets to establish if the sightlines are dominated by warm or cold clouds. We found that CN is detected in all cold (CH/CH$^+ >1$) clouds, but is frequently absent in warm clouds. We inspected the WISE ($3 - 22 \mu$m) observed emission morphology around our sightlines and excluded a circumstellar nature for the observed dust extinction. We found that most sightlines are dominated by cold clouds that are located far away from the heating source. For 132 stars, we derived the spectral type and the associated spectral type-luminosity distance. We also applied the interstellar Ca II distance scale, and compared both these distance estimates with GAIA parallaxes, finding that these distance estimates scatter by ~40%. By comparing spectrophotometric distances with GAIA we detected for nine sightlines a hidden dust component amounting to a few mag of extinction. This {\it Dark Dust} is populated by $>1 \mu$m large grains and appears predominantly in the field of the cold ISM.

The search for extrasolar planets has driven rapid advances in instrumentation, resulting in cameras such as SPHERE at the VLT, GPI at Gemini South and SCExAO at Subaru, capable of achieving very high contrast ($\sim10^{6}$) around bright stars with small inner working angles ($\sim 0\farcs{1}$). The optimal exploitation of data from these instruments depends on the availability of easy-to-use software to process and analyse their data products. We present a pure-python pipeline, {\sc precision}, which provides fast, memory-efficient reduction of data from the SPHERE/IRDIS near-infrared imager, and can be readily extended to other instruments. We apply {\sc precision} to observations of the extreme red supergiant VX~Sgr, the inner outflow of which is revealed to host complex, asymmetric structure in the near-IR. In addition, optical polarimetric imaging reveals clear extended polarised emission on $\sim0.5^{\prime\prime}$ scales which varies significantly with azimuth, confirming the asymmetry. While not conclusive, this could suggest that the ejecta are confined to a disc or torus, which we are viewing nearly face on, although other non-spherical or clumpy configurations remain possible. VX~Sgr has no known companions, making such a geometry difficult to explain, as there is no obvious source of angular momentum in the system.

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6". We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13", 16.3-18.3 mag at 0.5" and 18.8-20.4 mag at 1.5". The contrast limits closer in (<0.6") depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6"), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far.

Observations of circumstellar environments to look for the direct signal of exoplanets and the scattered light from disks has significant instrumental implications. In the past 15 years, major developments in adaptive optics, coronagraphy, optical manufacturing, wavefront sensing and data processing, together with a consistent global system analysis have enabled a new generation of high-contrast imagers and spectrographs on large ground-based telescopes with much better performance. One of the most productive is the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) designed and built for the ESO Very Large Telescope (VLT) in Chile. SPHERE includes an extreme adaptive optics system, a highly stable common path interface, several types of coronagraphs and three science instruments. Two of them, the Integral Field Spectrograph (IFS) and the Infra-Red Dual-band Imager and Spectrograph (IRDIS), are designed to efficiently cover the near-infrared (NIR) range in a single observation for efficient young planet search. The third one, ZIMPOL, is designed for visible (VIR) polarimetric observation to look for the reflected light of exoplanets and the light scattered by debris disks. This suite of three science instruments enables to study circumstellar environments at unprecedented angular resolution both in the visible and the near-infrared. In this work, we present the complete instrument and its on-sky performance after 4 years of operations at the VLT.

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 next interaction location. In polarised Monte Carlo radiative transfer with aligned non-spherical grains, the nature of dichroism complicates the concept of optical depth. Aims: We investigate in detail the relation between optical depth and the optical properties and density of the attenuating medium in polarised Monte Carlo radiative transfer codes that take into account dichroic extinction. Methods: Based on solutions for the radiative transfer equation, we discuss the optical depth scale in polarised radiative transfer with spheroidal grains. We compare the dichroic optical depth to the extinction and total optical depth scale. Results: In a dichroic medium, the optical depth is not equal to the usual extinction optical depth, nor to the total optical depth. For representative values of the optical properties of dust grains, the dichroic optical depth can differ from the extinction or total optical depth by several ten percent. A closed expression for the dichroic optical depth cannot be given, but it can be derived efficiently through an algorithm that is based on the analytical result corresponding to elongated grains with a uniform grain alignment. Conclusions: Optical depth is more complex in dichroic media than in systems without dichroic attenuation, and this complexity needs to be considered when generating random free path lengths in Monte Carlo radiative transfer simulations. There is no benefit in using approximations instead of the dichroic optical depth.

We present results of the analysis of a sample of 22 stars of spectral types from O7 to B5 and luminosity classes I-V for which Spitzer/IRS spectra are available. The IRS spectra of these stars are examined for signs of excess infrared (IR) emission by comparison with stellar atmospheric spectra. We find that the spectra of half of the studied stars are dominated by excess emission in the far-IR, including all six super- and bright giants. In order to examine the origin of the far-IR excess, we supplement the Spitzer data with optical high-resolution echelle spectroscopy ($\lambda/\Delta \lambda \sim 10^5$), near-IR high-contrast coronagraphic imaging taken with the SPHERE instrument at VLT with a spatial resolution of 0.05", and WISE and Herschel photometry. In the optical region, we detect various absorption and emission lines (H$\alpha$, CIII, and NIII) irrespective of the far-IR excess. Pfund($\alpha$) and Humphrey($\alpha$) lines are observed at the same time as the far-IR excess. These lines are stronger in stars with far-IR excess than in stars without excess. A scattered-light disk in the central r < 2.5" region of the far-IR excess stars HD149404, HD151804, and HD154368 can be excluded from H band imaging down to a 1$\sigma$ contrast of $F(r)/F_{*} \sim 10^{-6}$. The far-IR excess is fit either by a free-free component from ionized gas as for the winds of hot stars or a large (1pc) circumstellar dust shell. The putative dust envelopes required to explain the excess have a visual extinction as low as a few hundred $\mu$-mag.

We describe the Zurich Imaging Polarimeter (ZIMPOL), the visual focal plane subsystem of the SPHERE "VLT planet finder", which pushes the limits of current AO systems to shorter wavelengths, higher spatial resolution, and much improved polarimetric performance. We provide new benchmarks for the performance of high contrast instruments, in particular for polarimetric differential imaging. We have analyzed SPHERE/ZIMPOL point spread functions and measure the peak surface brightness, the encircled energy, and the full width half maximum (FWHM) for different wavelengths, atmospheric conditions, star brightness, and instrument modes. Coronagraphic images are described and analized and the performance for different coronagraphs is compared with tests for the binary alpha Hyi with a separation of 92 mas and a contrast of 6 mag. For the polarimetric mode we made the instrument calibrations using zero polarization and high polarization standard stars and here we give a recipe for the absolute calibration of polarimetric data. The data show a small <1 mas but disturbing differential polarimetric beam shifts, which can be explained as Goos-H\"ahnchen shifts from the inclined mirrors, and we discuss how to correct this effect. The polarimetric sensitivity is investigated with non-coronagraphic and deep, coronagraphic observations of the dust scattering around the symbiotic Mira variable R Aqr. SPHERE/ZIMPOL achieves imaging performances in the visual range with unprecedented characteristics, in particular very high spatial resolution and very high polarimetric contrast. This instrument opens up many new research opportunities for the detailed investigation of circumstellar dust, in scattered and therefore polarized light, for the investigation of faint companions, and for the mapping of circumstellar Halpha emission.

It is well known that the dust properties of the diffuse interstellar medium exhibit variations towards different sight-lines on a large scale. We have investigated the variability of the dust characteristics on a small scale, and from cloud-to-cloud. We use low-resolution spectro-polarimetric data obtained in the context of the Large Interstellar Polarisation Survey (LIPS) towards 59 sight-lines in the Southern Hemisphere, and we fit these data using a dust model composed of silicate and carbon particles with sizes from the molecular to the sub-micrometre domain. Large (> 6 nm) silicates of prolate shape account for the observed polarisation. For 32 sight-lines we complement our data set with UVES archive high-resolution spectra, which enable us to establish the presence of single-cloud or multiple-clouds towards individual sight-lines. We find that the majority of these 35 sight-lines intersect two or more clouds, while eight of them are dominated by a single absorbing cloud. We confirm several correlations between extinction and parameters of the Serkowski law with dust parameters, but we also find previously undetected correlations between these parameters that are valid only in single-cloud sight-lines. We find that interstellar polarisation from multiple-clouds is smaller than from single-cloud sight-lines, showing that the presence of a second or more clouds depolarises the incoming radiation. We find large variations of the dust characteristics from cloud-to-cloud. However, when we average a sufficiently large number of clouds in single-cloud or multiple-cloud sight-lines, we always retrieve similar mean dust parameters. The typical dust abundances of the single-cloud cases are [C]/[H] = 92 ppm and [Si]/[H] = 20 ppm.

Polarimetric studies of light transmitted through interstellar clouds may give constraints on the properties of the interstellar dust grains. Traditionally, broadband linear polarisation (BBLP) measurements have been considered an important diagnostic tool for the study of the interstellar dust, while comparatively less attention has been paid to spectropolarimetric measurements. However, spectropolarimetry may offer stronger constraints than BBLP, for example by revealing narrowband features, and by allowing us to distinguish the contribution of dust from the contribution of interstellar gas. Therefore, we have decided to carry out a Large Interstellar Polarisation Survey (LIPS) using spectropolarimetric facilities in both hemispheres. Here we present the results obtained in the Southern Hemisphere with the FORS2 instrument of the ESO Very Large Telescope. Our spectra cover the wavelength range 380--950\,nm at a spectral resolving power of about 880. We have produced a publicly available catalogue of 127 linear polarisation spectra of 101 targets. We also provide the Serkowski-curve parameters, as well as the wavelength gradient of the polarisation position angle for the interstellar polarisation along 76 different lines of sight. In agreement with previous literature, we found that the best-fit parameters of the Serkowski-curve are not independent of each other. However, the relationships that we obtained are not always consistent with what was found in previous studies.

We present observations of the known edge-on debris disk around HIP 79977 (HD 146897, F star in Upper Sco, 123 pc), taken with the ZIMPOL differential polarimeter of the SPHERE instrument in the Very Broad Band filter ($\lambda_c=735$ nm, $\Delta\lambda=290$ nm) with a spatial resolution of about 25 mas. We measure the polarization flux along and perpendicular to the disk spine of the highly inclined disk for projected separations between 0.2" (25 AU) and 1.6" (200 AU) and investigate the diagnostic potential of such data with model simulations. The polarized flux contrast ratio for the disk is $F_{pol}/F_\ast= (5.5 \pm 0.9) 10^{-4}$. The surface brightness reaches a maximum of 16.2 mag arcsec$^{-2}$ at a separation of $0.2"-0.5"$ along the disk spine with a maximum surface brightness contrast of 7.64 mag arcsec$^{-2}$. The polarized flux has a minimum near the star $<0.2"$ because no or only little polarization is produced by forward or backward scattering in the disk section lying in front of or behind the star. The data are modeled as a circular dust belt with an inclination $i=85(\pm 1.5)^\circ$ and a radius between $r_0$ = 60 AU and 90 AU. The radial density dependence is described by $(r/r_0)^{\alpha}$ with a steep power law index $\alpha=5$ inside $r_0$ and a more shallow index $\alpha=-2.5$ outside $r_0$. The scattering asymmetry factor lies between $g$ = 0.2 and 0.6 adopting a scattering angle-dependence for the fractional polarization as for Rayleigh scattering. Our data are qualitatively very similar to the case of AU Mic and they confirm that edge-on debris disks have a polarization minimum at a position near the star and a maximum near the projected separation of the main debris belt. The comparison of the polarized flux contrast ratio $F_{pol}/F_{\ast}$ with the fractional infrared excess provides strong constraints on the scattering albedo of the dust.

We present polarisation spectra of seven stars in the lines-of-sight towards the Sco OB1 association. Our spectra were obtained within the framework of the Large Interstellar Polarization Survey carried out with the FORS instrument of the ESO VLT. We have modelled the wavelength-dependence of extinction and linear polarisation with a dust model for the diffuse interstellar medium which consists of a mixture of particles with size ranging from the molecular domain of 0.5 nm up to 350 nm. We have included stochastically heated small dust grains with radii between 0.5 and 6 nm made of graphite and silicate, as well as polycyclic aromatic hydrocarbon molecules (PAHs), and we have assumed that larger particles are prolate spheroids made of amorphous carbon and silicate. Overall, a dust model with eight free parameters best reproduces the observations. Reducing the number of free parameters leads to results that are inconsistent with cosmic abundance constraints. We found that aligned silicates are the dominant contributor to the observed polarisation, and that the polarisation spectra are best-fit by a lower limit of the equivolume sphere radius of aligned grains of 70 - 200nm.

R Aqr is a symbiotic binary system consisting of a mira variable, a hot companion with a spectacular jet outflow, and an extended emission line nebula. We have used R Aqr as test target for the visual camera subsystem ZIMPOL, which is part of the new extreme adaptive optics (AO) instrument SPHERE at the Very Large Telescope (VLT). We compare our observations with data from the Hubble Space Telescope (HST) and illustrate the complementarity of the two instruments. We determine from the Halpha emission the position, size, geometric structure, and line fluxes of the jet source and the clouds in the innermost region (<2") of R Aqr and determine Halpha emissivities mean density, mass, recombination time scale, and other cloud parameters. Our data resolve for the first time the R Aqr binary and we measure for the jet source a relative position 46+/-1 mas West of the mira. The central jet source is the strongest Halpha component. North east and south west from the central source there are many clouds with very diverse structures. We see in the SW a string of bright clouds arranged in a zig-zag pattern and, further out, more extended bubbles. In the N and NE we see a bright, very elongated filamentary structure and faint perpendicular "wisps" further out. Some jet clouds are also detected in the ZIMPOL [OI] and He I filters, as well as in the HST line filters for Halpha, [OIII], [NII], and [OI]. We determine jet cloud parameters and find a very well defined anti-correlation between cloud density and distance to the central binary. Future Halpha observations will provide the orientation of the orbital plane of the binary and allow detailed hydrodynamical investigations of this jet outflow and its interaction with the wind of the red giant companion.

Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy.

We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $\mu$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $\mu$m, 10.3 $\mu$m, and 11.6 $\mu$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $\mu$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $\mu$m to 0.31 $\pm$ 0.15 per cent at 11.6 $\mu$m and follows a power-law function, $\lambda^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $\mu$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions.

Context: In the last years, there have been many studies on the omnipresence and structures of filaments in star-forming regions, as well as their role in the process of star formation. Those filaments are normally identified as elongated fibres across the plane of the sky. But how would we detect filaments that are inclined? Aims: We aim to learn more about whether, and how, total column density or dust temperature change with respect to the line of sight. Such variations would enable observers to use dust observations to identify and study filaments at any inclination and gain more insight on the distribution and orientations of filaments within the Galactic plane. Methods: As a first step, we perform numerical calculations on simple cylindrical models to evaluate the influence of filament geometry on the average flux density. After that, we apply our three-dimensional Monte Carlo dust radiative transfer code on two models of star-forming regions and derive maps of effective total column density and dust temperature at different viewing angles. Results: We see only slight changes of average flux density for all cylinders we study. For our more complex models, we find that the effective dust temperature is not sensitive to viewing angle, while the total column density is strongly influenced, with differences exceeding an order of magnitude. The variations are not injective with the viewing angle and depend on the structure of the object. Conclusions: We conclude that there is no single quantity in our analysis that can uniquely trace the inclination and three-dimensional structure of a filament based on dust observations alone. However, observing wide variations in total column density at a given effective dust temperature is indicative of inclined filaments.

Using the Herschel Space Observatory we have observed a representative sample of 87 powerful 3CR sources at redshift z < 1. The far-infrared (FIR, 70-500 micron) photometry is combined with mid-infrared (MIR) photometry from the Wide-Field Infrared Survey Explorer (WISE) and catalogued data to analyse the complete spectral energy distributions (SEDs) of each object from optical to radio wavelength. To disentangle the contributions of different components, the SEDs are fitted with a set of templates to derive the luminosities of host galaxy starlight, dust torus emission powered by active galactic nuclei (AGN) and cool dust heated by stars. The level of emission from relativistic jets is also estimated, in order to isolate the thermal host galaxy contribution. The new data are in line with the orientation-based unification of high-excitation radio-loud AGN, in that the dust torus becomes optically thin longwards of 30 micron. The low excitation radio galaxies and the MIR weak sources represent MIR- and FIR-faint AGN population different from the high-excitation MIR-bright objects; it remains an open question whether they are at a later evolutionary state or an intrinsically different population. The derived luminosities for host starlight and dust heated by star formation are converted to stellar masses and star formation rates (SFR). The host-normalized SFR of the bulk of the 3CR sources is low when compared to other galaxy populations at the same epoch. Estimates of the dust mass yield a 1--100 times lower dust/stellar mass ratio than for the Milky Way, indicating that these 3CR hosts have very low levels of interstellar matter explaining the low level of star formation. Less than 10% of the 3CR sources show levels of star formation above those of the main sequence of star forming galaxies.

Massive stars live short lives, losing large amounts of mass through their stellar wind. Their mass is a key factor determining how and when they explode as supernovae, enriching the interstellar medium with heavy elements and dust. During the red supergiant phase, mass-loss rates increase prodigiously, but the driving mechanism has proven elusive. Here we present high-contrast optical polarimetric-imaging observations of the extreme red supergiant VY Canis Majoris and its clumpy, dusty, mass-loss envelope, using the new extreme-adaptive-optics instrument SPHERE at the VLT. These observations allow us to make the first direct and unambiguous detection of submicron dust grains in the ejecta; we derive an average grain radius $\sim$ 0.5 $\mu$m, 50 times larger than in the diffuse ISM, large enough to receive significant radiation pressure by photon scattering. We find evidence for varying grain sizes throughout the ejecta, highlighting the dynamical nature of the envelope. Grains with 0.5 $\mu$m sizes are likely to reach a safe distance from the eventual explosion of VY Canis Majoris; hence it may inject upwards of 10$^{-2}$ M$_\odot$ of dust into the ISM.

(abridged) The dust content of the universe is primarily explored via its interaction with stellar photons, producing interstellar extinction. However, owing to the physical extension of the observing beam, observations may detect scattered photons, resulting in a change in the observed (or effective) extinction, depending on the spatial distribution of the dust and the resolution of the instrument. We investigate the influence of clumpy dust distributions on effective extinction toward embedded sources and those in the diffuse ISM. We use Monte Carlo radiative transfer to examine effective extinction for various geometries. By varying the number, optical depth and volume-filling factor of clumps in models of spherical shells and the diffuse ISM, we explore the evolution of extinction. Depending on the number of scattering events in the beam, the extinction curve steepens in homogeneous media and flattens in clumpy media. As a result, clumpy dust distributions can to reproduce extinction curves with arbitrary R_v, the effective ratio of total-to-selective extinction. The flattening `washes out' the 2175 \AA bump and shifts the peak to shorter wavelengths. The mean R_v of a shell correlates with the optical depth of an individual clump and the wavelength at which a clump becomes optically thick. Similar behaviour is seen for edge-on discs or tori. However, at grazing inclinations the combination of extinction and strong forward scattering results in chaotic behaviour. Caution is therefore advised when measuring extinction in, for example, AGN tori or toward SNIa or GRB afterglows. In face-on discs, the shape of the scattered continuum changes with clumpiness, however, unlike absorption features, individual features in the scattering cross-sections are preserved. Finally, we show that diffuse interstellar extinction is not modified by scattering on scales of a few kpc.

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 present Herschel (PACS and SPIRE) far-infrared (FIR) photometry of a complete sample of z>1 3CR sources, from the Herschel GT project The Herschel Legacy of distant radio-loud AGN (PI: Barthel). Combining these with existing Spitzer photometric data, we perform an infrared (IR) spectral energy distribution (SED) analysis of these landmark objects in extragalactic research to study the star formation in the hosts of some of the brightest active galactic nuclei (AGN) known at any epoch. Accounting for the contribution from an AGN-powered warm dust component to the IR SED, about 40% of our objects undergo episodes of prodigious, ULIRG-strength star formation, with rates of hundreds of solar masses per year, coeval with the growth of the central supermassive black hole. Median SEDs imply that the quasar and radio galaxy hosts have similar FIR properties, in agreement with the orientation-based unification for radio-loud AGN. The star-forming properties of the AGN hosts are similar to those of the general population of equally massive non-AGN galaxies at comparable redshifts, thus there is no strong evidence of universal quenching of star formation (negative feedback) within this sample. Massive galaxies at high redshift may be forming stars prodigiously, regardless of whether their supermassive black holes are accreting or not.

We present Gran Telescopio CANARIAS CanariCam 8.7$\mu$m imaging and 7.5-13$\mu$m spectroscopy of six local systems known to host an active galactic nucleus (AGN) and have nuclear star formation. Our main goal is to investigate whether the molecules responsible for the 11.3$\mu$m polyclyclic aromatic hydrocarbon (PAH) feature are destroyed in the close vicinity of an AGN. We detect 11.3$\mu$m PAH feature emission in the nuclear regions of the galaxies as well as extended PAH emission over a few hundred parsecs. The equivalent width (EW) of the feature shows a minimum at the nucleus but increases with increasing radial distances, reaching typical star-forming values a few hundred parsecs away from the nucleus. The reduced nuclear EW are interpreted as due to increased dilution from the AGN continuum rather than destruction of the PAH molecules. We conclude that at least those molecules responsible for the 11.3$\mu$m PAH feature survive in the nuclear environments as close as 10pc from the AGN and for Seyfert-like AGN luminosities. We propose that material in the dusty tori, nuclear gas disks, and/or host galaxies of AGN is likely to provide the column densities necessary to protect the PAH molecules from the AGN radiation field.

We present a revised model for the infrared emission of the hyperluminous infrared galaxy IRAS F10214+4724 which takes into account recent photometric data from Spitzer and Herschel that sample the peak of its spectral energy distribution. We first present and discuss a grid of smooth active galactic nucleus (AGN) torus models computed with the method of Efstathiou & Rowan-Robinson and demonstrate that the combination of these models and the starburst models of Efstathiou and coworkers, while able to give an excellent fit to the average spectrum of Seyfert 2s and spectra of individual type 2 quasars measured by Spitzer, fails to match the spectral energy distribution of IRAS F10214+4724. This is mainly due to the fact that the nuSnu distribution of the galaxy falls very steeply with increasing frequency (a characteristic that is usually indicative of heavy absorption by dust) but shows a silicate feature in emission. Such emission features are not expected in sources with optical/near-infrared type 2 AGN spectral signatures. The Herschel data show that there is more power emitted in the rest-frame 20-50um wavelength range compared with the model presented by Efstathiou which assumes three components of emission: an edge-on torus, clouds (at a temperature of 610 and 200K) that are associated with the narrow-line region (NLR) and a highly obscured starburst that dominates in the submillimetre. We present a revised version of that model that assumes an additional component of emission which we associate with NLR clouds at a temperature of 100K. The 100K dust component could also be explained by a highly obscured hot starburst. The model suggests that the NLR of IRAS F10214+4724 has an unusually high covering factor (>= 17%) or more likely the magnification of the emission from the NLR clouds is significantly higher than that of the emission from the torus.

We present a model for the diffuse interstellar dust that explains the observed wavelength-dependence of extinction, emission, linear and circular polarisation of light. The model is set-up with a small number of parameters. It consists of a mixture of amorphous carbon and silicate grains with sizes from the molecular domain of 0.5 up to about 500nm. Dust grains with radii larger than 6nm are spheroids. Spheroidal dust particles have a factor 1.5 - 3 larger absorption cross section in the far IR than spherical grains of the same volume. Mass estimates derived from submillimeter observations that ignore this effect are overestimated by the same amount. In the presence of a magnetic field, spheroids may be partly aligned and polarise light. We find that polarisation spectra help to determine the upper particle radius of the otherwise rather unconstrained dust size distribution. Stochastically heated small grains of graphite, silicates and polycyclic aromatic hydrocarbons (PAHs) are included. We tabulate parameters for PAH emission bands in various environments. They show a trend with the hardness of the radiation field that can be explained by the ionisation state or hydrogenation coverage of the molecules. For each dust component its relative weight is specified, so that absolute element abundances are not direct input parameters. The model is confronted with the average properties of the Milky Way, which seems to represent dust in the solar neighbourhood. It is then applied to four specific sight lines including the reflection nebula NGC2023. For these sight lines, we present linear and circular spectro-polarimetric observations obtained with FORS/VLT. Using prolate rather than oblate grains gives a better fit to observed spectra; the axial ratio of the spheroids is typically two and aligned silicates are the dominant contributor to the polarisation.

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 present a new model for the infrared emission of the high redshift hyperluminous infrared galaxy IRAS F10214+4724 which takes into account recent photometric data from Spitzer and Herschel that sample the peak of its spectral energy distribution. We first demonstrate that the combination of the AGN tapered disc and starburst models of Efstathiou and coworkers, while able to give an excellent fit to the average spectrum of type 2 AGN measured by Spitzer, fails to match the spectral energy distribution of IRAS F10214+4724. This is mainly due to the fact that the nuSnu distribution of the galaxy falls very steeply with increasing frequency (a characteristic of heavy absorption by dust) but shows a silicate feature in emission. We propose a model that assumes two components of emission: clouds that are associated with the narrow-line region and a highly obscured starburst. The emission from the clouds must suffer significantly stronger gravitational lensing compared to the emission from the torus to explain the observed spectral energy distribution.

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.

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.

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.

METIS is a mid-infrared instrument proposed for the European Extremely Large Telescope (E-ELT). It is designed to provide imaging and spectroscopic capabilities in the 3 - 14 micron region up to a spectral resolution of 100000. One of the novel concepts of METIS is that of a high-resolution integral field spectrograph (IFS) for a diffraction-limited mid-IR instrument. While this concept has many scientific and operational advantages over a long-slit spectrograph, one drawback is that the spectral resolution changes over the field of view. This has an impact on the procedures to correct for telluric absorption lines imprinted on the science spectra. They are a major obstacle in the quest to maximize spectral fidelity, the ability to distinguish a weak spectral feature from the continuum. The classical technique of division by a standard star spectrum, observed in a single IFS spaxel, cannot simply be applied to all spaxels, because the spectral resolution changes from spaxel to spaxel. Here we present and discuss possible techniques of telluric line correction of METIS IFS spectra, including the application of synthetic model spectra of telluric transmission, to maximize spectral fidelity.

We present observations and models for one of these MYSO candidates, NGC3603 IRS 9A. Our goal is to investigate with infrared interferometry the structure of IRS 9A on scales as small as 200AU, exploiting the fact that a cluster of O and B stars has blown away much of the obscuring foreground dust and gas. Observations in the N-band were carried out with the MIDI beam combiner attached to the VLTI. Additional interferometric observations which probe the structure of IRS 9A on larger scales were performed with an aperture mask installed in the T-ReCS instrument of Gemini South. The spectral energy distribution (SED) is constrained by the MIDI N-band spectrum and by data from the Spitzer Space Telescope. Our efforts to model the structure and SED of IRS 9A range from simple geometrical models of the brightness distribution to one- and two-dimensional radiative transfer computations. The target is resolved by T-ReCS, with an equivalent (elliptical) Gaussian width of 330mas by 280mas (2300 AU by 2000 AU). Despite this fact, a warm compact unresolved component was detected by MIDI which is possibly associated with the inner regions of a flattened dust distribution. Based on our interferometric data, no sign of multiplicity was found on scales between about 200AU and 700AU projected separation. A geometric model consisting of a warm (1000 K) ring (400 AU diameter) and a cool (140 K) large envelope provides a good fit to the data. No single model fitting all visibility and photometric data could be found, with disk models performing better than spherical models. While the data are clearly inconsistent with a spherical dust distribution they are insufficient to prove the existence of a disk but rather hint at a more complex dust distribution.

We present a technique to synthesise telluric absorption and emission features both for in-situ wavelength calibration and for their removal from astronomical spectra. While the presented technique is applicable for a wide variety of optical and infrared spectra, we concentrate in this paper on selected high-resolution near-infrared spectra obtained with the CRIRES spectrograph to demonstrate its performance and limitation. We find that synthetic spectra reproduce telluric absorption features to about 2%, even close to saturated line cores. Thus, synthetic telluric spectra could be used to replace the observation of telluric standard stars, saving valuable observing time. This technique also provides a precise in-situ wavelength calibration, especially useful for high-resolution near-infrared spectra in the absence of other calibration sources.

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.

We present results of performance modelling for METIS, the Mid-infrared European Extremely Large Telescope (E-ELT) Imager and Spectrograph. Designed by a consortium of NOVA (Netherlands), UK Astronomy Technology Centre (UK), MPIA Heidelberg (Germany), CEA Saclay (France) and KU Leuven (Belgium), METIS will cover the atmospheric windows in L, M and N-band and will offer imaging, medium-resolution slit spectroscopy (R~1000-3000) and high-resolution integral field spectroscopy (R~100,000). Our model uses a detailed set of input parameters for site characteristics and atmospheric profiles, optical design, thermal background and the most up-to-date IR detector specifications. We show that METIS will bring an orders-of-magnitude level improvement in sensitivity and resolution over current ground-based IR facilities, bringing mid-IR sensitivities to the micro-Jansky regime. As the only proposed E-ELT instrument to cover this entire spectral region, and the only mid-IR high-resolution integral field unit planned on the ground or in space, METIS will open up a huge discovery space in IR astronomy in the next decade.

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.

Mid-infrared imaging at resolutions of 300 mas of the central kpc region of 13 nearby, well-known active galaxies is presented. The bulk of the mid-IR emission is concentrated on an unresolved central source within a size of less than 5 to 130 pc, depending on the object distance. Further resolved emission is detected in 70% of the sample in the form of circumnuclear star-forming rings or diffuse nuclear extended emission. In the three cases with circumnuclear star formation, the stellar contribution is at least as important as that of the AGN. In those with extended nuclear emission -- a third of the sample -- this emission represents a few per cent of the total measured; however, this contribution may be underestimated because of the chopped nature of these observations. This extended emission is generally collimated in a preferential direction often coinciding with that of the extended ionized gas or the jet. In all cases, the nuclear fluxes measured at 11.8 and 18.7 micron represent a minor contribution of the flux levels measured by large aperture IRAS data at the nearest energy bands of 12 and 25 micron. This contribution ranges from 30% to less than 10%. In only three cases do the AGN fluxes agree with IRAS to within a factor of 2. In the AGNs with strong circumnuclear star formation, this component can well account for most of the IRAS flux measured in these objects. But in all other cases, either a low surface brightness component extending over galactic scales or strong extra-nuclear IR sources -- e.g. HII regions in spiral arms -- have to be the main source of the IRAS emission. In either case, the contribution of these components dwarfs that of the AGN at mid-IR wavelengths.

Polycyclic aromatic hydrocarbon (PAH) emission has been found in both starbursts and modestly starforming galaxies, but the relation between starforming activity and PAH luminosity is still a matter of debate. The different correlation degrees could be caused by the variety of optical and far-infrared sample selection criteria. In order to obtain a census of the typical properties of PAH emitting galaxies, we here study moderately distant galaxies which have been selected by their PAH emission. Combining the ISOCAM Parallel Survey at 6.7 micron with 2MASS we have colour-selected a sample of 120 candidates for strong PAH emission. We obtained optical and mid-infrared spectra of 75 and 19 sources, respectively, and analysed IRAS-ADDSCANs and available Spitzer 3.6-160 micron photometry. The Spitzer spectra exhibit clear PAH features and corroborate that our photometric selection criteria trace the PAH emission of galaxies fairly well. The optical spectra show emission lines, at median redshift z~0.1, as well as Hdelta and CaII absorption, indicating ongoing starformation as well as post-starbursts. The mid- and far-infrared spectral energy distributions (SEDs) provide evidence for a broad range of far-infrared (FIR) luminosities, but in general the dust is colder than for starburst galaxies like M82. For most galaxies the monocromatic luminosity peaks at about equal height at optical, 6.7 micron (PAH) and FIR wavelengths. In about 15% of the sources the FIR luminosity exceeds the optical and PAH energy output by about a factor 5-10 despite the cool dust temperature; in these galaxies a large dust mass of 10^8 - 10^9 M_sun is inferred.

We present radiative transfer models for submillimeter galaxies with spectroscopic redshifts and mid-infrared spectroscopy from Spitzer/IRS and analyze available Spitzer/MIPS 24, 70 and 160mu data. We use two types of starburst models, a cirrus model and a model for the emission of an AGN torus in order to investigate the nature of these objects. We find that for three of the objects (25%) cirrus emission alone can account for the mid-infrared spectrum and the MIPS and submillimeter data. For the remaining objects we find that we need a combination of starburst and cirrus in order to fit simultaneously the multi--wavelength data. We find that the typical submillimeter galaxy has comparable luminosity in the starburst (median L=10^12.5 Lo) and cirrus (median L=10^12.4 Lo) components. This could arise if the galaxies have been forming stars continuously for the last 250Myr with the star formation occurring in the last 5Myr being shrouded by high-optical-depth molecular cloud dust, whereas the rest of the starlight is attenuated by diffuse dust or cirrus with an Av of about 1mag.

The formation and evolution of the Milky Way bulge is not yet well understood and its classification is ambiguous. Constraints can, however, be obtained by studying the abundances of key elements in bulge stars. The aim of this study is to determine the chemical evolution of CNO, and a few other elements in stars in the Galactic bulge, and to discuss the sensitivities of the derived abundances from molecular lines. High-resolution, near-IR spectra in the H band were recorded using VLT/CRIRES. Due to the high and variable visual extinction in the line-of-sight towards the bulge, an analysis in the near-IR is preferred. The CNO abundances can all be determined simultaneously from the numerous molecular lines in the wavelength range observed. The three giant stars in Baade's window presented here are the first bulge stars observed with CRIRES. We have especially determined the CNO abundances, with uncertainties of less than 0.20 dex, from CO, CN, and OH lines. Since the systematic uncertainties in the derived CNO abundances due to uncertainties in the stellar fundamental parameters, notably Teff, are significant, a detailed discussion of the sensitivities of the derived abundances is included. We find good agreement between near-IR and optically determined O, Ti, Fe, and Si abundances. Two of our stars show a solar [C+N/Fe], suggesting that these giants have experienced the first dredge-up and that the oxygen abundance should reflect the original abundance of the giants. The two giants fit into the picture, in which there is no significant difference between the O abundance in bulge and thick-disk stars. Our determination of the S abundances is the first for bulge stars. The high [S/Fe] values for all the stars indicate a high star-formation rate in an early phase of the bulge evolution.

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.

METIS, the Mid-infrared ELT Imager and Spectrograph (formerly called MIDIR), is a proposed instrument for the European Extremely Large Telescope (E-ELT), currently undergoing a phase-A study. The study is carried out within the framework of the ESO-sponsored E-ELT instrumentation studies. METIS will be designed to cover the E-ELT science needs at wavelengths longward of 3um, where the thermal background requires different operating schemes. In this paper we discuss the main science drivers from which the instrument baseline has been derived. Specific emphasis has been given to observations that require very high spatial and spectral resolution, which can only be achieved with a ground-based ELT. We also discuss the challenging aspects of background suppression techniques, adaptive optics in the mid-IR, and telescope site considerations. The METIS instrument baseline includes imaging and spectroscopy at the atmospheric L, M, and N bands with a possible extension to Q band imaging. Both coronagraphy and polarimetry are also being considered. However, we note that the concept is still not yet fully consolidated. The METIS studies are being performed by an international consortium with institutes from the Netherlands, Germany, France, United Kingdom, and Belgium.

We present a refined orbital solution for the components A, B, and C of the nearby late-M type multiple system LHS 1070. By combining astrometric datapoints from NACO/VLT, CIAO/SUBARU, and PUEO/CFHT, as well as a radial velocity measurement from the newly commissioned near infrared high-resolution spectrograph CRIRES/VLT, we achieve a very precise orbital solution for the B and C components and a first realistic constraint on the much longer orbit of the A-BC system. Both orbits appear to be co-planar. Masses for the B and C components calculated from the new orbital solution (M_(B+C) = 0.157 +/- 0.009 M_sun) are in excellent agreement with theoretical models, but do not match empirical mass-luminosity tracks. The preliminary orbit of the A-BC system reveals no mass excess for the A component, giving no indication for a previously proposed fourth (D) component in LHS 1070.

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 present measurements of the fluorine abundance in a Galactic Bulge Asymptotic Giant Branch (AGB) star. The measurements were performed using high resolution K-band spectra obtained with the CRIRES spectrograph, which has been recently installed at ESO's VLT, together with state-of-the-art model atmospheres and synthetic spectra. This represents the first fluorine abundance measurement in a Galactic Bulge star, and one of few measurements of this kind in a third dredge-up oxygen-rich AGB star. The F abundance is found to be close to the solar value scaled down to the metallicity of the star, and in agreement with Disk giants that are comparable to the Bulge giant studied here. The measurement is of astrophysical interest also because the star's mass can be estimated rather accurately ($1.4 \lesssim M/\mathrm{M}_{\sun} \lesssim 2.0$). AGB nucleosynthesis models predict only a very mild enrichment of F in such low mass AGB stars. Thus, we suggest that the fluorine abundance found in the studied star is representative for the star's natal cloud, and that fluorine must have been produced at a similar level in the Bulge and in the Disk.

In the course of the NIR/MIR AGN search combining the 6.7 mu ISOCAM Parallel Survey and 2MASS we have discovered 24 type-1 quasars about a third of which are too red to be discriminated by optical/UV search techniques. Here we report on a detailed case study of the reddest type-1 quasar of our sample (J2341) at redshift z=0.236 with M_K=-25.8 and J-K=1.95. We performed spectroscopy in the optical with VLT/FORS1 and in the MIR with Spitzer as well as NIR imaging with ISPI at CTIO. The optical and NIR observations reveal a star forming emission-line galaxy at the same redshift as the quasar with a projected linear separation of 1.8 arcsec (6.7 kpc). The quasar and its companion are embedded in diffuse extended continuum emission. Compared with its companion the quasar exhibits redder optical-NIR colours, which we attribute to hot nuclear dust. The MIR spectrum shows only few emission lines superimposed on a power-law spectral energy distribution. However, the lack of strong FIR emission suggests that our potentially interacting object contains much less gas and dust and is in a stage different from dust reddened ULIRG-AGN like Mrk 231. The optical spectrum shows signatures for reddening in the emission-lines and no post-starburst stellar population is detected in the host galaxy of the quasar. The optical continuum emission of the active nucleus appears absorbed and diluted. Even the combination of absorption and host dilution is not able to match J2341 with standard quasar templates. While the BLR shows only a rather moderate absorption of E_(B-V)=0.3, the continuum shorter than 4500 AA requires strong obscuration with E_(B-V)=0.7, exceeding the constraints from the low upper limit on the 9.7 mu silicate absorption. This leads us to conclude that the continuum of J2341 is intrinsically redder than that of typical quasars.

Half of the Seyfert-2 galaxies escaped detection of broad lines in their polarised spectra observed so far. Some authors have suspected that these non-HBLRs contain real Sy2 nuclei without intrinsic broad line region hidden behind a dust torus. If this were true, then their nuclear structure would fundamentally differ from that of Sy2s with polarised broad lines: in particular, they would not be explained by orientation-based AGN unification. Further arguments for two physically different Sy2 populations have been derived from the warm and cool IRAS F25/F60 ratios. These ratios, however, refer to the entire host galaxies and are unsuitable to conclusively establish the absence of a nuclear dust torus. Instead, a study of the Seyfert-2 dichotomy should be performed on the basis of nuclear properties only. Here we present the first comparison between [OIII] 5007A and mid-infrared imaging at matching spatial resolution. Exploring the Seyfert-2 dichotomy we find that the distributions of nuclear mid-infrared/[OIII] luminosity ratios are indistinguishable for Sy1s and Sy2s with and without detected polarised broad lines and irrespective of having warm or cool IRAS F25/F60 ratios. We find no evidence for the existence of a population of real Sy2s with a deficit of nuclear dust emission. Our results suggest 1) that all Seyfert nuclei possess the same physical structure including the putative dust torus and 2) that the cool IRAS colours are caused by a low contrast of AGN to host galaxy. Then the Seyfert-2 dichotomy is explained in part by unification of non-HBLRs with narrow-line Sy1s and to a larger rate by observational biases caused by a low AGN/host contrast and/or an unfavourable scattering geometry.

We show individual high resolution spectra of components A, B, and C of the nearby late-M type multiple system LHS 1070. Component A is a mid-M star, B and C are known to have masses at the threshold to brown dwarfs. From our spectra we measure rotation velocities and the mean magnetic field for all three components individually. We find magnetic flux on the order of several kilo-Gauss in all components. The rotation velocities of the two late-M objects B and C are similar (vsini = 16km/s), the earlier A component is spinning only at about half that rate. This suggests weakening of net rotational braking at late-M spectral type, and that the lack of slowly rotating late-M and L dwarfs is real. Furthermore, we found that magnetic flux in the B component is about twice as strong as in component C at similar rotation rate. This indicates that rotational braking is not proportional to magnetic field strength in fully convective objects, and that a different field topology is the reason for the weak braking in low mass objects.

OO Serpentis is a deeply embedded pre-main sequence star that went into outburst in 1995 and gradually faded afterwards. Its eruption resembled the well-known FU Orionis-type or EX Lupi-type outbursts. Since very few such events have ever been documented at infrared wavelengths, our aim is to study the temporal evolution of OO Ser in the infrared. OO Ser was monitored with the Infrared Space Observatory starting 4 months after peak brightness and covering 20 months. In 2004-2006 we again observed OO Ser from the ground and complemented this dataset with archival Spitzer obsevations also from 2004. We analysed these data with special attention to source confusion and constructed light curves at 10 different wavelengths as well as spectral energy distributions. The outburst caused brightening in the whole infrared regime. According to the infrared light curves, OO Ser started a wavelength-independent fading after peak brightness. Later the flux decay became slower but stayed wavelength-independent. The fading is still ongoing, and current fading rates indicate that OO Ser will not return to quiescent state before 2011. The outburst timescale of OO Ser seems to be shorter than that of FUors, but longer than that of EXors. The outburst timescale and the moderate luminosity suggest that OO Ser is different from both FUors and EXors, and shows similarities to the recently erupted young star V1647 Ori. Based on its spectral energy distribution and bolometric temperature, OO Ser seems to be an early class I object, with an age of < 10^5 yr. The object is probably surrounded by an accretion disc and a dense envelope. Due to the shorter outburst timescales, the viscosity in the circumstellar disc of OO Ser is probably an order of magnitude higher than usual for FUors.

Context: A long-standing challenge of observational AGN research is to find type 2 quasars, the luminous analogues of Seyfert-2 galaxies. Aims: We search for luminous narrow-line type 2 AGN, characterise their properties, and compare them with broad-line type 1 AGN. Methods: Combining the ISOCAM parallel survey at 6.7 mu with 2MASS, we have selected AGN via near-mid-infrared colours caused by the hot nuclear dust emission. We performed spectroscopy in the optical and, for a subset of the sample, also in the mid-infrared with Spitzer. Results: We find nine type 2 AGN at redshift 0.1<z<0.5, three of them have even quasar-like [OIII] luminosities. At the given redshift and luminosity range the number of type 2 AGN is at least as high as that of type 1s. At z>0.5 we did not find type 2 AGN, probably because the hottest dust emission, still covered by the NIR filters, is obscured. The optical spectra of the type 2 host galaxies show young and old stellar populations. Only one object is an ultraluminous infrared galaxy with starburst. The 5-38 mu spectra of the two type 2 sources observed show a strong continuum with PAH emission in one case and silicate absorption in the other case. Conclusions: The near-mid-infrared selection is a successful strategy to find luminous type 2 AGN at low z. The objects exhibit a large range of properties so that it is difficult to infer details by means of popular SED fitting with simple average templates.

We investigate the morphology and the evolutionary stage of the hypercompact HII region M17-UC1 using observations at infrared wavelengths and NIR radiative transfer modeling. For the first time, this region is resolved into two emission areas separated by a dark lane reminiscent of an obscuring silhouette of a circumstellar disk. So far, the observational data as well as model calculations suggest that M17-UC1 is surrounded by a disk of cool dust. This direct direct detection of a circumstellar disk candidate around a hypercompact HII region is in agreement with the expectations of the disk accretion model for high-mass star formation.

We have observed seven powerful FR2 radiogalaxies and seven quasars with the Spitzer IRS. Both samples are comparable in both, redshift range and isotropic 178 Hz luminosity. Both samples are found to have similar distributions in the luminosity ratios of Mid-IR high- and low-excitation lines ([NeV]/[NeII]), and of Mid-IR high-excitation line to radio power ratio ([NeV]/P_178MHz). However, the MIR/FIR ratio is generally higher for quasars. We further observed Silicate features at 10 and 18micron in emission. In our sample only quasars show emission features, while silicate absorption is seen only in the radio galaxies. These observations are all in agreement with unification schemes that explain both groups as the same class of objects seen under different orientation angles.

We provide a library of some 7000 SEDs for starbursts and ultra luminous galaxies (http://www.eso.org/~rsiebenm/sb_models). Its purpose is to quickly obtain estimates of the basic parameters, such as luminosity, size and dust or gas mass 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. The objects may be in the local universe or at high z. We calculate the radiative transfer in spherical symmetry for a stellar cluster permeated by an interstellar medium with standard (Milky Way) dust properties. The cluster contains two stellar populations: old bulge stars and OB stars. Because the latter are young, a certain fraction of them will be embedded in compact clouds which constitute hot spots that determine the MIR fluxes. We present SEDs for a broad range of luminosities, sizes and obscurations. We argue that the assumption of spherical symmetry and the neglect of clumpiness of the medium are not severe shortcomings for computing the dust emission. The validity of the approach is demonstrated by matching the SED of the best studied galaxies, including M82 and Arp220, by library elements; one example is shown for a galaxy at high redshift (z ~ 3). Generally, one finds an element which fits the observed SED very well, and the parameters defining the element are in full accord with what is known about the galaxy from detailed studies.

We present a high-resolution (R~16,000) spectrum and a narrow-band image centered on the [NeII]12.8 micron line of the central kpc region of the starburst/Seyfert2 galaxy NGC 7582. The galaxy has a rotating circum-nuclear starburst disk, shown at great detail at a diffraction-limited resolution of 0.4 arcsec. The high spatial resolution allows us to probe the dynamics of the [NeII] gas in the nuclear regions, and to estimate the mass of the central black hole. We construct models of gas disks rotating in the combined gravitational potential from the stellar bulge and a central black hole, and derive a black hole mass of 5.5 x 10^7 solar masses with a 95% confidence interval of [3.6,8.1] x 10^7 solar masses. The black hole mass combined with stellar velocity dispersion measurements from the literature shows that the galaxy is consistent with the local M-sigma relation. This is the first time that a black hole mass in a galaxy except our own Milky Way system has been estimated from gas dynamics in the mid-infrared. We show that spatially resolved mid-infrared spectroscopy may be competitive with similar techniques in the optical and near-infrared, and may prove to be important for estimating black hole masses in galaxies with strong nuclear dust obscuration. The high spectral resolution allows us to determine the heliocentric systemic velocity of the galaxy to between 1614 and 1634 km/s. The mid-infrared image reveals several dense knots of dust-embedded star formation in the circum-nuclear disk, and we briefly discuss its morphology.

We provide a library of some 7000 SEDs (available at www.eso.org/~rsiebenm) for the nuclei of starburst and ultra luminous galaxies. Its purpose is to quickly obtain estimates of the basic parameters, such as luminosity, size and dust or gas mass 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. The objects may be in the local universe or at high z. We calculate the radiative transfer in spherical symmetry for a stellar cluster permeated by an interstellar medium with standard (Milky Way) dust properties. The cluster contains two stellar populations: old bulge stars and OB stars. Because the latter are young, a certain fraction of them will be embedded in compact clouds which constitute hot spots that determine the MIR fluxes. We present SEDs for a broad range of luminosities, sizes and obscurations. We argue that the assumption of spherical symmetry and the neglect of clumpiness of the medium are not severe shortcomings for computing the dust emission. The validity of the approach is demonstrated by matching the SED of seven of the best studied galaxies, including M82 and Arp220, by library elements. In all cases, one finds an element which fits the observed SED very well, and the parameters defining the element are in full accord with what is known about the galaxy from detailed studies. We also compare our method of computing SEDs with other techniques described in the literature.

With the Spitzer Space Telescope Infrared Spectrograph (IRS) we have observed seven powerful FR2 radio galaxies and seven quasars. Both samples, the galaxies and the quasars, are comparable in isotropic 178 MHz luminosity (10^26.5 W/Hz < P 178MHz < 10^29.5 W/Hz) and in redshift range (0.05 < z < 1.5). We find for both samples similar distributions in the luminosity ratios of high- to low-excitation lines ([NeV] 14.3 mu / [NeII] 12.8 mu) and of high-excitation line to radio power ([NeV] 24.3 \mu / P 178MHz). This solves the long debate about the apparent difference of quasars and radio galaxies in favor of the orientation-dependent unified schemes. Furthermore, the luminosity ratio [OIII] 500.7 nm / [OIV] 25.9 mu of most galaxies is by a factor of ten lower than that of the quasars. This suggests that the optical emission from the central NLR is essentially absorbed (A_V} > 3) in the powerful FR2 galaxies and that the [OIII] 500.7 nm luminosity does not serve as isotropic tracer for testing the unified schemes.

We combined the ISOCAM Parallel Mode Survey at 6.7 micron (LW2 filter) with the Two Micron All Sky Survey in order to obtain a powerful tool to search for AGN independent of dust extinction. Using moderate colour criteria H-K>0.5 and K-LW2>2.7 we have selected a sample of 77 AGN candidates in an effective area of about 10 square degrees. By means of optical spectroscopy we find 24 (30%) type-1 QSOs at redshifts 0.1<z<2.3; nine of them have z>0.8. About one third of the ISO-2MASS QSOs show so red optical colours, that they are missed in optical and UV AGN surveys like SDSS, 2DF, or HES. With a surface density of about 2 deg^(-2) down to R<18 mag the ISO-2MASS QSOs outnumber the 1.35 deg^(-2) of the SDSS quasar survey by 50%; we find a combined optical-IR QSO surface density of 2.7 deg^(-2). Since only two of the ISO-2MASS QSOs have also J-K>2, the inclusion of the ISO mid-infrared photometry significantly extends the capabilities of the pure 2MASS red AGN survey. We suggest that the newly found red AGN resemble young members of the quasar population, and that quasars spend much of their lifetime in a dust enshrouded phase.

Present AGN and starburst models aiming to account for the observed infrared SEDs consider a physical description of the dust and a solution of the radiative transfer problem. MIR spectra obtained at different spatial scales (SST-IRS, ISO and TIMMI2) are presented. They show that PAH bands are detected in starburst regions but significantly reduced near the centre of AGN. This is explained by examining the heating mechanism of PAHs after hard (FUV, X-ray) photon interactions. Economic radiative transfer models of starbursts and AGN are made available. The successful application of the starburst model is demonstrated by fitting broad band data and detailed Spitzer spectra of NGC7714. The AGN model is applied to ISO data of a sample of 68 radio galaxies and quasars of the 3CR catalogue. Radiative transfer models of galaxies with Hidden Broad Line Regions are shown. Their SED enable us to separate the contributions from the dusty disc of the AGN and the starbursts. The composite model is consistent with the unified scheme and the idea that the infrared emission of AGN is dominated in the MIR and starbursts in the FIR. According to the unified scheme, the observed diversity of AGN properties results from the different orientations relative to our line of sight. The strong resonance of silicate dust at 10mic is therefore, as expected, seen in absorption towards many type-2 AGN. In type-1 AGN, it should be seen in emission because the hot inner surface of the dust torus becomes visible. However, this has not been observed so far, thus challenging the unification scheme or leading to exotic modifications of the dust-torus model. Here the recent discovery of the 10mic silicate feature in emission in luminous quasar spectra is reported.

We present ISO mid-infrared spectrophotometry and far-infrared photometry of galaxies with Hidden Broad Line Regions (HBLR). We also present radiative transfer models of their spectral energy distributions which enable us to separate the contributions from the dusty disc of the AGN and the dusty starbursts. We find that the combination of tapered discs (discs whose thickness increases with distance from the central source in the inner part but stays constant in the outer part) and starbursts provide good fits to the data. The tapered discs dominate in the mid-infrared part of the spectrum and the starbursts in the far-infrared. After correcting the AGN luminosity for anisotropic emission we find that the ratio of the AGN luminosity to the starburst luminosity, L(AGN)/L(SB), ranges from about unity for IRAS14454-4343 to about 13 for IRAS01475-0740. Our results suggest that the warm IRAS colours of HBLR are due to the relatively high L(AGN)/L(SB). Our fits are consistent with the unified model and the idea that the infrared emission of AGN is dominated by a dusty disc in the mid-infrared and starbursts in the far-infrared.

We have imaged $\epsilon$ Indi B, the closest brown dwarf binary known, with VISIR at the VLT in three narrow-band mid-infrared bandpasses located around 8.6$\mu$m, 10.5$\mu$m and 11.3$\mu$m. We are able to spatially resolve both components, and determine accurate mid-infrared photometry for both components independently. In particular, our VISIR observations probe the NH$_3$ feature in the atmospheres of the cooler and warmer brown dwarfs. For the first time, we can disentangle the contributions of the two components, and find that % our photometry of $\epsilon$ IndiBb is in good agreement with recent ``cloud-free'' atmosphere models having an effective temperature of $T_\mathrm{eff}=800$ K. With an assumed age of 1 Gyr for the $\epsilon$ Indi system, component Ba agrees more with $T_\mathrm{eff} \approx 1100$ K rather than with $T_\mathrm{eff}=1200$ K, as suggested by SPITZER spectroscopic observations of the combined $\epsilon$ Indi B system (Roellig et al., 2004). Even higher effective temperatures appear inconsistent with our absolute photometry, as they would imply an unphysical small size of the brown dwarf $\epsilon$ IndiBa.

According to the unified scheme, AGN are surrounded by a dust-torus, and the observed diversity of AGN properties results from the different orientations relative to our line of sight. The strong resonance of silicate dust at 10 micron is therefore, as expected, seen in absorption towards many type-2 AGN. In type-1 AGN, it should be seen in emission because the hot inner surface of the dust torus becomes visible. However, this has not been observed so far, thus challenging the unification scheme or leading to exotic modifications of the dust-torus model. Here we report the discovery of the 10 micron silicate feature in emission in two luminous quasars with the Infrared Spectrograph of the Spitzer Space Telescope.

Present AGN models aiming to account for the observed infrared SEDs consider a physical description of the dust and a solution of the radiative transfer problem. Mid infrared spectra obtained at different spatial scales are presented. They show that PAH bands are detected in starburst regions but significantly reduced near the centre of AGN. This may be explained by examining the heating mechanism of PAHs after hard photon interactions. On the radiative transfer side first a most economic model is presented where three parameters, luminosity, effective size and extinction of the nucleus are varied to obtain SED fits. A full grid of model spectra are made available at: "http://www.eso.org/~rsiebenm/agn_models". This model is sufficient to account for ISO broad band data of a sample of 68 radio galaxies and quasars of the 3CR catalogue. The hot dust component detected is mainly due to small grains and PAHs. In such models, type 1 AGNs are represented by a compact dust distribution with warm grains and weak PAH emission. In AGNs of type 2, the dust appears to be more extended, relatively colder and PAH bands are strong. Realistic AGN models which are consistent with the unification need to explain the overall absence of the 9.7mic. silicate emission feature. This can be done by considering various geometries (tapered discs). Models which combine AGN and starburst activity are presented for galaxies with hidden broad line region. It is found that the AGN torus dominate the mid IR continuum emission and that the starbursts dominate the PAH band as well as the far infrared and submillimeter emission.

We compare three 3CR samples of 11 FR1 galaxies, 17 FR2 galaxies and 18 lobe-dominated quasars contained in the ISO Data Archive. In contrast to the powerful FR2 galaxies with edge-brightened lobes, the low radio power FR1 galaxies in our sample do not exhibit any high MIR or FIR dust luminosity, which is typical for a buried, intrinsically more luminous AGN. This consolidates the fact already inferred from optical studies that their AGNs have only a relatively low luminosity. Also the FR1 galaxies show a high FIR/MIR luminosity ratio, compared to quasars, suggesting that their FIR luminosity is substantially powered by the interstellar radiation field (ISRF) of the giant elliptical hosts. Finally, we discuss the FR1-FR2 morphological dichotomy. FR1 galaxies do not have more interstellar matter (ISM) than FR2s as traced -- on the large scale -- by the cool FIR emitting dust and -- in the nuclear region -- by the warm MIR emitting dust. Due to the lack of central gas we suggest that the black holes of our FR1 galaxies are fed at a lower accretion rate than those of the FR2 galaxies.

We obtained N-band observations of the Apollo asteroid 2002 NY40 during its close Earth fly-by in August 2002 with TIMMI2 at the ESO 3.6 m telescope. The photometric measurement allowed us to derive a radiometric diameter of 0.28+/-0.03 km and an albedo of 0.34+/-0.06 through the near-Earth asteroid thermal model (NEATM) and a thermophysical model (TPM). The values are in agreement with results from radar data, visual and near-IR observations. In this first comparison between these two model approaches we found that the empirical NEATM beaming parameter $\eta$=1.0 corresponds to a thermal inertia values of about 100 $\mathrm{J m^{-2} s^{-0.5} K^{-1}}$ for a typical range of surface roughness, assuming an equator-on viewing angle. Our TPM analysis indicated that the surface of 2002 NY40 consists of rocky material with a thin or no dust regolith. The asteroid very likely has a prograde sense of rotation with a cold terminator at the time of our observations. Although both model approaches can fit the thermal spectra taken at phase angles of 22$^{\circ}$ and 59$^{\circ}$, we did not find a consistent model solution that describes all pieces of photometric and spectroscopic data. In addition to the 2002 NY40 analysis, we discuss the possibilities to distinguish between different models with only very few photometric and/or spectroscopic measurements spread over a range of phase angles.

We present the complete set of ISOPHOT observations of 3CR radio galaxies and quasars, which are contained in the ISO Data Archive, providing 75 mid- and far-infrared spectral energy distributions (SEDs) between 5 and 200 micron. For 28 sources they are supplemented with MAMBO 1.2 mm observations and for 15 sources with new submillimetre data from the SCUBA archive. We check the orientation-dependent unified scheme, in which the powerful FR2 narrow line galaxies are quasars viewed at high inclination, so that their nuclei are hidden behind a dust torus intercepting the optical-ultraviolet AGN radiation and reemitting it in the infrared. We find that (1) both the quasars and the galaxies show a high mid- to far-infrared luminosity ratio typical for powerful AGNs and (2) -- when matched in 178 MHz luminosity -- both show the same ratio of isotropic far-infrared to isotropic 178 MHz lobe power. Therefore, from our large sample investigated here we find strong evidence for the orientation-dependent unification of the powerful FR2 galaxies with the quasars.

Since a large fraction of active galactic nuclei (AGN) is missed in common UV-excess surveys and is even hard to find in radio, near-IR and X-ray surveys, we have used a new AGN selection technique which is expected to be not affected by extinction. Within the scientific verification of the ISOCAM Parallel Survey at 6.7 micron we have discovered objects with exceptional mid-infrared (MIR) emission. They are essentially not detected on IRAS-ADDSCANs and only very few of them show up in the NVSS and FIRST radio surveys. Various colour criteria of the 6.7 micron data with 2MASS and optical wavebands show that the sources reach more extreme IR colours than the sources in the Hubble Deep Field-South and the ELAIS survey. The comparison with known object types suggests that we have found AGN with a pronounced MIR emission, probably due to circum-nuclear dust. First results from optical spectroscopy of ten candidates corroborate this interpretation showing four AGN, two reddened LINER and four extremely reddened emission-line galaxies with MIR/FIR flux ratios higher than for known pure starburst galaxies. The results will make a significant contribution to the debate on the entire AGN population.

We present a survey of all 3CR sources imaged with ISOCAM onboard the Infrared Space Observatory (ISO). For the source, we present spatially integrated mid--infrared (MIR, 5 - 18mic.) fluxes measured from newly calibrated ISOCAM images. In total, we detected 68 AGN of the 3CR catalogue, at redshifts z < 2.5. The one with the highest redshift is 4C+72.26 at z = 3.53. ISOCAM data are combined with other photometric measurements to construct the spectral energy distribution (SED) from optical to radio wavelengths. In order to describe dust emission we apply new radiative transfer models. By varying three parameters, luminosity, effective size and extinction, we obtain a fit to the SED for our objects. In the MIR the hot dust component is mainly due to small grains and PAHs. In the models, a type 1 AGN is represented by a compact dust distribution, the dust is therefore very warm and emission of PAHs is weak because of photo--destruction. In AGNs of type 2, the dust is relatively colder but PAH bands are strong.

We have measured mid-infrared (MIR) fluxes of 3C sources on images taken with ISOCAM onboard the Infrared Space Observatory (ISO). The photometric data were combined with existing photometry at other wavelengths to assemble the spectral energy distribution (SED) for each galaxy from infrared to radio wavelengths. In addition, we used ISOPHOT spectra to compute average MIR spectra for different types of Active Galactic Nuclei (AGNs). The MIR emission of 53 of our sources shows evidence of dust. We find a clear correlation between the SEDs, emission from polycyclic aromatic hydrocarbons (PAHs), and the type of AGN. Specifically, we find that in broad emission line radio galaxies (BLRGs), the dust emission peaks at ~40 microns and PAH bands are weak, whereas emission in comparable narrow line radio galaxies (NLRGs) peaks at a longer wavelength of \~100microns and the PAH bands are much stronger. Although less pronounced, a similar trend is seen when comparing quasars (QSOs) with high-luminosity NLRGs. We used dust radiative transfer models with a central heating source to describe the SEDs. The difference in the dust emission appears to be the effect of heating by radiation from the central engine. In type~1 AGNs (BLRGs and QSOs), the broad line region (BLR) is directly visible and therefore hot (T>~300K) dust dominates the emission. In this region, PAHs are destroyed which explains the weak PAH emission. On the other hand, the BLR is hidden in type~2 AGNs (NLRGs). Their SEDs are therefore dominated by cooler dust including PAHs.

We have re-analyzed the 6-12 micron ISO spectrum of the ultra-luminous infrared galaxy Arp220 with the conclusion that it is not consistent with that of a scaled up version of a typical starburst. Instead, both template fitting with spectra of the galaxies NGC4418 and M83 and with dust models suggest that it is best represented by combinations of a typical starburst component, exhibiting PAH emission features, and a heavily absorbed dust continuum which contributes ~40% of the 6-12 micron flux and likely dominates the luminosity. Of particular significance relative to previous studies of Arp220 is the fact that the emission feature at 7.7 micron comprises both PAH emission and a broader component resulting from ice and silicate absorption against a heavily absorbed continuum. Extinction to the PAH emitting source, however, appears to be relatively low. We tentatively associate the PAH emitting and heavily dust/ice absorbed components with the diffuse emission region and the two compact nuclei respectively identified by Soifer et al. (2002) in their higher spatial resolution 10 micron study. Both the similarity of the absorbed continuum with that of the embedded Galactic protostars and results of the dust models imply that the embedded source(s) in Arp220 could be powered by, albeit extremely dense, starburst activity. Due to the high extinction, it is not possible with the available data to exclude that AGN(s) also contribute some or all of the observed luminosity. In this case, however, the upper limit measured for its hard X-ray emission would require Arp220 to be the most highly obscured AGN known.

We investigate the mid-infrared radiation of galaxies that are powered by a starburst or by an AGN. For this end, we compare the spectra obtained at different spatial scales in a sample of infrared bright galaxies. ISO observations which include emission of the nucleus as well as most of the host galaxy are compared with TIMMI2 spectra of the nuclear region. We find that ISO spectra are generally dominated by strong PAH bands. However, this is no longer true when inspecting the mid-infrared emission of the pure nucleus. Here PAH emission is detected in starbursts whereas it is significantly reduced or completely absent in AGNs. A physical explanation of these new observational results is presented by examining the temperature fluctuation of a PAH after interaction with a photon. It turns out that the hardness of the radiation field is a key parameter for quantifying the photo-destruction of small grains. Our theoretical study predicts PAH evaporation in soft X-ray environments. Radiative transfer calculations of clumpy starbursts and AGN corroborate the observational fact that PAH emission is connected to starburst activity whereas PAHs are destroyed near an AGN. The radiative transfer models predict for starbursts a much larger mid-infrared size than for AGN. This is confirmed by our TIMMI2 acquisition images: We find that the mid-infrared emission of Seyferts is dominated by a compact core while most of the starbursts are spatially resolved.

We present polarimetric imaging of the protoplanetary nebula RAFGL 2688 obtained at 4.5 microns with the Infrared Space Observatory (ISO). We have deconvolved the images to remove the signature of the point spread function of the ISO telescope, to the extent possible. The deconvolved 4.5 micron image and polarimetric map reveal a bright point source with faint, surrounding reflection nebulosity. The reflection nebula is brightest to the north-northeast, in agreement with previous ground- and space-based infrared imaging. Comparison with previous near-infrared polarimetric imaging suggests that the polarization of starlight induced by the dust grains in RAFGL 2688 is more or less independent of wavelength between 2 microns and 4.5 microns. This, in turn, indicates that scattering dominates over thermal emission at wavelengths as long as ~5 microns, and that the dust grains have characteristic radii < 1 micron.

The mid infrared polarisation properties of four Ultraluminous Infrared Galaxies (ULIRGs) have been investigated by broad band filter observations with the ISOCAM instrument on board the Infrared Space Observatory (ISO). The wavelength region from 5 to 18 mic. was selected where the emission from the putative torus peaks. We report detection of polarisation for all ULIRGs studied. The fractional polarisation ranges from about 3% up to 8%. The highest polarisation is recorded in Mrk231 which has a clear AGN signature, whereas the lowest is for Arp220, which is generally thought to be powered predominantly by star formation. We discuss the various mechanisms that could give rise to the polarisation and conclude that the most likely interpretation is that it is due to magnetically aligned elongated dust grains. This is the same mechanism believed to be operating in a number of galactic sources. The position angle of polarisation could give the projected magnetic field direction and therefore constrain models for the formation of the tori.

The paper is devoted to the understanding of the infrared emission of nuclear regions in galaxies. a) ISO data of NGC1808 are presented: spectro--photometry from 5.1 to 16.4 mic., a 25"x25" map at 6 mic. and 170mic. photometry. b) The data are complemented by a polarization measurement at 170 mic. (2.5+/-0.4% at position angle 94+/-5deg) and a map at 6 mic. In the map, the degree of polarisation goes up to 20% in the outer regions. We argue that the polarisation is produced by emission of big grains and exclude very small grains and PAHs or scattering and extinction. c) The mid infrared spectrum shows, beside the main emission bands, a so far unknown plateau of PAH features in the >13 mic.region. d) The total spectrum can be fit under the assumption of optically thin emission. However, such a model fails to reproduce the 25mic. point and implies that the mid infrared is due to very small grains and PAHs. These particles would then also have to be responsible for the 6 mic. polarisation, which is unlikely. e) To avoid these difficulties, we successfully turn to a radiative transfer model whose major feature is the existence of (hot spots) produced by the dust clouds around OB stars. We demonstrate the decisive influence on the mid infrared spectrum of both the PAHs and the hot spots.

We present photometric ISO 60 and 170um measurements, complemented by some IRAS data at 60um, of a sample of 84 nearby main-sequence stars of spectral class A, F, G and K in order to determine the incidence of dust disks around such main-sequence stars. Of the stars younger than 400 Myr one in two has a disk; for the older stars this is true for only one in ten. We conclude that most stars arrive on the main sequence surrounded by a disk; this disk then decays in about 400 Myr. Because (i) the dust particles disappear and must be replenished on a much shorter time scale and (ii) the collision of planetesimals is a good source of new dust, we suggest that the rapid decay of the disks is caused by the destruction and escape of planetesimals. We suggest that the dissipation of the disk is related to the heavy bombardment phase in our Solar System. Whether all stars arrive on the main sequence surrounded by a disk cannot be established: some very young stars do not have a disk. And not all stars destroy their disk in a similar way: some stars as old as the Sun still have significant disks.

We probe the region around the protostar HH108MMS by deep mid infrared photometric and polarimetric imaging. The protostar is detected at 14mu in absorption against the diffuse background. Next to HH108MMS, we find a second absorbing core, named Q1, and the young stellar object IRAS18331-0035 which is more advanced in its evolution and already seen in emission at 12mu and 14mu. HH108MMS, Q1 and IRAS18331-0035 form a triplet along an extended filamentary absorption feature. From the variation of the surface brightness across the source, we derive for HH108MMS and Q1 the optical depth and density profile. As HH108MMS is also seen in the millimeter dust emission, we can derive the ratio of the dust extinction coefficients at 14mu and 1.3mm and obtain kappa(14mu) / kappa(1300mu) ~ 470. Because models for the dust in the diffuse interstellar medium predict a ratio of around 2000, our value points to fluffy composite grains which are expected to prevail in dense and cold environments. First mid infrared polarisation images of pre--stellar absorbing cores are presented. At 12mu and 14mu the polarisation of the region around HH108MMS is strong (> 15%) and tightly correlated with the source triplet. We demonstrate that the high degree of polarisation can be explained by extinction of rotationally aligned dust particles of moderate elongation.

The mid infrared spectra of the starbursts show the 9.7$\mu m$ silicate absorption feature and strong mid infrared emission bands centered at 6.2, 7.7, and 11.3 $\mu$m. Illustrative models of the {\it active} galaxies are presented: As the starbursts are most likely confined to the central region of the galaxy the radiative transfer in the nucleus is solved under the assumption of spherical symmetry. The distribution of stars and dust are adjusted until the complete infrared spectrum of the galaxies are modeled. The dust is described as a mixture of large grains, very small grains and PAH, which are undergoing temperature fluctuations. Although the galactic nuclei are deeply hidden in the dust its global structure can be estimated by the simple calculations presented. ISOPHOT samples of eleven {\it active} galaxies and seven {\it inactive} spirals are presented. The far infrared and submillimeter spectrum of the {\it active} galaxies can be described by a single modified black--body at a color temperature of $31.5 \pm 2.8$ K. This leads to a ratio of infrared luminosity to gas mass, $L_{\rm IR}/M_{\rm gas}$, of $\sim 90$ L$_\odot$/ M$_\odot$. In contrast, the spectral energy distributions of {\it inactive} spirals require, apart from warm dust of $31.8 \pm 2.8$ K, an additional very cold component of at most $12.9\pm 1.7$ K. This implies a $L_{\rm IR}/M_{\rm gas}$ ratio of $\sim 3$ L$_\odot$/ M$_\odot$ for the {\it inactive} spirals, a factor $\sim 30$ lower than for the {\it active} galaxies. The detection of such cold dust can be predicted by radiative transfer models.

We present mid IR spectro-photometric imaging of a sample of eight nearby ($D \leq 240$pc) Herbig Ae/Be stars. The spectra are dominated by photospheric emission (HR6000), featureless infrared excess emission (T~Cha), broad silicate emission feature (HR5999) and the infrared emission bands (HD 97048, HD 97300, TY~CrA, HD 176386). The spectrum of HD179218 shows both silicate emission and infrared emission bands (IEB). All stars of our sample where the spectrum is entirely dominated by IEB have an extended emission on scales of a few thousand AU ($\sim 10''$). We verify the derived source extension found with ISOCAM by multi--aperture photometry with ISOPHT and compare our ISOCAM spectral photometry with ISOSWS spectra.

We present illustrative models for the UV to millimeter emission of starburst galaxies which are treated as an ensemble of optically thick giant molecular clouds (GMCs) centrally illuminated by recently formed stars. The models follow the evolution of the GMCs due to the ionization-induced expansion of the HII regions and the evolution of the stellar population within the GMC according to the Bruzual & Charlot stellar population synthesis models. The effect of transiently heated dust grains/PAHs to the radiative transfer, as well as multiple scattering, is taken into account. The expansion of the HII regions and the formation of a narrow neutral shell naturally explains why the emission from PAHs dominates over that from hot dust in the near to mid-IR, an emerging characteristic of the infrared spectra of starburst galaxies. The models allow us to relate the observed properties of a galaxy to its age and star formation history. We find that exponentially decaying 10^7-10^8 yrs old bursts can explain the IRAS colours of starburst galaxies. The models are also shown to account satisfactorily for the multiwavelength data on the prototypical starburst galaxy M82 and NGC6090, a starburst galaxy recently observed by ISO. In M82 we find evidence for two bursts separated by 10^7yrs. In NGC6090 we find that at least part of the far-IR excess may be due to the age of the burst (6.4 x 10^7yrs). We also make predictions about the evolution of the luminosity of starbursts at different wavelengths which indicate that far-IR surveys may preferentially detect older starbursts than mid-IR surveys.

We present ISO photometric and spectroscopic observations of a sample of 57 bright Asymptotic Giant Branch stars and red supergiants in the Large Magellanic Cloud, selected on the basis of IRAS colours indicative of high mass-loss rates. PHOT-P and PHOT-C photometry at 12, 25 and 60 $\mu$m and CAM photometry at 12 $\mu$m are used in combination with quasi-simultaneous ground-based near-IR photometry to construct colour-colour diagrams for all stars in our sample. PHOT-S and CAM-CVF spectra in the 3 to 14 $\mu$m region are presented for 23 stars. From the colour-colour diagrams and the spectra, we establish the chemical types of the dust around 49 stars in this sample. Many stars have carbon-rich dust. The most luminous carbon star in the Magellanic Clouds has also a (minor) oxygen-rich component. OH/IR stars have silicate absorption with emission wings. The unique dataset presented here allows a detailed study of a representative sample of thermal-pulsing AGB stars with well-determined luminosities.

We describe ISO observations of the obscured Asymptotic Giant Branch (AGB) star IRAS04496-6958 in the Large Magellanic Cloud (LMC). This star has been classified as a carbon star. Our new ISOCAM CVF spectra show that it is the first carbon star with silicate dust known outside of the Milky Way. The existence of this object, and the fact that it is one of the highest luminosity AGB stars in the LMC, provide important information for theoretical models of AGB evolution and understanding the origin of silicate carbon stars.

We present a radiative-transfer calculation which reproduces the infrared spectrum of the planetary nebula BD~+30$^{\circ}$3639. We calculate the transfer process through absorption and scattering in a spherical-symmetric multi-grain dust shell. The emission of transiently heated particles is taken into account, as well as polycyclic aromatic hydrocarbons. We obtain an acceptable fit to most of the spectrum, including the PAH infrared bands. At submillimetre wavelengths the observed emission is larger than the model predicts, indicating that large dust conglomerates (``f{}luffy grains'') may be needed as an additional constituent. The fit favours a distance of $\ge 2 \,$kpc, which implies that BD~+30$^\circ$3639 has evolved from a massive progenitor of several solar masses. A low dust-to-gas mass ratio is found in the ionised region. The calculations yield an original mass-loss rate of $2\times10^{-4} \msolar \peryr$ on the Asymptotic Giant Branch. Using this mass-loss rate, we calculate how the infrared spectrum has evolved during the post-AGB evolution. We show in particular the evolution of the IRAS colours during the preceding post-AGB evolution.

We present a radiative-transfer calculation which reproduces the infrared spectrum of the planetary nebula BD+30$^{\circ}$3639, fitting both the spectral energy distribution and the spatial extent at infrared wavelengths. We obtain an acceptable fit to most of the spectrum, including the infrared bands. The fit requires a distance of $\ge 2 \kpc$, which implies that BD+30$^\circ$3639 has evolved from a massive progenitor of several solar masses. Two surprising results are (1) a very low dust-to-gas ratio, (2) an absence of the smallest PAH molecules. Extrapolating back in time, we calculate the previous infrared evolution of BD+30$^\circ$3639.