EX Lup is a well-studied T Tauri star that represents the prototype of young eruptive stars EXors. In this paper we analyze new adaptive optics imaging and spectroscopic observations of EX Lup and its circumstellar environment in near-infrared in its quiescent phase. We observed EX Lup with the high contrast imager SPHERE/IRDIS in the dual-beam polarimetric imaging mode to resolve the circumstellar environment in near-infrared scattered light. We complemented these data with earlier SINFONI spectroscopy. We resolve for the first time in scattered light a compact feature around EX Lup azimuthally extending from 280deg to 360deg, and radially extending from 0.3arcsec to 0.55arcsec in the plane of the disk. We explore two different scenarios for the detected emission. One accounting for the emission as coming from the brightened walls of the cavity excavated by the outflow whose presence was suggested by ALMA observations in the J=3-2 line of 12CO. The other accounts for the emission as coming from an inclined disk. We detect for the first time a more extended circumstellar disk in scattered light, which shows that a region between 10 and 30 au is depleted of mum-size grains. We compare the J-, H- and K-band spectra obtained with SINFONI in quiescence with the spectra taken during the outburst, showing that all the emission lines were due to the episodic accretion event. Conclusions. Based on the morphology analysis we favour the scenario in which the scattered light is coming from a circumstellar disk rather than the outflow around EX Lup. We analyze the origin of the observed feature either as coming from a continuous circumstellar disk with a cavity, or from the illuminated wall of the outer disk or from a shadowed disk. Moreover, we discuss what is the origin of the mum-size grains depleted region, exploring the possibility that a sub-stellar companion may be the cause of it.

A large fraction of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks, and trace evidence for binary-disk interactions. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyze the observed disk morphology and substructures. We run 2D hydrodynamical models simulating the evolution of the circumbinary ring over the lifetime of the disk. The disk, as well as the cavity and the inner region are highly structured with several shadowed regions, spiral structures, and streamer-like filaments, some of them detected for the first time. The streamer-like filaments appear to connect the outer ring with the northern arc. Their azimuthal spacing suggests that they may be generated by periodic perturbations by the binary, tearing off material from the inner edge of the outer disk once during each orbit. By comparing observations to hydrodynamical simulations we find that the main features, in particular the gap size, as well as the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with semi-major axis of $\sim$35 au on an orbit coplanar with the circumbinary ring.

We present observations of the young multiple system UX Tauri to look for circumstellar disks and for signs of dynamical interactions. We obtained SPHERE/IRDIS deep differential polarization images in the J and H bands. We also used ALMA archival CO data. Large extended spirals are well detected in scattered light coming out of the disk of UX Tau A. The southern spiral forms a bridge between UX Tau A and C. These spirals, including the bridge connecting the two stars, all have a CO (3-2) counterpart seen by ALMA. The disk of UX Tau C is detected in scattered light. It is much smaller than the disk of UX Tau A and has a major axis along a different position angle, suggesting a misalignment. We performed PHANTOM SPH hydrodynamical models to interpret the data. The scattered light spirals, CO emission spirals and velocity patterns of the rotating disks, and the compactness of the disk of UX Tau C all point to a scenario in which UX Tau A has been perturbed very recently (about 1000 years) by the close passage of UX Tau C.

Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11{\deg}. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk.

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.

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

Shadows in scattered light images of protoplanetary disks are a common feature and support the presence of warps or misalignments between disk regions. These warps are possibly due to an inclined (sub-)stellar companion embedded in the disk. We study the morphology of the protoplanetary disk around the Herbig Ae star HD 139614 based on the first scattered light observations of this disk, which we model with the radiative transfer code MCMax3D. We obtained J- and H-band observations in polarized scattered light with VLT/SPHERE that show strong azimuthal asymmetries. In the outer disk, beyond ~30 au, a broad shadow spans a range of ~240{\deg} in position angle, in the East. A bright ring at ~16 au also shows an azimuthally asymmetric brightness, with the faintest side roughly coincidental with the brightest region of the outer disk. Additionally, two arcs are detected at ~34 au and ~50 au. We created a simple 4-zone approximation to a warped disk model of HD 139614 in order to qualitatively reproduce these features. The location and misalignment of the disk components were constrained from the shape and location of the shadows they cast. We find that the shadow on the outer disk covers a range of position angle too wide to be explained by a single inner misaligned component. Our model requires a minimum of two separate misaligned zones -- or a continuously warped region -- to cast this broad shadow on the outer disk. A small misalignment of ~4{\deg} between adjacent components can reproduce most of the observed shadow features. Multiple misaligned disk zones, potentially mimicing a warp, can explain the observed broad shadows in the HD 139614 disk. A planetary mass companion in the disk, located on an inclined orbit, could be responsible for such a feature and for the dust depleted gap responsible for a dip in the SED.

Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.

Volatile molecules are critical to habitability, yet difficult to observe directly at the optically thick midplanes of protoplanetary disks, where planets form. We present a new technique to measure indirectly the bulk atomic composition and distribution of solids that have been removed from the gas accreting onto the young star TW Hydrae. Our results suggest that disk gaps and planetesimal formation in TW Hya efficiently prevent radial drift and effectively block replenishment of the inner disk solids from ice reservoirs. Volatile-enhancement of giant planet atmospheres cannot be explained by primary accretion of this gas. Half of the missing oxygen is locked in silicates inside the snowline, while carbon is predominantly locked in organics beyond the CO snowline, leaving terrestrial planet-forming region solids both 'dry' and carbon-poor. Earth's depletion of carbon could be common amongst terrestrial exoplanets.

The presence of rocky exoplanets with a large refractory carbon inventory is predicted by chemical evolution models of protoplanetary disks of stars with photospheric C/O >0.65, and by models studying the radial transport of refractory carbon. High-pressure high-temperature laboratory experiments show that most of the carbon in these exoplanets differentiates into a graphite outer shell. Our aim is to evaluate the effects of a graphite outer shell on the thermal evolution of rocky exoplanets containing a metallic core and a silicate mantle. We implement a parameterized model of mantle convection to determine the thermal evolution of rocky exoplanets with graphite layer thicknesses up to 1000 km. We find that, due to the high thermal conductivity of graphite, conduction is the dominant heat transport mechanism in a graphite layer for the long-term evolution (>200 Myr). The conductive graphite shell essentially behaves like a stagnant lid with a fixed thickness. Models of Kepler-37b (Mercury-size) and a Mars-sized exoplanet show that a planet with a graphite lid cools faster than a planet with a silicate lid, and a planet without a stagnant lid cools the fastest. A graphite lid needs to be approximately ten times thicker than a corresponding silicate lid in order to produce similar thermal evolution.

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

Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light.

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.

Most major planetary bodies in the solar system rotate in the same direction as their orbital motion: their spin is prograde. Theoretical studies to explain the direction as well as the magnitude of the spin vector have had mixed success. When the accreting building blocks are $\sim$ km-size planetesimals -- as predicted by the classical model -- the accretion process is so symmetric that it cancels out prograde with retrograde spin contributions, rendering the net spin minute. For this reason, the currently-favored model for the origin of planetary rotation is the giant impact model, in which a single collision suffices to deliver a spin, which magnitude is close to the breakup rotation rate. However, the giant impact model does not naturally explain the preference for prograde spin. Similarly, an increasing number of spin-vector measurement of asteroids also shows that the spin vector of large (primordial) asteroids is not isotropic. Here, we re-assess the viability of smaller particles to bestow planetary bodies with a net spin, focusing on the pebble accretion model in which gas drag and gravity join forces to accrete small particles at a large cross section. Similar to the classical calculation for planetesimals, we integrate the pebble equation of motion and measure the angular momentum transfer at impact. We consider a variety of disk conditions and pebble properties and conduct our calculations in the limits of 2D (planar) and 3D (homogeneous) pebble distributions. We find that in certain regions of the parameter space the angular momentum transfer is significant, much larger than with planetesimals and on par with or exceeding the current spin of planetary bodies.

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.

The structure of cometary dust is a tracer of growth processes in the formation of planetesimals. Instrumentation on board the Rosetta mission to comet 67P/Churyumov- Gerasimenko captured dust particles and analysed them in situ. However, these deposits are a product of a collision within the instrument. We conducted laboratory experiments with cometary dust analogues, simulating the collection process by Rosetta instruments (specifically COSIMA, MIDAS). In Paper I we reported that velocity is a key driver in determining the appearance of deposits. Here in Paper II we use materials with different monomer sizes, and study the effect of tensile strength on the appearance of deposits. We find that mass transfer efficiency increases from $\sim$1 up to $\sim$10% with increasing monomer diameter from 0.3 $\mu$m to 1.5 $\mu$m (i.e. tensile strength decreasing from $\sim$12 to $\sim$3 kPa), and velocities increasing from 0.5 to 6 m/s. Also, the relative abundance of small fragments after impact is higher for material with higher tensile strength. The degeneracy between the effects of velocity and material strength may be lifted by performing a closer study of the deposits. This experimental method makes it possible to estimate the mass transfer efficiency in the COSIMA instrument. Extrapolating these results implies that more than half of the dust collected during the Rosetta mission has not been imaged. We analysed two COSIMA targets containing deposits from single collisions. The collision that occurred closest to perihelion passage led to more small fragments on the target.

Several detections of wide-orbit planet-mass/sub-stellar companions around young solar-like stars were reported in the last decade. The origin of those possible planets is still unclear but accretion tracers and VLT/SPHERE observations indicate that they are surrounded by circumplanetary material or even a circumplanetary disk. We want to investigate if the gas component of disks around wide-orbit companions is detectable with current and future (sub)mm telescopes and what constraints such gas observations can provide on the nature of the circumplanetary material and on the mass of the companion. We applied the radiation thermo-chemical disk code ProDiMo to model the dust and gas component of passive circumplanetary disks and produced realistic synthetic observables. We considered different companion properties, disk parameters and radiative environments and compared the resulting synthetic observables to telescope sensitivities and to existing dust observations. The main criterion for a successful detection is the size of the circumplanetary disk. At a distance of about 150 pc, a circumplanetary disk with an outer radius of about 10 au is detectable with ALMA in about 6 hours in optically thick CO lines. Other aspects such as the companion's luminosity, disk inclination and background radiation fields are also relevant and should be considered to optimize the observing strategy for detection experiments. For most of the known wide-orbit planet-mass companions, their maximum theoretical disk size of one third of the Hill radius would be sufficient to allow detection of CO lines. It is therefore feasible to detect their gas disks and constrain the mass of the companion through the kinematic signature. Even in the case of non-detections such observations will provide stringent constraints on disk size and gas mass, information crucial for formation theories.

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.

Aims: We want to detect and quantify observables related to accretion processes occurring locally in circumstellar disks, which could be attributed to young forming planets. We focus on objects known to host protoplanet candidates and/or disk structures thought to be the result of interactions with planets. Methods: We analyzed observations of 6 young stars (age $3.5-10$ Myr) and their surrounding environments with the SPHERE/ZIMPOL instrument on the VLT in the H$\alpha$ filter (656 nm) and a nearby continuum filter (644.9 nm). Results: We re-detect the known accreting M-star companion HD142527 B with the highest published signal to noise to date in both H$\alpha$ and the continuum. We derive new astrometry ($r = 62.8^{+2.1}_{-2.7}$ mas and $\text{PA} = (98.7\,\pm1.8)^\circ$) and photometry ($\Delta$N_Ha=$6.3^{+0.2}_{-0.3}$ mag, $\Delta$B_Ha=$6.7\pm0.2$ mag and $\Delta$Cnt_Ha=$7.3^{+0.3}_{-0.2}$ mag) for the companion in agreement with previous studies, and estimate its mass accretion rate ($\dot{M}\approx1-2\,\times10^{-10}\,M_\odot\text{ yr}^{-1}$). A faint point-like source around HD135344 B (SAO206462) is also investigated, but a second deeper observation is required to reveal its nature. No other companions are detected. In the framework of our assumptions we estimate detection limits at the locations of companion candidates around HD100546, HD169142 and MWC758 and calculate that processes involving H$\alpha$ fluxes larger than $\sim8\times10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}>10^{-10}-10^{-12}\,M_\odot\text{ yr}^{-1}$) can be excluded. Furthermore, flux upper limits of $\sim10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}<10^{-11}-10^{-12}\,M_\odot \text{ yr}^{-1}$) are estimated within the gaps identified in the disks surrounding HD135344B and TW Hya.

Carbon-enriched rocky exoplanets have been proposed around dwarf stars as well as around binary stars, white dwarfs and pulsars. However, the mineralogical make up of such planets is poorly constrained. We performed high-pressure high-temperature laboratory experiments ($P$ = 1$-$2 GPa, $T$ = 1523$-$1823 K) on carbon-enriched chemical mixtures to investigate the deep interiors of Pluto- to Mars-size planets the upper mantles of larger planets. Our results show that these exoplanets, when fully-differentiated, comprise a metallic core, a silicate mantle and a graphite layer on top of the silicate mantle. The silicate mineralogy (olivine, orthopyroxene, clinopyroxene and spinel) is largely unaffected by the amount of carbon. Metals are either two immiscible iron-rich alloys (S-rich and S-poor) or a single iron-rich alloy in the Fe-C-S system with immiscibility depending on the S/Fe ratio and core pressure. Graphite is the dominant carbon-bearing phase at the conditions of our experiments with no traces of silicon carbide or carbonates. If the bulk carbon content is higher than needed to saturate the mantle and the core, graphite would be in the form of an additional layer on top of the silicate mantle assuming differentiation. For a thick enough graphite layer, diamonds would form at the bottom of this layer due to high pressures. We model the interior structure of Kepler-37b and show that a mere 10 wt% graphite layer would decrease its derived mass by 7%, suggesting future space missions that determine both radius and mass of rocky exoplanets with insignificant gaseous envelopes could provide quantitative limits on their carbon content. Future observations of rocky exoplanets with graphite-rich surfaces would show low albedos due to the low reflectance of graphite. The absence of life-bearing elements other than carbon on the surface likely makes them uninhabitable.

Context. The mechanisms governing the opening of cavities in transition disks are not fully understood. Several processes have been proposed but their occurrence rate is still unknown. Aims. We present spatially resolved observations of two transition disks and aim at constraining their vertical and radial structure using multiwavelength observations. Methods. We have obtained near-IR scattered light observations with VLT/SPHERE of the transition disks J1608 and J1852. We complement our datasets with ALMA observations and with unresolved photometric observations covering a wide range of wavelengths. We performed radiative transfer modeling to analyze the morphology of the disks and compare the results with a sample of 20 other transition disks observed with both SPHERE and ALMA. Results. The scattered light image of J1608 reveals a very inclined disk, with two bright lobes and a large cavity. J1852 shows an inner ring extending beyond the coronagraphic radius up to 15au, a gap and a second ring at 42au. Our radiative transfer model of J1608 indicates that the millimeter-sized grains are less extended vertically and radially than the micron-sized grains, indicating advanced settling and radial drift. We find good agreement with the observations of J1852 with a similar model, but due to the low inclination of the system, the model remains partly degenerate. The analysis of 22 transition disks shows that, in general, the cavities observed in scattered light are smaller than the ones detected at millimeter wavelengths. Conclusions. The analysis of a sample of transition disks indicates that the small grains can flow inward of the region where millimeter grains are trapped. While 15 out of the 22 cavities in our sample could be explained by a planet of less than 13 Jupiter masses, the others either require the presence of a more massive companion or of several low-mass planets.

Sub-Neptune planets are a very common type of planets. They are inferred to harbour a primordial (H/He) envelope, on top of a (rocky) core, which dominates the mass. Here, we investigate the long-term consequences of the core properties on the planet mass-radius relation. We consider the role of various core energy sources resulting from core formation, its differentiation, its solidification (latent heat), core contraction and radioactive decay. We divide the evolution of the rocky core into three phases: the formation phase, which sets the initial conditions, the magma ocean phase, characterized by rapid heat transport, and the solid state phase, where cooling is inefficient. We find that for typical sub-Neptune planets of ~2-10 Earth masses and envelope mass fractions of 0.5-10% the magma ocean phase lasts several Gyrs, much longer than for terrestrial planets. The magma ocean phase effectively erases any signs of the initial core thermodynamic state. After solidification, the reduced heat flux from the rocky core causes a significant drop in the rocky core surface temperature, but its effect on the planet radius is limited. In the long run, radioactive heating is the most significant core energy source in our model. Overall, the long term radius uncertainty by core thermal effects is up to 15%.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

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

Low brightness dips have been recently observed in images of protoplanetary disks, and they are believed to be shadows by the inner disk. We present VLT/SPHERE polarimetric differential imaging of the transition disk around the dipper star RX J1604.3-2130. We gathered 11 epochs that cover a large temporal baseline, to search for variability over timescales of years, months, weeks, and days. Our observations unambiguously reveal two dips along an almost face-on narrow ring (with a width of $\sim$20 au), and the location of the peak of this ring is at $\sim$65 au. The ring lies inside the ring-like structure observed with ALMA, which peaks at $\sim$83 au. This segregation can result from particle trapping in pressure bumps, potentially due to planet(s). We find that the dips are variable, both in morphology and in position. The eastern dip, at a position angle (PA) of $\sim$83.7$\pm$13.7$^\circ$, has an amplitude that varies between 40% to 90%, and its angular width varies from 10$^{\circ}$ to 34$^{\circ}$. The western dip, at a PA of $\sim$265.90$\pm$13.0$^\circ$, is more variable, with amplitude and width variations of 31% to 95% and 12$^{\circ}$ to 53$^{\circ}$, respectively. The separation between the dips is 178.3$^{\circ}\pm$14.5$^{\circ}$, corresponding to a large misalignment between the inner and outer disk, supporting the classification of J1604 as an aperiodic dipper. The variability indicates that the innermost regions are highly dynamic, possibly due to a massive companion or to a complex magnetic field topology.

We present the first part of our DARTTS-S (Disks ARound TTauri Stars with SPHERE) survey: Observations of 8 TTauri stars which were selected based on their strong (sub-)mm excesses using SPHERE / IRDIS polarimetric differential imaging (PDI) in the J and H bands. All observations successfully detect the disks, which appear vastly different in size, from $\approx$80 au in scattered light to $>$400 au, and display total polarized disk fluxes between 0.06% and 0.89% of the stellar flux. For five of these disks, we are able to determine the three-dimensional structure and the flaring of the disk surface, which appears to be relatively consistent across the different disks, with flaring exponents $\alpha$ between $\approx$1.1 and $\approx$1.6. We also confirm literature results w.r.t. the inclination and position angle of several of our disk, and are able to determine which side is the near side of the disk in most cases. While there is a clear trend of disk mass with stellar ages ($\approx$1 Myr to $>$10 Myr), no correlations of disk structures with age were found. There are also no correlations with either stellar mass or sub-mm flux. We do not detect significant differences between the J and H bands. However, we note that while a high fraction (7/8) of the disks in our sample show ring-shaped sub-structures, none of them display spirals, in contrast to the disks around more massive Herbig Ae/Be stars, where spiral features are common.

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.

The variety of observed protoplanetary disks in polarimetric light motivates a taxonomical study to constrain their evolution and establish the current framework of this type of observations. We classified 58 disks with available polarimetric observations into six major categories (Ring, Spiral, Giant, Rim, Faint, and Small disks) based on their appearance in scattered light. We re-calculated the stellar and disk properties from the newly available GAIA DR2 and related these properties with the disk categories. More than a half of our sample shows disk sub-structures. For the remaining sources, the absence of detected features is due to their faintness, to their small size, or to the disk geometry. Faint disks are typically found around young stars and typically host no cavity. There is a possible dichotomy in the near-IR excess of sources with spiral-disks (high) and ring-disks (low). Like spirals, shadows are associated with a high near-IR excess. If we account for the pre-main sequence evolutionary timescale of stars with different mass, spiral arms are likely associated to old disks. We also found a loose, shallow declining trend for the disk dust mass with time. Protoplanetary disks may form sub-structures like rings very early in their evolution but their detectability in scattered light is limited to relatively old sources (more than 5 Myr) where the recurrently detected disk cavities allow to illuminate the outer disk. The shallow decrease of disk mass with time might be due to a selection effect, where disks observed thus far in scattered light are typically massive, bright transition disks with longer lifetime than most disks. Our study points toward spirals and shadows being generated by planets of fraction-to-few Jupiter masses that leave their (observed) imprint on both the inner disk near the star and the outer disk cavity.

Theoretical models predict the condensation of silicon carbide around host stars with C/O ratios higher than 0.65 (cf. C/O$_{\mathrm{Sun}}$ = 0.54), in addition to its observations in meteorites, interstellar medium and protoplanetary disks. Consequently, the interiors of rocky exoplanets born from carbon-enriched refractory material are often assumed to contain large amounts of silicon carbide. Here we aim to investigate the stability of silicon carbide in the interior of carbon-enriched rocky exoplanets and to derive the reaction leading to its transformation. We performed a high-pressure high-temperature experiment to investigate the reaction between a silicon carbide layer and a layer representative of the bulk composition of a carbon-enriched rocky exoplanet. We report the reaction leading to oxidation of silicon carbide producing quartz, graphite, and molten iron silicide. Combined with previous studies, we show that in order to stabilize silicon carbide, carbon saturation is not sufficient, and a complete reduction of Fe$^{2+}$ to Fe$^{0}$ in a planetary mantle is required, suggesting that future spectroscopic detection of Fe$^{2+}$ or Fe$^{3+}$ on the surface of rocky exoplanets would imply the absence of silicon carbide in their interiors.

The Earth is strongly depleted in carbon compared to the dust in the ISM, implying efficient removal of refractory carbon before parent body formation. It has been argued that grains get rid of their carbon through oxidation and photolysis in the exposed upper disk layers. We assess the efficacy of these C-removal mechanisms accounting for the vertical and radial transport of grains. We obtain the carbon and carbon free mass budget of solids by solving two 1D advection-diffusion equations, accounting for the dust grain size distribution and radial transport. The carbon removal acts on the fraction of the grains that are in the exposed layer and requires efficient vertical transport. In models without radial transport, oxidation and photolysis can destroy most of the refractory carbon in terrestrial planet formation region. But it only reaches the observed depletion levels for extreme parameter combinations and requires that parent body formation was delayed by 1 Myr. Adding radial transport of solids prevents the depletion entirely, leaving refractory carbon equally distributed throughout the disk. It is unlikely that the observed carbon depletion can ultimately be attributed to mechanisms operating on small grains in the disk surface layers. Other mechanisms need to be studied, for example flash heating events or FU Ori outbursts in order to remove carbon quickly and deeply. However, a sustained drift barrier or strongly reduced radial grain mobility are necessary to prevent replenishment of carbon from the outer disk.

While planet formation is thought to occur early in the history of a protoplanetary disk, the presence of planets embedded in disks, or of other processes driving disk evolution, might be traced from their imprints on the disk structure. We observed the T Tauri star HD 143006, located in the 5-11 Myr-old Upper Sco region, in polarized scattered light with VLT/SPHERE at near-infrared wavelengths, reaching an angular resolution of ~0.037" (~6 au). We obtained two datasets, one with a 145 mas diameter coronagraph, and the other without, enabling us to probe the disk structure down to an angular separation of ~0.06" (~10 au). In our observations, the disk of HD 143006 is clearly resolved up to ~0.5" and shows a clear large-scale asymmetry with the eastern side brighter than the western side. We detect a number of additional features, including two gaps and a ring. The ring shows an overbrightness at a position angle (PA) of ~140 deg, extending over a range in position angle of ~60 deg, and two narrow dark regions. The two narrow dark lanes and the overall large-scale asymmetry are indicative of shadowing effects, likely due to a misaligned inner disk. We demonstrate the remarkable resemblance between the scattered light image of HD 143006 and a model prediction of a warped disk due to an inclined binary companion. The warped disk model, based on the hydrodynamic simulations combined with 3D radiative transfer calculations, reproduces all major morphological features. However, it does not account for the observed overbrightness at PA~140 deg. Shadows have been detected in several protoplanetary disks, suggesting that misalignment in disks is not uncommon. However, the origin of the misalignment is not clear. As-yet-undetected stellar or massive planetary companions could be responsible for them, and naturally account for the presence of depleted inner cavities.

The nearby Herbig Be star HD100546 is known to be a laboratory for the study of protoplanets and their relation with the circumstellar disk that is carved by at least 2 gaps. We observed the HD100546 environment with high contrast imaging exploiting several different observing modes of SPHERE, including datasets with/without coronagraphs, dual band imaging, integral field spectroscopy and polarimetry. The picture emerging from these different data sets is complex. Flux-conservative algorithms images clearly show the disk up to 200au. More aggressive algorithms reveal several rings and warped arms overlapping the main disk. The bright parts of this ring lie at considerable height over the disk mid-plane at about 30au. Our images demonstrate that the brightest wings close to the star in the near side of the disk are a unique structure, corresponding to the outer edge of the intermediate disk at ~40au. Modeling of the scattered light from the disk with a geometrical algorithm reveals that a moderately thin structure can well reproduce the light distribution in the flux-conservative images. We suggest that the gap between 44 and 113 au span between the 1:2 and 3:2 resonance orbits of a massive body located at ~70au that might coincide with the candidate planet HD100546b detected with previous thermal IR observations. In this picture, the two wings can be the near side of a ring formed by disk material brought out of the disk at the 1:2 resonance with the same massive object. While we find no clear evidence confirming detection of the planet candidate HD100546c in our data, we find a diffuse emission close to the expected position of HD100546b. This source can be described as an extremely reddened substellar object surrounded by a dust cloud or its circumplanetary disk. Its astrometry is broadly consistent with a circular orbital motion on the disk plane.

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

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

The nearby and young M star AU Mic is surrounded by a debris disk in which we previously identified a series of large-scale arch-like structures that have never been seen before in any other debris disk and that move outward at high velocities. We initiated a monitoring program with the following objectives: 1) track the location of the structures and better constrain their projected speeds, 2) search for new features emerging closer in, and ultimately 3) understand the mechanism responsible for the motion and production of the disk features. AU Mic was observed at 11 different epochs between August 2014 and October 2017 with the IR camera and spectrograph of SPHERE. These high-contrast imaging data were processed with a variety of angular, spectral, and polarimetric differential imaging techniques to reveal the faintest structures in the disk. We measured the projected separations of the features in a systematic way for all epochs. We also applied the very same measurements to older observations from the Hubble Space Telescope (HST) with the visible cameras STIS and ACS. The main outcomes of this work are 1) the recovery of the five southeastern broad arch-like structures we identified in our first study, and confirmation of their fast motion (projected speed in the range 4-12 km/s); 2) the confirmation that the very first structures observed in 2004 with ACS are indeed connected to those observed later with STIS and now SPHERE; 3) the discovery of two new very compact structures at the northwest side of the disk (at 0.40" and 0.55" in May 2015) that move to the southeast at low speed; and 4) the identification of a new arch-like structure that might be emerging at the southeast side at about 0.4" from the star (as of May 2016). Abridged.

We report the discovery of a bright, brown dwarf companion to the star HIP 64892, imaged with VLT/SPHERE during the SHINE exoplanet survey. The host is a B9.5V member of the Lower-Centaurus-Crux subgroup of the Scorpius Centaurus OB association. The measured angular separation of the companion ($1.2705\pm0.0023$") corresponds to a projected distance of $159\pm12$ AU. We observed the target with the dual-band imaging and long-slit spectroscopy modes of the IRDIS imager to obtain its SED and astrometry. In addition, we reprocessed archival NACO L-band data, from which we also recover the companion. Its SED is consistent with a young (<30 Myr), low surface gravity object with a spectral type of M9$_{\gamma}\pm1$. From comparison with the BT-Settl atmospheric models we estimate an effective temperature of $T_{\textrm{eff}}=2600 \pm 100$ K, and comparison of the companion photometry to the COND evolutionary models yields a mass of $\sim29-37$ M$_{\text{J}}$ at the estimated age of $16^{+15}_{-7}$ Myr for the system. HIP 64892 is a rare example of an extreme-mass ratio system ($q\sim0.01$) and will be useful for testing models relating to the formation and evolution of such low-mass objects.

Transition disks are considered sites of ongoing planet formation, and their dust and gas distributions could be signposts of embedded planets. The transition disk around the T Tauri star RY Lup has an inner dust cavity and displays a strong silicate emission feature. Using high-resolution imaging we study the disk geometry, including non-axisymmetric features, and its surface dust grain, to gain a better understanding of the disk evolutionary process. Moreover, we search for companion candidates, possibly connected to the disk. We obtained high-contrast and high angular resolution data in the near-infrared with the VLT/SPHERE extreme adaptive optics instrument whose goal is to study the planet formation by detecting and characterizing these planets and their formation environments through direct imaging. We performed polarimetric imaging of the RY~Lup disk with IRDIS (at 1.6 microns), and obtained intensity images with the IRDIS dual-band imaging camera simultaneously with the IFS spectro-imager (0.9-1.3 microns). We resolved for the first time the scattered light from the nearly edge-on circumstellar disk around RY~Lup, at projected separations in the 100 \,au range. The shape of the disk and its sharp features are clearly detectable at wavelengths ranging from 0.9 to 1.6 microns. We show that the observed morphology can be interpreted as spiral arms in the disk. This interpretation is supported by in-depth numerical simulations. We also demonstrate that these features can be produced by one planet interacting with the disk. We also detect several point sources which are classified as probable background objects.

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.

Sub-Neptune planets are very common in our galaxy and show a large diversity in their mass-radius relation. In sub-Neptunes most of the planet mass is in the rocky part (hereafter core) which is surrounded by a modest hydrogen-helium envelope. As a result, the total initial heat content of such a planet is dominated by that of the core. Nonetheless, most studies contend that the core cooling will only have a minor effect on the radius evolution of the gaseous envelope, because the core's cooling is in sync with the envelope, i.e., most of the initial heat is released early on timescales of about 10-100 Myr. In this Letter we examine the importance of the core cooling rate for the thermal evolution of the envelope. Thus, we relax the early core cooling assumption and present a model where the core is characterized by two parameters: the initial temperature and the cooling time. We find that core cooling can significantly enhance the radius of the planet when it operates on a timescale similar to the observed age, i.e. several Gyr. Consequently, the interpretation of sub-Neptunes' mass-radius observations depends on the assumed core thermal properties and the uncertainty therein. The degeneracy of composition and core thermal properties can be reduced by obtaining better estimates of the planet ages (in addition to their radii and masses) as envisioned by future observations.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

HD$\,$169142 is an excellent target to investigate signs of planet-disk interaction due to the previous evidence of gap structures. We performed J-band (~1.2{\mu}m) polarized intensity imaging of HD169142 with VLT/SPHERE. We observe polarized scattered light down to 0.16" (~19 au) and find an inner gap with a significantly reduced scattered light flux. We confirm the previously detected double ring structure peaking at 0.18" (~21 au) and 0.56" (~66 au), and marginally detect a faint third gap at 0.70"-0.73" (~82-85 au). We explore dust evolution models in a disk perturbed by two giant planets, as well as models with a parameterized dust size distribution. The dust evolution model is able to reproduce the ring locations and gap widths in polarized intensity, but fails to reproduce their depths. It, however, gives a good match with the ALMA dust continuum image at 1.3 mm. Models with a parameterized dust size distribution better reproduce the gap depth in scattered light, suggesting that dust filtration at the outer edges of the gaps is less effective. The pile-up of millimeter grains in a dust trap and the continuous distribution of small grains throughout the gap likely require a more efficient dust fragmentation and dust diffusion in the dust trap. Alternatively, turbulence or charging effects might lead to a reservoir of small grains at the surface layer that is not affected by the dust growth and fragmentation cycle dominating the dense disk midplane. The exploration of models shows that extracting planet properties such as mass from observed gap profiles is highly degenerate.

Spatially resolving the immediate surroundings of young stars is a key challenge for the planet formation community. SPHERE on the VLT represents an important step forward by increasing the opportunities offered by optical or near-infrared imaging instruments to image protoplanetary discs. The Guaranteed Time Observation Disc team has concentrated much of its efforts on polarimetric differential imaging, a technique that enables the efficient removal of stellar light and thus facilitates the detection of light scattered by the disc within a few au from the central star. These images reveal intriguing complex disc structures and diverse morphological features that are possibly caused by ongoing planet formation in the disc. An overview of the recent advances enabled by SPHERE is presented.

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.

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

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.

The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErschel Surveys of Star forming regions (CHESS), a very low HCO$^{+}$/N$_{2}$H$^{+}$ abundance ratio of about 3-4, has been reported toward the protocluster OMC-2 FIR4. This finding suggests a cosmic-ray ionization rate in excess of 10$^{-14}$ s$^{-1}$, much higher than the canonical value of $\zeta$ = 3$\times$10$^{-17}$ s$^{-1}$ (value expected in quiescent dense clouds). To assess the specificity of OMC-2 FIR4, we have extended this study to a sample of sources in low- and intermediate mass. More specifically, we seek to measure the HCO$^{+}$/N$_2$H$^{+}$ abundance ratio from high energy lines (J $\ge$ 6) toward this source sample in order to infer the flux of energetic particles in the warm and dense gas surrounding the protostars. We use observations performed with the Heterodyne Instrument for the FarInfrared spectrometer on board the Herschel Space Observatory toward a sample of 9 protostars. We report HCO$^{+}$/N$_2$H$^{+}$ abundance ratios in the range of 5 up to 73 toward our source sample. The large error bars do not allow us to conclude whether OMC-2~FIR4 is a peculiar source. Nonetheless, an important result is that the measured HCO$^{+}$/N$_2$H$^{+}$ ratio does not vary with the source luminosity. At the present time, OMC-2 FIR4 remains the only source where a high flux of energetic particles is clearly evident. More sensitive and higher angular resolution observations are required to further investigate this process.

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.

High-resolution imaging reveals a large morphological variety of protoplanetary disks. To date, no constraints on their global evolution have been found from this census. An evolutionary classification of disks was proposed based on their IR spectral energy distribution, with the Group I sources showing a prominent cold component ascribed to an earlier stage of evolution than Group II. Disk evolution can be constrained from the comparison of disks with different properties. A first attempt of disk taxonomy is now possible thanks to the increasing number of high-resolution images of Herbig Ae/Be stars becoming available. Near-IR images of six Group II disks in scattered light were obtained with VLT/NACO in Polarimetric Differential Imaging, which is the most efficient technique to image the light scattered by the disk material close to the stars. We compare the stellar/disk properties of this sample with those of well-studied Group I sources available from the literature. Three Group II disks are detected. The brightness distribution in the disk of HD163296 indicates the presence of a persistent ring-like structure with a possible connection with the CO snowline. A rather compact (less than 100 AU) disk is detected around HD142666 and AK Sco. A taxonomic analysis of 17 Herbig Ae/Be sources reveals that the difference between Group I and Group II is due to the presence or absence of a large disk cavity (larger than 5 AU). There is no evidence supporting the evolution from Group I to Group II. Group II are not evolved version of the Group I. Within the Group II disks, very different geometries (both self-shadowed and compact) exist. HD163296 could be the primordial version of a typical Group I. Other Group II, like AK Sco and HD142666, could be smaller counterpart of Group I unable to open cavities as large as those of Group I.

We have observed the protoplanetary disk of the well-known young Herbig star HD 142527 using ZIMPOL Polarimetric Differential Imaging with the VBB (Very Broad Band, ~600-900nm) filter. We obtained two datasets in May 2015 and March 2016. Our data allow us to explore dust scattering around the star down to a radius of ~0.025" (~4au). The well-known outer disk is clearly detected, at higher resolution than before, and shows previously unknown sub-structures, including spirals going inwards into the cavity. Close to the star, dust scattering is detected at high signal-to-noise ratio, but it is unclear whether the signal represents the inner disk, which has been linked to the two prominent local minima in the scattering of the outer disk, interpreted as shadows. An interpretation of an inclined inner disk combined with a dust halo is compatible with both our and previous observations, but other arrangements of the dust cannot be ruled out. Dust scattering is also present within the large gap between ~30 and ~140au. The comparison of the two datasets suggests rapid evolution of the inner regions of the disk, potentially driven by the interaction with the close-in M-dwarf companion, around which no polarimetric signal is detected.

The transition disk around the T Tauri star T Cha possesses a large gap, making it a prime target for high-resolution imaging in the context of planet formation. We aim to find signs of disk evolutionary processes by studying the disk geometry and the dust grain properties at its surface, and to search for companion candidates. We analyze a set of VLT/SPHERE data at near-infrared and optical wavelengths. We performed polarimetric imaging of T Cha with IRDIS (1.6 $\mu$m) and ZIMPOL (0.5-0.9 $\mu$m), and obtained intensity images from IRDIS dual-band imaging with simultaneous spectro-imaging with IFS (0.9-1.3 $\mu$m). The disk around T Cha is detected in all observing modes and its outer disk is resolved in scattered light with unprecedented angular resolution and signal-to-noise. The images reveal a highly inclined disk with a noticeable east-west brightness asymmetry. The significant amount of non-azimuthal polarization signal in the $U_{\phi}$ images, with a $U_{\phi}$/$Q_{\phi}$ peak-to-peak value of 14%, is in accordance with theoretical studies on multiple scattering. Our optimal axisymmetric radiative transfer model considers two coplanar inner and outer disks, separated by a gap of 0.28" (~30au) in size. We derive a disk inclination of ~69 deg and PA of ~114 deg. In order to self-consistently reproduce the intensity and polarimetric images, the dust grains, responsible for the scattered light, need to be dominated by sizes of around ten microns. A point source is detected at an angular distance of 3.5" from the central star. It is, however, found not to be co-moving. We confirm that the dominant source of emission is forward scattered light from the near edge of the outer disk. Our point source analysis rules out the presence of a companion heavier than ~8.5 $M_{\mathrm{jup}}$ between 0.1" and 0.3". The detection limit decreases to ~2 $M_{\mathrm{jup}}$ for 0.3" to 4.0".

Cometary dust provides a unique window on dust growth mechanisms during the onset of planet formation. Measurements by the Rosetta spacecraft show that the dust in the coma of comet 67P/Churyumov-Gerasimenko has a granular structure at size scales from sub-um up to several hundreds of um, indicating hierarchical growth took place across these size scales. However, these dust particles may have been modified during their collection by the spacecraft instruments. Here we present the results of laboratory experiments that simulate the impact of dust on the collection surfaces of COSIMA and MIDAS, instruments onboard the Rosetta spacecraft. We map the size and structure of the footprints left by the dust particles as a function of their initial size (up to several hundred um) and velocity (up to 6 m/s). We find that in most collisions, only part of the dust particle is left on the target; velocity is the main driver of the appearance of these deposits. A boundary between sticking/bouncing and fragmentation as an outcome of the particle-target collision is found at v ~ 2 m/s. For velocities below this value, particles either stick and leave a single deposit on the target plate, or bounce, leaving a shallow footprint of monomers. At velocities > 2 m/s and sizes > 80 um, particles fragment upon collision, transferring up to 50 per cent of their mass in a rubble-pile-like deposit on the target plate. The amount of mass transferred increases with the impact velocity. The morphologies of the deposits are qualitatively similar to those found by the COSIMA instrument.

Context. The SAO 206462 (HD 135344B) disk is one of the few known transitional disks showing asymmetric features in scattered light and thermal emission. Near-infrared scattered-light images revealed two bright outer spiral arms and an inner cavity depleted in dust. Giant protoplanets have been proposed to account for the disk morphology. Aims. We aim to search for giant planets responsible for the disk features and, in the case of non-detection, to constrain recent planet predictions using the data detection limits. Methods. We obtained new high-contrast and high-resolution total intensity images of the target spanning the Y to the K bands (0.95-2.3 mic) using the VLT/SPHERE near-infrared camera and integral field spectrometer. Results. The spiral arms and the outer cavity edge are revealed at high resolutions and sensitivities without the need for image post-processing techniques, which introduce photometric biases. We do not detect any close-in companions. For the derivation of the detection limits on putative giant planets embedded in the disk, we show that the knowledge of the disk aspect ratio and viscosity is critical for the estimation of the attenuation of a planet signal by the protoplanetary dust because of the gaps that these putative planets may open. Given assumptions on these parameters, the mass limits can vary from ~2-5 to ~4-7 Jupiter masses at separations beyond the disk spiral arms. The SPHERE detection limits are more stringent than those derived from archival NaCo/L' data and provide new constraints on a few recent predictions of massive planets (4-15 MJ) based on the spiral density wave theory. The SPHERE and ALMA data do not favor the hypotheses on massive giant planets in the outer disk (beyond 0.6). There could still be low-mass planets in the outer disk and/or planets inside the cavity.

We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are generally much lower than predictions of thermo-chemical models with canonical inputs. One ground-state transition is newly detected in the stacked spectrum of AA Tau, DM Tau, LkCa 15, and MWC 480. We run a grid of models to show that the abundance of gas-phase oxygen needs to be reduced by a factor of at least ~100 to be consistent with the observational upper limits (and positive detections) if a dust-to-gas mass ratio of 0.01 were to be assumed. As a continuation of previous ideas, we propose that the underlying reason for the depletion of oxygen (hence the low detection rate) is the freeze-out of volatiles such as water and CO onto dust grains followed by grain growth and settling/migration, which permanently removes these gas-phase molecules from the emissive upper layers of the outer disk. Such depletion of volatiles is likely ubiquitous among different disks, though not necessarily to the same degree. The volatiles might be returned back to the gas phase in the inner disk (within about 15 AU), which is consistent with current constraints. Comparison with studies on disk dispersal due to photoevaporation indicates that the timescale for volatile depletion is shorter than that of photoevaporation.

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

Context. Ice lines are suggested to play a significant role in grain growth and planetesimal formation in protoplanetary disks. Evaporation fronts directly influence the gas and ice abundances of volatile species in the disk and therefore the coagulation physics and efficiency and the chemical composition of the resulting planetesimals. Aims. In this work we investigate the influence of the existence of the CO ice line on the particle growth and on the distribution of CO in the disk. Methods. We include the possibility of tracking the CO content and/or other volatiles in particles and in the gas in our existing dust coagulation and disk evolution model and developed a method for evaporation and condensation of CO using the Hertz-Knudsen equation. Our model does not include fragmentation, yet, which will be part of further investigations. Results. We find no enhanced grain growth just outside the ice line where the particle size is limited by radial drift. Instead we find a depletion of solid material inside the ice line which is solely due to evaporation of the CO. Such a depression inside the ice line may be observable and may help to quantify the processes described in this work. Furthermore, we find that the viscosity and diffusivity of the gas heavily influence the re-distribution of vaporized CO at the ice line and can lead to an increase in the CO abundance by up to a factors of a few in the region just inside the ice line. Depending on the strength of the gaseous transport mechanisms the position of the ice line in our model can change by up to 10 AU and consequently, the temperature at that location can range from 21 K to 23 K.

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

Dust coagulation in interstellar space and protoplanetary disks is usually treated as one of 2 extreme cases: Particle-Cluster Aggregation and Cluster-Cluster Aggregation. In this paper we study the process of hierarchical growth, where aggregates are built from significantly smaller aggregates (but not monomers). We show that this process can be understood as a modified, PCA-like process that produces porous, but non-fractal particles whose filling factor is chiefly determined by the porosity of the building blocks. We also show that in a coagulation environment where relative velocities are driven by turbulence, a logarithmically flat mass distribution (equal mass per mass decade) as it is typically found in environments where fragmentation replenishes small grains, leads to a situation where small particles and aggregates dominate the growth of large ones. Therefore, in such environments, hierarchical growth should be seen as the norm. Consequently, we predict that the aggregates in such environments are not fractals with extremely low densities as they would result from extrapolation fractal laws to large sizes. The compactification of aggregates does not only result from collisions with enough energy to restructure aggregates - it starts already earlier by filling voids in particles with smaller grains that contribute to the growth.

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.

We search for signs of ongoing planet-disk interaction and study the distribution of small grains at the surface of the transition disk around RXJ1615.3-3255 (RX J1615). We observed RXJ1615 with VLT/SPHERE. We image the disk for the first time in scattered light and detect two arcs, two rings, a gap and an inner disk with marginal evidence for an inner cavity. The shapes of the arcs suggest that they probably are segments of full rings. Ellipse fitting for the two rings and inner disk yield a disk inclination i = 47 \pm 2 degrees and find semi-major axes of 1.50 \pm 0.01" (278 au), 1.06 \pm 0.01" (196 au) and 0.30 \pm 0.01" (56 au), respectively. We determine the scattering surface height above the midplane, based on the projected ring center offsets. Nine point sources are detected between 2.1" and 8.0" separation and considered as companion candidates. With NACO data we recover four of the nine point sources, which we determine not to be co-moving, and therefore unbound to the system. We present the first detection of the transition disk of RXJ1615 in scattered light. The height of the rings indicate limited flaring of the disk surface, which enables partial self-shadowing in the disk. The outermost arc either traces the bottom of the disk or it is another ring with semi-major axis > 2.35" (435 au). We explore both scenarios, extrapolating the complete shape of the feature, which will allow to distinguish between the two in future observations. The most interesting scenario, where the arc traces the bottom of the outer ring, requires the disk truncated at r ~ 360 au. The closest companion candidate, if indeed orbiting the disk at 540 au, would then be the most likely cause for such truncation. This companion candidate, as well as the remaining four, require follow up observations to determine if they are bound to the system.

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 present an overview of the scientific potential of MATISSE, the Multi Aperture mid-Infrared SpectroScopic Experiment for the Very Large Telescope Interferometer. For this purpose we outline selected case studies from various areas, such as star and planet formation, active galactic nuclei, evolved stars, extrasolar planets, and solar system minor bodies and discuss strategies for the planning and analysis of future MATISSE observations. Moreover, the importance of MATISSE observations in the context of complementary high-angular resolution observations at near-infrared and submillimeter/millimeter wavelengths is highlighted.

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.

LkCa 15 hosts a pre-transitional disk as well as at least one accreting protoplanet orbiting in its gap. Previous disk observations have focused mainly on the outer disk, which is cleared inward of ~50 au. The planet candidates, on the other hand, reside at orbital radii around 15 au, where disk observations have been unreliable until recently. Here we present new J-band imaging polarimetry of LkCa 15 with SPHERE IRDIS, yielding the most accurate and detailed scattered-light images of the disk to date down to the planet-hosting inner regions. We find what appear to be persistent asymmetric structures in the scattering material at the location of the planet candidates, which could be responsible at least for parts of the signals measured with sparse-aperture masking. These images further allow us to trace the gap edge in scattered light at all position angles and search the inner and outer disks for morphological substructure. The outer disk appears smooth with slight azimuthal variations in polarized surface brightness, which may be due to shadowing from the inner disk or a two-peaked polarized phase function. We find that the near-side gap edge revealed by polarimetry matches the sharp crescent seen in previous ADI imaging very well. Finally, the ratio of polarized disk to stellar flux is more than six times larger in J-band than in the RI bands.

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.

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances by opening new avenues in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 um, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with up to four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we present one of the main science objectives, the study of protoplanetary disks, that has driven the instrument design and motivated several VLTI upgrades (GRA4MAT and NAOMI). We introduce the physical concept of MATISSE including a description of the signal on the detectors and an evaluation of the expected performances. We also discuss the current status of the MATISSE instrument, which is entering its testing phase, and the foreseen schedule for the next two years that will lead to the first light at Paranal.

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.

Debris disks offer valuable insights into the latest stages of circumstellar disk evolution, and can possibly help us to trace the outcomes of planetary formation processes. In the age range 10 to 100\,Myr, most of the gas is expected to have been removed from the system, giant planets (if any) must have already been formed, and the formation of terrestrial planets may be on-going. Pluto-sized planetesimals, and their debris released in a collisional cascade, are under their mutual gravitational influence, which may result into non-axisymmetric structures in the debris disk. High angular resolution observations are required to investigate these effects and constrain the dynamical evolution of debris disks. Furthermore, multi-wavelength observations can provide information about the dust dynamics by probing different grain sizes. Here we present new VLT/SPHERE and ALMA observations of the debris disk around the 40\,Myr-old solar-type star HD\,61005. We resolve the disk at unprecedented resolution both in the near-infrared (in scattered and polarized light) and at millimeter wavelengths. Thanks to the new observations, we propose a solution for both the radial and azimuthal distribution of the dust grains in the debris disk. We find that the disk has a moderate eccentricity ($e \sim 0.1$) and that the dust density is two times larger at the pericenter compared to the apocenter. With no giant planets detected in our observations, we investigate alternative explanations besides planet-disk interactions to interpret the inferred disk morphology. We postulate that the morphology of the disk could be the consequence of a massive collision between $\sim$\,1000\,km-sized bodies at $\sim$\,61\,au. If this interpretation holds, it would put stringent constraints on the formation of massive planetesimals at large distances from the star.

The young star beta Pictoris is well known for its dusty debris disk, produced through the grinding down by collisions of planetesimals, kilometre-sized bodies in orbit around the star. In addition to dust, small amounts of gas are also known to orbit the star, likely the result from vaporisation of violently colliding dust grains. The disk is seen edge on and from previous absorption spectroscopy we know that the gas is very rich in carbon relative to other elements. The oxygen content has been more difficult to assess, however, with early estimates finding very little oxygen in the gas at a C/O ratio 20x higher than the cosmic value. A C/O ratio that high is difficult to explain and would have far-reaching consequences for planet formation. Here we report on observations by the far-infrared space telescope Herschel, using PACS, of emission lines from ionised carbon and neutral oxygen. The detected emission from C+ is consistent with that previously reported being observed by the HIFI instrument on Herschel, while the emission from O is hard to explain without assuming a higher-density region in the disk, perhaps in the shape of a clump or a dense torus, required to sufficiently excite the O atoms. A possible scenario is that the C/O gas is produced by the same process responsible for the CO clump recently observed by ALMA in the disk, and that the re-distribution of the gas takes longer than previously assumed. A more detailed estimate of the C/O ratio and the mass of O will have to await better constraints on the C/O gas spatial distribution.

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.

We image with unprecedented spatial resolution and sensitivity disk features that could be potential signs of planet-disk interaction. Two companion candidates have been claimed in the disk around the young Herbig Ae/Be star HD100546. Thus, this object serves as an excellent target for our investigation of the natal environment of giant planets. We exploit the power of extreme adaptive optics operating in conjunction with the new high-contrast imager SPHERE to image HD100546 in scattered light. We obtain the first polarized light observations of this source in the visible (with resolution as fine as 2 AU) and new H and K band total intensity images that we analyze with the Pynpoint package. The disk shows a complex azimuthal morphology, where multiple scattering of photons most likely plays an important role. High brightness contrasts and arm-like structures are ubiquitous in the disk. A double-wing structure (partly due to ADI processing) resembles a morphology newly observed in inclined disks. Given the cavity size in the visible (11 AU), the CO emission associated to the planet candidate 'c' might arise from within the circumstellar disk. We find an extended emission in the K band at the expected location of 'b'. The surrounding large-scale region is the brightest in scattered light. There is no sign of any disk gap associated to 'b'.

The journey from dust particle to planetesimal involves physical processes acting on scales ranging from micrometers (the sticking and restructuring of aggregates) to hundreds of astronomical units (the size of the turbulent protoplanetary nebula). Considering these processes simultaneously is essential when studying planetesimal formation. We develop a novel, global, semi-analytical model for the evolution of the mass-dominating dust particles in a turbulent protoplanetary disk that takes into account the evolution of the dust surface density while preserving the essential characteristics of the porous coagulation process. This panoptic model is used to study the growth from submicron to planetesimal sizes in disks around Sun-like stars. For highly porous ices, unaffected by collisional fragmentation and erosion, rapid growth to planetesimal sizes is possible in a zone stretching out to ${\sim}10\mathrm{~AU}$ for massive disks. When porous coagulation is limited by erosive collisions, the formation of planetesimals through direct coagulation is not possible, but the creation of a large population of aggregates with Stokes numbers close to unity might trigger the streaming instability (SI). However, we find that reaching conditions necessary for SI is difficult and limited to dust-rich disks, (very) cold disks, or disks with weak turbulence. Behind the snow-line, porosity-driven aggregation of icy grains results in rapid (${\sim}10^{4}\mathrm{~yr}$) formation of planetesimals. If erosive collisions prevent this, SI might be triggered for specific disk conditions. The numerical approach introduced in this work is ideally suited for studying planetesimal formation and pebble delivery simultaneously and will help build a coherent picture of the start of the planet formation process.

We combined hydrodynamical simulations of planet-disk interactions with dust evolution models that include coagulation and fragmentation of dust grains over a large range of radii and derived observational properties using radiative transfer calculations. We studied the role of the snow line in the survival of the inner disk of transition disks. Inside the snow line, the lack of ice mantles in dust particles decreases the sticking efficiency between grains. As a consequence, particles fragment at lower collision velocities than in regions beyond the snow line. This effect allows small particles to be maintained for up to a few Myrs within the first astronomical unit. These particles are closely coupled to the gas and do not drift significantly with respect to the gas. For lower mass planets (1$M_{\rm{Jup}}$), the pre-transition appearance can be maintained even longer because dust still trickles through the gap created by the planet, moves invisibly and quickly in the form of relatively large grains through the gap, and becomes visible again as it fragments and gets slowed down inside of the snow line. The global study of dust evolution of a disk with an embedded planet, including the changes of the dust aerodynamics near the snow line, can explain the concentration of millimetre-sized particles in the outer disk and the survival of the dust in the inner disk if a large dust trap is present in the outer disk. This behaviour solves the conundrum of the combination of both near-infrared excess and ring-like millimetre emission observed in several transition disks.

[Abridged] Context. The young systems PZ Tel and HD 1160, hosting known low-mass companions, were observed during the commissioning of the new planet finder SPHERE with several imaging and spectroscopic modes. Aims. We aim to refine the physical properties and architecture of both systems. Methods. We use SPHERE commissioning data and REM observations, as well as literature and unpublished data from VLT/SINFONI, VLT/NaCo, Gemini/NICI, and Keck/NIRC2. Results. We derive new photometry and confirm the nearly daily photometric variability of PZ Tel A. Using literature data spanning 38 yr, we show that the star also exhibits a long-term variability trend. The 0.63-3.8 mic SED of PZ Tel B allows us to revise its properties: spectral type M7+/-1, Teff=2700+/-100 K, log(g)<4.5 dex, log(L/L_Sun)=-2.51+/-0.10 dex, and mass 38-72 MJ. The 1-3.8 mic SED of HD 1160 B suggests a massive brown dwarf or a low-mass star with spectral type M5.5-7.0, Teff=3000+/-100 K, [M/H]=-0.5-0.0 dex, log(L/L_Sun)=-2.81+/-0.10 dex, and mass 39-168 MJ. We confirm the deceleration and high eccentricity (e>0.66) of PZ Tel B. For e<0.9, the inclination, longitude of the ascending node, and time of periastron passage are well constrained. The system is seen close to an edge-on geometry. We reject other brown dwarf candidates outside 0.25" for both systems, and massive giant planets (>4 MJ) outside 0.5" for the PZ Tel system. We also show that K1-K2 color can be used with YJH low-resolution spectra to identify young L-type companions, provided high photometric accuracy (<0.05 mag) is achieved. Conclusions. SPHERE opens new horizons in the study of young brown dwarfs and giant exoplanets thanks to high-contrast imaging capabilities at optical and near-infrared wavelengths, as well as high signal-to-noise spectroscopy in the near-infrared from low (R~30-50) to medium resolutions (R~350).

The system of four planets around HR8799 offers a unique opportunity to probe the physics and chemistry at play in the atmospheres of self-luminous young (~30 Myr) planets. We recently obtained new photometry of the four planets and low-resolution (R~30) spectra of HR8799 d and e with the SPHERE instrument (paper III). In this paper (paper IV), we compare the available spectra and photometry of the planets to known objects and atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to characterize the atmospheric properties of the planets. We find that HR8799d and e properties are well reproduced by those of L6-L8 dusty dwarfs discovered in the field, among which some are candidate members of young nearby associations. No known object reproduces well the properties of planets b and c. Nevertheless, we find that the spectra and WISE photometry of peculiar and/or young early-T dwarfs reddened by submicron grains made of corundum, iron, enstatite, or forsterite successfully reproduce the SED of these two planets. Our analysis confirms that only the Exo-REM models with thick clouds fit (within 2{\sigma}) the whole set of spectrophotometric datapoints available for HR8799 d and e for Teff = 1200 K, log g in the range 3.0-4.5, and M/H=+0.5. The models still fail to reproduce the SED of HR8799c and b. The determination of the metallicity, log g, and cloud thickness are degenerate. We conclude that an enhanced content in dust and decreased CIA of H2 is certainly responsible for the deviation of the properties of the planet with respect to field dwarfs. The analysis suggests in addition that HR8799c and b have later spectral types than the two other planets, and therefore could both have lower masses.

GJ758 B is a brown dwarf companion to a nearby (15.76 pc) solar-type, metal-rich (M/H = +0.2 dex) main-sequence star (G9V) that was discovered with Subaru/HiCIAO in 2009. From previous studies, it has drawn attention as being the coldest (~600K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the VLT. The data was obtained in Y-, J-, H-, and Ks-bands with the dual-band imaging (DBI) mode of IRDIS, providing a broad coverage of the full near-infrared (near-IR) range at higher contrast and better spectral sampling than previously reported. In this new set of high-quality data, we report the re-detection of the companion, as well as the first detection of a new candidate closer-in to the star. We use the new 8 photometric points for an extended comparison of GJ758 B with empirical objects and 4 families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison object can accurately represent the observed near-IR fluxes of GJ758 B. From comparison to atmospheric models, we attribute a Teff = 600K $\pm$ 100K, but we find that no atmospheric model can adequately fit all the fluxes of GJ758 B. The photometry of the new candidate companion is broadly consistent with L-type objects, but a second epoch with improved photometry is necessary to clarify its status. The new astrometry of GJ758 B shows a significant proper motion since the last epoch. We use this result to improve the determination of the orbital characteristics using two fitting approaches, Least-Square Monte Carlo and Markov Chain Monte Carlo. Finally, we analyze the sensitivity of our data to additional closer-in companions and reject the possibility of other massive brown dwarf companions down to 4-5 AU. [abridged]

The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $\mu$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $\mu$m, 1.667 $\mu$m), $K1K2$ (2.110 $\mu$m, 2.251 $\mu$m), and broadband $J$ (1.245 $\mu$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $\mu$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Lo\`eve image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.

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.

Context. Some of transition disks show asymmetric structures in thermal sub-millimetre emission and optical scattered light. These structures can be the result of planet(s) or companions embedded in the disk. Aims. We aim to detect and analyse the scattered light of the transition disk J160421.7-213028, identify disk structures, and compare the results with previous observations of this disk at other wavelengths. Methods. We obtained and analysed new polarised intensity observations of the transition disk J160421.7-213028 with VLT/SPHERE using the visible light instrument ZIMPOL at $R'$ band (0.626$\mu$m). We probe the disk gap down to a radius of confidence of 0.1'' (${\sim}15$ AU at 145 pc). We interpret the results in the context of dust evolution when planets interact with the parental disk. Results. We observe a gap from 0.1 to 0.3'' (${\sim}15$ to 40AU) and a bright annulus as previously detected by HiCIAO $H$ band observations at $1.65\mu$m. The radial width of the annulus is around $40$AU, and its centre is at ${\sim}61$AU from the central star. The peak of the reflected light at 0.626$\mu$m is located 20 AU inward of the cavity detected in the submillimetre. In addition, we detect a dip at a position angle of ${\sim}46.2 \pm 5.4^\circ$. A dip was also detected with HiCIAO, but located at ${\sim}85^\circ$. If the dip observed with HiCIAO is the same, this suggests an average dip rotation of ${\sim}12^\circ/$year, which is inconsistent with the local Keplerian angular velocity of $\sim$0.8$^\circ$/yr at $\sim$61AU. Conclusions. The spatial discrepancy in the radial emission in J160421.7-213028 at different wavelengths is consistent with dust filtration at the outer edge of a gap carved by a massive planet. The dip rotation can be interpreted by fast variability of the inner disk and/or the presence of a warp or circumplanetary material of a planet at ${\sim}9.6$ AU.

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 present the first optical (590--890 nm) imaging polarimetry observations of the pre-transitional protoplanetary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous studies. We detect the previously unseen far side of the disk gap, confirm the highly eccentric scattered-light gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from millimeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This new source of scattered light may contribute to the near-infrared interferometric signal attributed to the protoplanet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present a new model for the system architecture of LkCa 15 that ties these new findings together. These observations were taken during science verification of SPHERE ZIMPOL and demonstrate this facility's performance for faint guide stars under adverse observing conditions.

One of the striking discoveries of protoplanetary disc research in recent years are the spiral arms seen in several transitional discs in polarised scattered light. An interesting interpretation of the observed spiral features is that they are density waves launched by one or more embedded (proto-)planets in the disc. In this paper we investigate whether planets can be held responsible for the excitation mechanism of the observed spirals. We use locally isothermal hydrodynamic simulations as well as analytic formulae to model the spiral waves launched by planets. Then H-band scattered light images are calculated using a 3D continuum radiative transfer code to study the effect of surface density and pressure scale height perturbation on the detectability of the spirals. We find that a relative change of about 3.5 in the surface density is required for the spirals to be detected with current telescopes in the near-infrared for sources at the distance of typical star-forming regions (140pc). This value is a factor of eight higher than what is seen in hydrodynamic simulations. We also find that a relative change of only 0.2 in pressure scale height is sufficient to create detectable signatures under the same conditions. Therefore, we suggest that the spiral arms observed to date in protoplanetary discs are the results of changes in the vertical structure of the disc (e.g. pressure scale height perturbation) instead of surface density perturbations.

The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r~0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

The coagulation of microscopic dust into planetesimals is the first step towards planet formation. The size and shape of the growing aggregates determine the efficiency of this early growth. It has been proposed that fluffy ice aggregates can grow very efficiently, suffering less from the bouncing and radial drift barriers. While the collision velocity between icy aggregates of similar size is thought to stay below the fragmentation threshold, they may nonetheless lose mass from collisions with much smaller projectiles. We investigate the effect of these erosive collisions on the ability of porous ice aggregates to cross the radial drift barrier. We develop a Monte Carlo code that calculates the evolution of the growing aggregates, while resolving the entire mass distribution at all times. The aggregate's porosity is treated independently of its mass, and is determined by collisions, gas compaction, and eventually self-gravity compaction. For erosion threshold velocities of 20-40 m/s, high-velocity collisions with small projectiles prevent the largest aggregates from growing when they start to drift. In these cases, our local simulations result in a steady-state distribution, with the majority of the dust mass in particles with Stokes numbers close to unity. Only for the highest erosion threshold considered (60 m/s), do porous aggregates manage to cross the radial drift barrier in the inner 10 AU of MMSN-like disks. Erosive collisions are more effective in limiting the growth than fragmentary collisions between similar-size particles. Conceivably, erosion limits the growth before the radial drift barrier, although the robustness of this statement depends on (uncertain) material properties of icy aggregates. If erosion inhibits planetesimal formation through direct sticking, the sea of ${\sim}10^9$ g, highly porous particles appears well-suited for triggering streaming instability.

We study the emission and absorption of CO ro-vibrational lines in the spectra of intermediate mass pre-main-sequence stars with the aim to determine both the spatial distribution of the CO gas and its physical properties. We also aim to correlate CO emission properties with disk geometry. Using high-resolution spectra containing fundamental and first overtone CO ro-vibrational emission, observed with CRIRES on the VLT, we probe the physical properties of the circumstellar gas by studying its kinematics and excitation conditions. We detect and spectrally resolve CO fundamental ro-vibrational emission in 12 of the 13 stars observed, and in two cases in absorption. Keeping in mind that we studied a limited sample, we find that the physical properties and spatial distribution of the CO gas correlate with disk geometry. Flaring disks show highly excited CO fundamental emission up to v$_u$ = 5, while self-shadowed disks show CO emission that is not as highly excited. Rotational temperatures range between ~250-2000 K. The $^{13}$CO rotational temperatures are lower than those of $^{12}$CO. The vibrational temperatures in self-shadowed disks are similar to or slightly below the rotational temperatures, suggesting that thermal excitation or IR pumping is important in these lines. In flaring disks the vibrational temperatures reach as high as 6000 K, suggesting fluorescent pumping. Using a simple kinematic model we show that the CO inner radius of the emitting region is $\approx$10 au for flaring disks and $\leq$ 1 au for self-shadowed disks. (abstract truncated due to 1920 character limit for arXiv)

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

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.

$\tau$ Ceti is a nearby, mature G-type star very similar to our Sun, with a massive Kuiper Belt analogue (Greaves et al. 2004) and possible multiplanet system (Tuomi et al. 2013) that has been compared to our Solar System. We present Herschel Space Observatory images of the debris disk, finding the disk is resolved at 70 and 160 microns, and marginally resolved at 250 microns. The Herschel images and infrared photometry from the literature are best modelled using a wide dust annulus with an inner edge between 1-10 AU and an outer edge at ~55 AU, inclined from face-on by 35$\pm$10 degrees, and with no significant azimuthal structure. We model the proposed tightly-packed planetary system of five super-Earths and find that the innermost dynamically stable disk orbits are consistent with the inner edge found by the observations. The photometric modelling, however, cannot rule out a disk inner edge as close to the star as 1 AU, though larger distances produce a better fit to the data. Dynamical modelling shows that the 5 planet system is stable with the addition of a Neptune or smaller mass planet on an orbit outside 5 AU, where the Tuomi et al. analysis would not have detected a planet of this mass.

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.

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

(abridged) In the core accretion scenario for the formation of planetary rocky cores, the first step toward planet formation is the growth of dust grains into larger and larger aggregates and eventually planetesimals. Although dust grains are thought to grow from the submicron sizes typical of interstellar dust to micron size particles in the dense regions of molecular clouds and cores, the growth from micron size particles to pebbles and kilometre size bodies must occur in protoplanetary disks. This step in the formation of planetary systems is the last stage of solids evolution that can be observed directly in young extrasolar systems. In this chapter we review the constraints on the physics of grain-grain collisions as they have emerged from laboratory experiments and numerical computations. We then review the current theoretical understanding of the global processes governing the evolution of solids in protoplanetary disks, including dust settling, growth, and radial transport. The predicted observational signatures are summarized. We discuss recent developments in the study of grain growth in molecular cloud cores and in collapsing envelopes of protostars as these provide the initial conditions for the dust in disks. We discuss the observational evidence for the growth of grains in young disks from mm surveys, as well as the recent evidence of radial variations of the dust properties in disks. We include a brief discussion of the constraints on the small end of the grain size distribution and on dust settling as derived from optical and IR observations. The observations are discussed in the context of global dust evolution models, in particular we focus on the emerging evidence for a very efficient early growth of grains and the radial distribution of grain sizes in disks. We also highlight the limits of current models, including the need to slow the radial drift of grains.

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.

We present H- and Ks-band polarized differential images (PDI) of the Herbig Ae/Be star HD142527, revealing its optically thick outer disk and the nearly empty gap. The very small inner working angle (~0.1") and high resolution achievable with an 8m-class telescope, together with a careful polarimetric calibration strategy, allow us to achieve images that surpass the quality of previous scattered light images. Previously known substructures are resolved more clearly and new structures are seen. Specifically, we are able to resolve 1) half a dozen spiral structures in the disk, including previously known outer-disk spirals as well as new spiral arms and arcs close to the inner rim of the disk; 2) peculiar holes in the polarized surface brightness at position angles of ~0{\deg} and ~160{\deg}; 3) the inner rim on the eastern side of the disk; 4) the gap between the outer and inner disk, ranging from the inner working angle of 0.1" out to between 0.7 and 1.0", which is nearly devoid of dust. We then use a Markov-chain Monte-Carlo algorithm to determine several structural parameters of the disk, using very simple assumptions, including its inclination, eccentricity, and the scale height of the inner rim. We compare our results to previous work on this object, and try to produce a consistent picture of the system and its transition disk.

Transitional disks represent a short stage of the evolution of circumstellar material. Studies of dust grains in these objects can provide pivotal information on the mechanisms of planet formation. Dissimilarities in the spatial distribution of small (micron-size) and large (millimeter-size) dust grains have recently been pointed out. Constraints on the small dust grains can be obtained by imaging the distribution of scattered light at near-infrared wavelengths. We aim at resolving structures in the surface layer of transitional disks (with particular emphasis on the inner 10 - 50 AU), thus increasing the scarce sample of high resolution images of these objects. We obtained VLT/NACO near-IR high-resolution polarimetric differential imaging observations of SAO 206462 (HD135344B). This technique allows one to image the polarized scattered light from the disk without any occulting mask and to reach an inner working angle of 0.1''. A face-on disk is detected in H and Ks bands between 0.1'' and 0.9''. No significant differences are seen between the H and Ks images. In addition to the spiral arms, these new data allow us to resolve for the first time an inner cavity for small dust grains. The cavity size (about 28 AU) is much smaller than what is inferred for large dust grains from (sub)mm observations (39 to 50 AU). The interaction between the disk and potential orbiting companion(s) can explain both the spiral arm structure and the discrepant cavity sizes for small and large dust grains. One planet may be carving out the gas (and, thus, the small grains) at 28 AU, and generating a pressure bump at larger radii (39 AU), which holds back the large grains. We analytically estimate that, in this scenario, a single giant planet (with a mass between 5 and 15 Jupiter masses) at 17 to 20 AU from the star is consistent with the observed cavity sizes.

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

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.

CO is an important component of a protoplanetary disc as it is one of the most abundant gas phase species. Furthermore, observations of the CO transitions can be used as a diagnostic of the gas, tracing conditions in both the inner and outer disc. We present Herschel/PACS spectroscopy of a sample of 22 Herbig Ae/Be (HAEBEs) and 8 T Tauri stars (TTS), covering the pure rotational CO transitions from J = 14-13 up to J = 49-48. CO is detected in only 5 HAEBEs, AB Aur, HD 36112, HD 97048, HD 100546 and IRS 48, and in 4 TTS, AS 205, S CrA, RU Lup, and DG Tau. The highest transition seen is J = 36-35, with Eup of 3669 K, detected in HD 100546 and DG Tau. We construct rotational diagrams for the discs with at least 3 CO detections to derive Trot, and find average temperatures of 270 K for the HAEBEs and 485 K for the TTS. HD 100546 requires an extra temperature component with Trot ~ 900-1000 K, suggesting a range of temperatures in its disc atmosphere, consistent with thermo-chemical disc models. In HAEBEs, the objects with CO detections all have flared discs in which the gas and dust are thermally decoupled. We use a small model grid to analyse our observations and find that an increased amount of flaring means higher line flux, as it increases the mass in warm gas. CO is not detected in our flat discs as the emission is below the detection limit. We find that HAEBE sources with CO detections have high Luv and strong PAH emission, again connected to the heating of the gas. In TTS, the objects with CO detections are all sources with evidence for a disc wind or outflow. For both groups of objects, sources with CO detections generally have high UV luminosity (either stellar in HAEBEs or due to accretion in TTS), but this is not a sufficient condition for the detection of the far-IR CO lines.

Giant planets form in protoplanetary disks while these disks are still gas-rich, and can reveal their presence through the annular gaps they carve out. HD 100546 is a gas-rich disk with a wide gap between between a radius of ~1 and 13 AU, possibly cleared out by a planetary companion or planetary system. We want to identify the nature of the unseen companion near the far end of the disk gap. We use mid-infrared interferometry at multiple baselines to constrain the curvature of the disk wall at the far end of the gap. We use 2D hydrodynamical simulations of embedded planets and brown dwarfs to estimate viscosity of the disk and the mass of a companion close to the disk wall. We find that the disk wall at the far end of the gap is not vertical, but rounded-off by a gradient in the surface density. Such a gradient can be reproduced in hydrodynamical simulations with a single, heavy companion (>=30...80 M_Jup) while the disk has viscosity of at least \alpha{} >~ 5 10^-3. Taking into account the changes in the temperature structure after gap opening reduces the lower limit on the planet mass and disk viscosity to 20 M_Jup and \alpha{} = 2 10^-3. The object in the disk gap of HD 100546 that shapes the disk wall is most likely a 60(+20)(-40) M_Jup brown dwarf, while the disk viscosity is estimated to be at least \alpha{} = 2 10^-3. The disk viscosity is an important factor in estimating planetary masses from disk morphologies: more viscous disks need heavier planets to open an equally deep gap.

[Abridged] Debris disks are extrasolar analogs to the solar system planetesimal belts. The star Fomalhaut harbors a cold debris belt at 140 AU as well as evidence of a warm dust component, which is suspected of being a bright analog to the solar system's zodiacal dust. Interferometric observations obtained with the VLTI and the KIN have identified near- and mid-infrared excesses attributed to hot and warm exozodiacal dust in the inner few AU of the star. We performed parametric modeling of the exozodiacal disk using the GRaTeR radiative transfer code to reproduce the interferometric data, complemented by mid- to far-infrared measurements. A detailed treatment of sublimation temperatures was introduced to explore the hot population at the sublimation rim. We then used an analytical approach to successively testing several source mechanisms. A good fit to the data is found by two distinct dust populations: (1) very small, hence unbound, hot dust grains confined in a narrow region at the sublimation rim of carbonaceous material; (2) bound grains at 2 AU that are protected from sublimation and have a higher mass despite their fainter flux level. We propose that the hot dust is produced by the release of small carbon grains following the disruption of aggregates that originate from the warm component. A mechanism, such as gas braking, is required to further confine the small grains for a long enough time. In situ dust production could hardly be ensured for the age of the star, so the observed amount of dust must be triggered by intense dynamical activity. Fomalhaut may be representative of exozodis that are currently being surveyed worldwide. We propose a framework for reconciling the hot exozodi phenomenon with theoretical constraints: the hot component of Fomalhaut is likely the tip of the iceberg since it could originate from a warm counterpart residing near the ice line.

Broadband spectral surveys of protostars offer a rich view of the physical, chemical and dynamical structure and evolution of star-forming regions. The Herschel Space Observatory opened up the terahertz regime to such surveys, giving access to the fundamental transitions of many hydrides and to the high-energy transitions of many other species. A comparative analysis of the chemical inventories and physical processes and properties of protostars of various masses and evolutionary states is the goal of the Herschel CHEmical Surveys of Star forming regions (CHESS) key program. This paper focusses on the intermediate-mass protostar, OMC-2 FIR 4. We obtained a spectrum of OMC-2 FIR 4 in the 480 to 1902 GHz range with the HIFI spectrometer onboard Herschel and carried out the reduction, line identification, and a broad analysis of the line profile components, excitation, and cooling. We detect 719 spectral lines from 40 species and isotopologs. The line flux is dominated by CO, H2O, and CH3OH. The line profiles are complex and vary with species and upper level energy, but clearly contain signatures from quiescent gas, a broad component likely due to an outflow, and a foreground cloud. We find abundant evidence for warm, dense gas, as well as for an outflow in the field of view. Line flux represents 2% of the 7 L_Sol luminosity detected with HIFI in the 480 to 1250 GHz range. Of the total line flux, 60% is from CO, 13% from H2O and 9% from CH3OH. A comparison with similar HIFI spectra of other sources is set to provide much new insight into star formation regions, a case in point being a difference of two orders of magnitude in the relative contribution of sulphur oxides to the line cooling of Orion KL and OMC-2 FIR 4.

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

Planets have been detected around a variety of stars, including low-mass objects, such as brown dwarfs. However, such extreme cases are challenging for planet formation models. Recent sub-millimeter observations of disks around brown dwarf measured low spectral indices of the continuum emission that suggest that dust grains grow to mm-sizes even in these very low mass environments. To understand the first steps of planet formation in scaled-down versions of T-Tauri disks, we investigate the physical conditions that can theoretically explain the growth from interstellar dust to millimeter-sized grains in disks around brown dwarf. We modeled the evolution of dust particles under conditions of low-mass disks around brown dwarfs. We used coagulation, fragmentation and disk-structure models to simulate the evolution of dust, with zero and non-zero radial drift. For the non-zero radial drift, we considered strong inhomogeneities in the gas surface density profile that mimic long-lived pressure bumps in the disk. We studied different scenarios that could lead to an agreement between theoretical models and the spectral slope found by millimeter observations. We find that fragmentation is less likely and rapid inward drift is more significant for particles in brown dwarf disks than in T-Tauri disks. We present different scenarios that can nevertheless explain millimeter-sized grains. As an example, a model that combines the following parameters can fit the millimeter fluxes measured for brown dwarf disks: strong pressure inhomogeneities of $\sim$ 40% of amplitude, a small radial extent $\sim$ 15 AU, a moderate turbulence strength $\alpha_{\mathrm{turb}}= 10^{-3}$, and average fragmentation velocities for ices $v_f = 10 m s^{-1}$.

OMC-2 FIR 4 is one of the closest known young intermediate-mass protoclusters, located at a distance of 420 pc in Orion. This region is one of the few where the complete 500-2000 GHz spectrum has been observed with the heterodyne spectrometer HIFI on board the Herschel satellite, and unbiased spectral surveys at 0.8, 1, 2 and 3 mm have been obtained with the JCMT and IRAM 30-m telescopes. In order to investigate the morphology of this region, we used the IRAM Plateau de Bure Interferometer to image OMC-2 FIR 4 in the 2-mm continuum emission, as well as in DCO+(2-1), DCN(2-1), C34S(3-2), and several CH3OH lines. In addition, we analysed observations of the NH3(1,1) and (2,2) inversion transitions made with the Very Large Array of the NRAO. The resulting maps have an angular resolution which allows us to resolve structures of 5", equivalent to 2000 AU. Our observations reveal three spatially resolved sources within OMC-2 FIR 4, of one or several solar masses each, with hints of further unresolved substructure within them. Two of these sources have elongated shapes and are associated with dust continuum emission peaks, thus likely containing at least one molecular core each. One of them also displays radio continuum emission, which may be attributed to a young B3-B4 star that dominates the overall luminosity output of the region. The third source identified displays a DCO+(2-1) emission peak, and weak dust continuum emission. Its higher abundance of DCO+ relative to the other two regions suggests a lower temperature and therefore its possible association with either a younger low-mass protostar or a starless core. It may alternatively be part of the colder envelope of OMC-2 FIR 4. Our interferometric observations evidence the complexity of this region, where multiple cores, chemical differentiation and an ionised region all coexist within an area of only 10000 AU.

Context. Transition disks typically appear in resolved millimeter observations as giant dust rings surrounding their young host stars. More accurate observations with ALMA have shown several of these rings to be in fact asymmetric: they have lopsided shapes. It has been speculated that these rings act as dust traps, which would make them important laboratories for studying planet formation. It has been shown that an elongated giant vortex produced in a disk with a strong viscosity jump strikingly resembles the observed asymmetric rings. Aims. We aim to study a similar behavior for a disk in which a giant planet is embedded. However, a giant planet can induce two kinds of asymmetries: (1) a giant vortex, and (2) an eccentric disk. We studied under which conditions each of these can appear, and how one can observationally distinguish between them. This is important because only a vortex can trap particles both radially and azimuthally, while the eccentric ring can only trap particles in radial direction. Methods. We used the FARGO code to conduct the hydro-simulations. We set up a disk with an embedded giant planet and took a radial grid spanning from 0.1 to 7 times the planet semi-major axis. We ran the simulations with various viscosity values and planet masses for 1000 planet orbits to allow a fully developed vortex or disk eccentricity. Afterwards, we compared the dust distribution in a vortex-holding disk with an eccentric disk using dust simulations. Results. We find that vorticity and eccentricity are distinguishable by looking at the azimuthal contrast of the dust density. While vortices, as particle traps, produce very pronounced azimuthal asymmetries, eccentric features are not able to accumulate millimeter dust particles in azimuthal direction, and therefore the asymmetries are expected to be modest.

Context: We have analysed Herschel-PACS spectra of 32 circumstellar disks around Herbig Ae/Be and T-Tauri stars obtained within the Herschel key programme DIGIT. In this paper we focus on the 69mu emission band of the crystalline silicate forsterite. Aims: This work provides an overview of the 69mu forsterite bands in the DIGIT sample. We aim to derive the temperature and composition of the forsterite grains. With this information, constraints can be placed on the spatial distribution of the forsterite in the disk and its formation history. Methods: Position and shape of the 69mu band are used to derive the temperature and composition of the dust by comparison to laboratory spectra of that band. We combine our data with existing Spitzer IRS spectra to compare the presence and strength of the 69mu band to the forsterite bands at shorter wavelengths. Results: A total of 32 sources have been observed, 8 of them show a 69mu emission band that can be attributed to forsterite. With the exception of the T-Tauri star AS205, all of the detections are for disks associated with Herbig Ae/Be stars. Most of the forsterite grains that give rise to the 69mu bands are warm (~100-200 K) and iron-poor (less than ~2% iron). Only AB-Aur requires approximately 3-4% of iron. Conclusions: Our findings support the hypothesis that the forsterite grains form through an equilibrium condensation process at high temperatures. The connection between the strength of the 69 and 33mu bands shows that at least part of the emission in these bands originates from the same dust grains. Further, any model that explains the PACS and the Spitzer IRS observations must take the effects of a wavelength dependent optical depth into account. We find indications of a correlation of the detection rate of the 69mu band with the spectral type of the host stars. However, our sample is too small to obtain a definitive result.

We present Subaru/HiCIAO H-band high-contrast images of the debris disk around HIP 79977, whose pres- ence was recently inferred from an infrared excess. Our images resolve the disk for the first time, allowing characterization of its shape, size, and dust grain properties. We use angular differential imaging (ADI) to reveal the disk geometry in unpolarized light out to a radius of ~2", as well as polarized differential imaging (PDI) to measure the degree of scattering polarization out to ~1.5". In order to strike a favorable balance between suppression of the stellar halo and conservation of disk flux, we explore the application of principal component analysis (PCA) to both ADI and reference star subtraction. This allows accurate forward modeling of the effects of data reduction on simulated disk images, and thus direct comparison with the imaged disk. The resulting best-fit values and well-fitting intervals for the model parameters are a surface brightness power-law slope of S_out = -3.2 [-3.6,-2.9], an inclination of i = 84{\deg} [81{\deg},86{\deg}], a high Henyey-Greenstein forward-scattering parameter of g = 0.45 [0.35, 0.60], and a non-significant disk-star offset of u = 3.0 [-1.5, 7.5] AU = 24 [-13, 61] mas along the line of nodes. Furthermore, the tangential linear polarization along the disk rises from ~10% at 0.5" to ~45% at 1.5". These measurements paint a consistent picture of a disk of dust grains produced by collisional cascades and blown out to larger radii by stellar radiation pressure.

As part of the Herschel Guaranteed Time Key Programme CHESS, we present the discovery of a diffuse gas component in the foreground of the intermediate-mass protostar OMC-2 FIR 4, located in the Orion A region. Making use of the full HIFI spectrum of OMC-2 FIR 4 obtained in CHESS, we detected several ground-state lines from OH+, H2O+, HF, and CH+, all of them seen in absorption against the dust continuum emission of the protostar's envelope. The lines peak at a velocity of 9 km/s, which is blue-shifted by 2 km/s with respect to the systemic velocity of OMC-2 FIR 4 (V = 11.4 km/s). We derived column densities for each species, as well as an upper limit to the column density of the undetected H3O+. In order to model and characterise the foreground cloud, we used the Meudon PDR code to run a homogeneous grid of models that spans a reasonable range of densities, visual extinctions, cosmic ray ionisation rates and far-ultraviolet (FUV) radiation fields. The results of our modelling indicate that the foreground cloud is composed of predominantly neutral diffuse gas (n(H) = 100 cm-3) and is heavily irradiated by an external source of FUV that most likely arises from the nearby Trapezium OB association. The cloud is 6 pc thick and bears many similarities with the so-called C+ interface between Orion-KL and the Trapezium cluster, 2 pc south of OMC-2 FIR 4. We conclude that the foreground cloud we detected is an extension of the C+ interface seen in the direction of Orion KL, and interpret it to be the remains of the parental cloud of OMC-1, which extends from OMC-1 up to OMC-2.

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.

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

The progenitors of supernovae (SNe) type Ia are usually assumed to be either a single white dwarf (WD) accreting from a non-degenerate companion (the SD channel) or the result of two merging WDs (DD channel). However, no consensus currently exists as to which progenitor scenario is the correct one, or whether the observed SN Ia rate is produced by a combination of both channels. Unlike a DD progenitor a SD progenitor is expected to emit supersoft X-rays for a prolonged period of time (~1 Myr) as a result of the burning of accreted matter on the surface of the WD. An argument against the SD channel as a significant producer of SNe type Ia has been the lack of observed supersoft X-ray sources (SSS) and the lower-than-expected integrated soft X-ray flux from elliptical galaxies. We wish to determine if it is possible to obscure the supersoft X-ray emission from a nuclear burning white dwarf in an accreting single degenerate binary system. In case of obscured systems we wish to determine their general observational characteristics. We examine the emergent X-ray emission from a canonical SSS system surrounded by a spherically symmetric configuration of material, assuming a black body spectrum with T_BB=50 eV and L=10^38 erg/s. The circumbinary material is assumed to be of solar chemical abundances, and we leave the mechanism behind the mass loss into the circumbinary region unspecified. If steadily accreting, nuclear burning WDs are canonical SSS our analysis suggests that they can be obscured by relatively modest circumbinary mass loss rates. This may explain the discrepancy of SSS compared to the SN Ia rate inferred from observations if the SD progenitor scenario contributes significantly to the SN Ia rate. Recycled emissions from obscured systems may be visible in other wavebands than X-rays. It may also explain the lack of observed SSS in symbiotic binary systems.

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.

Dust settling and grain growth are the first steps in the planet-formation process in protoplanetary disks. These disks are observed around stars with different spectral types, and there are indications that the disks around lower mass stars are significantly flatter, which could indicate that they settle and evolve faster, or in a different way. We aim to test this assumption by modeling the median spectral energy distributions (SEDs) of three samples of protoplanetary disks: around Herbig stars, T Tauri stars and brown dwarfs. We focus on the turbulent mixing strength to avoid a strong observational bias from disk and stellar properties that depend on stellar mass. We generated SEDs with the radiative transfer code MCMax, using a hydrostatic disk structure and settling the dust in a self-consistent way with the alpha-prescription to probe the turbulent mixing strength. We are able to fit all three samples with a disk with the same input parameters, scaling the inner edge to the dust evaporation radius and disk mass to millimeter photometry. The Herbig stars require a special treatment for the inner rim regions, while the T-Tauri stars require viscous heating, and the brown dwarfs lack a good estimate of the disk mass because only few millimeter detections exist. We find that the turbulent mixing strength does not vary across the stellar mass range for a fixed grain size distribution and gas-to-dust ratio. Regions with the same temperature have a self-similar vertical structure independent of stellar mass, but regions at the same distance from the central star appear more settled in disks around lower mass stars. We find a relatively low turbulent mixing strength of alpha = 10^(-4) for a standard grain size distribution, but our results are also consistent with alpha = 0.01 for a grain size distribution with fewer small grains or a lower gas-to-dust ratio.

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (~9 Jupiter masses at 0.5", ~3 Jupiter masses at 1").

We investigate the possibility of obscuring supersoft X-ray sources in the winds of companion stars. We derive limits on the amount of circumstellar material needed to fully obscure a 'canonical' supersoft X-ray source in the Large Magellanic Cloud, as observed with the Chandra X-ray Observatory.

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.

A small number of young stellar objects show signs of a halo-like structure of optically thin dust. This halo or torus is located within a few AU of the star, but its origin has not yet been understood. A dynamically excited cloud of planetesimals colliding to eventually form dust could produce such a structure. The cause of the dynamical excitation could be one or more planets. This work investigates an inwardly migrating planet that is dynamically scattering planetesimals as a possible cause for the observed structures. If this mechanism is responsible, the observed halo-like structure could be used to infer the existence of planets in these systems. We present analytical estimates on the maximum inclination reached owing to dynamical interactions between planetesimals and a migrating planet. A symplectic integrator is used to simulate the effect of a migrating planet on a population of planetesimals. It is found that an inwardly migrating planet is only able to scatter the material it encounters to highly-inclined orbits if the material is on an eccentric orbit. Such eccentric orbits can be the result of resonance trapping and eccentricity pumping. Simulations show that for a certain range of migration rates and planet masses, resonance capture combined with planetary migration indeed causes the planetesimals to reach eccentric orbits and subsequently get scattered to highly-inclined orbits. The size distribution of the resulting dust is calculated determined to find the total mass and optical depth, which are found to compare reasonably well with the observed structures. Dynamical scattering of planetesimals caused by a planet migrating in, followed by the grinding down of these planetesimals to dust grains, appears to be a promising explanation for the inferred circumstellar dust clouds. Further study is needed to see if the haloes can be used to infer the presence of planets.

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.

We study the effect of radiation pressure on the dust in the inner rim of transitional disks with large inner holes. In particular, we evaluate whether radiation pressure can be responsible for keeping the inner holes dust-free, while allowing gas accretion to proceed. This has been proposed in a paper by Chiang and Murray-Clay (2007, Nature Physics 3, p. 604) who explain the formation of these holes as an inside-out evacuation due to X- ray-triggered accretion of the innermost layer of the disk rim outside of the hole. We show that radiation pressure is clearly incapable of stopping dust from flowing into the hole because of dust pile-up and optical depth effects, and also because of viscous mixing. Other mechanisms need to be found to explain the persistence of the opacity hole in the presence of accretion, and we speculate on possible solutions.

Astrometric monitoring of the Sirius binary system over the past century has yielded several predictions for an unseen third system component, the most recent one suggesting a \leq50 MJup object in a ~6.3-year orbit around Sirius A. Here we present two epochs of high-contrast imaging observations performed with Subaru IRCS and AO188 in the 4.05 \mum narrow-band Br alpha filter. These data surpass previous observations by an order of magnitude in detectable companion mass, allowing us to probe the relevant separation range down to the planetary mass regime (6-12 M_Jup at 1", 2-4 M_Jup at 2", and 1.6 M_Jup beyond 4"). We complement these data with one epoch of M-band observations from MMT/AO Clio, which reach comparable performance. No dataset reveals any companion candidates above the 5-sigma level, allowing us to refute the existence of Sirius C as suggested by the previous astrometric analysis. Furthermore, our Br alpha photometry of Sirius B confirms the lack of an infrared excess beyond the white dwarf's blackbody spectrum.

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

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

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 Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

CONTEXT: Water vapour maser emission from evolved oxygen-rich stars remains poorly understood. Additional observations, including polarisation studies and simultaneous observation of different maser transitions may ultimately lead to greater insight. AIMS: We have aimed to elucidate the nature and structure of the VY CMa water vapour masers in part by observationally testing a theoretical prediction of the relative strengths of the 620.701 GHz and the 22.235 GHz maser components of ortho water vapour. METHODS: In its high-resolution mode (HRS) the Herschel Heterodyne Instrument for the Infrared (HIFI) offers a frequency resolution of 0.125 MHz, corresponding to a line-of-sight velocity of 0.06 km/s, which we employed to obtain the strength and linear polarisation of maser spikes in the spectrum of VY CMa at 620.701 GHz. Simultaneous ground based observations of the 22.235 GHz maser with the Max-Planck-Institut f\"ur Radioastronomie 100-meter telescope at Effelsberg, provided a ratio of 620.701 GHz to 22.235 GHz emission. RESULTS:We report the first astronomical detection to date of water vapour maser emission at 620.701 GHz. In VY CMa both the 620.701 and the 22.235 GHz polarisation are weak. At 620.701 GHz the maser peaks are superposed on what appears to be a broad emission component, jointly ejected asymmetrically from the star. We observed the 620.701 GHz emission at two epochs 21 days apart, both to measure the potential direction of linearly polarised maser components and to obtain a measure of the longevity of these components. Although we do not detect significant polarisation levels in the core of the line, they rise up to approximately 6% in its wings.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

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.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predictions from conventional chemical models by many orders of magnitude. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen-rich stars using the NH3 ortho J_K = 1_0 - 0_0 ground-state line. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 J_K = 1_0 - 0_0 transition near 572.5 GHz, simultaneously with the ortho-H2O J_Ka,Kc = 1_1,0 -1_0,1 transition, toward VY CMa, OH 26.5+0.6, IRC+10420, and IK Tau. We conducted non-LTE radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the latter two stars, Very Large Array imaging of NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those obtained for the ground state H2O line. The NH3 abundances relative to H2 are very high and range from 2 x 10-7 to 3 x 10-6 for the objects we have studied. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield satisfactory fits to the observed line profiles, it leads to abundance estimates that confirm the very high values found in earlier studies. New ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

The Heterodyne Instrument for the Far Infrared (HIFI) onboard the Herschel Space Observatory allows the first observations of light diatomic molecules at high spectral resolution and in multiple transitions. Here, we report deep integrations using HIFI in different lines of hydrides towards the high-mass star forming region AFGL 2591. Detected are CH, CH+, NH, OH+, H2O+, while NH+ and SH+ have not been detected. All molecules except for CH and CH+ are seen in absorption with low excitation temperatures and at velocities different from the systemic velocity of the protostellar envelope. Surprisingly, the CH(JF,P = 3/2_2,- - 1/2_1,+) and CH+(J = 1 - 0, J = 2 - 1) lines are detected in emission at the systemic velocity. We can assign the absorption features to a foreground cloud and an outflow lobe, while the CH and CH+ emission stems from the envelope. The observed abundance and excitation of CH and CH+ can be explained in the scenario of FUV irradiated outflow walls, where a cavity etched out by the outflow allows protostellar FUV photons to irradiate and heat the envelope at larger distances driving the chemical reactions that produce these molecules.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

Water in Star-forming regions with Herschel (WISH) is a Herschel Key Program investigating the water chemistry in young stellar objects (YSOs) during protostellar evolution. Hydroxyl (OH) is one of the reactants in the chemical network most closely linked to the formation and destruction of H2O. High-temperature chemistry connects OH and H2O through the OH + H2 <-> H2O + H reactions. Formation of H2O from OH is efficient in the high-temperature regime found in shocks and the innermost part of protostellar envelopes. Moreover, in the presence of UV photons, OH can be produced from the photo-dissociation of H2O. High-resolution spectroscopy of the OH 163.12 micron triplet towards HH 46 and NGC 1333 IRAS 2A was carried out with the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel. The low- and intermediate-mass YSOs HH 46, TMR 1, IRAS 15398-3359, DK Cha, NGC 7129 FIRS 2, and NGC 1333 IRAS 2A were observed with the Photodetector Array Camera and Spectrometer (PACS) in four transitions of OH and two [OI] lines. The OH transitions at 79, 84, 119, and 163 micron and [OI] emission at 63 and 145 micron were detected with PACS towards the class I low-mass YSOs as well as the intermediate-mass and class I Herbig Ae sources. No OH emission was detected from the class 0 YSO NGC 1333 IRAS 2A, though the 119 micron was detected in absorption. With HIFI, the 163.12 micron was not detected from HH 46 and only tentatively detected from NGC 1333 IRAS 2A. The combination of the PACS and HIFI results for HH 46 constrains the line width (FWHM > 11 km/s) and indicates that the OH emission likely originates from shocked gas. This scenario is supported by trends of the OH flux increasing with the [OI] flux and the bolometric luminosity. Similar OH line ratios for most sources suggest that OH has comparable excitation temperatures despite the different physical properties of the sources.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

We performed Herschel/HIFI observations of several CO lines in the far-infrared/sub-mm in the protoplanetary nebula CRL618. The high spectral resolution provided by HIFI allows measurement of the line profiles. Since the dynamics and structure of the nebula is well known from mm-wave interferometric maps, it is possible to identify the contributions of the different nebular components (fast bipolar outflows, double shells, compact slow shell) to the line profiles. The observation of these relatively high-energy transitions allows an accurate study of the excitation conditions in these components, particularly in the warm ones, which cannot be properly studied from the low-energy lines. The 12CO J=16-15, 10-9, and 6-5 lines are easily detected in this source. 13CO J=10-9 and 6-5 are also detected. Wide profiles showing spectacular line wings have been found, particularly in 12CO 16-15. Other lines observed simultaneously with CO are also shown. Our analysis of the CO high-J transitions, when compared with the existing models, confirms the very low expansion velocity of the central, dense component, which probably indicates that the shells ejected during the last AGB phases were driven by radiation pressure under a regime of maximum transfer of momentum. No contribution of the diffuse halo found from mm-wave data is identified in our spectra, because of its low temperature. We find that the fast bipolar outflow is quite hot, much hotter than previously estimated; for instance, gas flowing at 100 km/s must have a temperature higher than ~ 200 K. Probably, this very fast outflow, with a kinematic age < 100 yr, has been accelerated by a shock and has not yet cooled down. The double empty shell found from mm-wave mapping must also be relatively hot, in agreement with the previous estimate.

A set of new, sensitive, and spectrally resolved, sub-millimeter line observations are used to probe the warm circumstellar gas around the S-type AGB star chi Cyg. The observed lines involve high rotational quantum numbers, which, combined with previously obtained lower-frequency data, make it possible to study in detail the chemical and physical properties of, essentially, the entire circumstellar envelope of chi Cyg. The data were obtained using the HIFI instrument aboard Herschel, whose high spectral resolution provides valuable information about the line profiles. Detailed, non-LTE, radiative transfer modelling, including dust radiative transfer coupled with a dynamical model, has been performed to derive the temperature, density, and velocity structure of the circumstellar envelope. We report the first detection of circumstellar H2O rotational emission lines in an S-star. Using the high-J CO lines to derive the parameters for the circumstellar envelope, we modelled both the ortho- and para-H2O lines. Our modelling results are consistent with the velocity structure expected for a dust-driven wind. The derived total H2O abundance (relative to H2) is (1.1 +/- 0.2)E-5, much lower than that in O-rich stars. The derived ortho-to-para ratio of 2.1 +/- 0.6 is close to the high-temperature equilibrium limit, consistent with H2O being formed in the photosphere.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

We report the discovery of water vapour toward the carbon star V Cygni. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1(11) - 0(00) para-water transition at 1113.3430 GHz in the upper sideband of the Band 4b receiver. The observed spectral line profile is nearly parabolic, but with a slight asymmetry associated with blueshifted absorption, and the integrated antenna temperature is 1.69 \pm 0.17 K km/s. This detection of thermal water vapour emission, carried out as part of a small survey of water in carbon-rich stars, is only the second such detection toward a carbon-rich AGB star, the first having been obtained by the Submillimeter Wave Astronomy Satellite toward IRC+10216. For an assumed ortho-to-para ratio of 3 for water, the observed line intensity implies a water outflow rate ~ (3 - 6) E-5 Earth masses per year and a water abundance relative to H2 of ~ (2-5) E-6. This value is a factor of at least 1E+4 larger than the expected photospheric abundance in a carbon-rich environment, and - as in IRC+10216 - raises the intriguing possibility that the observed water is produced by the vapourisation of orbiting comets or dwarf planets. However, observations of the single line observed to date do not permit us to place strong constraints upon the spatial distribution or origin of the observed water, but future observations of additional transitions will allow us to determine the inner radius of the H2O-emitting zone, and the H2O ortho-to-para ratio, and thereby to place important constraints upon the origin of the observed water emission.

We are carrying out a physical and chemical study of the protostellar envelopes in a representative sample of IM Class 0 protostars. In our first paper (Crimier et al. 2010), we determined the physical structure (density-temperature radial profiles) of the protostellar envelopes. Here, we study the CO depletion and N2H+ deuteration. We observed the millimeter lines of C18O, C17O, N2H+ and N2D+ toward the protostars using the IRAM 30m telescope. Based on these observations, we derived the C18O, N2H+ and N2D+ radial abundance profiles across their envelopes using a radiative transfer code. In addition, we modeled the chemistry of the protostellar envelopes. All the C18O 1-0 maps are well fit assuming that the C18O abundance decreases inwards within the protostellar envelope until the gas and dust reach the CO evaporation temperature, 20-25K, where the CO is released back to the gas phase. The N2H+ deuterium fractionation in Class 0 IMs is [N2D+]/[N2H+]=0.005-0.014, two orders of magnitude higher than the elemental [D/H] value in the interstellar medium, but a factor of 10 lower than in pre-stellar clumps. Chemical models account for the C18O and N2H+ observations if we assume the CO abundance is 2 times lower than the canonical value in the inner envelope. This could be the consequence of the CO being converted into CH3OH on the grain surfaces prior to the evaporation and/or the photodissociation of CO by the stellar UV radiation. The deuterium fractionation is not fitted by chemical models. This discrepancy is very likely caused by the simplicity of our model that assumes spherical geometry and neglects important phenomena like the effect of bipolar outflows and UV radiation from the star. More important, the deuterium fractionation is dependent on the ortho-to-para H2 ratio, which is not likely to reach the steady-state value in the dynamical time scales of these protostars.

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.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

We obtained Herschel PACS and SPIRE images of the thermal emission of the debris disk around the A5V star {\beta} Pic. The disk is well resolved in the PACS filters at 70, 100, and 160 {\mu}m. The surface brightness profiles between 70 and 160 {\mu}m show no significant asymmetries along the disk, and are compatible with 90% of the emission between 70 and 160 {\mu}m originating in a region closer than 200 AU to the star. Although only marginally resolving the debris disk, the maps obtained in the SPIRE 250 - 500 {\mu}m filters provide full-disk photometry, completing the SED over a few octaves in wavelength that had been previously inaccessible. The small far-infrared spectral index ({\beta} = 0.34) indicates that the grain size distribution in the inner disk (<200AU) is inconsistent with a local collisional equilibrium. The size distribution is either modified by non-equilibrium effects, or exhibits a wavy pattern, caused by an under-abundance of impactors which have been removed by radiation pressure.

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

We present five band imaging of the Vega debris disc obtained using the Herschel Space Observatory. These data span a wavelength range of 70-500 um with full-width half-maximum angular resolutions of 5.6-36.9". The disc is well resolved in all bands, with the ring structure visible at 70 and 160 um. Radial profiles of the disc surface brightness are produced, and a disc radius of 11" (~ 85 AU) is determined. The disc is seen to have a smooth structure thoughout the entire wavelength range, suggesting that the disc is in a steady state, rather than being an ephemeral structure caused by the recent collision of two large planetesimals.

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

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.

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 coagulation of dust aggregates occurs in various astrophysical environments. Each one is characterized by different conditions that influence the growth, e.g. relative velocities, composition, and size of the smallest constituents (monomers). Here we study the microphysics of collisions of dust aggregates in a four-dimensional parameter space. The parameters are the collision energy, the initial compactness of agglomerates, the mass ratio of collision partners, and the impact parameter. For this purpose we employ a state of the art molecular dynamics type of model that has been extensively and successfully tested against laboratory experiments. It simulates the motion of individual monomers interacting dynamically via van der Waals surface forces. The structure of aggregates is quantified by the filling factor that provides information about the internal structure, the packing density of monomers, and the projected surface area of aggregates. Our results show the importance of the impact parameter that causes formation of elongated particles, due to tensile forces acting in offset collisions. We also describe in detail the internal structure of the resulting aggregates. Our findings are summarized in the form of a simple collision recipe. The recipe specifies the outcome of a collision, averaged over the impact parameter. It is provided in tabular form for a range of physical parameters such as impact energy and pre-collision filling factor. The dependence on the mass ratio of impactor and target is taken into account by providing both a local and a global branch of the recipe.

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.

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.

In dense molecular clouds collisions between dust grains alter the ISM-dust size distribution. We study this process by inserting the results from detailed numerical simulations of two colliding dust aggregates into a coagulation model that computes the dust size distribution with time. All collisional outcomes -- sticking, fragmentation (shattering, breakage, and erosion) -- are included and the effects on the internal structure of the aggregates are also tabulated. The dust aggregate evolution model is applied to an homogeneous and static cloud of temperature 10 K and gas densities between 10^3 and 10^7 cm^-3. The coagulation is followed locally on timescales of ~10^7 yr. We find that the growth can be divided into two stages: a growth dominated phase and a fragmentation dominated phase. Initially, the mass distribution is relatively narrow and shifts to larger sizes with time. At a certain point, dependent on the material properties of the grains as well as on the gas density, collision velocities will become sufficiently energetic to fragment particles, halting the growth and replenishing particles of lower mass. Eventually, a steady state is reached, where the mass distribution is characterized by a mass spectrum of approximately equal amount of mass per logarithmic size bin. The amount of growth that is achieved depends on the cloud's lifetime. If clouds exist on free-fall timescales the effects of coagulation on the dust size distribution are very minor. On the other hand, if clouds have long-term support mechanisms, the impact of coagulation is important, resulting in a significant decrease of the opacity on timescales longer than the initial collision timescale between big grains.

We investigate the physical properties and spatial distribution of Carbon Monoxide (CO) gas in the disks around the Herbig Ae/Be stars HD 97048 and HD 100546. Using high-spectral-resolution 4.588-4.715 $\mu$m spectra containing fundamental CO emission taken with CRIRES on the VLT, we probe the circumstellar gas and model the kinematics of the emission lines. By using spectro-astrometry on the spatially resolved targets, we constrain the physical size of the emitting regions in the disks. We resolve, spectrally and spatially, the emission of the $^{13}$CO v(1-0) vibrational band and the $^{12}$CO $v=1-0, v=2-1, v=3-2$ and $v=4-3$ vibrational bands in both targets, as well as the $^{12}$CO $v=5-4$ band in HD 100546. Modeling of the CO emission with a homogeneous disk in Keplerian motion, yields a best fit with an inner and outer radius of the CO emitting region of 11 and $\geq$ 100 AU for HD 97048. HD 100546 is not fit well with our model, but we derive a lower limit on the inner radius of 8 AU. The fact that gaseous [OI] emission was previously detected in both targets at significantly smaller radii suggests that CO may be effectively destroyed at small radii in the surface layers of these disks

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.

We present a parameter study of self-consistent models of protoplanetary disks around Herbig AeBe stars. We use the code developed by Dullemond and Dominik, which solves the 2D radiative transfer problem including an iteration for the vertical hydrostatic structure of the disk. This grid of models will be used for several studies on disk emission and mineralogy in followup papers. In this paper we take a first look on the new models, compare them with previous modeling attempts and focus on the effects of various parameters on the overall structure of the SED that leads to the classification of Herbig AeBe stars into two groups, with a flaring (group I) or self-shadowed (group II) SED. We find that the parameter of overriding importance to the SED is the total mass in grains smaller than 25um, confirming the earlier results by Dullemond and Dominik. All other parameters studied have only minor influences, and will alter the SED type only in borderline cases. We find that there is no natural dichotomy between group I and II. From a modeling point of view, the transition is a continuous function of the small dust mass. We also show that moderate grain growth produces spectra with weak or no 1um feature, both for flaring (Group I) and non-flaring (Group II) sources. The fact that sources with weak features have been found mostly in Group I sources is therefore surprising and must be due to observational biases or evolutionary effects.

Turbulent coagulation in protoplanetary disks is known to operate on timescale far shorter than the lifetime of the disk. In the absence of mechanisms that replenish the small dust grain population, protoplanetary disks would rapidly lose their continuum opacity-bearing dust. This is inconsistent with infrared observations of disks around T Tauri stars and Herbig Ae/Be stars, which are usually optically thick at visual wavelengths and show signatures of small (a<~ 3um) grains. A plausible replenishing mechanism of small grains is collisional fragmentation or erosion of large dust aggregates, which model calculations predict to play an important role in protoplanetary disks. If optically thick disks are to be seen as proof for ongoing fragmentation or erosion, then alternative explanations for the existence of optically thick disks must be studied carefully. In this study we explore two scenarios. First, we study the effect of residual, low-level infall of matter onto the disk surface. We find that infall rates as low as 10^{-11} Msun/yr can, in principle, replenish the small grain population to a level that keeps the disk marginally optically thick. However, it remains to be seen if the assumption of such inflow is realistic for star+disk systems at the age of several Myrs, at which winds and jets are expected to have removed any residual envelope. In summary, fragmentation or erosion still appear to be the most promising processes to explain the abundant presence of small grains in old disks.

Aims: We wish to investigate what the effect of dust sedimentation is on the observed 10 mum feature of protoplanetary disks and how this may affect the interpretation of the observations. Methods: Using a combination of modeling tools, we simulated the sedimentation of a dust grain size distribution in an axisymmetric 2-D model of a turbulent protoplanetary disk, and we used a radiative transfer program to compute the resulting spectra. Results: We find that the sedimentation can turn a flat feature into a pointy one, but only to a limited degree and for a very limited set of particle size distributions. Only if we have a bimodal size distribution, i.e. a very small grain population and a bigger grain population, do we find that the transformation from a flat to a pointy feature upon dust sedimentation is strong. However, our model shows that, if sedimentation is the sole reason for the variety of silicate feature strengths observed in protoplanetary disks, then we would expect to find a correlation such that disks with weak mid- to far-infrared excess have a stronger 10 mum silicate feature than disks with a strong mid- to far-infrared excess. If this is contrary to what is observed, then this would indicate that sedimentation cannot be the main reason for the variety of 10 mum silicate features observed in protoplanetary disks.

Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy.

The formation of planetesimals requires the growth of dust particles through collisions. Micron-sized particles must grow by many orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical properties, and the interaction cross sections of aggregates with surrounding gas must be well understood. Recent advances in experimental (laboratory) studies now provide the background for pushing numerical aggregate models onto a new level. We present the calibration of a previously tested model of aggregate dynamics. We use plastic deformation of surface asperities as the physical model to bring critical velocities for sticking into accordance with experimental results. The modified code is then used to compute compression strength and the velocity of sound in the aggregate at different densities. We compare these predictions with experimental results and conclude that the new code is capable of studying the properties of small aggregates.

We present high spectral resolution optical spectra of three young intermediate mass stars, in all of which we spectrally resolve the 6300 Angstrom [OI] emission line. Two of these have a double peaked line profile. We fit these data with a simple model of the [OI] emission caused by photo-dissociation of OH molecules in the upper layer of a circumstellar disk by stellar UV radiation and thus translate the Doppler broadened [OI] emission profile into an amount of emission as a function of distance from the central star. The resulting spectra are in agreement with the expected disk shapes as derived from their spectral energy distribution. We find evidence for shadowing by an inner rim in the disk surrounding HD101412 and see a flaring disk structure in HD179218 while the [OI] spectrum of HD135344 is more complex. The [OI] emission starts for all three targets at velocities corresponding to their dust sublimation radius and extends up to radii of 10 -- 90 AU. This shows that this method can be a valuable tool in the future investigation of circumstellar disks.

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.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

The formation of planets is directly linked to the evolution of the circumstellar (CS) disk from which they are born. The dissipation timescales of CS disks are, therefore, of direct astrophysical importance in evaluating the time available for planet formation. We employ Spitzer Space Telescope spectra to complete the CS disk census for the late-type members of the ~8 Myr-old eta Chamaeleontis star cluster. Of the 15 K- and M-type members, eight show excess emission. We find that the presence of a CS disk is anti-correlated with binarity, with all but one disk associated with single stars. With nine single stars in total, about 80% retain a CS disk. Of the six known or suspected close binaries the only CS disk is associated with the primary of RECX 9. No circumbinary disks have been detected. We also find that stars with disks are slow rotators with surface values of specific angular momentum j = 2-15 j_sun. All high specific angular momentum systems with j = 20-30 j_sun are confined to the primary stars of binaries. This provides novel empirical evidence for rotational disk locking and again demonstrates the much shorter disk lifetimes in close binary systems compared to single star systems. We estimate the characteristic mean disk dissipation timescale to be ~5 Myr and ~9 Myr for the binary and single star systems, respectively.

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.

Prestellar cores are unique laboratories for studies of the chemical and physical conditions preceding star formation. We observed the prestellar core L1544 in the fundamental transition of ortho-H2D+ (1_1,0-1_1,1) at different positions over 100", and found a strong correlation between its abundance and the CO depletion factor. We also present a tentative detection of the fundamental transition of para-D2H+ (1_1,0-1_0,1) at the dust emission peak. Maps in N2H+, N2D+, HCO+ and DCO+ are used, and interpreted with the aid of a spherically symmetric chemical model that predicts the column densities and abundances of these species as a function of radius. The correlation between the observed deuterium fractionation of H3+, N2H+ and HCO+ and the observed integrated CO depletion factor across the core can be reproduced by this chemical model. In addition a simpler model is used to study the H2D+ ortho-to-para ratio. We conclude that, in order to reproduce the observed ortho-H2D+ observations, the grain radius should be larger than 0.3 microns.

We discuss the results of laboratory measurements and theoretical models concerning the aggregation of dust in protoplanetary disks, as the initial step toward planet formation. Small particles easily stick when they collide and form aggregates with an open, often fractal structure, depending on the growth process. Larger particles are still expected to grow at collision velocities of about 1m/s. Experiments also show that, after an intermezzo of destructive velocities, high collision velocities above 10m/s on porous materials again lead to net growth of the target. Considerations of dust-gas interactions show that collision velocities for particles not too different in surface-to-mass ratio remain limited up to sizes about 1m, and growth seems to be guaranteed to reach these sizes quickly and easily. For meter sizes, coupling to nebula turbulence makes destructive processes more likely. Global aggregation models show that in a turbulent nebula, small particles are swept up too fast to be consistent with observations of disks. An extended phase may therefore exist in the nebula during which the small particle component is kept alive through collisions driven by turbulence which frustrates growth to planetesimals until conditions are more favorable for one or more reasons.

It has been suggested that the dark halos of galaxies are constituted by cloudlets of cold (< 10 K) H_2 and dense (> 10^7$ cm^{-3}) molecular gas. Such gas is extremely difficult to detect, because the classical tracers of molecular gas, CO and/or dust grains, have very low abundances and their emission is exceedingly weak. For this reason, the cloudlet hypothesis remains so far substantially unproven. In this Letter we propose a new method to probe the presence of cold H_2 clouds in galactic halos: the ground transition of ortho-H_2D^+ at 372 GHz. We discuss why the H_2D^+ is abundant under the physical conditions appropriate for the cloudlets, and present a chemical model that predicts the H_2D^+ abundance as function of four key parameters: gas density and metallicity, cosmic ray ionization rate and dust grain size. We conclude that current ground-based instruments might detect the ortho-H_2D^+ line emitted by the cloudlets halo, and prove, therefore, the existence of large quantities of dark baryonic matter around galaxies.

We study the effect of rotation during the collision between dust aggregates, in order to address a mismatch between previous model calculations of Brownian motion driven aggregation and experiments. We show that rotation during the collision does influence the shape and internal structure of the aggregates formed. The effect is limited in the ballistic regime when aggregates can be considered to move on straight lines during a collision. However, if the stopping length of an aggregate becomes smaller than its physical size, extremely elongated aggregates can be produced. We show that this effect may have played a role in the inner regions of the solar nebula where densities were high.

Recent observations of the ground state transition of HDO at 464 GHz towards the protoplanetary disk of DM Tau have detected the presence of water vapor in the regions just above the outer disk midplane (Ceccarelli et al 2005). In the absence of non-thermal desorption processes, water should be almost entirely frozen onto the grain mantles and HDO undetectable. In this Letter we present a chemical model that explores the possibility that the icy mantles are photo-desorbed by FUV (6eV < h nu < 13.6eV) photons. We show that the average Interstellar FUV field is enough to create a layer of water vapor above the disk midplane over the entire disk. Assuming a photo-desorption yield of 10^{-3}, the water abundance in this layer is predicted to be ~ 3 x 10^{-7} and the average H2O column density is ~ 1.6x 10^{15} cm^{-2}. The predictions are very weakly dependent on the details of the model, like the incident FUV radiation field, and the gas density in the disk. Based on this model, we predict a gaseous HDO/H2O ratio in DM Tau of ~1%. In addition, we predict the ground state transition of water at 557 GHz to be undetectable with ODIN and/or HSO-HIFI.

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.

We report the first detection of the ground transition of the deuterated water at 464 GHz in the young proto-planetary disk surrounding the solar type protostar DM Tau. The line is observed in absorption against the continuum from the cold dust in the disk midplane, with a line to continuum ratio close to unity. The observation implies that deuterated gaseous water is present, with a relatively large abundance ($\sim 3\times10^{-9}$), in the outer disk above the midplane, where the density is, within a factor ten, $\sim 10^6$ cm$^{-3}$ and the temperature is lower than about 25 K. In these conditions, the H$_2$O condensation timescale is much smaller than the DM Tau disk age, and, therefore, water should be fully frozen onto the grain mantles. We suggest that UV photons and/or X-rays sublimate part of the mantles re-injecting the ices into the gas phase. Even though there is currently no measurement of H$_2$O, we provide arguments that the HDO/H$_2$O ratio should be about 0.01 or larger, which would be hundreds of times larger than the values measured in Solar System objects. This suggests the need of strong caution in comparing and linking the HDO/H$_2$O in Solar System and star forming environments.

We present simulations of the interferometric visibilities of Herbig Ae star disks. We investigate whether interferometric measurements in the 10 micrometer atmospheric window are sensitive to the presence of an increased scale height at the inner disk edge, predicted by recent models. Furthermore, we investigate whether such measurements can discriminate between disks with a ``flaring'' geometry and disks with a ``flat'' geometry. We show that both these questions can be addressed, using measurements at a small number of appropriately chosen baselines. The classification of Herbig Ae stars in two groups, based on the appearance of the spectral energy distribution (SED), has been attributed to a difference in disk geometry. Sources with a group I SED would have a flaring outer disk geometry, whereas the disk of group II sources is proposed to be flat (or ``self-shadowed''). We show that this hypothesis can be tested using long-baseline interferometric measurements in the micrometer atmospheric window.

Probing the gas and dust in proto-planetary disks is central for understanding the process of planet formation. In disks surrounding solar type protostars, the bulk of the disk mass resides in the outer midplane, which is cold ($\leq$20 K), dense ($\geq 10^7$ cm$^{-3}$) and depleted of CO. Observing the disk midplane has proved, therefore, to be a formidable challenge. Ceccarelli et al. (2004) detected H$_2$D$^+$ emission in a proto-planetary disk and claimed that it probes the midplane gas. Indeed, since all heavy-elements bearing molecules condense out onto the grain mantles, the most abundant ions in the disk midplane are predicted to be H$_3^+$ and its isotopomers. In this article, we carry out a theoretical study of the chemical structure of the outer midplane of proto-planetary disks. Using a self-consistent physical model for the flaring disk structure, we compute the abundances of H$_3^+$ and its deuterated forms across the disk midplane. We also provide the average column densities across the disk of H$_3^+$, H$_2$D$^+$, HD$_2^+$ and D$_3^+$, and line intensities of the ground transitions of the ortho and para forms of H$_2$D$^+$ and HD$_2^+$ respectively. We discuss how the results depend on the cosmic ray ionization rate, dust-to-gas ratio and average grain radius, and general stellar/disk parameters. An important factor is the poorly understood freeze-out of N$_2$ molecules onto grains, which we investigate in depth. We finally summarize the diagnostic values of observations of the H$_3^+$ isotopomers.

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.

We present extended Spectral Energy Distributions (SEDs) of seven classical RV Tauri stars, using newly obtained submillimetre continuum measurements and Geneva optical photometry supplemented with literature data. The broad-band SEDs show a large IR excess with a black-body slope at long wavelengths in six of the seven stars, R Sct being the noticeable exception. This long wavelength slope is best explained assuming the presence of a dust component of large grains in the circumstellar material. We show that the most likely distribution of the circumstellar dust around the six systems is that the dust resides in a disc. Moreover, very small outflow velocities are needed to explain the presence of dust near the sublimation temperature and we speculate that the discs are Keplerian. The structure and evolution of these compact discs are as yet not understood but a likely prerequisite for their formation is that the dusty RV Tauri stars are binaries.

We model the process of dust coagulation in protoplanetary disks and calculate how it affects their observational appearance. Our model involves the detailed solution of the coagulation equation at every location in the disk. At regular time intervals we feed the resulting 3-D dust distribution functions into a continuum radiative transfer code to obtain spectral energy distributions. We find that, even if only the very basic -- and well understood -- coagulation mechanisms are included, the process of grain growth is much too quick to be consistent with infrared observations of T Tauri disks. Small grains are removed so efficiently that, long before the disk reaches an age of 1E6 years typical of T Tauri stars, the SED shows only very weak infrared excess. This is inconsistent with observed SEDs of most classical T Tauri stars. Small grains somehow need to be replenished, for instance by aggregate fragmentation through high-speed collisions. A very simplified calculation shows that when aggregate fragmentation is included, a quasi-stationary grain size distribution is obtained in which growth and fragmentation are in equilibrium. This quasi-stationary state may last 1E6 years or even longer, dependent on the circumstances in the disk, and may bring the time scales into the right regime. We conclude that a simple evolutionary scenario in which grains slowly grow from pristine 0.1 micron grains to larger grains over a period of a few Myr is most likely incorrect.

We analyze how the process of dust settling affects the spectral energy distribution and optical appearance of protoplanetary disks. Using simple analytic estimates on the one hand, and detailed 1+1-D models on the other hand, we show that, while the time scale for settling down to the equator may exceed the life time of the disk, it takes much less time for even small grains of 0.1 micron to settle down to a few pressure scale heights. This is often well below the original location of the disk's photosphere, and the disk therefore becomes effectively 'flatter'. If turbulent stirring is included, a steady state solution can be found, which is typically reached after a few times 1E5 years. Dependent on the strength of the turbulence, the shape of the disk in such a steady state can be either fully flaring, or flaring only up to a certain radius and self-shadowed beyond that radius. We show that these partly self-shadowed disks have a much weaker mid- to far-infrared flux than the fully flaring ones. We also show that these self-shadowed regions of the disk are very weak in resolved images of scattered light, in contrast to the fully flaring disks. From the calculations with compact grains it follows that, after about 1E6 years, most disks should be self-shadowed. The fact that some older disks are still observed with the characteristics of flaring disks therefore seems somewhat inconsistent with the time scales predicted by the settling model based on compact grains. This suggests that perhaps even the small grains (lesssim 0.1 micron) have a porous or fractal structure, or that the different geometries of observed disks is merely a reflection of the turbulent state of these disks.

Isolated Herbig Ae stars can be divided into two groups (Meeus et al. 2001): those with an almost flat spectral energy distribution in the mid-infrared (`group I'), and those with a strong decline towards the far-infrared (`group II'). In this paper we show that the group I vs. II distinction can be understood as arising from flaring vs. self-shadowed disks. We show that these two types of disks are natural solutions of the 2-D radiation-hydrostatic structure equations. Disks with high optical depth turn out to be flaring and have a strong far-IR emission, while disks with an optical depth below a certain threshold drop into the shadow of their own puffed-up inner rim and are weak in the far-IR. In spite of not having a directly irradiated surface layer, self-shadowed disks still display dust features in emission, in agreement with observations of group II sources. We propose an evolutionary scenario in which a disk starts out with a flaring shape (group I source), and then goes through the process of grain growth, causing the optical depth of the disk to drop and the disk to become self-shadowed (group II source). We show that this scenario predicts that the (sub-)millimeter slope of the disk changes from steep (small grains) to Rayleigh-Jeans-like (large grains) in the early stages of evolution, so that all group II sources are expected to have Rayleigh-Jeans-like slopes, while some group I sources may still have steep (sub-)millimeter slopes.

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.

In a separate paper (Decin et al 2003), we have re-examined the observations of IR excess obtained with the ISO satellite and discussed the ages of stars with excess. The amount of dust (measured by the luminosity fraction \fdust=L_\mathrm{IR}/L_{\star}) seen around main-sequence stars of different ages shows several interesting trends. To discuss these results in the context of a physical model, we develop in this paper an analytical model for the dust production in Vega-type systems. Previously it has been claimed that a powerlaw slope of about -2 in the diagram plotting amount of dust versus time could be explained by a simple collisional cascade. We show that such a cascade in fact results in a powerlaw \fdust\propto t^{-1} if the dust removal processes are dominated by collisions. A powerlaw \fdust\propto t^{-2} only results when the dust removal processes become dominated by Pointing-Robertson drag. This may be the case in the Kuiper Belt of our own solar system, but it is certainly not the case in any of the observed disks. A steeper slope can, however, be created by including continuous stirring into the models. We show that the existence of both young and old Vega-like systems with large amounts of dust (\fdust\simeq 10^{-3}) can be explained qualitatively by Kuiper-Belt-like structures with \emph{delayed stirring}. Finally, the absence of young stars with intermediate amounts of dust may be due to the fact that stirring due to planet formation may not be active in young low-mass disks. The considerations in this paper support the picture of simultaneous stirring and dust production proposed by Kenyon and Bromley (2002).

We study the time dependency of Vega-like excesses using infrared studies obtained with the imaging photopolarimeter ISOPHOT on board of ISO. We review the different studies published on this issue, and critically check and revise ages and fractional luminosities in the different samples. The conclusions of our study differ significantly from those obtained by other authors (e.g. Holland et al. 1998; Spangler et al. 2001) who suggested that there is a global power-law governing the amount of dust seen in debris disks as a function of time. Our investigations lead us to conclude that (i) for stars at most ages, a large spread in fractional luminosity occurs, but (ii) there are few very young stars with intermediate or small excesses; (iii) the maximum excess seen in stars of a given age is about a fractional luminosity of approximately 10^{-3}, independent of time; and (iv) Vega-like excess is more common in young stars than in old stars.

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.

We have investigated the spatial distribution, and the properties and chemical composition of the dust orbiting HD 100546. This system is remarkably different from other isolated Herbig Ae/Be stars in both the strength of the mid-IR excess and the composition of the circumstellar dust. We speculate that the formation and spatial distribution of the crystalline dust observed in the HD 100546 system may be linked to the formation of a proto-Jupiter in the disk around HD 100546. Such a proto-Jupiter could gravitationally stir the disk leading to a collisional cascade of asteroidal sized objects producing small crystalline grains, or it could cause shocks by tidal interaction with the disk which might produce crystalline dust grains through flash heating. As shown by Malfait et al. (1998), the infrared spectrum of HD 100546 is very similar to that of C/1995 O1 Hale-Bopp (Crovisier et al. 1997). Using an identical methodology, we have therefore also studied this solar system comet. Both objects have an almost identical grain composition, but with the important difference that the individual dust species in Hale-Bopp are in thermal contact with each other, while this is not the case in HD 100546. This suggests that if similar processes leading to the dust composition as seen in HD 100546 also occurred in our own solar system, that Hale-Bopp formed after the formation of one or more proto-gas giants.

We discuss spectral energy distributions of a sample of Herbig Ae/Be stars in the context of a passive irradiated disk model. The data have been presented earlier by Meeus et al (2001), and preliminary interpretations of these data were given in that paper. While the spectra of Herbig Ae stars all show similarities, there is significant variation between the spectra, in particular in the shape of the mid-IR rise and in the presence or absence of a silicate feature. We explore the hypothesis that all these different spectra can be interpreted as pure disk spectra without additional components. Using the model of Dullemond, Dominik and Natta (2001) we deduce the disk parameters of a number of the sources, and find that for a large fraction of investigated sources, satisfactory fits can be obtained. The derived model parameters show that some group Ia sources can only be fit with radially increasing surface densities, indicating the presence of depleted inner disk regions. The steep-sloped SEDs of group IIa sources can be fit with very compact disks, probably representing disks with collapsed outer regions. The largest difficulties arise from sources that do not show significant silicate emission features. Our attempts to explain these objects with a pure geometric effect are only partially successful. It seems that these stars indeed require a strong depletion of small silicate grains.

We present the Spectral Energy Distribution of HR 4049 based on literature data and new continuum measurements at 850 micron. The SED shows variable absorption in the UV, and a large IR excess, both caused by circumstellar dust. The shape of the IR excess from 1 micron all the way down to 850 micron can be nearly perfectly fitted with a single blackbody function at T ~ 1150 K or alternatively with a sum of blackbodies in a narrow temperature range. The energy emitted in this IR continuum radiation is about one-third of the stellar luminosity. We show that this blackbody radiation must be due to the presence of a circumbinary disk with a large height. This disk must also be gas-rich, in agreement with the observations of molecular bands in the ISO-SWS spectrum. We present two possible scenario's for explaining the shape and the intensity of the IR excess. The first scenario involves large grains (a >1 mm) that each radiate like a blackbody. The second scenario argues that the blackbody radiation is due to a very optically thick circumbinary disk. We investigate if such a disk would indeed produce blackbody radiation by presenting results from radiative transfer calculations. We further quantify the properties of such a disk and its stability in the framework of (hydro)dynamics, grain settling, radiation pressure and grain drift. The virtues and shortcomings of both models for the origin of the IR blackbody are discussed by contrasting them with other observations and assessing them in the framework of (binary) (post-)AGB evolution.

We analyse the optical and IR spectra, as well as the spectral energy distribution (UV to mm) of the candidate Herbig Ae star HD100453. This star is particular, as it shows an energy distribution similar to that of other isolated Herbig Ae/Be stars (HAEBEs), but unlike most of them, it does not have a silicate emission feature at 10 um, as is shown in Meeus (2001). We confirm the HAEBE nature of HD100453 through an analysis of its optical spectrum and derived location in the H-R diagram. The IR spectrum of HD100453 is modelled by an optically thin radiative transfer code, from which we derive constraints on the composition, grain-size and temperature distribution of the circumstellar dust. We show that it is both possible to explain the lack of the silicate feature as (1) a grain-size effect - lack of small silicate grains, and (2) a temperature effect - lack of small, hot silicates, as proposed by Dullemond (2001), and discuss both possibilities.

Focussing on preplanetary grains growth, we discuss the properties of dust aggregation driven by magnetic dipole forces. While there is no direct evidence for the existence of magnetic grains present in the solar nebula, there are reasons to assume they may have been present. We derive analytical expressions for the cross-section of two interacting dipoles. The effective cross section depends upon the strength of the magnetic dipoles and the initial velocities. For typical conditions the magnetic cross section is between 2 and 3 orders of magnitude larger than the geometric cross section. We study the growth dynamics of magnetic grains and find that the mass of the aggregates should increase with time like t^3.2 whereas Brownian motion growth behaves like t^2. A numerical tool is introduced which can be used to model dust aggregation in great detail, including the treatment of contact forces, aggregate restructuring processes and long-range forces. This tool is used to simulate collisions between magnetic grains or clusters and to validate the analytical cross-sections. The numerically derived cross section is in excellent agreement with the analytical expression. The numerical tool is also used to demonstrate that structural changes in the aggregates during collisions can be significant.

We report the discovery of a double ring structure in the waist of the nebula surrounding eta Carinae. The rings are detected in the mid-IR dust continuum at wavelengths of 7.9, 11.9, 12.9 and 20 micronmeter. The dust in the rings has a temperature of about 300 K. The orientation of the rings is inclined with respect to the axis of the homunculus by either 37 or 58 degrees. The central star is not in the projected centre of the structure defined by the two rings. This geometry is reminiscent of that seen in SN1987A and some planetary nebulae. We discuss several possible origins for this remarkable geometry and its orientation.

A model for irradiated dust disks around Herbig Ae stars is proposed. The model is based on the flaring disk model of Chiang & Goldreich (1997), but with the central regions of the disk removed. The inner rim of the disk is puffed up and is much hotter than the rest of the disk, because it is directly exposed to the stellar flux. If located at the dust evaporation radius, its reemitted flux produces a conspicuous bump in the SED which peaks at 2-3 micron. We propose that this emission is the explanation for the near-infrared bump observed in the SEDs of Herbig Ae stars. We study for which stellar parameters this bump would be observable, and find that it is the case for Herbig Ae stellar parameters but not for T-Tauri stars, confirming what is found from the observations. We also study the effects of the shadow cast by the inner rim over the rest of the flaring disk. The shadowed region can be quite large, and under some circumstances the entire disk may lie in the shadow. This shadowed region will be much cooler than an unshadowed flaring disk, since its only heating sources are radial radiative diffusion and possible indirect sources of irradiation. Under certain special circumstances the shadowing effect can suppress, or even completely eliminate, the 10 micron emission feature from the spectrum, which might explain the anomalous SEDs of some isolated Herbig Ae stars in the sample of Meeus et al. (2001). At much larger radii the disk emerges from the shadow, and continues as a flaring disk towards the outer edge. The complete model, including structure of the inner edge, shadowed region and the flared outer part, is described in detail in this paper, and we show examples of the general behavior of the model for varying parameters.

We have combined time dependent hydrodynamics with a two-fluid model for dust driven AGB winds. Our calculations include self-consistent gas chemistry, grain formation and growth, and a new implementation of the viscous momentum transfer between grains and gas. This allows us to perform calculations in which no assumptions about the completeness of momentum coupling are made. We derive new expressions to treat time dependent and non-equilibrium drift in a hydro code. Using a stationary state calculation for IRC +10216 as initial model, the time dependent integration leads to a quasi-periodic mass loss in the case where dust drift is taken into account. The time scale of the variation is of the order of a few hundred years, which corresponds to the time scale needed to explain the shell structure of the envelope of IRC +10216 and other AGB and post-AGB stars, which has been a puzzle since its discovery. No such periodicity is observed in comparison models without drift between dust and gas.

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

We have developed an implementation for the momentum transfer force in numerical two fluid hydrodynamics. This form of the frictional coupling between gas and grains is consistent with the microscopic interactions between the two components. The coupling force gives rise to a drift velocity of the grains with respect to the gas. We apply this mechanism to the outflow of (Post-) AGB objects. Our numerical hydrodynamics code calculates self consistently the dynamics of these outflows, as well as the nucleation and growth of grains and equilibrium chemistry of the gas. Grain nucleation and growth are processes that depend strongly on the rate of gas-grain collisions. Hence, the drift velocity becomes an important variable. The tight connection between grain chemistry and drift causes the system to become extremely sensitive to small changes in almost any parameter. This may be a cause for deviation from (spherical) symmetry and structure.