Aims: ARCiS, a novel code for the analysis of exoplanet transmission and emission spectra is presented. The aim of the modelling framework is to provide a tool able to link observations to physical models of exoplanet atmospheres. Methods: The modelling philosophy chosen in this paper is to use physical and chemical models to constrain certain parameters while keeping free the parts where our physical understanding is still more limited. This approach, in between full physical modelling and full parameterisation, allows us to use the processes we understand well and parameterise those less understood. A Bayesian retrieval framework is implemented and applied to the transit spectra of a set of 10 hot Jupiters. The code contains chemistry and cloud formation and has the option for self consistent temperature structure computations. Results: The code presented is fast and flexible enough to be used for retrieval and for target list simulations for e.g. JWST or the ESA Ariel missions. We present results for the retrieval of elemental abundance ratios using the physical retrieval framework and compare this to results obtained using a parameterised retrieval setup. Conclusions: We conclude that for most of the targets considered the current dataset is not constraining enough to reliably pin down the elemental abundance ratios. We find no significant correlations between different physical parameters. We confirm that planets in our sample with a strong slope in the optical transmission spectrum are the planets where we find cloud formation to be most active. Finally, we conclude that with ARCiS we have a computationally efficient tool to analyse exoplanet observations in the context of physical and chemical models.

Having a short orbital period and being tidally locked makes WASP-43b an ideal candidate for JWST observations. Phase curve observations of an entire orbit will enable the mapping of the atmospheric structure across the planet, with different wavelengths of observation allowing different atmospheric depths to be seen. We provide insight into the details of the clouds that may form on WASP-43b in order to prepare the forthcoming interpretation of the JWST and follow-up data. We utilize 3D GCM results as input for a kinetic, non-equilibrium model for mineral cloud particles, and for a kinetic model to study a photochemicaly-driven hydrocarbon haze component. Mineral condensation seeds form throughout the atmosphere of WASP-43b. This is in stark contrast to the ultra-hot Jupiters, like WASP-18b and HAT-P-7b. The dayside is loaded with few but large mineral cloud particles in addition to hydrocarbon haze particles of comparable abundance. Photochemically driven hydrocarbon haze appears on the dayside, but does not contribute to the cloud formation on the nightside. The geometrical cloud extension differs across the globe due to the changing thermodynamic conditions. Day and night differ by 6000km in pressure scale height. As reported for other planets, the C/O is not constant throughout the atmosphere. The mean molecular weight is approximately constant in a H2-dominated WASP-43b. WASP-43b is expected to be fully covered in clouds which are not homogeneously distributed throughout the atmosphere. The dayside and the terminator clouds will be a combination of mineral particles of locally varying size and composition, and of hydrocarbon hazes. The optical depth of hydrocarbon hazes is considerably lower than that of mineral cloud particles such that a wavelength-dependent radius measurement of WASP-43b would be determined by the mineral cloud particles but not by hazes.

Exoplanet atmosphere characterisation has become an important tool in understanding exoplanet formation, evolution. However, clouds remain a key challenge for characterisation: upcoming space telescopes (e.g. JWST, ARIEL) and ground-based high-resolution spectrographs (e.g. CRIRES+) will produce data requiring detailed understanding of cloud formation and cloud effects. We aim to understand how the micro-porosity of cloud particles affects the cloud structure, particle size, and material composition. We examine the spectroscopic effects of micro-porous particles, the particle size distribution, and non-spherical cloud particles. We expanded our kinetic non-equilibrium cloud formation model and use a grid of prescribed 1D (T_gas-p_gas) DRIFT-PHOENIX profiles. We applied the effective medium theory and the Mie theory to model the spectroscopic properties of cloud particles with micro-porosity and a derived particle size distribution. We used a statistical distribution of hollow spheres to represent the effects of non-spherical cloud particles. Highly micro-porous cloud particles (90% vacuum) have a larger surface area, enabling efficient bulk growth higher in the atmosphere than for compact particles. Increases in single-scattering albedo and cross-sectional area for these mineral snowflakes cause the cloud deck to become optically thin only at a wavelength of ~100 ${\rm \mu m}$ instead of at the ~20 ${\rm \mu m}$ for compact cloud particles. A significant enhancement in albedo is also seen when cloud particles occur with a locally changing Gaussian size distribution. Non-spherical particles increase the opacity of silicate spectral features, which further increases the wavelength at which the clouds become optically thin. JWST MIRI will be sensitive to signatures of micro-porous and non-spherical cloud particles based on the wavelength at which clouds are optically thin.

We have conducted a re-analysis of publicly available Hubble Space Telescope Wide Field Camera 3 (HST WFC3) transmission data for the hot-Jupiter exoplanet WASP-43b, using the Bayesian retrieval package Tau-REx. We report evidence of AlO in transmission to a high level of statistical significance (> 5-sigma in comparison to a flat model, and 3.4-sigma in comparison to a model with H2O only). We find no evidence of the presence of CO, CO2, or CH4 based on the available HST WFC3 data or on Spitzer IRAC data. We demonstrate that AlO is the molecule that fits the data to the highest level of confidence out of all molecules for which high-temperature opacity data currently exists in the infrared region covered by the HST WFC3 instrument, and that the subsequent inclusion of Spitzer IRAC data points in our retrieval further supports the presence of AlO. H2O is the only other molecule we find to be statistically significant in this region. AlO is not expected from the equilibrium chemistry at the temperatures and pressures of the atmospheric layer that is being probed by the observed data. Its presence therefore implies direct evidence of some disequilibrium processes with links to atmospheric dynamics. Implications for future study using instruments such as the James Webb Space Telescope (JWST) are discussed, along with future opacity needs. Comparisons are made with previous studies into WASP-43b.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

Our aim was to directly detect the thermal emission of the putative exoring system responsible for the complex deep transits observed in the light curve for the young Sco-Cen star 1SWASP J140747.93-394542.6 (V1400 Cen, hereafter J1407), confirming it as the occulter seen in May 2007, and to determine its orbital parameters with respect to the star. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the field centred on J1407 in the 340 GHz (Band 7) continuum in order to determine the flux and astrometric location of the ring system relative to the star. We used the VLT/NACO camera to observe the J1407 system in March 2019 and to search for the central planetary mass object at thermal infrared wavelengths. We detect no point source at the expected location of J1407, and derive an upper limit $3\sigma$ level of $57.6~\mu\rm{Jy}$. There is a point source detected at an angular separation consistent with the expected location for a free-floating ring system that occulted J1407 in May 2007, with a flux of $89~\mu\rm{Jy}$ consistent with optically thin dust surrounding a massive substellar companion. At 3.8 microns with the NACO camera, we detect the star J1407 but no other additional point sources within 1.3 arcseconds of the star, with a lower bound on the sensitivity of $6M_{Jup}$ at the location of the ALMA source, and down to $4M_{Jup}$ in the sky background limit. The ALMA upper limit at the location of J1407 implies that a hypothesised bound ring system is composed of dust smaller than $1\rm{~mm}$ in size, implying a young ring structure. The detected ALMA source has multiple interpretations, including: (i) it is an unbound substellar object surrounded by warm dust in Sco-Cen with an upper mass limit of $6M_{Jup}$, or (ii) it is a background galaxy.

Dust grains are the building {blocks} of future planets. They evolve in size, shape and composition during the life cycle of the interstellar medium. We seek to understand the process which leads from diffuse medium grains to dust grains in the vicinity of protostars inside disks. As a first step, we propose to characterize the dust evolution inside pre-stellar cores thanks to multi-wavelength observations. We will present how NIKA2 maps are crucial to better constrain dust properties and {we will} introduce SIGMA: a new flexible dust model in open access.

Post-AGB binaries are surrounded by circumbinary disks of gas and dust that are similar to protoplanetary disks found around young stars. We aim to understand the structure of these disks and identify the physical phenomena at play in their very inner regions. We want to understand the disk-binary interaction and to further investigate the comparison with protoplanetary disks. We have conducted an interferometric snapshot survey of 23 post-AGB binaries in the near-infrared (H-band) using VLTI/PIONIER. We have fitted the multiwavelength visibilities and closure phases with purely geometrical models with an increasing complexity in order to retrieve the sizes, temperatures and flux ratios of the different components All sources are resolved and the different components contributing to the H-band flux are dissected. The environment of these targets is very complex: 13/23 targets need models with thirteen or more parameters to fit the data. We find that the inner disk rims follow and extend the size-luminosity relation established for disks around young stars with an offset toward larger sizes. The measured temperature of the near-infrared circumstellar emission of post-AGB binaries is lower (Tsub~1200K) than for young stars, probably due to a different dust mineralogy and/or gas density in the dust sublimation region. The dusty inner rims of the circumbinary disks around post-AGB binaries are ruled by dust sublimation physics. Additionally, a significant amount of the circumstellar H-band flux is over-resolved (14 targets have more than 10% of their non-stellar flux over-resolved) hinting for more structure from a yet unknown origin (disk structure or outflow). The amount of over-resolved flux is larger than around young stars. Due to the complexity of these targets, interferometric imaging is a necessary tool to reveal the interacting inner regions in a model-independent way.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

The scattering properties of the dust originating from debris discs are still poorly known. The analysis of scattered light is however a powerful remote-sensing tool to understand the physical properties of dust particles orbiting other stars. Scattered light is indeed widely used to characterise the properties of cometary dust in the solar system. We aim to measure the morphology and scattering properties of the dust from the debris ring around HR4796A in polarised optical light. We obtained high-contrast polarimetric images of HR4796A in the wavelength range 600-900nm with the SPHERE / ZIMPOL instrument on the Very Large Telescope. We measured for the first time the polarised phase function of the dust in a debris system over a wide range of scattering angles in the optical. We confirm that it is incompatible with dust particles being compact spheres under the assumption of the Mie theory, and propose alternative scenarios compatible with the observations, such as particles with irregular surface roughness or aggregate particles.

Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating/cooling balance between dust and gas and, of course, chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines we present a homogeneous, panchromatic collection of data on a sample of 85 T Tauri and Herbig Ae objects for which data cover a range from X-rays to centimeter wavelengths. Datasets consist of photometric measurements, spectra, along with results from the data analysis such as line fluxes from atomic and molecular transitions. Additional properties resulting from modeling of the sources such as disc mass and shape parameters. dust size and PAH properties are also provided for completeness. Targets were selected based on their properties data availability. Data from more than 50 different telescopes and facilities were retrieved and combined in homogeneous datasets directly from public data archives or after being extracted from more than 100 published articles. X-ray data for a subset of 56 sources represent an exception as they were reduced from scratch and are presented here for the first time. Compiled datasets along with a subset of continuum and emission-line models are stored in a dedicated database and distributed through a publicly accessible online system. All datasets contain metadata descriptors that allow to backtrack them to their original resources. The graphical user interface of the online system allows the user to visually inspect individual objects but also compare between datasets and models. It also offers to the user the possibility to download any of the stored data and metadata for further processing.

Understanding of clouds is instrumental in interpreting current and future spectroscopic observations of exoplanets. Modelling clouds consistently is complex, since it involves many facets of chemistry, nucleation theory, condensation physics, coagulation, and particle transport. We develop a simple physical model for cloud formation and transport, efficient and versatile enough that it can be used in modular fashion for parameter optimization searches of exoplanet atmosphere spectra. The transport equations are formulated in 1D, accounting for sedimentation and diffusion. The grain size is obtained through a moment method. For simplicity, only one cloud species is considered and the nucleation rate is parametrized. From the resulting physical profiles we simulate transmission spectra covering the visual to mid-IR wavelength range. We apply our models towards KCl clouds in the atmosphere of GJ1214 b and towards MgSiO3 clouds of a canonical hot-Jupiter. We find that larger cloud diffusivity $K_{zz}$ increases the thickness of the cloud, pushing the $\tau=1$ surface to a lower pressure layer higher in the atmosphere. A larger nucleation rate also increases the cloud thickness while it suppresses the grain size. Coagulation is most important at high nuclei injection rates ($\dot\Sigma_n$) and low $K_{zz}$. We find that the investigated combinations of $K_{zz}$ and $\dot\Sigma_n$ greatly affect the transmission spectra in terms of the slope at near-IR wavelength (a proxy for grain size), the molecular features seen at ~1\micr (which disappear for thick clouds, high in the atmosphere), and the 10\micr silicate feature, which becomes prominent for small grains high in the atmosphere. The result of our hybrid approach -- aimed to provide a good balance between physical consistency and computational efficiency -- is ideal towards interpreting (future) spectroscopic observations of exoplanets.

The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana).

Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal cloud variations, circumplanetary disks, atmospheric oblateness, and cloud particle properties on the integrated degree and direction of polarization in the near-infrared. We have developed a three-dimensional Monte Carlo radiative transfer code (ARTES) for scattered light simulations in (exo)planetary atmospheres. The code is applicable to calculations of reflected light and thermal radiation in a spherical grid with a parameterized distribution of gas, clouds, hazes, and circumplanetary material. The disk-integrated degree of polarization of a horizontally-inhomogeneous atmosphere is maximal when the planet is flattened, the optical thickness of the equatorial clouds is large compared to the polar clouds, and the clouds are located at high altitude. For a flattened planet, the integrated polarization can both increase or decrease with respect to a spherical planet which depends on the horizontal distribution and optical thickness of the clouds. The direction of polarization can be either parallel or perpendicular to the projected direction of the rotation axis when clouds are zonally distributed. Rayleigh scattering by submicron-sized cloud particles will maximize the polarimetric signal whereas the integrated degree of polarization is significantly reduced with micron-sized cloud particles as a result of forward scattering. The presence of a cold or hot circumplanetary disk may also produce a detectable degree of polarization ($\lesssim$1%) even with a uniform cloud layer in the atmosphere.

We introduce a new modelling framework called FLiTs to simulate infrared line emission spectra from protoplanetary discs. This paper focuses on the mid-IR spectral region between 9.7 um to 40 um for T Tauri stars. The generated spectra contain several tens of thousands of molecular emission lines of H2O, OH, CO, CO2, HCN, C2H2, H2 and a few other molecules, as well as the forbidden atomic emission lines of SI, SII, SIII, SiII, FeII, NeII, NeIII, ArII and ArIII. In contrast to previously published works, we do not treat the abundances of the molecules nor the temperature in the disc as free parameters, but use the complex results of detailed 2D ProDiMo disc models concerning gas and dust temperature structure, and molecular concentrations. FLiTs computes the line emission spectra by ray tracing in an efficient, fast and reliable way. The results are broadly consistent with R=600 Spitzer/IRS observational data of T Tauri stars concerning line strengths, colour, and line ratios. In order to achieve that agreement, however, we need to assume either a high gas/dust mass ratio of order 1000, or the presence of illuminated disc walls at distances of a few au. The molecules in these walls cannot be photo-dissociated easily by UV because of the large densities in the walls which favour their re-formation. Most observable molecular emission lines are found to be optically thick, rendering a standard analysis with column densities difficult. We find that the difference between gas and dust temperatures in the disc surface is important for the line formation. We briefly discuss the effects of C/O ratio and choice of chemical rate network on these results. Our analysis offers new ways to infer the chemical and temperature structure of T Tauri discs from future JWST/MIRI observations, and to possibly detect secondary illuminated disc walls based on their specific mid-IR molecular signature.

Context. Kepler observations have revealed a class of short period exoplanets, of which Kepler-1520 b is the prototype, which have comet-like dust tails thought to be the result of small, rocky planets losing mass. The shape and chromaticity of the transits constrain the properties of the dust particles originating from the planet's surface, offering a unique opportunity to probe the composition and geophysics of rocky exoplanets. Aims. We aim to approximate the average Kepler long-cadence light curve of Kepler-1520 b and investigate how the optical thickness and transit cross-section of a general dust tail can affect the observed wavelength dependence and depth of transit light curves. Methods. We developed a new 3D model that ejects sublimating particles from the planet surface to build up a dust tail, assuming it to be optically thin, and used 3D radiative transfer computations that fully treat scattering using the distribution of hollow spheres (DHS) method, to generate transit light curves between 0.45 and 2.5 $\mu$m. Results. We show that the transit depth is wavelength independent for optically thick tails, potentially explaining why only some observations indicate a wavelength dependence. From the 3D nature of our simulated tails, we show that their transit cross-sections are related to the component of particle ejection velocity perpendicular to the planet's orbital plane and use this to derive a minimum ejection velocity of 1.2 kms$^{-1}$. To fit the average transit depth of Kepler-1520 b of 0.87%, we require a high dust mas-loss rate of 7 $-$ 80 M$_\oplus$ Gyr$^{-1}$ which implies planet lifetimes that may be inconsistent with the observed sample. Therefore, these mass-loss rates should be considered to be upper limits.

Context: Post-Asymptotic Giant Branch (AGB) binaries are surrounded by stable dusty and gaseous disks similar to the ones around young stellar objects. Whereas significant effort is spent on modeling observations of disks around young stellar objects, the disks around post-AGB binaries receive significantly less attention, even though they pose significant constraints on theories of disk physics and binary evolution. Aims: We want to examine the structure of and phenomena at play in circumbinary disks around post-AGB stars. We continue the analysis of our near-infrared interferometric image of the inner rim of the circumbinary disk around IRAS08544-4431. We want to understand the physics governing this inner disk rim. Methods: We use a radiative transfer model of a dusty disk to reproduce simultaneously the photometry as well as the near-infrared interferometric dataset on IRAS08544-4431. The model assumes hydrostatic equilibrium and takes dust settling self-consistently into account. Results: The best-fit radiative transfer model shows excellent agreement with the spectral energy distribution up to mm wavelengths as well as with the PIONIER visibility data. It requires a rounded inner rim structure, starting at a radius of 8.25 au. However, the model does not fully reproduce the detected over-resolved flux nor the azimuthal flux distribution of the inner rim. While the asymmetric inner disk rim structure is likely to be the consequence of disk-binary interactions, the origin of the additional over-resolved flux remains unclear. Conclusions: As in young stellar objects, the disk inner rim of IRAS08544-4431 is ruled by dust sublimation physics. Additional observations are needed to understand the origin of the extended flux and the azimuthal perturbation at the inner rim of the disk.

The study of extrasolar planets and of the Solar System provides complementary pieces of the mosaic represented by the process of planetary formation. Exoplanets are essential to fully grasp the huge diversity of outcomes that planetary formation and the subsequent evolution of the planetary systems can produce. The orbital and basic physical data we currently possess for the bulk of the exoplanetary population, however, do not provide enough information to break the intrinsic degeneracy of their histories, as different evolutionary tracks can result in the same final configurations. The lessons learned from the Solar System indicate us that the solution to this problem lies in the information contained in the composition of planets. The goal of the Atmospheric Remote-Sensing Infrared Exoplanet Large-survey (ARIEL), one of the three candidates as ESA M4 space mission, is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. ARIEL will focus on warm and hot planets to take advantage of their well-mixed atmospheres, which should show minimal condensation and sequestration of high-Z materials and thus reveal their bulk composition across all main cosmochemical elements. In this work we will review the most outstanding open questions concerning the way planets form and the mechanisms that contribute to create habitable environments that the compositional information gathered by ARIEL will allow to tackle

This paper investigates how the far-IR water ice features can be used to infer properties of disks around T Tauri stars and the water ice thermal history. We explore the power of future observations with SOFIA/HIRMES and SPICA's proposed far-IR instrument SAFARI. A series of detailed radiative transfer disk models around a representative T Tauri star are used to investigate how the far-IR water ice features at 45 and 63 micron change with key disk properties: disk size, grain sizes, disk dust mass, dust settling, and ice thickness. In addition, a series of models is devised to calculate the water ice emission features from warmup, direct deposit and cooldown scenarios of the water ice in disks. Photodesorption from icy grains in disk surfaces weakens the mid-IR water ice features by factors 4-5. The far-IR water ice emission features originate from small grains at the surface snow line in disks at distance of 10-100 au. Unless this reservoir is missing in disks (e.g. transitional disks with large cavities), the feature strength is not changing. Grains larger than 10 micron do not contribute to the features. Grain settling (using turbulent description) is affecting the strength of the ice features by at most 15%. The strength of the ice feature scales with the disk dust mass and water ice fraction on the grains, but saturates for dust masses larger than 1.e-4 Msun and for ice mantles that increase the dust mass by more than 50%. The various thermal histories of water ice leave an imprint on the shape of the features (crystalline/amorphous) as well as on the peak strength and position of the 45 micron feature. SOFIA/HIRMES can only detect crystalline ice features much stronger than simulated in our standard T Tauri disk model in deep exposures (1 hr). SPICA/SAFARI can detect the typical ice features in our standard T Tauri disk model in short exposures (10 min). (abbreviated)

Context. Multi-wavelength observations are indispensable in studying disk geometry and dust evolution processes in protoplanetary disks. Aims. We aimed to construct a 3-dimensional model of HD 163296 capable of reproducing simultaneously new observations of the disk surface in scattered light with the SPHERE instrument and thermal emission continuum observations of the disk midplane with ALMA. We want to determine why the SED of HD 163296 is intermediary between the otherwise well-separated group I and group II Herbig stars. Methods. The disk was modelled using the Monte Carlo radiative transfer code MCMax3D. The radial dust surface density profile was modelled after the ALMA observations, while the polarized scattered light observations were used to constrain the inclination of the inner disk component and turbulence and grain growth in the outer disk. Results. While three rings are observed in the disk midplane in millimeter thermal emission at $\sim$80, 124 and 200 AU, only the innermost of these is observed in polarized scattered light, indicating a lack of small dust grains on the surface of the outer disk. We provide two models capable of explaining this difference. The first model uses increased settling in the outer disk as a mechanism to bring the small dust grains on the surface of the disk closer to the midplane, and into the shadow cast by the first ring. The second model uses depletion of the smallest dust grains in the outer disk as a mechanism for decreasing the optical depth at optical and NIR wavelengths. In the region outside the fragmentation-dominated regime, such depletion is expected from state-of-the-art dust evolution models. We studied the effect of creating an artificial inner cavity in our models, and conclude that HD 163296 might be a precursor to typical group I sources.

Meteorites contain minerals from Solar System asteroids with different properties (like size, presence of water, core formation). We provide new mid-IR transmission spectra of powdered meteorites to obtain templates of how mid-IR spectra of asteroidal debris would look like. This is essential for interpreting mid-IR spectra of past and future space observatories, like the James Webb Space Telescope. We show that the transmission spectra of wet and dry chondrites, carbonaceous and ordinary chondrites and achondrite and chondrite meteorites are distinctly different in a way one can distinguish in astronomical mid-IR spectra. The two observables that spectroscopically separate the different meteorites groups (and thus the different types of parent bodies) are the pyroxene-olivine feature strength ratio and the peak shift of the olivine spectral features due to an increase in the iron concentration of the olivine.

The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure.

We present multi-epoch VLT/SPHERE observations of the protoplanetary disk around HD 135344B (SAO 206462). The $J$-band scattered light imagery reveal, with high spatial resolution ($\sim$41 mas, 6.4 au), the disk surface beyond $\sim$20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution in the inner disk. These shadows manifest themselves as narrow lanes, cast by localized density enhancements, and broader features which possibly trace the larger scale dynamics of the inner disk. We acquired visible and near-infrared photometry which shows variations up to 10% in the $JHK$ bands, possibly correlated with the presence of the shadows. Analysis of archival VLTI/PIONIER $H$-band visibilities constrain the orientation of the inner disk to $i = 18.2\deg^{+3.4}_{-4.1}$ and ${\rm PA} = 57.3\deg \pm 5.7\deg$, consistent with an alignment with the outer disk or a minor disk warp of several degrees. The latter scenario could explain the broad, quasi-stationary shadowing in N-NW direction in case the inclination of the outer disk is slightly larger. The correlation between the shadowing and the near-infrared excess is quantified with a grid of radiative transfer models. The variability of the scattered light contrast requires extended variations in the inner disk atmosphere ($H/r \lesssim 0.2$). Possible mechanisms that may cause asymmetric variations in the optical depth ($\Delta\tau\lesssim1$) through the atmosphere of the inner disk include turbulent fluctuations, planetesimal collisions, or a dusty disk wind, possibly enhanced by a minor disk warp. A fine temporal sampling is required to follow day-to-day changes of the shadow patterns which may be a face-on variant of the UX Orionis phenomenon.

Opacities of molecules in exoplanet atmospheres rely on increasingly detailed line-lists for these molecules. The line lists available today contain for many species up to several billions of lines. Computation of the spectral line profile created by pressure and temperature broadening, the Voigt profile, of all of these lines is becoming a computational challenge. We aim to create a method to compute the Voigt profile in a way that automatically focusses the computation time into the strongest lines, while still maintaining the continuum contribution of the high number of weaker lines. Here, we outline a statistical line sampling technique that samples the Voigt profile quickly and with high accuracy. The number of samples is adjusted to the strength of the line and the local spectral line density. This automatically provides high accuracy line shapes for strong lines or lines that are spectrally isolated. The line sampling technique automatically preserves the integrated line opacity for all lines, thereby also providing the continuum opacity created by the large number of weak lines at very low computational cost. The line sampling technique is tested for accuracy when computing line spectra and correlated-k tables. Extremely fast computations (~3.5e5 lines per second per core on a standard current day desktop computer) with high accuracy (< 1 % almost everywhere) are obtained. A detailed recipe on how to perform the computations is given.

Shadows in transitional disks are generally interpreted as signs of a misaligned inner disk. This disk is usually beyond the reach of current day high contrast imaging facilities. However, the location and morphology of the shadow features allow us to reconstruct the inner disk geometry. We derive analytic equations of the locations of the shadow features as a function of the orientation of the inner and outer disk and the height of the outer disk wall. In contrast to previous claims in the literature, we show that the position angle of the line connecting the shadows cannot be directly related to the position angle of the inner disk. We show how the analytic framework derived here can be applied to transitional disks with shadow features. We use estimates of the outer disk height to put constraints on the inner disk orientation. In contrast with the results from Long et al. (2017), we derive that for the disk surrounding HD100453 the analytic estimates and interferometric observations result in a consistent picture of the orientation of the inner disk. The elegant consistency in our analytic framework between observation and theory strongly support both the interpretation of the shadow features as coming from a misaligned inner disk as well as the diagnostic value of near infrared interferometry for inner disk geometry.

Exoplanet Proxima b will be an important laboratory for the search for extraterrestrial life for the decades ahead. Here we discuss the prospects of detecting carbon dioxide at 15 um using a spectral filtering technique with the Medium Resolution Spectrograph (MRS) mode of the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST). At superior conjunction, the planet is expected to show a contrast of up to 100 ppm with respect to the star. At a spectral resolving power of R=1790-2640, about 100 spectral CO2 features are visible within the 13.2-15.8 um (3B) band, which can be combined to boost the planet atmospheric signal by a factor 3-4, depending on the atmospheric temperature structure and CO2 abundance. If atmospheric conditions are favorable (assuming an Earth-like atmosphere), with this new application to the cross-correlation technique carbon dioxide can be detected within a few days of JWST observations. However, this can only be achieved if both the instrumental spectral response and the stellar spectrum can be determined to a relative precision of <= 1e-4 between adjacent spectral channels. Absolute flux calibration is not required, and the method is insensitive to the strong broadband variability of the host star. Precise calibration of the spectral features of the host star may only be attainable by obtaining deep observations of the system during inferior conjunction that serve as a reference. The high-pass filter spectroscopic technique with the MIRI MRS can be tested on warm Jupiters, Neptunes, and super-Earths with significantly higher planet/star contrast ratios than the Proxima system.

We present experimental phase functions of three types of millimeter-sized dust grains consisting of enstatite, quartz and volcanic material from Mount Etna, respectively. The three grains present similar sizes but different absorbing properties. The measurements are performed at 527 nm covering the scattering angle range from 3 to 170 degrees. The measured phase functions show two well defined regions i) soft forward peaks and ii) a continuous increase with the scattering angle at side- and back-scattering regions. This behavior at side- and back-scattering regions are in agreement with the observed phase functions for the Fomalhaut and HR 4796A dust rings. Further computations and measurements (including polarization) for millimeter sized-grains are needed to draw some conclusions about the fluffy or compact structure of the dust grains.

Temperature changes in the planet forming disc midplanes carry important physico-chemical consequences, such as the effect on the locations of the condensation fronts of molecules - the snowlines. Snowlines impose major chemical gradients and possibly foster grain growth. The aim of this paper is to understand how disc midplane temperature changes with gas and dust evolution, and identify trends that may influence planet formation or allow to constrain disc evolution observationally. We calculate disc temperature, hydrostatic equilibrium and dust settling in a mutually consistent way from a grid of disc models at different stages of gas loss, grain growth and hole opening. We find that the CO snowline location depends very strongly on disc properties. The CO snowline location migrates closer to the star for increasing degrees of gas dispersal and dust growth. Around a typical A type star, the snowline can be anywhere between several tens and a few hundred au, depending on the disc properties such as gas mass and grain size. In fact, gas loss is as efficient as dust evolution in settling discs, and flat discs may be gas-poor counterparts of flared discs. Our results, in the context of different pre-main sequence evolution of the luminosity in low- and intermediate-mass stars suggests very different thermal (and hence chemical) histories in these two types of discs. Discs of T Tauri stars settle and cool down while discs of Herbig Ae stars may remain rather warm throughout the pre-main sequence.

HR4796A is surrounded by a debris disc, observed in scattered light as an inclined ring. Past observations raised several questions. First, a strong brightness asymmetry detected in polarized reflected light recently challenged our understanding of scattering by the dust particles in this system. Secondly, the morphology of the ring strongly suggests the presence of planets, although no planets have been detected to date. We obtained high-angular resolution coronagraphic images of the circumstellar environment around HR4796A with VLT/SPHERE during the commissioning of the instrument in May 2014 and during guaranteed-time observations in February 2015. The observations reveal for the first time the entire ring of dust, including the semi-minor axis that was previously hidden either behind the coronagraphic spot or in the speckle noise. We determine empirically the scattering phase function of the dust in the H band from 13.6deg to 166.6deg. It shows a prominent peak of forward scattering, never detected before, for scattering angles below 30deg. We analyse the reflectance spectra of the disc from the 0.95 to 1.6 microns, confirming the red colour of the dust, and derive detection limits on the presence of planetary mass objects. We confirm which side of the disc is inclined towards the Earth. The analysis of the phase function suggests that the dust population is dominated by particles much larger than the observation wavelength, of about 20 microns. Compact Mie grains of this size are incompatible with the spectral energy distribution of the disc, however the observed rise in scattering efficiency beyond 50deg points towards aggregates which could reconcile both observables. We do not detect companions orbiting the star but our high-contrast observations provide the most stringent constraints yet on the presence of planets responsible for the morphology of the dust.

To understand the chemical composition of planets, it is important to know the chemical composition of the region where they form in protoplanetary disks. Due to its fundamental role in chemical and biological processes, carbon is a key element to trace. We aim to identify the carriers and processes behind the extended NIR flux observed around several Herbig stars. We compare the extended NIR flux from objects in the PIONIER Herbig Ae/Be survey with their flux in the PAH features. HD 100453 is used as a benchmark case to investigate the influence of quantum heated particles, like PAHs or very small carbonaceous grains, in more detail. We use the Monte Carlo radiative transfer code MCMax to do a parameter study of the QHP size and scale- height and examine the influence of quantum heating on the amount of extended flux in the NIR visibilities. There is a correlation between the PAH feature flux of a disk and the amount of its extended NIR flux. We find that very small carbonaceous grains create the observed extended NIR flux around HD 100453 and still lead to a realistic SED. These results can not be achieved without using quantum heating effects, e.g. only with scattered light and grains in thermal equilibrium. It is possible to explain the extended NIR emission around Herbig stars with the presence of carbonaceous, quantum heated particles. Interferometric observations can be used to constrain the spatial distribution and typical size of carbonaceous material in the terrestrial planet forming region.

The James Webb Space Telescope will enable astronomers to obtain exoplanet spectra of unprecedented precision. Especially the MIRI instrument may shed light on the nature of the cloud particles obscuring planetary transmission spectra in the optical and near-infrared. We provide self-consistent atmospheric models and synthetic JWST observations for prime exoplanet targets in order to identify spectral regions of interest and estimate the number of transits needed to distinguish between model setups. We select targets which span a wide range in planetary temperature and surface gravity, ranging from super-Earths to giant planets, and have a high expected SNR. For all targets we vary the enrichment, C/O ratio, presence of optical absorbers (TiO/VO) and cloud treatment. We calculate atmospheric structures and emission and transmission spectra for all targets and use a radiometric model to obtain simulated observations. We analyze JWST's ability to distinguish between various scenarios. We find that in very cloudy planets such as GJ 1214b less than 10 transits with NIRSpec may be enough to reveal molecular features. Further, the presence of small silicate grains in atmospheres of hot Jupiters may be detectable with a single JWST MIRI transit. For a more detailed characterization of such particles less than 10 transits are necessary. Finally, we find that some of the hottest hot Jupiters are well fitted by models which neglect the redistribution of the insolation and harbor inversions, and that 1-4 eclipse measurements with NIRSpec are needed to distinguish between the inversion models. Wet thus demonstrate the capabilities of JWST for solving some of the most intriguing puzzles in current exoplanet atmospheric research. Further, by publishing all models calculated for this study we enable the community to carry out similar or retrieval analyses for all planets included in our target list.

The protoplanetary disk around the F-type star HD 135344B (SAO 206462) is in a transition stage and shows many intriguing structures both in scattered light and thermal (sub-)millimeter emission which are possibly related to planet formation processes. We study the morphology and surface brightness of the disk in scattered light to gain insight into the innermost disk regions, the formation of protoplanets, planet-disk interactions traced in the surface and midplane layers, and the dust grain properties of the disk surface. We have carried out high-contrast polarimetric differential imaging (PDI) observations with VLT/SPHERE and obtained polarized scattered light images with ZIMPOL in R- and I-band and with IRDIS in Y- and J-band. The scattered light images reveal with unprecedented angular resolution and sensitivity the spiral arms as well as the 25 au cavity of the disk. Multiple shadow features are discovered on the outer disk with one shadow only being present during the second observation epoch. A positive surface brightness gradient is observed in the stellar irradiation corrected images in southwest direction possibly due to an azimuthally asymmetric perturbation of the temperature and/or surface density by the passing spiral arms. The disk integrated polarized flux, normalized to the stellar flux, shows a positive trend towards longer wavelengths which we attribute to large aggregate dust grains in the disk surface. Part of the the non-azimuthal polarization signal in the Uphi image of the J-band observation could be the result of multiple scattering in the disk. The detected shadow features and their possible variability have the potential to provide insight into the structure of and processes occurring in the innermost disk regions.

Understanding the diversity of planets requires to study the morphology and the physical conditions in the protoplanetary disks in which they form. We observed and spatially resolved the disk around the ~10 Myr old protoplanetary disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and near-infrared wavelengths, reaching an angular resolution of ~0.02", and an inner working angle of ~0.09". We detect polarized scattered light up to ~0.42" (~48 au) and detect a cavity, a rim with azimuthal brightness variations at an inclination of 38 degrees, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint spiral-like feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of ~119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by ~72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim, induces an azimuthal variation of the scale height that can contribute to the brightness variations. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. The origin of such a misalignment in HD 100453, and of the spirals, is unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit.

We present a mid-IR interferometric survey of the circumstellar environment of a specific class of post-Asymptotic Giant Branch (post-AGB) binaries. For this class the presence of a compact dusty disk has been postulated on the basis of various spatially unresolved measurements. Our interferometric survey was performed with the MIDI instrument on the VLTI. In total 19 different systems were observed using variable baseline configurations. Combining all the visibilities at a single wavelength at 10.7 micron, we fitted two parametric models to the data: a uniform disk (UD) and a ring model mimicking a temperature gradient. We compared our observables of the whole sample, with synthetic data computed from a grid of radiative transfer models of passively irradiated disks in hydrostatic equilibrium. These models are computed with a Monte Carlo code that has been widely applied to describe the structure of protoplanetary disks around young stellar objects (YSO). The spatially resolved observations show that the majority of our targets cluster closely together in the distance-independent size-colour diagram, and have extremely compact N-band emission regions. The typical uniform disk diameter of the N-band emission region is about 40 mass which corresponds to a typical brightness temperature of 400-600~K. The resolved objects display very similar characteristics in the interferometric observables and in the spectral energy distributions. Therefore, the physical properties of the disks around our targets must be similar. The grid of protoplanetary disk models covers very well the observed objects. Much like for young stars, the spatially resolved N-band emission region is determined by the hot inner rim of the disk. Continued comparisons between post-AGB and protoplanetary disks will help to understand grain growth and disk evolution processes,

We present scattered light images of the TW Hya disk performed with SPHERE in PDI mode at 0.63, 0.79, 1.24 and 1.62 micron. We also present H2/H3-band ADI observations. Three distinct radial depressions in the polarized intensity distribution are seen, around 85, 21, and 6~au. The overall intensity distribution has a high degree of azimuthal symmetry; the disk is somewhat brighter than average towards the South and darker towards the North-West. The ADI observations yielded no signifiant detection of point sources in the disk. Our observations have a linear spatial resolution of 1 to 2au, similar to that of recent ALMA dust continuum observations. The sub-micron sized dust grains that dominate the light scattering in the disk surface are strongly coupled to the gas. We created a radiative transfer disk model with self-consistent temperature and vertical structure iteration and including grain size-dependent dust settling. This method may provide independent constraints on the gas distribution at higher spatial resolution than is feasible with ALMA gas line observations. We find that the gas surface density in the "gaps" is reduced by 50% to 80% relative to an unperturbed model. Should embedded planets be responsible for carving the gaps then their masses are at most a few 10 Mearth. The observed gaps are wider, with shallower flanks, than expected for planet-disk interaction with such low-mass planets. If forming planetary bodies have undergone collapse and are in the "detachted phase" then they may be directly observable with future facilities such as METIS at the E-ELT.

Whether BP Piscium (BP Psc) is either a pre-main sequence T Tauri star at d ~ 80 pc, or a post-main sequence G giant at d ~ 300 pc is still not clear. As a first-ascent giant, it is the first to be observed with a molecular and dust disk. Alternatively, BP Psc would be among the nearest T Tauri stars with a protoplanetary disk (PPD). We investigate whether the disk geometry resembles typical PPDs, by comparing polarimetric images with radiative transfer models. Our VLT/SPHERE/ZIMPOL observations allow us to perform Polarimetric Differential Imaging; Reference Star Differential Imaging; and Richardson-Lucy deconvolution. We present the first visible light polarization and intensity images of the disk of BP Psc. Our deconvolution confirms the disk shape as detected before, mainly showing the southern side of the disk. In polarized intensity the disk is imaged at larger detail and also shows the northern side, giving it the typical shape of high inclination flared disks. We explain the observed disk features by retrieving the large-scale geometry with MCMax radiative transfer modeling, which yields a strongly flared model, atypical for disks of T Tauri stars.

High-contrast scattered light observations have revealed the surface morphology of several dozens of protoplanetary disks at optical and near-infrared wavelengths. Inclined disks offer the opportunity to measure part of the phase function of the dust grains that reside in the disk surface which is essential for our understanding of protoplanetary dust properties and the early stages of planet formation. We aim to construct a method which takes into account how the flaring shape of the scattering surface of an (optically thick) protoplanetary disk projects onto the image plane of the observer. This allows us to map physical quantities (scattering radius and scattering angle) onto scattered light images and retrieve stellar irradiation corrected (r^2-scaled) images and dust phase functions. We apply the method on archival polarized intensity images of the protoplanetary disk around HD 100546 that were obtained with VLT/SPHERE in R'-band and VLT/NACO in H- and Ks-band. The brightest side of the r^2-scaled R'-band polarized intensity image of HD 100546 changes from the far to the near side of the disk when a flaring instead of a geometrically flat disk surface is used for the r^2-scaling. The decrease in polarized surface brightness in the scattering angle range of ~40-70 deg is likely a result of the dust phase function and degree of polarization which peak in different scattering angle regimes. The derived phase functions show part of a forward scattering peak which indicates that large, aggregate dust grains dominate the scattering opacity in the disk surface. Projection effects of a protoplanetary disk surface need to be taken into account to correctly interpret scattered light images. Applying the correct scaling for the correction of stellar irradiation is crucial for the interpretation of the images and the derivation of the dust properties in the disk surface layer.

We studied the well known circumstellar disk around the Herbig Ae/Be star HD 97048 with high angular resolution to reveal undetected structures in the disk, which may be indicative of disk evolutionary processes such as planet formation. We used the IRDIS near-IR subsystem of the extreme adaptive optics imager SPHERE at the ESO/VLT to study the scattered light from the circumstellar disk via high resolution polarimetry and angular differential imaging. We imaged the disk in unprecedented detail and revealed four ring-like brightness enhancements and corresponding gaps in the scattered light from the disk surface with radii between 39 au and 341 au. We derived the inclination and position angle as well as the height of the scattering surface of the disk from our observational data. We found that the surface height profile can be described by a single power law up to a separation ~270 au. Using the surface height profile we measured the scattering phase function of the disk and found that it is well consistent with theoretical models of compact dust aggregates. We discuss the origin of the detected features and find that low mass (< 1 M_Jup) nascent planets are a possible explanation.

Evaporating rocky exoplanets, such as KIC 12557548b, eject large amounts of dust grains, which can trail the planet in a comet-like tail. When such objects occult their host star, the resulting transit signal contains information about the dust in the tail. We aim to use the detailed shape of the Kepler light curve of KIC 12557548b to constrain the size and composition of the dust grains that make up the tail, as well as the mass loss rate of the planet. Using a self-consistent numerical model of the dust dynamics and sublimation, we calculate the shape of the tail by following dust grains from their ejection from the planet to their destruction due to sublimation. From this dust cloud shape, we generate synthetic light curves (incorporating the effects of extinction and angle-dependent scattering), which are then compared with the phase-folded Kepler light curve. We explore the free-parameter space thoroughly using a Markov chain Monte Carlo method. Our physics-based model is capable of reproducing the observed light curve in detail. Good fits are found for initial grain sizes between 0.2 and 5.6 micron and dust mass loss rates of 0.6 to 15.6 M_earth/Gyr (2-sigma ranges). We find that only certain combinations of material parameters yield the correct tail length. These constraints are consistent with dust made of corundum (Al2O3), but do not agree with a range of carbonaceous, silicate, or iron compositions. Using a detailed, physically motivated model, it is possible to constrain the composition of the dust in the tails of evaporating rocky exoplanets. This provides a unique opportunity to probe to interior composition of the smallest known exoplanets.

The mass of gas in protoplanetary discs is a quantity of great interest for assessing their planet formation potential. Disc gas masses are, however, traditionally inferred from measured dust masses by applying an assumed standard gas-to-dust ratio of $g/d=100$. Furthermore, measuring gas masses based on CO observations has been hindered by the effects of CO freeze-out. Here we present a novel approach to study the mid-plane gas by combining C$^{18}$O line modelling, CO snowline observations and the spectral energy distribution (SED) and selectively study the inner tens of au where freeze-out is not relevant. We apply the modelling technique to the disc around the Herbig Ae star HD 163296 with particular focus on the regions within the CO snowline radius, measured to be at 90 au in this disc. Our models yield the mass of C$^{18}$O in this inner disc region of $M_{\text{C}^{18}\text{O}}(<90\,\text{au})\sim 2\times10^{-8}$ M$_\odot$. We find that most of our models yield a notably low $g/d<20$, especially in the disc mid-plane ($g/d<1$). Our only models with a more interstellar medium (ISM)-like $g/d$ require C$^{18}$O to be underabundant with respect to the ISM abundances and a significant depletion of sub-micron grains, which is not supported by scattered light observations. Our technique can be applied to a range of discs and opens up a possibility of measuring gas and dust masses in discs within the CO snowline location without making assumptions about the gas-to-dust ratio.

The presence or absence of ice in protoplanetary disks is of great importance for the formation of planets. By enhancing the solid surface density and increasing the sticking efficiency, ice catalyzes the rapid formation of planetesimals and decreases the time scale for giant planet core accretion. Aims: In this paper we analyse the composition of the outer disk around the Herbig star HD~142527. We focus on the composition of the water ice, but also analyse the abundances of previously proposed minerals. Methods: We present new Herschel far infrared spectra and a re-reduction of archival data from the Infrared Space Observatory (ISO). We model the disk using full 3D radiative transfer to obtain the disk structure. Also, we use an optically thin analysis of the outer disk spectrum to obtain firm constraints on the composition of the dust component. Results: The water ice in the disk around HD~142527 contains a large reservoir of crystalline water ice. We determine the local abundance of water ice in the outer disk (i.e. beyond 130\,AU). The re-reduced ISO spectrum differs significantly from that previously published, but matches the new Herschel spectrum at their common wavelength range. In particular, we do not detect any significant contribution from carbonates or hydrous silicates, in contrast to earlier claims. Conclusions: The amount of water ice detected in the outer disk requires $\sim80\,$\% of the oxygen atoms. This is comparable to the water ice abundance in the outer solar system, in comets and in dense interstellar clouds. The water ice is highly crystalline while the temperatures where we detect it are too low to crystallize the water on relevant time scales. We discuss the implications of this finding.

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$\alpha$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$\alpha$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$\alpha$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$\alpha$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $\rho(r)=\rho_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.

We investigate the effect that the turbulent mixing strength parameter $\alpha_{\rm{turb}}$ plays on near-infrared polarimetric and sub-millimetre interferometric imaging observations of transitional discs (TDs) with a gap carved by a planet. We generate synthetic observations of these objects with ALMA and VLT/SPHERE-ZIMPOL by combining hydrodynamical, dust evolution, radiative transfer and instrument models for values of $\alpha_{\rm{turb}}=[10^{-4}, 10^{-3}, 10^{-2}]$. We find that, through a combination of effects on the viscosity of the gas, the turbulent mixing and dust evolution processes, $\alpha_{\rm{turb}}$ strongly affects the morphology of the dust distribution that can be traced with these observations. We constrain the value of $\alpha_{\rm{turb}}$ to be within an order of magnitude of $10^{-3}$ in TD sources that show cavities in sub-mm continuum images while featuring continuous distribution of dust or smaller cavities in NIR-polarimetric images.

Context. Extinction and emission of dust models need for observational constraints to be validated. The coreshine phenomenon has already shown the importance of scattering in the 3 to 5 micron range and its ability to validate dust properties for dense cores. Aims. We want to investigate whether scattering can also play a role at longer wavelengths and to place even tighter constraints on the dust properties. Methods. We analyze the inversion of the Spitzer 8 micron map of the dense molecular cloud L183, to examine the importance of scattering as a potential contributor to the line-of-sight extinction. Results. The column density deduced from the inversion of the 8 micron map, when we neglect scattering, disagrees with all the other column density measurements of the same region. Modeling confirms that scattering at 8 microns is not negligible with an intensity of several hundred kJy per sr. This demonstrates the need of efficiently scattering dust grains at MIR wavelengths up to 8 microns. Coagulated aggregates are good candidates and might also explain the discrepancy at high extinction between E(J-K) et tau(9.7) toward dense molecular clouds. Further investigation requires considering efficiently scattering dust grains including ices as realistic dust models.

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

In protoplanetary disks micron-size dust grains coagulate to form larger structures with complex shapes and compositions. The coagulation process changes the absorption and scattering properties of particles in the disk in significant ways. To properly interpret observations of protoplanetary disks and to place these observations in the context of the first steps of planet formation, it is crucial to understand the optical properties of these complex structures. We derive the optical properties of dust aggregates using detailed computations of aggregate structures and compare these computa- tionally demanding results with approximate methods that are cheaper to compute in practice. In this way we wish to understand the merits and problems of approximate methods and define the context in which they can or cannot be used to analyze observations of objects where significant grain growth is taking place. For the detailed computations we used the discrete dipole approximation (DDA), a method able to compute the interaction of light with a complexly shaped, inhomogeneous particle. We compared the results to those obtained using spherical and irregular, homogeneous and inhomogeneous particles. While no approximate method properly reproduces all characteristics of large dust aggregates, the thermal properties of dust can be analyzed using irregularly shaped, porous, inhomogeneous grains. The asymmetry of the scattering phase function is a good indicator of aggregate size, while the degree of polarization is probably determined by the size of the constituent particles. Optical properties derived from aggregates significantly differ from the most frequently used standard ("astronomical silicate" in spherical grains). We outline a computationally fast and relatively accurate method that can be used for a multiwavelength analysis of aggregate dust in protoplanetary disks.

We aim to constrain the structure of the circumstellar material around the post-AGB binary and RV Tauri pulsator AC Her. We want to constrain the spatial distribution of the amorphous as well as of the crystalline dust. We present very high-quality mid-IR interferometric data that were obtained with MIDI/VLTI. We analyse the MIDI data and the full SED, using the MCMax radiative transfer code, to find a good structure model of AC Her's circumbinary disk. We include a grain size distribution and midplane settling of dust self-consistently. The spatial distribution of crystalline forsterite in the disk is investigated with the mid-IR features, the 69~$\mu$m band and the 11.3~$\mu$m signatures in the interferometric data. All the data are well fitted. The inclination and position angle of the disk are well determined at i=50+-8 and PA=305+-10. We firmly establish that the inner disk radius is about an order of magnitude larger than the dust sublimation radius. Significant grain growth has occurred, with mm-sized grains being settled to the midplane of the disk. A large dust mass is needed to fit the sub-mm fluxes. By assuming {\alpha}=0.01, a good fit is obtained with a small grain size power law index of 3.25, combined with a small gas/dust ratio <10. The resulting gas mass is compatible with recent estimates employing direct gas diagnostics. The spatial distribution of the forsterite is different from the amorphous dust, as more warm forsterite is needed in the surface layers of the inner disk. The disk in AC Her is very evolved, with its small gas/dust ratio and large inner hole. Mid-IR interferometry offers unique constraints, complementary to mid-IR features, for studying the mineralogy in disks. A better uv coverage is needed to constrain in detail the distribution of the crystalline forsterite in AC Her, but we find strong similarities with the protoplanetary disk HD100546.

Context: AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H2O. Aims: The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8+0.0, as traced by Herschel observations, is studied, with a focus on the H2O content and the dust-to-gas ratio. We report detecting a large number of H2O vapor emission lines up to J = 9 in the Herschel data, for which we present the measured line strengths. Methods: The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is illustrated for both low and high mass-loss-rate sources. Specifically, we discuss different ways of assessing the dust-to-gas ratio: 1) from the dust thermal emission spectrum and the CO molecular gas line strengths; 2) from the momentum transfer from dust to gas and the measured gas terminal velocity; and 3) from the determination of the required amount of dust to reproduce H2O lines for a given H2O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation of a possible radial dependence of the dust-to-gas ratio. Results: ... Continued in article.

The water ice or snow line is one of the key properties of protoplanetary disks that determines the water content of terrestrial planets in the habitable zone. Its location is determined by the properties of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in the disk. We calculate the snow line location from recent observations of mass accretion rates and as a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 0.4 dex naturally leads to a dispersion in snow line locations of 0.2 dex. For ISM-like dust sizes, the one-sigma snow line location among solar mass stars of the same age ranges from 2 to 5 au. For more realistic dust opacities that include larger grains, the snow line is located up to two times closer to the star. We use these locations and the outcome of N-body simulations to predict the amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find that the dispersion in snow line locations leads to a large range in water content. For ISM-like dust sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and no water is delivered to the habitable zones of low-mass M stars (less than half a solar mass) as in previous works. The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a significant fraction of M stars and all FGK stars. Considering their larger numbers and higher planet occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets are able to hold on to their water in their subsequent evolution.

Low- and intermediate-mass stars go through a period of intense mass-loss at the end of their lives in a phase known as the asymptotic giant branch (AGB). During the AGB a significant fraction of their initial mass is expelled in a stellar wind. This process controls the final stages of their evolution and contributes to the chemical evolution of galaxies. However, the wind-driving mechanism of AGB stars is not yet well understood, especially so for oxygen-rich sources. Characterizing both the present-day mass-loss and wind structure and the evolution of the mass-loss rate of such stars is paramount to advancing our understanding of this processes. We modelled the dust envelope of W Hya using an advanced radiative transfer code. The dust model was analysed in the light of a previously calculated gas-phase wind model and compared to measurements available in the literature, such as infrared spectra, infrared images, and optical scattered light fractions. We find that the dust spectrum of W Hya can partly be explained by a gravitationally bound dust shell that probably is responsible for most of the amorphous Al$_2$O$_3$ emission. The composition of the large ($\sim$\,0.3\,$\mu$m) grains needed to explain the scattered light cannot be constrained, but probably is dominated by silicates. Silicate emission in the thermal infrared was found to originate from beyond 40 AU from the star and we find that they need to have substantial near-infrared opacities to be visible at such large distances. The increase in near-infrared opacity of the dust at these distances roughly coincides with a sudden increase in expansion velocity as deduced from the gas-phase CO lines. Finally, the recent mass loss of W Hya is confirmed to be highly variable and we identify a strong peak in the mass-loss rate that occurred about 3500 years ago and lasted for a few hundred years.

We study the grain properties and location of the forsterite crystals in the circumstellar environment of the pre-planetary nebula (PPN) IRAS 17150-3224 in order to learn more about the as yet poorly understood evolutionary phase prior to the PPN. We use the best-fit model for IRAS 17150-3224 of Meixner et al. (2002) and add forsterite to this model. We investigate different spatial distributions and grain sizes of the forsterite crystals in the circumstellar environment. We compare the spectral bands of forsterite in the mid-infrared and at 69 micrometre in radiative transport models to those in ISO-SWS and Herschel/PACS observations. We can reproduce the non-detection of the mid-infrared bands and the detection of the 69 micrometre feature with models where the forsterite is distributed in the whole outflow, in the superwind region, or in the AGB-wind region emitted previous to the superwind, but we cannot discriminate between these three models. To reproduce the observed spectral bands with these three models, the forsterite crystals need to be dominated by a grain size population of 2 micrometre up to 6 micrometre. We hypothesise that the large forsterite crystals were formed after the superwind phase of IRAS 17150-3224, where the star developed an as yet unknown hyperwind with an extremely high mass-loss rate (10^-3 Msol/yr). The high densities of such a hyperwind could be responsible for the efficient grain growth of both amorphous and crystalline dust in the outflow. Several mechanisms are discussed that might explain the lower-limit of 2 micrometre found for the forsterite grains, but none are satisfactory. Among the mechanisms explored is a possible selection effect due to radiation pressure based on photon scattering on micron-sized grains.

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of particle trapping and grain growth in pressure bumps originating from the disk interaction with a planetary companion. The multiple features observed in some transition disks, such as SR 21, suggest the presence of more than one planet. We study the gas and dust distributions of a disk hosting two massive planets as function of different disk and dust parameters. Observational signatures, such as the spectral energy distribution, sub-millimetre, and polarised images are simulated for the various parameters. We confirm that planets can lead to particle trapping, although for a disk with high viscosity ($\alpha_{\rm{turb}}=10^{-2}$), the planet should be more massive than $5 M_{\rm{Jup}}$ and dust fragmentation should occur with low efficiency ($v_{f}\sim30\rm{m s}^{-1}$). This will lead to a ring-like feature as observed in transition disks in the millimetre. When trapping occurs, we find that a smooth distribution of micron sized grains throughout the disk, sometimes observed in scattered light, can only happen if the combination of planet mass and turbulence is such that small grains are not fully filtered out. A high disk viscosity ($\alpha_{\rm{turb}}=10^{-2}$) ensures a replenishment of the cavity in micron-sized dust, while for lower viscosity ($\alpha_{\rm{turb}}=10^{-3}$), the planet mass is constrained to be less than $5 M_{\rm{Jup}}$. In these cases, the gas distribution is likely to show low-amplitude azimuthal asymmetries caused by disk eccentricity rather than by long-lived vortices.

Binaries with circumbinary disks are commonly found among optically bright post-AGB stars. Although clearly linked to binary interaction processes, the formation, evolution and fate of these disks are still badly understood. Due to their compactness, interferometric techniques are required to resolve them. Here, we discuss our high-quality multiwavelength interferometric data of two prototypical yet very different post-AGB binaries, AC and 89 Herculis, as well as the modeling thereof with radiative transfer models. A detailed account of the data and models of both objects is published in three separate papers elsewhere; here we focus on comparing the modeling results for the two objects. In particular we discuss the successes and limitations of the models which were developed for protoplanetary disks around young stars. We conclude that multiwavelength high-angular-resolution observations and radiative transfer disk models are indispensible to understand these complex interacting objects and their place in the grand scheme of the (binary) evolution of low and intermediate mass stars.

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$\mu$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

Recently, two exoplanet candidates have been discovered, KIC 12557548b and KOI-2700b, whose transit profiles show evidence for a comet-like tail of dust trailing the planet, thought to be fed by the evaporation of the planet's surface. We aim to put constraints on the composition of the dust ejected by these objects from the shape of their transit light curves. We derive a semi-analytical expression for the attenuation of dust cross-section in the tail, incorporating the sublimation of dust grains as well as their drift away from the planet. This expression shows that the length of the tail is highly sensitive to the sublimation properties of the dust material. We compute tail lengths for several possible dust compositions, and compare these to observational estimates of the tail lengths of KIC 12557548b and KOI-2700b, inferred from their light curves. The observed tail lengths are consistent with dust grains composed of corundum (Al2O3) or iron-rich silicate minerals (e.g., fayalite, Fe2SiO4). Pure iron and carbonaceous compositions are disfavoured. In addition, we estimate dust mass loss rates of 1.7 +/- 0.5 M_earth/Gyr for KIC 12557548b, and > 0.007 M_earth/Gyr (1-sigma lower limit) for KOI-2700b.

Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient medium when measured hours to days after the burst. High values of optical polarization detected minutes after burst in GRB 120308A indicate the presence of large-scale ordered magnetic fields originating from the central engine (the power source of the GRB). Theoretical models predict low degrees of linear polarization and negligable circular polarization at late times, when the energy in the original ejecta is quickly transferred to the ambient medium and propagates farther into the medium as a blastwave. Here we report the detection of circularly polarized optical light in the afterglow of GRB 121024A, measured 0.15 days after the burst. We show that the circular polarization is intrinsic to the afterglow and unlikely to be produced by dust scattering or plasma propagation effects. A possible explanation is to invoke anisotropic (rather than the commonly assumed isotropic) electron pitch angle distributions, and we suggest that new models are required to produce the complex microphysics of realistic shocks in relativistic jets.

We present an analysis of mid-and far-infrared (IR) spectrum and spectral energy distribution (SED) of the LBV R71 in the LMC.This work aims to understand the overall contribution of high-mass LBVs to the total dust-mass budget of the interstellar medium (ISM) of the LMC and compare this with the contribution from low-mass asymptotic giant branch (AGB) stars. As a case study, we analyze the SED of R71. We compiled all the available photometric and spectroscopic observational fluxes from various telescopes for a wide wavelength range (0.36 -- 250\,$\mu$m). We determined the dust composition from the spectroscopic data, and derived the ejected dust mass, dust mass-loss rate, and other dust shell properties by modeling the SED of R71. We noted nine spectral features in the dust shell of R71 by analyzing Spitzer spectroscopic data. Among these, we identified three new crystalline silicate features. We computed our model spectrum by using 3D radiative transfer code MCMax. Our model calculation shows that dust is dominated by amorphous silicates, with some crystalline silicates, metallic iron, and a very tiny amount of polycyclic aromatic hydrocarbon (PAH) molecules. The presence of both silicates and PAHs indicates that the dust has a mixed chemistry. We derived a dust mass of 0.01 M$_\odot$, from which we arrive at a total ejected mass of $\approx$ 5 M$_\odot$. This implies a time-averaged dust mass-loss rate of 2.5$\times$10$^{-6}$ M$_\odot$\,yr$^{-1}$ with an explosion about 4000 years ago. We assume that the other five confirmed dusty LBVs in the LMC loose mass at a similar rate, and estimate the total contribution to the mass budget of the LMC to be $\approx$ 10$^{-5}$ M$_\odot$\,yr$^{-1}$, which is comparable to the contribution by all the AGB stars in the LMC. Based on our analysis on R71, we speculate that LBVs as a class may be an important dust source in the ISM of the LMC.

The spectra of protoplanetary disks contain mid- and far- infrared emission features produced by forsterite dust grains. The spectral features contain information about the forsterite temperature, chemical composition and grain size. We aim to characterize how the 23 and 69 micron features can be used to constrain the physical locations of forsterite in disks. We check for consistency between two independent forsterite temperature measurements: the 23/69 feature strength ratio and the shape of the 69 micron band. We performed radiative transfer modeling to study the effect of disk properties to the forsterite spectral features. Temperature-dependent forsterite opacities were considered in self-consistent models to compute forsterite emission from protoplanetary disks. Modelling grids are presented to study the effects of grain size, disk gaps, radial mixing and optical depth to the forsterite features. Independent temperature estimates derived from the 23/69 feature strength ratio and the 69 micron band shape are most inconsistent for HD141569 and Oph IRS 48. A case study of the disk of HD141569 shows two solutions to fit the forsterite spectrum. A model with T ~ 40 K, iron-rich (~0-1 % Fe) and 1 micron forsterite grains, and a model with warmer (T ~ 100 K), iron-free, and larger (10 micron) grains. We find that for disks with low upper limits of the 69 micron feature (most notably in flat, self-shadowed disks), the forsterite must be hot, and thus close to the star. We find no correlation between disk gaps and the presence or absence of forsterite features. We argue that the 69 micron feature of the evolved transitional disks HD141569 and Oph IRS 48 is most likely a tracer of larger (i.e. ~10 micron) forsterite grains.

The presence of disks and outflows is widespread among post-AGB binaries. In the first paper of this series, a surprisingly large fraction of optical light was found to be resolved in the 89 Her post-AGB system. The data showed this flux to arise from close to the central binary. Scattering off the inner rim of the circumbinary disk, or in a dusty outflow were suggested as two possible origins. With detailed dust radiative transfer models of the disk we aim to discriminate between these two configurations. By including Herschel/SPIRE photometry, we extend the SED such that it now fully covers UV to sub-mm wavelengths. The MCMax radiative transfer code is used to create a large grid of disk models. Our models include a self-consistent treatment of dust settling as well as of scattering. A Si-rich composition with two additional opacity sources, metallic Fe or amorphous C, are tested. The SED is fit together with mid-IR (MIDI) visibilities as well as the optical and near-IR visibilities of Paper I, to constrain the structure of the disk and in particular of its inner rim. The near-IR visibility data require a smooth inner rim, here obtained with a two-power-law parameterization of the radial surface density distribution. A model can be found that fits all the IR photometric and interferometric data well, with either of the two continuum opacity sources. Our best-fit passive models are characterized by a significant amount of mm-sized grains, which are settled to the midplane of the disk. Not a single disk model fits our data at optical wavelengths though, the reason being the opposing constraints imposed by the optical and near-IR interferometric data. A geometry in which a passive, dusty, and puffed-up circumbinary disk is present, can reproduce all the IR but not the optical observations of 89 Her. Another dusty, outflow or halo, component therefore needs to be added to the system.

Hot exozodiacal dust is thought to be responsible for excess near-infrared (NIR) emission emanating from the innermost parts of some debris disks. The origin of this dust, however, is still a matter of debate. We test whether hot exozodiacal dust can be supplied from an exterior parent belt by Poynting-Robertson (P-R) drag, paying special attention to the pile-up of dust that occurs due to the interplay of P-R drag and dust sublimation. Specifically, we investigate whether pile-ups still occur when collisions are taken into account, and if they can explain the observed NIR excess. We compute the steady-state distribution of dust in the inner disk by solving the continuity equation. First, we derive an analytical solution under a number of simplifying assumptions. Second, we develop a numerical debris disk model that for the first time treats the complex interaction of collisions, P-R drag, and sublimation in a self-consistent way. From the resulting dust distributions we generate thermal emission spectra and compare these to observed excess NIR fluxes. We confirm that P-R drag always supplies a small amount of dust to the sublimation zone, but find that a fully consistent treatment yields a maximum amount of dust that is about 7 times lower than that given by analytical estimates. The NIR excess due this material is much smaller (<10^-3 for A-type stars with parent belts at >1 AU) than the values derived from interferometric observations (~10^-2). Pile-up of dust still occurs when collisions are considered, but its effect on the NIR flux is insignificant. Finally, the cross-section in the innermost regions is clearly dominated by barely bound grains.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15, and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67 [-0.11,+0.18]) and a significantly tapered gap edge ('round wall'), but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images. We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27 [-20,+19] mas) is notably lower than that found in our earlier H-band images (Thalmann et al. 2010). Intriguingly, we also find, at high significance, an offset of y = 69 [-25, +49] mas perpendicular to the line of nodes. If confirmed by future observations, this would imply a highly elliptical -- or otherwise asymmetric -- disk gap with an effective eccentricity of e = ~0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk, and appears brighter than the far side because of strong forward scattering.

Planet-forming disks of gas and dust around young stars contain polycyclic aromatic hydrocarbons (PAHs). We aim to characterize how the charge state of PAHs can be used as a probe of flows of gas through protoplanetary gaps. In this context, our goal is to understand the PAH spectra of four transitional disks. In addition, we want to explain the observed correlation between PAH ionization (traced by the 6.2/11.3 feature ratio) and the disk mass (traced by the 1.3 mm luminosity). We implement a model to calculate the charge state of PAHs in the radiative transfer code MCMax. The emission spectra and ionization balance are calculated. A benchmark modeling grid is presented that shows how PAH ionization and luminosity behave as a function of star and disk properties. The PAH ionization is most sensitive to ultraviolet (UV) radiation and the electron density. In optically thick disks, where the UV field is low and the electron density is high, PAHs are predominantly neutral. Ionized PAHs trace low-density optically thin disk regions where the UV field is high and the electron density is low. Such regions are characteristic of gas flows through the gaps of transitional disks. We demonstrate that fitting the PAH spectra of four transitional disks requires a contribution of ionized PAHs in gas flows through the gap. The PAH spectra of transitional disks can be understood as superpositions of neutral and ionized PAHs. For HD97048, neutral PAHs in the optically thick disk dominate the spectrum. In the cases of HD169142, HD135344B and Oph IRS 48, small amounts of ionized PAHs located in the gas flows through the gap are strong contributors to the total PAH luminosity. The observed trend between the disk mass and PAH ionization may imply that lower-mass disks have larger gaps. Ionized PAHs in gas flows through these gaps contribute strongly to their spectra.

A region of roughly half of the solar system scale around the star HD 100546 is largely cleared of gas and dust, in contrast to the bright outer disc. However, some material is observed in the immediate vicinity of the star. We investigate how the dust is distributed within and outside the gap, and constrain the disc geometry with mid-infrared interferometric observations using VLTI/MIDI. With baseline lengths of 40m, our long baseline observations are sensitive to the inner few AU from the star, and we combined them with observations at shorter, 15m baselines, to probe emission beyond the gap at up to 20AU from the star. We modelled the mid-infrared emission using radial temperature profiles. Our model is composed of infinitesimal concentric annuli emitting as black bodies, and it has distinct inner and outer disc components. We derived an upper limit of 0.7AU for the radial size of the inner disc, from our longest baseline data. This small dusty disc is separated from the edge of the outer disc by a large, roughly 10AU wide gap. Our short baseline data place a bright ring of emission at 11+-1AU, consistent with prior observations of the transition region between the gap and the outer disc, known as the disc wall. The inclination and position angle are constrained by our data to i=53+-8deg and PA=145+-5deg. Compared to the rim and outer disc geometry this suggests co-planarity. Brightness asymmetry is evident in both short and long baseline data, and it is unequivocally discernible from any atmospheric or instrumental effects. The origin of the asymmetry is consistent with the bright disc wall, which we find to be 1-2AU wide. The gap is cleared of micron-sized dust, but we cannot rule out the presence of larger particles and/or perturbing bodies.

Spectra of OH/IR stars show prominent spectral bands of crystalline olivine (Mg$_{(2-2x)}$Fe$_{(2x)}$SiO$_{4}$). To learn more about the timescale of the outflows of OH/IR stars, we study the spectral band of crystalline olivine at 69 {\mu}m. The 69 {\mu}m band is of interest because its width and peak wavelength position are sensitive to the grain temperature and to the exact composition of the crystalline olivine. With Herschel/PACS, we observed the 69 {\mu}m band in the outflow of 14 OH/IR stars. By comparing the crystalline olivine features of our sample with those of model spectra, we determined the size of the outflow and its crystalline olivine abundance. The temperature indicated by the observed 69 {\mu}m bands can only be reproduced by models with a geometrically compact superwind ($R_{\rm{SW}}\lesssim$ 2500 AU = 1400 R$_{*}$).This means that the superwind started less than 1200 years ago (assuming an outflow velocity of 10 km/s). The small amount of mass lost in one superwind and the high progenitor mass of the OH/IR stars introduce a mass loss and thus evolutionary problem for these objects, which has not yet been understood.

The circumstellar environments of classical T Tauri stars are challenging to directly image because of their high star-to-disk contrast ratio. One method to overcome this is by using imaging polarimetry where scattered and consequently polarised starlight from the star's circumstellar disk can be separated from the unpolarised light of the central star. We present images of the circumstellar environment of SU Aur, a classical T Tauri star at the transition of T Tauri to Herbig stars. The images directly show that the disk extends out to ~500 au with an inclination angle of $\sim$ 50$^\circ$. Using interpretive models, we derived very small grains in the surface layers of its disk, with a very steep size- and surface-density distribution. Additionally, we resolved a large and extended nebulosity in our images that is most likely a remnant of the prenatal molecular cloud. The position angle of the disk, determined directly from our images, rules out a polar outflow or jet as the cause of this large-scale nebulosity.

Transitional discs are a special type of protoplanetary discs where planet formation is thought to be taking place. These objects feature characteristic inner cavities and/or gaps of a few tens of AUs in the sub-millimitre images of the disc. This signature suggests a localised depletion of matter in the disc that could be caused by planet formation processes. However, recent observations have revealed differences in the structures imaged at different wavelengths in some of these discs. In this paper, we aim to explain these observational differences using self-consistent physical 2-D hydrodynamical and dust evolution models of such objects, assuming their morphology is indeed generated by the presence of a planet. We use these models to derive the distribution of gas and dust in a theoretical planet-hosting disc, for various planet masses and orbital separations. We then simulate observations of the emitted and scattered light from these models with VLT/SPHERE ZIMPOL, Subaru/HiCIAO, VLT/VISIR and ALMA. We do this by first computing the full resolution images of the models at different wavelengths, and then simulating the observations accounting for the characteristics of each particular instrument. The presence of the planet generates pressure bumps in the gas distribution of the disc whose characteristics strongly depend on the planet mass and position. These bumps cause large grains to accumulate while small grains are allowed into inner regions. This spatial differentiation of the grain sizes explains the differences in the observations since different wavelengths and observing techniques trace different parts of the dust size distribution. Based on this effect, we conclude that the combination of visible/near-infrared polarimetric and sub-mm images is the best strategy to constrain the properties of the unseen planet responsible for the disc structure.

*Context The evolution of young massive protoplanetary disks toward planetary systems is expected to include the formation of gaps and the depletion of dust and gas. *Aims A special group of flaring disks around Herbig Ae/Be stars do not show prominent silicate emission features. We focus our attention on four key Herbig Ae/Be stars to understand the structural properties responsible for the absence of silicate feature emission. *Methods We investigate Q- and N-band images taken with Subaru/COMICS, Gemini South/T-ReCS and VLT/VISIR. Our radiative transfer modeling solutions require a separation of inner- and outer- disks by a large gap. From this we characterize the radial density structure of dust and PAHs in the disk. *Results The inner edge of the outer disk has a high surface brightness and a typical temperature between ~100-150 K and therefore dominates the emission in the Q-band. We derive radii of the inner edge of the outer disk of 34, 23, 30 and 63 AU for HD97048, HD169142, HD135344B and Oph IRS 48 respectively. For HD97048 this is the first detection of a disk gap. The continuum emission in the N-band is not due to emission in the wings of PAHs. This continuum emission can be due to VSGs or to thermal emission from the inner disk. We find that PAH emission is not always dominated by PAHs on the surface of the outer disk. *Conclusions. The absence of silicate emission features is due to the presence of large gaps in the critical temperature regime. Many, if not all Herbig disks with Spectral Energy Distribution (SED) classification `group I' are disks with large gaps and can be characterized as (pre-) transitional. An evolutionary path from the observed group I to the observed group II sources seems no longer likely. Instead, both might derive from a common ancestor.

Context. Many protoplanetary nebulae show strong asymmetries in their surrounding shell, pointing to asymmetries during the mass loss phase. Questions concerning the origin and the onset of deviations from spherical symmetry are important for our understanding of the evolution of these objects. Here we focus on the circumstellar shell of the post-AGB star HD 161796. Aims. We aim at detecting signatures of an aspherical outflow, as well as to derive the properties of it. Methods. We use the imaging polarimeter ExPo (the extreme polarimeter), a visitor instrument at the William Herschel Telescope, to accurately image the dust shell surrounding HD 161796 in various wavelength filters. Imaging polarimetry allows us to separate the faint, polarized, light from circumstellar material from the bright, unpolarized, light from the central star. Results. The shell around HD 161796 is highly aspherical. A clear signature of an equatorial density enhancement can be seen. This structure is optically thick at short wavelengths and changes its appearance to optically thin at longer wavelengths. In the classification of the two different appearances of planetary nebulae from HST images it changes from being classified as DUPLEX at short wavelengths to SOLE at longer wavelengths. This strengthens the interpretation that these two appearances are manifestations of the same physical structure. Furthermore, we find that the central star is hotter than often assumed and the relatively high observed reddening is due to circumstellar rather than interstellar extinction.

Well over 700 exoplanets have been detected to date. Only a handful of these have been observed directly. Direct observation is extremely challenging due to the small separation and very large contrast involved. Imaging polarimetry offers a way to decrease the contrast between the unpolarized starlight and the light that has become linearly polarized after scattering by circumstellar material. This material can be the dust and debris found in circumstellar disks, but also the atmosphere or surface of an exoplanet. We present the design, calibration approach, polarimetric performance and sample observation results of the Extreme Polarimeter, an imaging polarimeter for the study of circumstellar environments in scattered light at visible wavelengths. The polarimeter uses the beam-exchange technique, in which the two orthogonal polarization states are imaged simultaneously and a polarization modulator swaps the polarization states of the two beams before the next image is taken. The instrument currently operates without the aid of Adaptive Optics. To reduce the effects of atmospheric seeing on the polarimetry, the images are taken at a frame rate of 35 fps, and large numbers of frames are combined to obtain the polarization images. Four successful observing runs have been performed using this instrument at the 4.2 m William Herschel Telescope on La Palma, targeting young stars with protoplanetary disks as well as evolved stars surrounded by dusty envelopes. In terms of fractional polarization, the instrument sensitivity is better than 10^-4. The contrast achieved between the central star and the circumstellar source is of the order 10^-6. We show that our calibration approach yields absolute polarization errors below 1%.

Some planetary systems harbour debris disks containing planetesimals such as asteroids and comets. Collisions between such bodies produce small dust particles, the spectral features of which reveal their composition and, hence, that of their parent bodies. A measurement of the composition of olivine crystals has been done for the protoplanetary disk HD100546 and for olivine crystals in the warm inner parts of planetary systems. The latter compares well with the iron-rich olivine in asteroids (x<0.29). In the cold outskirts of the Beta Pictoris system, an analogue to the young Solar System, olivine crystals were detected but their composition remained undetermined, leaving unknown how the composition of the bulk of Solar System cometary olivine grains compares with that of extrasolar comets. Here we report the detection of the 69-micrometre-wavelength band of olivine crystals in the spectrum of Beta Pictoris. Because the disk is optically thin, we can associate the crystals with an extrasolar proto-Kuiper belt a distance of 15-45 astronomical units from the star (one astronomical unit is the Sun-Earth distance), determine their magnesium-rich composition (x=0.01+/-0.001) and show that they make up 3.6+/-1.0 per cent of the total dust mass. These values are strikingly similar to those for the dust emitted by the most primitive comets in the Solar System, even though Beta Pictoris is more massive and more luminous and has a different planetary system architecture.

Scattered light images of circumstellar disks play an important role in characterizing the planet forming environments around young stars. The characteristic size of the scattering dust grains can be estimated from the observed brightness asymmetry between the front and back side of the disk, for example using standard Mie theory. However such models often overpredict their brightness by one or two orders of magnitude, and have difficulty explaining very red disk colors. We aim to develop a dust model that explains simultaneously the observed disk surface brightness, colors and asymmetry in scattered light, focusing on constraining grain sizes. We use the 2D radiative transfer code MCMax with anisotropic scattering to explore the effects of grain size on synthetic scattered light images of circumstellar disks. We compare the results with surface brightness profiles of the protoplanetary disk HD 100546 in scattered light at wavelengths from 0.4 to 2.2 micron. We find that extreme forward scattering by micron sized particles lowers the /effective/ dust albedo and creates a faint and red disk that /appears/ only slightly forward scattering. For the outer (>~100 AU) disk of HD 100546 we derive a minimum grain size of 2.5 micron, likely present in the form of aggregates. Intermediate sized grains are too bright, whereas smaller grains are faint and scatter more isotropically, but also produce disk colors that are too blue. Observed surface brightness asymmetries alone are not sufficient to constrain the grain size in circumstellar disks. Additional information, such as the brightness and colors of the disk are needed to provide additional constraints.

Context. The Kepler object KIC 12557548 b is peculiar. It exhibits transit-like features every 15.7 hours that vary in depth between 0.2% and 1.2%. Rappaport et al. (2012) explain the observations in terms of a disintegrating, rocky planet that has a trailing cloud of dust created and constantly replenished by thermal surface erosion. The variability of the transit depth is then a consequence of changes in the cloud optical depth. Aims. We aim to validate the disintegrating-planet scenario by modeling the detailed shape of the observed light curve, and thereby constrain the cloud particle properties to better understand the nature of this intriguing object. Methods. We analysed the six publicly-available quarters of raw Kepler data, phase-folded the light curve and fitted it to a model for the trailing dust cloud. Constraints on the particle properties were investigated with a light-scattering code. Results. The light curve exhibits clear signatures of light scattering and absorption by dust, including a brightening in flux just before ingress correlated with the transit depth and explained by forward scattering, and an asymmetry in the transit light curve shape, which is easily reproduced by an exponentially decaying distribution of optically thin dust, with a typical grain size of 0.1 micron. Conclusions. Our quantitative analysis supports the hypothesis that the transit signal of KIC 12557548 b is due to a variable cloud of dust, most likely originating from a disintegrating object.

The broad 30 \mu m feature in carbon stars is commonly attributed to MgS dust particles. However, reproducing the 30 \mu m feature with homogeneous MgS grains would require much more sulfur relative to the solar abundance. Direct gas-phase condensation of MgS occurs at a low efficiency. Precipitation of MgS on SiC precursor grains provides a more efficient formation mechanism, such that the assumption of homogeneous MgS grains may not be correct. Using a Monte Carlo-based radiative transfer code, we aim to model the 30 \mu m feature of the extreme carbon star LL Peg with MgS dust particles. We find that for LL Peg this modeling is insensitive to the unknown MgS optical properties at \lambda < 10 \mu m. When MgS is allowed to be in thermal contact with amorphous carbon and SiC, the amount of MgS required to reproduce the strength of 30 \mu m feature agrees with the solar abundance of sulfur, thereby resolving the reported MgS mass problem. We conclude that MgS is a valid candidate to be the carrier of the 30 \mu m feature when it is part of a composite grain population that has optical properties representative of an ensemble of particle shapes.

Z CMa is a complex binary system, composed of a Herbig Be and an FU Ori star. The Herbig star is surrounded by a dust cocoon of variable geometry, and the whole system is surrounded by an infalling envelope. Previous spectropolarimetric observations have reported a preferred orientation of the polarization angle, perpendicular to the direction of a large, parsec-sized jet associated with the Herbig star. The variability in the amount of polarized light has been associated to changes in the geometry of the dust cocoon that surrounds the Herbig star. We aim to constrain the properties of Z CMa by means of imaging polarimetry at optical wavelengths. Using ExPo, a dual-beam imaging polarimeter which operates at optical wavelengths, we have obtained imaging (linear) polarimetric data of Z CMa. Our observations were secured during the return to quiescence after the 2008 outburst. We detect three polarized features over Z CMa. Two of these features are related to the two jets reported in this system: the large jet associated to the Herbig star, and the micro-jet associated to the FU Ori star. Our results suggest that the micro-jet extends to a distance ten times larger than reported in previous studies. The third feature suggests the presence of a hole in the dust cocoon that surrounds the Herbig star of this system. According to our simulations, this hole can produce a pencil beam of light that we see scattered off the low-density envelope surrounding the system.

One of the key predictions of modeling from the IR excess of Herbig Ae stars is that for protoplanetary disks, where significant grain growth and settling has occurred, the dust disk has flattened to the point that it can be partially or largely shadowed by the innermost material at or near the dust sublimation radius. When the self-shadowing has already started, the outer disk is expected to be detected in scattered light only in the exceptional cases that the scale height of the dust disk at the sublimation radius is smaller than usual. High-contrast imaging combined with the IR spectral energy distribution allow us to measure the degree of flattening of the disk, as well as to determine the properties of the outer disk. We present polarimetric differential imaging in $H$ band obtained with Subaru/HiCIAO of one such system, MWC 480. The HiCIAO data were obtained at a historic minimum of the NIR excess. The disk is detected in scattered light from 0\farcs2-1\farcs0 (27.4-137AU). Together with the marginal detection of the disk from 1998 February 24 by HST/NICMOS, our data constrain the opening half angle for the disk to lie between 1.3$\leq\theta\leq 2.2^\circ$. When compared with similar measures in CO for the gas disk from the literature, the dust disk subtends only $\sim$30% of the gas disk scale height (H/R$\sim$0.03). Such a dust disk is a factor of 5-7 flatter than transitional disks, which have structural signatures that giant planets have formed.

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

Fomalhaut is a young, nearby star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and 500 micrometer. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

We present recent polarimetric images of the highly variable star R CrB using ExPo and archival WFPC2 images from the HST. We observed R CrB during its current dramatic minimum where it decreased more than 9 mag due to the formation of an obscuring dust cloud. Since the dust cloud is only in the line-of-sight, it mimics a coronograph allowing the imaging of the star's circumstellar environment. Our polarimetric observations surprisingly show another scattering dust cloud at approximately 1.3" or 2000 AU from the star. We find that to obtain a decrease in the stellar light of 9 mag and with 30% of the light being reemitted at infrared wavelengths (from R CrB's SED) the grains in R CrB's circumstellar environment must have a very low albedo of approximately 0.07%. We show that the properties of the dust clouds formed around R CrB are best fitted using a combination of two distinct populations of grains size. The first are the extremely small 5 nm grains, formed in the low density continuous wind, and the second population of large grains (~0.14 {\mu}m) which are found in the ejected dust clouds. The observed scattering cloud, not only contains such large grains, but is exceptionally massive compared to the average cloud.

SPICES (Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems) is a five-year M-class mission proposed to ESA Cosmic Vision. Its purpose is to image and characterize long-period extrasolar planets and circumstellar disks in the visible (450 - 900 nm) at a spectral resolution of about 40 using both spectroscopy and polarimetry. By 2020/22, present and near-term instruments will have found several tens of planets that SPICES will be able to observe and study in detail. Equipped with a 1.5 m telescope, SPICES can preferentially access exoplanets located at several AUs (0.5-10 AU) from nearby stars ($<$25 pc) with masses ranging from a few Jupiter masses to Super Earths ($\sim$2 Earth radii, $\sim$10 M$_{\oplus}$) as well as circumstellar disks as faint as a few times the zodiacal light in the Solar System.

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

Imaging polarimetry is a powerful tool for imaging faint circumstellar material. For a correct analysis of observations we need to fully understand the effects of dust particle parameters, as well as the effects of the telescope, atmospheric seeing, and assumptions about the data reduction and processing of the observed signal. Here we study the major effects of dust particle structure, size-dependent grain settling, and instrumental properties. We performed radiative transfer modeling using different dust particle models and disk structures. To study the influence of seeing and telescope diffraction we ran the models through an instrument simulator for the ExPo dual-beam imaging polarimeter mounted at the 4.2m William Herschel Telescope (WHT). Particle shape and size have a strong influence on the brightness and detectability of the disks. In the simulated observations, the central resolution element also contains contributions from the inner regions of the protoplanetary disk besides the unpolarized central star. This causes the central resolution element to be polarized, making simple corrections for instrumental polarization difficult. This effect strongly depends on the spatial resolution, so adaptive optics systems are needed for proper polarization calibration. We find that the commonly employed homogeneous sphere model gives results that differ significantly from more realistic models. For a proper analysis of the wealth of data available now or in the near future, one must properly take the effects of particle types and disk structure into account. The observed signal depends strongly on the properties of these more realistic models, thus providing a potentially powerful diagnostic. We conclude that it is important to correctly understand telescope depolarization and calibration effects for a correct interpretation of the degree of polarization.

Context. The circumstellar environment of evolved stars is generally rich in molecular gas and dust. Typically, the entire environment is either oxygen-rich or carbon-rich, depending on the evolution of the central star. Aims. In this paper we discuss three evolved disc sources with evidence of atypical emission lines in their infrared spectra. The stars were taken from a larger sample of post-AGB binaries for which we have Spitzer infrared spectra, characterised by the presence of a stable oxygen-rich circumbinary disc. Our previous studies have shown that the infrared spectra of post-AGB disc sources are dominated by silicate dust emission, often with an extremely high crystallinity fraction. However, the three sources described here are selected because they show a peculiar molecular chemistry. Methods. Using Spitzer infrared spectroscopy, we study in detail the peculiar mineralogy of the three sample stars. Using the observed emission features, we identify the different observed dust, molecular and gas species. Results. The infrared spectra show emission features due to various oxygen-rich dust components, as well as CO2 gas. All three sources show the strong infrared bands generally ascribed to polycyclic aromatic hydrocarbons. Furthermore, two sample sources show C60 fullerene bands. Conclusions. Even though the majority of post-AGB disc sources are dominated by silicate dust in their circumstellar environment, we do find evidence that, for some sources at least, additional processing must occur to explain the presence of large carbonaceous molecules. There is evidence that some of these sources are still oxygen-rich, which makes the detection of these molecules even more surprising.

The dust size distribution in molecular clouds can be strongly affected by ice-mantle formation and (subsequent) grain coagulation. Following previous work where the dust size distribution has been calculated from a state-of-the art collision model for dust aggregates that involves both coagulation and fragmentation (Paper I), the corresponding opacities are presented in this study. The opacities are calculated by applying the effective medium theory assuming that the dust aggregates are a mix of 0.1{\mu}m silicate and graphite grains and vacuum. In particular, we explore how the coagulation affects the near-IR opacities and the opacity in the 9.7{\mu}m silicate feature. We find that as dust aggregates grow to {\mu}m-sizes both the near-IR color excess and the opacity in the 9.7 {\mu}m feature increases. Despite their coagulation, porous aggregates help to prolong the presence of the 9.7{\mu}m feature. We find that the ratio between the opacity in the silicate feature and the near-IR color excess becomes lower with respect to the ISM, in accordance with many observations of dark clouds. However, this trend is primarily a result of ice mantle formation and the mixed material composition of the aggregates, rather than being driven by coagulation. With stronger growth, when most of the dust mass resides in particles of size 10{\mu}m or larger, both the near-IR color excess and the 9.7{\mu}m silicate feature significantly diminish. Observations at additional wavelengths, in particular in the sub-mm range, are essential to provide quantitative constraints on the dust size distribution within dense cores. Our results indicate that the sub-mm index {\beta} will increase appreciably, if aggregates grow to ~100{\mu}m in size.

Forsterite is one of the crystalline dust species that is often observed in protoplanetary disks and solar system comets. Being absent in the interstellar medium, it must be produced during the disk lifetime. It can therefore serve as a tracer of dust processing and disk evolution, which can lead to a better understanding of the physical processes occurring in the disk, and possibly planet formation. However, the connection of these processes with the overall disk crystallinity remains unclear. We aim to characterize the forsterite abundance and spatial distribution in the disk of the Herbig Be star HD 100546, to investigate if a connection exists with the large disk gap. We use a 2D radiative transfer code, MCMax, to model the circumstellar dust around HD 100546. We use VISIR Q-band imaging to probe the outer disk geometry and mid-infrared features to model the spatial distribution of forsterite. The temperature-dependent shape of the 69 micron feature observed with Herschel PACS is used as a critical tool to constrain this distribution. We find a crystalline mass fraction of 40 - 60 %, located close to the disk wall between 13 and 20 AU, and possibly farther out at the disk surface. The forsterite is in thermal contact with the other dust species. We put an upper limit on the iron content of forsterite of 0.3 %. Optical depth effects play a key role in explaining the observed forsterite features, hiding warm forsterite from view at short wavelengths. The disk wall acts as a showcase: it displays a localized high abundance of forsterite, which gives rise to a high observed crystallinity, while the overall mass fraction of forsterite is a factor of ten lower.

Imaging polarimetry is a powerful tool for detecting and characterizing exoplanets and circumstellar environments. Polarimetry allows a separation of the light coming from an unpolarized source such as a star and the polarized source such as a planet or a protoplanetary disk. Future facilities like SPHERE at the VLT or EPICS at the E-ELT will incorporate imaging polarimetry to detect exoplanets. The Extreme Polarimeter (ExPo) is a dual-beam imaging polarimeter that currently can reach contrast ratios of 10^5, enough to characterize circumstellar environments. We present the data reduction steps for a dual-beam imaging polarimeter that can reach contrast ratios of 10^5. The data obtained with ExPo at the William Herschel Telescope (WHT) are analyzed. Instrumental artifacts and noise sources are discussed for an unpolarized star and for a protoplanetary disk (AB Aurigae). The combination of fast modulation and dual-beam techniques allow us to minimize instrumental artifacts. A proper data processing and alignment of the images is fundamental when dealing with large contrasts. Imaging polarimetry proves to be a powerful method to resolve circumstellar environments even without a coronagraph mask or an Adaptive Optics system.

The recent findings of gas giant planets around young A-type stars suggest that disks surrounding Herbig Ae/Be stars will develop planetary systems. An interesting case is HD142527, for which previous observations revealed a complex circumstellar environment and an unusually high ratio of infrared to stellar luminosity. Its properties differ considerably from other Herbig Ae/Be stars. This suggests that the disk surrounding HD142527 is in an uncommon evolutionary stage. We aim for a better understanding of the geometry and evolutionary status of the circumstellar material around the Herbig Ae/Be star HD142527. We map the composition and spatial distribution of the dust around HD142527. We analyze SEST and ATCA millimeter data, VISIR N and Q-band imaging and spectroscopy. We gather additional relevant data from the literature. We use the radiative transfer code MCMax to construct a model of the geometry and density structure of the circumstellar matter, which fits all of the observables satisfactorily. We find that the disk of HD142527 has three geometrically distinct components separated by a disk gap running from 30 to 130 AU. There is a geometrically flat inner disk running from 0.3 AU up to 30 AU; an optically thin halo-like component of dust in the inner disk regions; and a massive self-shadowed outer disk running from 130 AU up to 200 AU. We derived a total dust mass in small grains of 10^-3 M0 and a vertical height of the inner wall of the outer disk of h = 60 AU. ....

Studying the composition of dust in the interstellar medium (ISM) is crucial in understanding the cycle of dust in our galaxy. The mid-infrared spectral signature of amorphous silicates, the most abundant dust species in the ISM, is studied in different lines-of-sight through the Galactic plane, thus probing different conditions in the ISM. We have analysed 10 spectra from the Spitzer archive, of which 6 lines-of-sight probe diffuse interstellar medium material and 4 probe molecular cloud material. The 9.7 um silicate absorption features in 7 of these spectra were studied in terms of their shape and strength. In addition, the shape of the 18 um silicate absorption features in 4 of the diffuse sightline spectra were analysed. The 9.7 um silicate absorption bands in the diffuse sightlines show a strikingly similar band shape. This is also the case for all but one of the 18 um silicate absorption bands observed in diffuse lines-of-sight. The 9.7 um bands in the 4 molecular sightlines show small variations in shape. These modest variations in the band shape are inconsistent with the interpretation of the large variations in {\tau}_9.7/E(J-K) between diffuse and molecular sightlines in terms of silicate grain growth. Instead, we suggest that the large changes in {\tau}_9.7 / E(J-K) must be due to changes in E(J-K).

The precise location of the water ice condensation front ('snow line') in the protosolar nebula has been a debate for a long time. Its importance stems from the expected substantial jump in the abundance of solids beyond the snow line, which is conducive to planet formation, and from the higher stickiness in collisions of ice-coated dust grains, which may help the process of coagulation of dust and the formation of planetesimals. In an optically thin nebula, the location of the snow line is easily calculated to be around 3 AU. However, in its first 5 to 10 million years, the solar nebula was optically thick, implying a smaller snow line radius due to shielding from direct sunlight, but also a larger radius because of viscous heating. Several models have attempted to treat these opposing effects. However, until recently treatments beyond an approximate 1+1D radiative transfer were unfeasible. We revisit the problem with a fully self-consistent 3D treatment in an axisymmetric disk model, including a density-dependent treatment of the dust and ice sublimation. We find that the location of the snow line is very sensitive to the opacities of the dust grains and the mass accretion rate of the disk. We show that previous approximate treatments are quite efficient at determining the location of the snow line if the energy budget is locally dominated by viscous accretion. We derive an analytic estimate of the location of the snow line that compares very well with results from this and previous studies. Using solar abundances of the elements we compute the abundance of dust and ice and find that the expected jump in solid surface density at the snow line is smaller than previously assumed. We further show that in the inner few AU the refractory species are partly evaporated, leading to a significantly smaller solid state surface density in the regions where the rocky planets were formed.

We present a diagnostic tool to determine the forsterite abundance of the dust ejected by AGB stars. Our method is based on a comparison between the observed strength of spectral bands of forsterite and model calculations. We show that the 11.3 {\mu}m forsterite band is a robust indicator of the forsterite abundance of the current mass-loss period for AGB stars with an optically thick dust shell. The 33.6 {\mu}m band of forsterite is sensitive to changes in the density and the geometry of the emitting dust shell, and so a less robust indicator. We apply this method to six high mass-loss rate AGB stars, showing that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum abundance of 5%.

In this paper we present mid-infrared spectra of a comprehensive set of Herbig Ae/Be stars observed with the Spitzer Space Telescope. The signal-to-noise ratio of these spectra is very high, ranging between about a hundred and several hundreds. During the analysis of these data we tested the validity of standard protoplanetary dust models and studied grain growth and crystal formation. On the basis of the analyzed spectra, the major constituents of protoplanetary dust around Herbig Ae/Be stars are amorphous silicates with olivine and pyroxene stoichiometry, crystalline forsterite and enstatite and silica. No other solid state features, indicating other abundant dust species, are present in the Spitzer spectra. Deviations of the synthetic spectra from the observations are most likely related to grain shape effects and uncertainties in the iron content of the dust grains. Our analysis revealed that larger grains are more abundant in the disk atmosphere of flatter disks than in that of flared disks, indicating that grain growth and sedimentation decrease the disk flaring. We did not find, however, correlations between the value of crystallinity and any of the investigated system parameters. Our analysis shows that enstatite is more concentrated toward the warm inner disk than forsterite, in contrast to predictions of equilibrium condensation models. None of the three crystal formation mechanisms proposed so far can alone explain all our findings. It is very likely that all three play at least some role in the formation of crystalline silicates.

We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

With the legacy of Spitzer and current advances in (sub)mm astronomy, a large number of 'transitional' disks has been identified which are believed to contain gaps or have developped large inner holes, some filled with dust. This may indicate that complex geometries may be a key feature in disk evolution that has to be understood and modeled correctly. The disk around LkCa 15 is such a disk, with a gap ranging from ~5 - 50 AU, as identified by Espaillat et al. (2007) using 1+1D radiative transfer modelling. To fit the SED, they propose 2 possible scenarios for the inner (<5 AU) disk - optically thick or optically thin - and one scenario for the outer disk. We use the gapped disk of LkCa 15 as a showcase to illustrate the importance of 2D radiative transfer in transitional disks, by showing how the vertical dust distribution in dust-filled inner holes determines not only the radial optical depth but also the outer disk geometry. We use MCMax, a 2D radiative transfer code with a self-consistent vertical structure, to model the SED. We identify two possible geometries for the inner and outer disk, that are both different from those in Espaillat et al. (2007). An inner disk in hydrostatic equilibrium reprocesses enough starlight to fit the near infrared flux, but also casts a shadow on the inner rim of the outer disk. This requires the outer disk scale height to be large enough to rise out of the shadow. An optically thin inner disk does not cast such a shadow, and the SED can be fitted with a smaller outer disk scale height. For the dust in the inner regions to become optically thin however, the scale height would have to be so much larger than its hydrostatic equilibrium value that it effectively becomes a dust shell. It is currently unclear if a physical mechanism exists which could provide for such a configuration.

Context. Based on the far infrared excess the Herbig class of stars is divided into a group with flaring circumstellar disks (group I) and a group with flat circumstellar disks (group II). Dust sedimentation is generally proposed as an evolution mechanism to transform flaring disks into flat disks. Theory predicts that during this process the disks preserve their gas content, however observations of group II Herbig Ae stars demonstrate a lack of gas. Aims. We map the spatial distribution of the gas and dust around the group II Herbig Ae star HD 95881. Methods. We analyze optical photometry, Q-band imaging, infrared spectroscopy, and K and N-band interferometric spectroscopy. We use a Monte Carlo radiative transfer code to create a model for the density and temperature structure which quite accurately reproduces all the observables. Results. We derive a consistent picture in which the disk consists of a thick puffed up inner rim and an outer region which has a flaring gas surface and is relatively void of 'visible' dust grains. Conclusions. HD 95881 is in a transition phase from a gas rich flaring disk to a gas poor self-shadowed disk.

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.

Aims. We present a diagnostic tool to determine the abundance of the crystalline silicate forsterite in AGB stars surrounded by a thick shell of silicate dust. Using six infrared spectra of high mass-loss oxygen rich AGB stars we obtain the forsterite abundance of their dust shells. Methods. We use a monte carlo radiative transfer code to calculate infrared spectra of dust enshrouded AGB stars. We vary the dust composition, mass-loss rate and outer radius. We focus on the strength of the 11.3 and the 33.6 \mu m forsterite bands, that probe the most recent (11.3 \mu m) and older (33.6 \mu m) mass-loss history of the star. Simple diagnostic diagrams are derived, allowing direct comparison to observed band strengths. Results. Our analysis shows that the 11.3 \mu m forsterite band is a robust indicator for the forsterite abundance of the current mass-loss period for AGB stars with an optically thick dust shell. The 33.6 \mu m band of forsterite is sensitive to changes in the density and the geometry of the emitting dust shell, and so a less robust indicator. Applying our method to six high mass-loss rate AGB stars shows that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum abundance of 5%.

Optical images of the nearby star Fomalhaut show a ring of dust orbiting the central star. This dust is in many respects expected to be similar to the zodiacal dust in the solar system. The ring displays a clear brightness asymmetry, attributed to asymmetric scattering of the central starlight by the circumstellar dust grains. Recent measurements show that the bright side of the Fomalhaut ring is oriented away from us. This implies that the grains in this system scatter most of the light in the backward direction, in sharp contrast to the forward-scattering nature of the grains in the solar system. In this letter, we show that grains considerably larger than those dominating the solar system zodiacal dust cloud provide a natural explanation for the apparent backward scattering behavior. In fact, we see the phases of the dust grains in the same way as we can observe the phases of the Moon and other large solar system bodies. We outline how the theory of the scattering behavior of planetesimals can be used to explain the Fomalhaut dust properties. This indicates that the Fomalhaut dust ring is dominated by very large grains. The material orbiting Fomalhaut, which is at the transition between dust and planetesimals, can, with respect to their optical behavior, best be described as micro-asteroids.

The inner boundary of protoplanetary discs is structured by the dramatic opacity changes at the transition from the dust-containing to a dust-free zone. This paper explores the variety and limits of inner rim structures in passively heated dusty discs. For this study, we implemented detailed sublimation physics in a fast Monte Carlo radiative transfer code. We show that the inner rim in dusty discs is not an infinitely sharp wall but a diffuse region which may be narrow or wide. Furthermore, high surface densities and large silicate grains as well as iron and corundum grains decrease the rim radius, from a 2.2AU radius for small silicates around a 47 Solar luminosity Herbig Ae star typically to 0.4AU and as close as 0.2AU. A passive disc with grain growth and a diverse dust composition must thus have a small inner rim radius. Finally, an analytical expression is presented for the rim location as a function of dust, disc and stellar properties.

Analysis of thermal dust emission spectra for dust mineralogy and physical grain properties depends on laboratory-measured or calculated comparison spectra. Often, the agreement between these two kinds of spectra is not satisfactory because of the strong influence of the grain morphology on the spectra. We investigate the ability of the statistical light-scattering model with a distribution of form factors (DFF model) to reproduce experimentally measured infrared extinction spectra for particles that are small compared to the wavelength. We take advantage of new experimental spectra measured for free particles dispersed in air with accompanying information on the grain morphology. For the calculations, we used DFFs that were derived for aggregates of spherical grains, as well as for compact grain shapes corresponding to Gaussian random spheres. Irregular particle shapes require a DFF similar to that of a Gaussian random sphere with sigma=0.3, whereas roundish grain shapes are best fitted with that of a fractal aggregate of a fractal dimension 2.4-1.8. In addition we used a fitting algorithm to obtain the best-fit DFFs for the various laboratory samples. In this way we can independently derive information on the shape of the grains from their infrared spectra. For anisotropic materials, different DFFs are needed for the different crystallographic axes. This is due to a theoretical problem, which is inherent to all models that are simply averaging the contributions of the crystallographic directions.

Based on the SED, Herbig stars have been categorized into two observational groups, reflecting their overall disk structure: group I members have disks with a higher degree of flaring than their group II counterparts. We investigate the 5-35 um Spitzer IRS spectra of a sample of 13 group I sources and 20 group II sources. We focus on the continuum emission to study the underlying disk geometry. We have determined the [30/13.5] and [13.5/7] continuum flux ratios. The 7-um flux excess with respect to the stellar photosphere is measured, as a marker for the strength of the near-IR emission produced by the inner disk. We have compared our data to self-consistent passive-disk model spectra, for which the same quantities were derived. We confirm the literature result that the difference in continuum emission between group I and II sources can largely be explained by a different amount of small dust grains. However, we report a strong correlation between the [30/13.5] and [13.5/7] flux ratios for Meeus group II sources. Moreover, the [30/13.5] flux ratio decreases with increasing 7-um excess for all targets in the sample. To explain these correlations with the models, we need to introduce an artificial scaling factor for the inner disk height. In roughly 50% of the Herbig Ae/Be stars in our sample, the inner disk must be inflated by a factor 2 to 3 beyond what hydrostatic calculations predict. The total disk mass in small dust grains determines the degree of flaring. We conclude, however, that for any given disk mass in small dust grains, the shadowing of the outer (tens of AU) disk is determined by the scale height of the inner disk (1 AU). The inner disk partially obscures the outer disk, reducing the disk surface temperature. Here, for the first time, we prove these effects observationally.

Aim: We aim to study in detail the peculiar mineralogy and structure of the circumstellar environment of two binary post-AGB stars, EPLyr and HD52961. Both stars were selected from a larger sample of evolved disc sources observed with Spitzer and show unique solid-state and gas features in their infrared spectra. Moreover, they show a very small infrared excess in comparison with the other sample stars. Methods: The different dust and gas species are identified on the basis of high-resolution Spitzer-IRS spectra. We fit the full spectrum to constrain grain sizes and temperature distributions in the discs. This, combined with our broad-band spectral energy distribution and interferometric measurements, allows us to study the physical structure of the disc, using a self-consistent 2D radiative-transfer disc model. Results: We find that both stars have strong emission features due to CO_2 gas, dominated by ^{12}C^{16}O_2, but with clear ^{13}C^{16}O_2 and even ^{16}O^{12}C^{18}O isotopic signatures. Crystalline silicates are apparent in both sources but proved very hard to model. EP Lyr also shows evidence of mixed chemistry, with emission features of the rare class-C PAHs. Whether these PAHs reside in the oxygen-rich disc or in a carbon-rich outflow is still unclear. With the strongly processed silicates, the mixed chemistry and the low ^{12}C/^{13}C ratio, EP Lyr resembles some silicate J-type stars, although the depleted photosphere makes nucleosynthetic signatures difficult to probe. We find that the disc environment of both sources is, to a first approximation, well modelled with a passive disc, but additional physics such as grain settling, radial dust distributions, and an outflow component must be included to explain the details of the observed spectral energy distributions in both stars.

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.

In this paper we present two efficient implementations of the diffusion approximation to be employed in Monte Carlo computations of radiative transfer in dusty media of massive circumstellar disks. The aim is to improve the accuracy of the computed temperature structure and to decrease the computation time. The accuracy, efficiency and applicability of the methods in various corners of parameter space are investigated. The effects of using these methods on the vertical structure of the circumstellar disk as obtained from hydrostatic equilibrium computations are also addressed. Two methods are presented. First, an energy diffusion approximation is used to improve the accuracy of the temperature structure in highly obscured regions of the disk, where photon counts are low. Second, a modified random walk approximation is employed to decrease the computation time. This modified random walk ensures that the photons that end up in the high-density regions can quickly escape to the lower density regions, while the energy deposited by these photons in the disk is still computed accurately. A new radiative transfer code, MCMax, is presented in which both these diffusion approximations are implemented. These can be used simultaneously to increase both computational speed and decrease statistical noise. We conclude that the diffusion approximations allow for fast and accurate computations of the temperature structure, vertical disk structure and observables of very optically thick circumstellar disks.

Aims: We investigate the mineralogy and dust processing in the circumbinary discs of binary post-AGB stars using high-resolution TIMMI2 and SPITZER infrared spectra. Methods: We perform a full spectral fitting to the infrared spectra using the most recent opacities of amorphous and crystalline dust species. This allows for the identification of the carriers of the different emission bands. Our fits also constrain the physical properties of different dust species and grain sizes responsible for the observed emission features. Results: In all stars the dust is oxygen-rich: amorphous and crystalline silicate dust species prevail and no features of a carbon-rich component can be found, the exception being EPLyr, where a mixed chemistry of both oxygen- and carbon-rich species is found. Our full spectral fitting indicates a high degree of dust grain processing. The mineralogy of our sample stars shows that the dust is constituted of irregularly shaped and relatively large grains, with typical grain sizes larger than 2 micron. The spectra of nearly all stars show a high degree of crystallinity, where magnesium-rich end members of olivine and pyroxene silicates dominate. Other dust features of e.g. silica or alumina are not present at detectable levels. Temperature estimates from our fitting routine show that a significant fraction of grains must be cool, significantly cooler than the glass temperature. This shows that radial mixing is very efficient is these discs and/or indicates different thermal conditions at grain formation. Our results show that strong grain processing is not limited to young stellar objects and that the physical processes occurring in the discs are very similar to those in protoplanetary discs.

In order to deduce properties of dust in astrophysical environments where dust growth through aggregation is important, knowledge of the way aggregated particles interact with radiation, and what information is encoded in the thermal radiation they emit, is needed. The emission characteristics are determined by the size and structure of the aggregate and the composition and shape of the constituents. We thus aim at performing computations of compositionally inhomogeneous aggregates composed of irregularly shaped constituents. In addition we aim at developing an empirical recipe to compute the optical properties of such aggregates in a fast and accurate manner. We performed CDA computations for aggregates of irregularly shaped particles with various compositions. The constituents of the aggregate are assumed to be in the Rayleigh regime, and in addition we assume that the dominant interaction of the aggregate constituents is through dipole-dipole interactions. We computed the spectral structure of the emission efficiency in the 10 micron region for aggregates with 30% amorphous carbon and 70% silicates by volume with various fractions of crystalline and amorphous components. We find that the spectral appearance of the various components of the aggregate are very different and depend on their abundances. Most notably, materials that have a very low abundance appear spectroscopically as if they were in very small grains, while more abundant materials appear, spectroscopically to reside in larger grains. We construct a fast empirical approximate method, based on the idea of an effective medium approximation, to construct the spectra for these aggregates which almost perfectly reproduces the more exact computations. This new method is fast enough to be easily implemented in fitting procedures.

Context: We study the structure and evolution of circumstellar discs around evolved binaries and their impact on the evolution of the central system. Aims: To study in detail the binary nature of RUCen and ACHer, as well as the structure and mineralogy of the circumstellar environment. Methods: We combine multi-wavelength observations with a 2D radiative transfer study. Our radial velocity program studies the central stars, while our Spitzer spectra and broad-band SEDs are used to constrain mineralogy, grain sizes and physical structure of the circumstellar environment. Results: We determine the orbital elements of RUCen showing that the orbit is highly eccentric with a rather long period of 1500 days. The infrared spectra of both objects are very similar and the spectral dust features are dominated by Mg-rich crystalline silicates. The small peak-to-continuum ratios are interpreted as being due to large grains. Our model contains two components with a cold midplain dominated by large grains, and the near- and mid-IR which is dominated by the emission of smaller silicates. The infrared excess is well modelled assuming a hydrostatic passive irradiated disc. The profile-fitting of the dust resonances shows that the grains must be very irregular. Conclusions: These two prototypical RVTauri pulsators with circumstellar dust are binaries where the dust is trapped in a stable disc. The mineralogy and grain sizes show that the dust is highly processed, both in crystallinity and grain size. The cool crystals show that either radial mixing is very efficient and/or that the thermal history at grain formation has been very different from that in outflows. The physical processes governing the structure of these discs are similar to those observed in protoplanetary discs around young stellar objects.

In the generally accepted, but poorly documented model, silicate J-type C-stars are binary objects for which the silicate emission originates from a circumbinary or a circumcompanion disc. We aim at testing this hypothesis by a thorough spectral and spatial observational study of one object: IRAS18006-3213. We obtained, analysed and modeled high spatial resolution interferometric VLTI/MIDI observations on multiple baselines ranging from 45 m to 100 m. All observations resolved the object and show the very compact nature of the N-band emission (~30 mas). In addition, the highest spatial resolution data show a significant differential phase jump around 8.3 micron. This demonstrates the asymmetric nature of the N-band emission. Moreover, the single telescope N-band spectrum shows the signature of highly processed silicate grains. These data are used to confirm the model on silicate J-type C-stars for IRAS18006-3213. We show that the most favourable model of the dust geometry is a stable circumbinary disc around the system, seen under an intermediate inclination. The data presented on the silicate J-type C-star IRAS18006-3213 provide evidence that the oxygen rich dust is trapped in a circumbinary disc. The formation of this disc is probably linked to the binary nature of the central star.

We investigate the composition and shape distribution of silicate dust grains in the interstellar medium. The effect of the amount of magnesium in the silicate lattice is studied. We fit the spectral shape of the interstellar 10 mu extinction feature as observed towards the galactic center. We use very irregularly shaped coated and non-coated porous Gaussian Random Field particles as well as a statistical approach to model shape effects. For the dust materials we use amorphous and crystalline silicates with various composition and SiC. The results of our analysis of the 10 mu feature are used to compute the shape of the 20 mu silicate feature and to compare this with observations. By using realistic particle shapes we are, for the first time, able to derive the magnesium fraction in interstellar silicates. We find that the interstellar silicates are highly magnesium rich (Mg/(Fe+Mg)>0.9) and that the stoichiometry lies between pyroxene and olivine type silicates. This composition is not consistent with that of the glassy material found in GEMS in interplanetary dust particles indicating that these are, in general, not unprocessed remnants from the interstellar medium. Also, we find a significant fraction of SiC (~3%). We discuss the implications of our results for the formation and evolutionary history of cometary and circumstellar dust. We argue that the fact that crystalline silicates in cometary and circumstellar grains are almost purely magnesium silicates is a natural consequence of our findings that the amorphous silicates from which they were formed were already magnesium rich.

We report on observations of circumstellar disks around young stars that have been obtained with the MIDI instrument, which is mounted on the VLT Interferometer and operates in the 10 micrometer atmospheric window. The maximum spatial resolution of 5 milli-arcsec corresponds to sub-AU scales at the distance to nearby star formation regions. Thus, we can study the disks on the spatial scales at which important processes occur, such as accretion, dust processing, and planet formation. The main results obtained so far can be summarized as follows: 1. The measured interferometric visibilities are in good qualitative agreement with those predicted by models of circumstellar disks. In particular, a predicted correlation between the strength of the far-infrared excess and the spatial structure of the disk is confirmed by direct measurements; 2. In several objects strong evidence for deviations from circular symmetry is present, indicating that an inclined disk is indeed the dominant component seen in the mid-infrared; 3. The dust properties are not uniform over the disk, but are instead a strong function of distance to the central star. The dust in the innermost disk regions is observed to be more ``processed'' than the dust further out, both in Herbig Ae star disks and in those around T-Tauri stars.

We present the first mid-IR long baseline interferometric observations of the circumstellar matter around binary post-AGB stars. Two objects, SX Cen and HD 52961, were observed using the VLTI/MIDI instrument during Science Demonstration Time. Both objects are known binaries for which a stable circumbinary disc is proposed to explain the SED characteristics. This is corroborated by our N-band spectrum showing a crystallinity fraction of more than 50 % for both objects, pointing to a stable environment where dust processing can occur. Surprisingly, the dust surrounding SX Cen is not resolved in the interferometric observations providing an upper limit of 11 mas (or 18 AU at the distance of this object) on the diameter of the dust emission. This confirms the very compact nature of its circumstellar environment. The dust emission around HD 52961 originates from a very small but resolved region, estimated to be ~ 35 mas at 8 micron and ~ 55 mas at 13 micron. These results confirm the disc interpretation of the SED of both stars. In HD 52961, the dust is not homogeneous in its chemical composition: the crystallinity is clearly concentrated in the hotter inner region. Whether this is a result of the formation process of the disc, or due to annealing during the long storage time in the disc is not clear.

Particle shape and aggregation have a strong influence on the spectral profiles of infrared phonon bands of solid dust grains. Calculating these effects is difficult due to the often extreme refractive index values in these bands. In this paper, we use the Discrete Dipole Approximation (DDA) and the T-matrix method to compute the absorption band profiles for simple clusters of touching spherical grains. We invest reasonable amounts of computation time in order to reach high dipole grid resolutions and take high multi-polar orders into account, respectively. The infrared phonon bands of three different refractory materials of astrophysical relevance are considered - Silicon Carbide (SiC), Wustite (FeO) and Silicon Dioxide (SiO2). We demonstrate that even though these materials display a range of material properties and therefore different strengths of the surface resonances, a complete convergence is obtained with none of the approaches. For the DDA, we find a strong dependence of the calculated band profiles on the exact dipole distribution within the aggregates, especially in the vicinity of the contact points between their spherical constituents. By applying a recently developed method to separate the material optical constants from the geometrical parameters in the DDA approach, we are able to demonstrate that the most critical material properties are those where the real part of the refractive index is much smaller than unity.

We model the 10 micron absorption spectra of nonspherical particles composed of amorphous silicate. We consider two classes of particles, compact ones and fractal aggregates composed of homogeneous spheres. For the compact particles we consider Gaussian random spheres with various degrees of non-sphericity. For the fractal aggregates we compute the absorption spectra for various fractal dimensions. The 10 micron spectra are computed for ensembles of these particles in random orientation using the well-known Discrete Dipole Approximation. We compare our results to spectra obtained when using volume equivalent homogeneous spheres and to those computed using a porous sphere approximation. We conclude that, in general, nonspherical particles show a spectral signature that is similar to that of homogeneous spheres with a smaller material volume. This effect is overestimated when approximating the particles by porous spheres with the same volume filling fraction. For aggregates with fractal dimensions typically predicted for cosmic dust, we show that the spectral signature characteristic of very small homogeneous spheres (with a volume equivalent radius r_V<0.5 micron) can be detected even in very large particles. We conclude that particle sizes are underestimated when using homogeneous spheres to model the emission spectra of astronomical sources. In contrast, the particle sizes are severely overestimated when using equivalent porous spheres to fit observations of 10 micron silicate emission.

The core of the nebula surrounding Eta Carinae has recently been observed with VLT/NACO, VLTI/VINCI, VLTI/MIDI and VLTI/AMBER in order to spatially and spectrally constrain the warm dusty environment and the central object. Narrow-band images at 3.74 and 4.05 micron reveal the structured butterfly-shaped dusty environment close to the central star with an unprecedented spatial resolution of about 60 mas. VINCI has resolved the present-day stellar wind of Eta Carinae on a scale of several stellar radii owing to the spatial resolution of the order of 5 mas (11 AU). The VINCI observations show that the object is elongated with a de-projected axis ratio of approximately 1.5. Moreover the major axis is aligned with that of the large bipolar nebula that was ejected in the 19th century. Fringes have also been obtained in the Mid-IR with MIDI using baselines of 75m. A peak of correlated flux of 100 Jy is detected 0.3" south-east from the photocenter of the nebula at 8.7 micron is detected. This correlated flux is partly attributed to the central object but it is worth noting that at these wavelengths, virtually all the 0.5" x 0.5" central area can generate detectable fringes witnessing the large clumping of the dusty ejecta. These observations provide an upper limit for the SED of the central source from 3.8 to 13.5 micron and constrain some parameters of the stellar wind which can be compared to Hillier's model. Lastly, we present the great potential of the AMBER instrument to study the numerous near-IR emissive lines from the star and its close vicinity. In particular, we discuss its ability to detect and follow the faint companion.

We discuss the composition and size distribution of the dust in the coma of comet Hale-Bopp. We do this by fitting simultaneously the infrared emission spectrum measured by the infrared space observatory (ISO) and the measured degree of linear polarization of scattered light at various phase angles and 12 different wavelengths. The effects of particle shape on the modeled optical properties of the dust grains are taken into account. We constrain our fit by forcing the abundances of the major rock forming chemical elements to be solar. The infrared spectrum at long wavelengths reveals that large grains are needed in order to fit the spectral slope. The size and shape distribution we employ allows us to estimate the sizes of the crystalline silicates. The ratios of the strength of various forsterite features show that the crystalline silicate grains in Hale-Bopp must be submicron sized. We exclude the presence of large crystalline silicate grains in the coma. Because of this lack of large crystalline grains combined with the fact that we do need large amorphous grains to fit the emission spectrum at long wavelengths, we need only approximately 4% of crystalline silicates by mass. After correcting for possible hidden crystalline material included in large amorphous grains, our best estimate of the total mass fraction of crystalline material is approximately 7.5%, significantly lower than deduced in previous studies in which the typical derived crystallinity is 20-30%. The implications of this on the possible origin and evolution of the comet are discussed. The crystallinity we observe in Hale-Bopp is consistent with the production of crystalline silicates in the inner solar system by thermal annealing and subsequent radial mixing to the comet forming region.

We present spectroscopic observations of a large sample of Herbig Ae stars in the 10 micrometer spectral region. We perform compositional fits of the spectra based on properties of homogeneous as well as inhomogeneous spherical particles, and derive the mineralogy and typical grain sizes of the dust responsible for the 10 $\mu$m emission. Several trends are reported that can constrain theoretical models of dust processing in these systems: i) none of the sources consists of fully pristine dust comparable to that found in the interstellar medium, ii) all sources with a high fraction of crystalline silicates are dominated by large grains, iii) the disks around more massive stars (M >~ 2.5 M_sun, L >~ 60 L_sun) have a higher fraction of crystalline silicates than those around lower mass stars, iv) in the subset of lower mass stars (M <~ 2.5 M_sun) there is no correlation between stellar parameters and the derived crystallinity of the dust. The correlation between the shape and strength of the 10 micron silicate feature reported by van Boekel et al. (2003) is reconfirmed with this larger sample. The evidence presented in this paper is combined with that of other studies to present a likely scenario of dust processing in Herbig Ae systems. We conclude that the present data favour a scenario in which the crystalline silicates are produced in the innermost regions of the disk, close to the star, and transported outward to the regions where they can be detected by means of 10 micron spectroscopy. Additionally, we conclude that the final crystallinity of these disks is reached very soon after active accretion has stopped.

We provide a theoretical foundation for the statistical approach for computing the absorption properties of particles in the Rayleigh domain. We present a general method based on the Discrete Dipole Approximation (DDA) to compute the absorption and scattering properties of particles in the Rayleigh domain. The method allows to separate the geometrical aspects of a particle from its material properties. Doing the computation of the optical properties of a particle once, provides them for any set of refractive indices, wavelengths and orientations. This allows for fast computations of e.g. absorption spectra of arbitrarily shaped particles. Other practical applications of the method are in the interpretation of atmospheric and radar measurements as well as computations of the scattering matrix of small particles as a function of the scattering angle. In the statistical approach, the optical properties of irregularly shaped particles are represented by the average properties of an ensemble of particles with simple shapes. We show that the absorption cross section of an ensemble of arbitrarily shaped particles with arbitrary orientations can always be uniquely represented by the average absorption cross section of an ensemble of spheroidal particles with the same composition and fixed orientation. This proves for the first time that the statistical approach is generally viable in the Rayleigh domain.

In this paper we study the combined effects of size and shape of small solid state particles on the absorption, emission and scattering characteristics. We use the statistical approach to calculate these optical properties. In this approach the average optical properties of an ensemble of particles in random orientation are represented by the average optical properties of an ensemble of simple shapes. The validity of this approach is studied in detail for a uniform distribution of hollow spheres where the fractional volume of the central inclusion is varied. We apply the results to two different areas of interest, i) infrared spectroscopy and ii) polarization of scattered light. The effects of particle size and shape on the optical characteristics are discussed. We compare the results using the distribution of hollow spheres with those obtained by using randomly oriented spheroids. Also we compare the results with observations and laboratory measurements. The distribution of hollow spheres is very successful in reproducing laboratory measurements of the scattering angle distribution of the degree of linear polarization for incident unpolarized light of randomly oriented irregular quartz particles. Furthermore, we show that we are able to derive the size distribution of dust grains by fitting the measured degree of linear polarization using computational result for hollow spheres. It is shown that the distribution of hollow spheres is a powerful tool for studying light scattering, absorption and emission by cosmic dust grains and in particular when large numbers of particle parameters need to be considered since the computational demand of the distribution of hollow spheres is small.

The core of the nebula surrounding Eta Carinae has been observed with the VLT Adaptive Optics system NACO and with the interferometer VLTI/MIDI to constrain spatially and spectrally the warm dusty environment and the central object. In particular, narrow-band images at 3.74 and 4.05 micron reveal the butterfly shaped dusty environment close to the central star with unprecedented spatial resolution. A void whose radius corresponds to the expected sublimation radius has been discovered around the central source. Fringes have been obtained in the Mid-IR which reveal a correlated flux of about 100Jy situated 0.3" south-east of the photocenter of the nebula at 8.7 micron, which corresponds with the location of the star as seen in other wavelengths. This correlated flux is partly attributed to the central object, and these observations provide an upper limit for the SED of the central source from 2.2 to 13.5 micron. Moreover, we have been able to spectrally disperse the signal from the nebula itself at PA=318 degree, i.e. in the direction of the bipolar nebula 310 degree) within the MIDI field of view of 3". A large amount of corundum (Al2O3) is discovered, peaking at 0.6-1.2" south-east from the star, whereas the dust content of the Weigelt blobs is dominated b silicates. We discuss the mechanisms of dust formation which are closely related to the geometry of this Butterfly nebulae.

We examine the thermal infrared spectra of large dust grains of different chemical composition and mineralogy. Strong resonances in the optical properties result in detectable spectral structure even when the grain is much larger than the wavelength at which it radiates. We apply this to the thermal infrared spectra of compact amorphous and crystalline silicates. The weak resonances of amorphous silicates at 9.7 and 18 micron virtually disappear for grains larger than about 10 micron. In contrast, the strong resonances of crystalline silicates produce emission dips in the infrared spectra of large grains; these emission dips are shifted in wavelength compared to the emission peaks commonly seen in small crystalline silicate grains. We discuss the effect of a fluffy or compact grain structure on the infrared emission spectra of large grains, and apply our theory to the dust shell surrounding Vega.

We present the 2-200 um Infrared Space Observatory (ISO) spectrum and 3.8-20 um ISAAC and TIMMI2 images of the extreme OH/IR star IRAS 16342-3814. Amorphous silicate absorption features are seen, together with crystalline silicate absorption features up to almost 45 um. No other OH/IR star is known to have crystalline silicate features in absorption up to these wavelengths. This suggests that IRAS 16342-3814 must have, or recently had, an extremely high mass-loss rate. Preliminary radiative transfer calculations suggest that the mass-loss rate may be as large as 10^{-3} Msun/yr. The 3.8 um ISAAC image shows a bipolar reflection nebula with a dark equatorial waist or torus, similar to that seen in optical Hubble Space Telescope (HST) images. The position angle of the nebula decreases significantly with increasing wavelength, suggesting that the dominant source of emission changes from scattering to thermal emission. Still, even up to 20 um the nebula is oriented approximately along the major axis of the nebula seen in the HST and ISAAC images, suggesting that the torus must be very cold, in agreement with the very red ISO spectrum. The 20 um image shows a roughly spherically symmetric extended halo, approximately 6'' in diameter, which is probably due to a previous phase of mass-loss on the AGB, suggesting a transition from a (more) spherically symmetric to a (more) axial symmetric form of mass-loss at the end of the AGB. We estimate the maximum dust particle sizes in the torus and in the reflection nebula to be 1.3 and 0.09 um respectively. The size of the particles in the torus is large compared to typical ISM values, but in agreement with high mass-loss rate objects like AFGL 4106 and HD161796. We discuss the possible reason for the difference in particle size between the torus and the reflection nebula.

We present a model of the growth of water ice on silicate grains in the circumstellar envelopes of Asymptotic Giant Branch (AGB) stars and Red Super Giants. We consider the growth of ice by gas grain collisions, the thermal evaporation of ice from a grain, and sputtering. Our model contains several improvements compared to earlier models, including a detailed treatment of the effects of sputtering, a detailed calculation of the radiation pressure on the grain, and the treatment of subsonic drift velocities. In terms of drift velocity between the grains and gas in the envelope, we find that the ice formation process can be divided into three regimes: (i) a sputtering dominated regime where ice growth is heavily suppressed, (ii) an intermediate regime with moderately efficient condensation and (iii) a thermally dominated regime where ice formation is highly efficient. Sputtering is the critical factor which determines if ice formation can occur at all. We find that in Red Supergiants, ice formation is suppressed, while the winds of OH/IR stars allow for efficient condensation and can convert significant fractions of the available water vapor (tens of percent) into ice mantles on grains. Population II AGB stars hardly form ice due to their low dust to gas ratios. We also modify an analytical equation describing condensation and depletion (Jura and Morris, 1985) in order to give reasonable results for high and low drift velocities. Initially, ice will condense in crystalline form, but continuing condensation at low temperatures, and damage caused by interstellar UV photons favor the production of amorphous ice as well. We predict that a significant fraction of the ice formed will be amorphous.