The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared interferometry. We provide reconstructed images of 15 objects selected from the Herbig AeBe survey carried out with PIONIER at the VLTI, using SPARCO. We find that 40% of the systems are centrosymmetric at the angular resolution of the observations. For the rest of the objects, we find evidence for asymmetric emission due to moderate-to-strong inclination of a disk-like structure for 30% of the objects and noncentrosymmetric morphology due to a nonaxisymmetric and possibly variable environment (30%). Among the systems with a disk-like structure, 20% show a resolved dust-free cavity. The image reconstruction process is a powerful tool to reveal complex disk inner rim morphologies. At the angular resolution reached by near-infrared interferometric observations, most of the images are compatible with a centrally peaked emission (no cavity). For the most resolved targets, image reconstruction reveals morphologies that cannot be reproduced by generic parametric models. Moreover, the nonaxisymmetric disks show that the spatial resolution probed by optical interferometers makes the observations of the near-infrared emission sensitive to temporal evolution with a time-scale down to a few weeks. The evidence of nonaxisymmetric emission that cannot be explained by simple inclination and radiative transfer effects requires alternative explanations, such as a warping of the inner disks. Interferometric observations can, therefore, be used to follow the evolution of the asymmetry of those disks at a sub-au scale.

Annular structures in proto-planetary discs, such as gaps and rings, are now ubiquitously found by high-resolution ALMA observations. Under the hypothesis that they are opened by planets, in this paper we investigate how the minimum planet mass needed to open a gap varies across different stellar host masses and distances from the star. The dependence on the stellar host mass is particularly interesting because, at least in principle, gap opening around low mass stars should be possible for lower mass planets, giving us a look into the young, low mass planet population. Using dusty hydrodynamical simulations, we find however the opposite behaviour, as a result of the fact that discs around low mass stars are geometrically thicker: gap opening around low mass stars can require more massive planets. Depending on the theoretical isochrone employed to predict the relationship between stellar mass and luminosity, the gap opening planet mass could also be independent of stellar mass, but in no case we find that gap opening becomes easier around low mass stars. This would lead to the expectation of a {\it lower} incidence of such structures in lower mass stars, since exoplanet surveys show that low mass stars have a lower fraction of giant planets. More generally, our study enables future imaging observations as a function of stellar mass to be interpreted using information on the mass vs. luminosity relations of the observed samples.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

Recent observations show that rings and gaps are ubiquitous in protoplanetary discs. These features are often interpreted as being due to the presence of planets; however, the effect of planetary migration on the observed morphology has not been investigated hitherto. In this work we investigate whether multiwavelength mm/submm observations can detect signatures of planet migration, using 2D dusty hydrodynamic simulations to model the structures generated by migrating planets and synthesising ALMA continuum observations at 0.85 and 3 mm. We identify three possible morphologies for a migrating planet: a slowly migrating planet is associated with a single ring outside the planet's orbit, a rapidly migrating planet is associated with a single ring inside the planet's orbit while a planet migrating at intermediate speed generates one ring on each side of the planet's orbit. We argue that multiwavelength data can distinguish multiple rings produced by a migrating planet from other scenarios for creating multiple rings, such as multiple planets or discs with low viscosity. The signature of migration is that the outer ring has a lower spectral index, due to larger dust grains being trapped there. Of the recent ALMA observations revealing protoplanetary discs with multiple rings and gaps, we suggest that Elias 24 is the best candidate for a planet migrating in the intermediate speed regime.

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.

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 Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.

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.

We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions.

We present spatially resolved ALMA images of CO J=3-2 emission from the protoplanetary disk around HD100546. We model the spatially-resolved kinematic structure of the CO emission. Assuming a velocity profile which prescribes a flat or flared emitting surface in Keplerian rotation, we uncover significant residuals with a peak of $\approx7\delta v$, where $\delta v = 0.21$ km s$^{-1}$ is the width of a spectral resolution element. The residuals reveal the possible presence of a severely warped and twisted inner disk extending to at most 100au. Adapting the model to include a misaligned inner gas disk with (i) an inclination almost edge-on to the line of sight, and (ii) a position angle almost orthogonal to that of the outer disk reduces the residuals to $< 3\delta v$. However, these findings are contrasted by recent VLT/SPHERE, MagAO/GPI, and VLTI/PIONIER observations of HD100546 that show no evidence of a severely misaligned inner dust disk down to spatial scales of $\sim 1$au. An alternative explanation for the observed kinematics are fast radial flows mediated by (proto)planets. Inclusion of a radial velocity component at close to free-fall speeds and inwards of $\approx 50$au results in residuals of $\approx 4 \delta v$. Hence, the model including a radial velocity component only does not reproduce the data as well as that including a twisted and misaligned inner gas disk. Molecular emission data at a higher spatial resolution (of order 10au) are required to further constrain the kinematics within $\lesssim 100$au. HD100546 joins several other protoplanetary disks for which high spectral resolution molecular emission shows that the gas velocity structure cannot be described by a purely Keplerian velocity profile with a universal inclination and position angle. Regardless of the process, the most likely cause is the presence of an unseen planetary companion. (Abridged)

Recent sub-millimeter observations show non-axisymmetric brightness distributions with a horseshoe-like morphology for more than a dozen transition disks. The most accepted explanation for the observed asymmetries is the accumulation of dust in large-scale vortices. Protoplanetary disks vortices can form by the excitation of Rossby-wave instability in the vicinity of a steep pressure gradient, which can develop at the edges of a giant planet carved gap or at the edges of an accretionally inactive zone. We studied the formation and evolution of vortices formed in these two distinct scenarios by means of two-dimensional locally isothermal hydrodynamic simulations. We found that the vortex formed at the edge of a planetary gap is short-lived, unless the disk is nearly inviscid. In contrast, the vortex formed at the outer edge of a dead zone is long-lived. The vortex morphology can be significantly different in the two scenarios: the vortex radial and azimuthal extensions are ~1.5 and ~3.5 times larger for the dead zone edge compared to gap models. In some particular cases, the vortex aspect ratios can be similar in the two scenarios, however, the vortex azimuthal extensions can be used to distinguish the vortex formation mechanisms. We calculate predictions for vortex observability in the sub-millimeter continuum with ALMA. We found that the azimuthal and radial extent of brightness asymmetry correlates with vortex formation process, within the limitations of alpha-viscosity prescription.

Spiral arms have been observed in nearly a dozen protoplanetary discs in near-infrared scattered light and recently also in the sub-millimetre continuum. While one of the most compelling explanations is that they are driven by planetary or stellar companions, in all but one cases such companions have not yet been detected and there is even ambiguity on whether the planet should be located inside or outside the spirals. Here we use 3D hydrodynamic simulations to study the morphology of spiral density waves launched by embedded planets taking into account the vertical temperature gradient, a natural consequence of stellar irradiation. Our simulations show that the pitch angle of the spirals in thermally stratified discs is the lowest in the disc mid-plane and increases towards the disc surface. We combine the hydrodynamic simulations with 3D radiative transfer calculations to predict that the pitch-angle of planetary spirals observed in the near-infrared is higher than in the sub-millimetre. We also find that in both cases the spirals converge towards the planet. This provides a new powerful observational method to determine if the perturbing planet is inside our outside the spirals, as well as map the thermal stratification of the disc.

The young star Elias 2-27 has recently been observed to posses a massive circumstellar disc with two prominent large-scale spiral arms. In this Letter we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modelling, synthetic ALMA imaging and an unsharped masking technique to explore three possibilities for the origin of the observed structures -- an undetected companion either internal or external to the spirals, and a self-gravitating disc. We find that a gravitationally unstable disc and a disc with an external companion can produce morphology that is consistent with the observations. In addition, for the latter, we find that the companion could be a relatively massive planetary mass companion (less than approximately 10 - 13 MJup) and located at large radial distances (between approximately 300 - 700 au). We therefore suggest that Elias 2-27 may be one of the first detections of a disc undergoing gravitational instabilities, or a disc that has recently undergone fragmentation to produce a massive companion.

In recent years an increasing number of observational studies have hinted at the presence of warps in protoplanetary discs, however a general comprehensive description of observational diagnostics of warped discs was missing. We performed a series of 3D SPH hydrodynamic simulations and combined them with 3D radiative transfer calculations to study the observability of warps in circumbinary discs, whose plane is misaligned with respect to the orbital plane of the central binary. Our numerical hydrodynamic simulations confirm previous analytical results on the dependence of the warp structure on the viscosity and the initial misalignment between the binary and the disc. To study the observational signatures of warps we calculate images in the continuum at near-infrared and sub-millimetre wavelengths and in the pure rotational transition of CO in the sub-millimetre. Warped circumbinary discs show surface brightness asymmetry in near-infrared scattered light images as well as in optically thick gas lines at sub-millimetre wavelengths. The asymmetry is caused by self-shadowing of the disc by the inner warped regions, thus the strength of the asymmetry depends on the strength of the warp. The projected velocity field, derived from line observations, shows characteristic deviations, twists and a change in the slope of the rotation curve, from that of an unperturbed disc. In extreme cases even the direction of rotation appears to change in the disc inwards of a characteristic radius. The strength of the kinematical signatures of warps decreases with increasing inclination. The strength of all warp signatures decreases with decreasing viscosity.

Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_r>=0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r~0.2, flaring to z/r~0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class.

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

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

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

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.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

This paper presents a self-consistent model for the evolution of gas produced in the debris disc of $\beta$ Pictoris. Our model proposes that atomic carbon and oxygen are created from the photodissociation of CO, which is itself released from volatile-rich bodies in the debris disc due to grain-grain collisions or photodesorption. While the CO lasts less than one orbit, the atomic gas evolves by viscous spreading resulting in an accretion disc inside the parent belt and a decretion disc outside. The temperature, ionisation fraction and population levels of carbon and oxygen are followed with the photodissociation region model Cloudy, which is coupled to a dynamical viscous $\alpha$ model. We present new gas observations of $\beta$ Pic, of C I observed with APEX and O I observed with Herschel, and show that these along with published C II and CO observations can all be explained with this new model. Our model requires a viscosity $\alpha$ > 0.1, similar to that found in sufficiently ionised discs of other astronomical objects; we propose that the magnetorotational instability is at play in this highly ionised and dilute medium. This new model can be tested from its predictions for high resolution ALMA observations of C I. We also constrain the water content of the planetesimals in $\beta$ Pic. The scenario proposed here might be at play in all debris discs and this model could be used more generally on all discs with C, O or CO detections.

Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself. Using the ALMA Science Verification data of HL Tau, we investigate the high velocity molecular gas being removed from the system as a result of the star formation process. We aim to place these motions in context with the optically detected jet, and the disk. With these high resolution ($\sim 1"$) ALMA observations of CO (J=1-0), we quantify the outwards motions of the molecular gas. We find evidence for a bipolar outwards flow, with an opening angle, as measured in the red-shifted lobe, starting off at 90$^\circ$, and narrowing to 60$^\circ$ further from the disk, likely because of magnetic collimation. Its outwards velocity, corrected for inclination angle is of order 2.4 km s$^{-1}$.

We investigate the minimum planet mass that produces observable signatures in infrared scattered light and submm continuum images and demonstrate how these images can be used to measure planet masses to within a factor of about two. To this end we perform multi-fluid gas and dust simulations of discs containing low mass planets, generating simulated observations at $1.65 \mu$m, $10 \mu$m and $850 \mu$m. We show that the minimum planet mass that produces a detectable signature is $\sim 15 M_\oplus$: this value is strongly dependent on disc temperature and changes slightly with wavelength (favouring the submm). We also confirm previous results that there is a minimum planet mass of $\sim 20 M_\oplus$ that produces a pressure maximum in the disc: only planets above this threshold mass generate a dust trap that can eventually create a hole in the submm dust. Below this mass, planets produce annular enhancements in dust outward of the planet and a reduction in the vicinity of the planet. These features are in steady state and can be understood in terms of variations in the dust radial velocity, imposed by the perturbed gas pressure radial profile, analogous to a traffic jam. We also show how planet masses can be derived from structure in scattered light and sub-mm images. We emphasise that simulations with dust need to be run over thousands of planetary orbits so as to allow the gas profile to achieve a steady state and caution against the estimation of planet masses using gas only simulations.

Scattered light images of transition discs in the near-infrared often show non-axisymmetric structures in the form of wide-open spiral arms in addition to their characteristic low-opacity inner gap region. We study self-gravitating discs and investigate the influence of gravitational instability on the shape and contrast of spiral arms induced by planet-disc interactions. Two-dimensional non-isothermal hydrodynamical simulations including viscous heating and a cooling prescription are combined with three-dimensional dust continuum radiative transfer models for direct comparison to observations. We find that the resulting contrast between the spirals and the surrounding disc in scattered light is by far higher for pressure scale height variations, i.e. thermal perturbations, than for pure surface density variations. Self-gravity effects suppress any vortex modes and tend to reduce the opening angle of planet-induced spirals, making them more tightly wound. If the disc is only marginally gravitationally stable with a Toomre parameter around unity, an embedded massive planet (planet-to-star mass ratio of $10^{-2}$) can trigger gravitational instability in the outer disc. The spirals created by this instability and the density waves launched by the planet can overlap resulting in large-scale, more open spiral arms in the outer disc. The contrast of these spirals is well above the detection limit of current telescopes.

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.

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

We report Herschel/PACS photometric observations at 70 {\mu}m and 160 {\mu}m of LRLL54361 - a suspected binary protostar that exhibits periodic (P=25.34 days) flux variations at shorter wavelengths (3.6 {\mu}m and 4.5 {\mu}m) thought to be due to pulsed accretion caused by binary motion. The PACS observations show unprecedented flux variation at these far-infrared wavelengths that are well cor- related with the variations at shorter wavelengths. At 70 {\mu}m the object increases its flux by a factor of six while at 160{\mu}m the change is about a factor of two, consistent with the wavelength dependence seen in the far-infrared spectra. The source is marginally resolved at 70 {\mu}m with varying FWHM. Deconvolved images of the sources show elongations exactly matching the outflow cavities traced by the scattered light observations. The spatial variations are anti-correlated with the flux variation indicating that a light echo is responsible for the changes in FWHM. The observed far-infrared flux variability indicates that the disk and en- velope of this source is periodically heated by the accretion pulses of the central source, and suggests that such long wavelength variability in general may provide a reasonable proxy for accretion variations in protostars.

The major physical processes responsible for shaping and sculpting pillars in the clouds surrounding massive stars (i.e. the `Pillars of Creation') are now being robustly incorporated into models quantifying the ionizing radiation from massive stars. The detailed gas dynamics within these pillars can now be compared with observations. Our goal is to quantify the gas dynamics in a pillar being sculpted by a nearby massive star. To do this, we use the CO, $^{13}$CO, and C$^{18}$O J=1-0 emission towards a pillar in the Vulpecula Rift. These data are a combination of CARMA and FCRAO observations providing high resolution ($\sim5''$) imaging of large scale pillar structures ($>100''$). We find that this cold ($\sim18$ K), low density material ($8\times10^3$ cm$^{-3}$) material is fragmenting on Jeans scales, has very low velocity dispersions ($\sim0.5$ km s$^{-1}$), and appears to be moving away from the ionizing source. We are able to draw direct comparisons with three models from the literature, and find that those with lower velocity dispersions best fit our data, although the dynamics of any one model do not completely agree with our observations. We do however, find that our observed pillar exhibits many of the characteristics expected from simulations.

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

The signatures of dust processing and grain growth - related to the formation of rocky planets - are easily seen in mid-infrared spectral features. One important diagnostic tool in this context is the silicate feature in the spectra of young stellar objects (YSO). The low-mass YSO, DG Tau shows unique temporal variations in its silicate feature. We conducted multi-epoch observations of DG Tau with the MID-Infrared Interferometric instrument of the Very Large Telescope Interferometer to obtain the spectra of the inner and outer disk regions in order to investigate where the previously reported variations of the silicate feature originate from. Here we present the preliminary results of the first two epochs of observations. We found that on a time-scale of two months, the source showed significant brightening. At the same time the mid-infrared emitting region expanded. While the identification of the silicate feature is difficult, our results qualitatively agree with the scenario explaining the varying silicate feature with dust lifted up above the disk.

Aims: The aim of this work is to study the structure of the protoplanetary disk surrounding the Herbig Ae star HD 163296. Methods: We have used high-resolution and high-sensitivity ALMA observations of the CO(3-2) emission line and the continuum at 850 microns, as well as the 3- dimensional radiative transfer code MCFOST to model the data presented in this work. Results: The CO(3-2) emission unveils for the first time at sub-millimeter frequencies the vertical structure details of a gaseous disk in Keplerian rotation, showing the back- and the front-side of a flared disk. Continuum emission at 850 microns reveals a compact dust disk with a 240 AU outer radius and a surface brightness profile that shows a very steep decline at radius larger than 125 AU. The gaseous disk is more than two times larger than the dust disk, with a similar critical radius but with a shallower radial profile. Radiative transfer models of the continuum data confirms the need for a sharp outer edge to the dust disk. The models for the CO(3-2) channel map require the disk to be slightly more geometrically thick than previous models suggested, and that the temperature at which CO gas becomes depleted (frozen-out) from the outer regions of the disk midplane is T < 20 K, in agreement with previous studies.

Disk winds have been postulated as a mechanism for angular momentum release in protostellar systems for decades. HD 163296 is a Herbig Ae star surrounded by a disk and has been shown to host a series of HH knots (HH 409) with bow shocks associated with the farthest knots. Here we present ALMA Science Verification data of CO J=2-1 and J=3-2 emission which are spatially coincident with the blue shifted jet of HH knots, and offset from the disk by -18.6 km/s. The emission has a double corkscrew morphology and extends more than 10" from the disk with embedded emission clumps coincident with jet knots. We interpret this double corkscrew as emission from material in a molecular disk wind, and that the compact emission near the jet knots is being heated by the jet which is moving at much higher velocities. We show that the J=3-2 emission is likely heavily filtered by the interferometer, but the J=2-1 emission suffers less due to the larger beam and measurable angular scales. Excitation analysis suggests temperatures exceeding 900 K in these compact features, with the wind mass, momentum and energy being of order 10^{-5} M_sun, 10^{-4} M_sun km/s and 10^{40} erg respectively. The high mass loss rate suggests that this star is dispersing the disk faster than it is funneling mass onto the star.

Context: It is hypothesized that low-mass young stellar objects undergo eruptive phases during their early evolution. The outburst of V1647 Ori between 2003 and 2006 offered a rare opportunity to investigate such an accretion event. Aims: By means of our interferometry observing campaign during this outburst, supplemented by other observations, we investigate the temporal evolution of the inner circumstellar structure of V1647 Ori We also study the role of the changing extinction in the brightening of the object and separate it from the accretional brightening. Methods: We observed V1647 Ori with MIDI/VLTI at two epochs in this outburst. First, during the slowly fading plateau phase (2005 March) and second, just before the rapid fading of the object (2005 September), which ended the outburst. We used the radiative transfer code MC3D to fit the interferometry data and the spectral energy distributions from five different epochs at different stages of the outburst. The comparison of these models allowed us to trace structural changes in the system on AU-scales. We also considered qualitative alternatives for the interpretation of our data. Results: We found that the disk and the envelope are similar to those of non-eruptive young stars and that the accretion rate varied during the outburst. We also found evidence for the increase of the inner radii of the circumstellar disk and envelope at the beginning of the outburst. Furthermore, the change of the interferometric visibilities indicates structural changes in the circumstellar material. We test a few scenarios to interpret these data. We also speculate that the changes are caused by the fading of the central source, which is not immediately followed by the fading of the outer regions. However, if the delay in the fading of the disk is responsible for the changes seen in the MIDI data, the effect should be confirmed by dynamical modeling.

Warm debris disks are a sub-sample of the large population of debris disks, and display excess emission in the mid-IR. Around solar-type stars, very few objects show emission features in mid-IR spectroscopic observations, that are attributed to small, warm silicate dust grains. The origin of this warm dust can possibly be explained either by a collision between several bodies or by transport from an outer belt. We present and analyse new far-IR Herschel/Pacs observations, supplemented by ground-based data in the mid-IR (VLTI/Midi and VLT/Visir), for one of these rare systems: the 10-16 Myr old debris disk around HD 113766 A. We improve an existing model to account for these new observations, and better constrain the spatial distribution of the dust and its composition. We underline the limitations of SED modelling and the need for spatially resolved observations. We find that the system is best described by an inner disk located within the first AU, well constrained by the Midi data, and an outer disk located between 9-13 AU. In the inner dust belt, our previous finding of Fe-rich crystalline olivine grains still holds. We do not observe time variability of the emission features over at least a 8 years time span, in a environment subjected to strong radiation pressure. The time stability of the emission features indicates that {\mu}m-sized dust grains are constantly replenished from the same reservoir, with a possible depletion of sub-{\mu}m-sized grains. We suggest that the emission features may arise from multi-composition aggregates. We discuss possible scenarios concerning the origin of the warm dust. The compactness of the innermost regions as probed by Midi, as well as the dust composition, suggest that we are witnessing the outcomes of (at least) one collision between partially differentiated bodies, in an environment possibly rendered unstable by terrestrial planetary formation.

We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micron HST/NICMOS data. While sub-millimeter studies suggested there is a dust-depleted cavity with r=0.35, we find scattered light as close as 0.1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks. We find two small-scaled spiral structures which asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h ~ 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5+3,-4 Mj, in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks, we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0.5. We reach 5-sigma contrasts limiting companions to planetary masses, 3-4 MJ at 1.0 and 2 MJ at 1.55 using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

(Abridged) Debris disks trace remnant reservoirs of leftover planetesimals in planetary systems. A handful of "warm" debris disks have been discovered in the last years, where emission in excess starts in the mid-infrared. An interesting subset within these warm debris disks are those where emission features are detected in mid-IR spectra, which points towards the presence of warm micron-sized dust grains. Given the ages of the host stars, the presence of these grains is puzzling, and questions their origin and survival in time. This study focuses on determining the mineralogy of the dust around 7 debris disks with evidence for warm dust, based on Spitzer/IRS spectroscopic data, in order to provide new insights into the origin of the dust grains. We present a new radiative transfer code dedicated to SED modeling of optically thin disks. We make use of this code on the SEDs of seven warm debris disks, in combination with recent laboratory experiments on dust optical properties. We find that most, if not all, debris disks in our sample are experiencing a transient phase, suggesting a production of small dust grains on relatively short timescales. From a mineralogical point of view, we find that enstatite grains have small abundances compared to crystalline olivine grains. The main result of our study is that we find evidences for Fe-rich crystalline olivine grains (Fe / [Mg + Fe] ~ 0.2) for several debris disks. This finding contrasts with studies of gas-rich protoplanetary disks. The presence of Fe-rich olivine grains, and the overall differences between the mineralogy of dust in Class II disks compared to debris disks suggest that the transient crystalline dust is of a new generation. We discuss possible crystallization routes to explain our results, and comment on the mechanisms that may be responsible for the production of small dust grains.

We present a comparison of archival Spitzer spectra of the strongly variable T Tauri EX Lupi, observed before and during its 2008 outburst. We analyze the mid-infrared emission from gas-phase molecules thought to originate in a circumstellar disk. In quiescence the emission shows a forest of H2O lines, highly excited OH lines, and the Q branches of the organics C2H2, HCN, and CO2, similar to the emission observed toward several T Tauri systems. The outburst emission shows instead remarkable changes: H2O and OH line fluxes increase, new OH, H2, and HI transitions are detected, and organics are no longer seen. We adopt a simple model of a single-temperature slab of gas in local thermal equilibrium, a common approach for molecular analyses of Spitzer spectra, and derive the excitation temperature, column density, and emitting area of H2O and OH. We show how model results strongly depend on the selection of emission lines fitted, and that spectrally-resolved observations are essential for a correct interpretation of the molecular emission from disks, particularly in the case of water. Using H2O lines that can be approximated as thermalized to a single temperature, our results are consistent with a column density decrease in outburst while the emitting area of warm gas increases. A rotation diagram analysis suggests that the OH emission can be explained with two temperature components, which remarkably increase in column density in outburst. The relative change of H2O and OH emission suggests a key role for UV radiation in the disk surface chemistry.

EX Lup is the prototype of the EXor class of eruptive young stars. These objects show optical outbursts which are thought to be related to runaway accretion onto the star. In a previous study we observed in-situ crystal formation in the disk of EX Lup during its latest outburst in 2008, making the object an ideal laboratory to investigate circumstellar crystal formation and transport. This outburst was monitored by a campaign of ground-based and Spitzer Space Telescope observations. Here we modeled the spectral energy distribution of EX Lup in the outburst from optical to millimeter wavelengths with a 2D radiative transfer code. Our results showed that the shape of the SED at optical wavelengths was more consistent with a single temperature blackbody than a temperature distribution. We also found that this single temperature component emitted 80-100 % of the total accretion luminosity. We concluded that a thermal instability, the most widely accepted model of EXor outbursts, was likely not the triggering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra revealed that the strength of all crystalline bands between 8 and 30 um increased right after the end of the outburst. Six months later, however, the crystallinity in the 10 um silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution of EXLup showed that, although vertical mixing should be stronger during the outburst than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) was required to reproduce the observed mid-infrared spectra.

Submillimetre images of transition discs are expected to reflect the distribution of the optically thin dust. Former observation of three transition discs LkHa330, SR21N, and HD1353444B at submillimetre wavelengths revealed images which cannot be modelled by a simple axisymmetric disc. We show that a large-scale anticyclonic vortex that develops where the viscosity has a large gradient (e.g., at the edge of the disc dead zone), might be accountable for these large-scale asymmetries. We modelled the long-term evolution of vortices being triggered by the Rossby wave instability. We found that a horseshoe-shaped (azimuthal wavenumber m=1) large-scale vortex forms by coalescing of smaller vortices within 5x10^4 yr, and can survive on the disc life-time (~5x10^6 yr), depending on the magnitude of global viscosity and the thickness of the viscosity gradient. The two-dimensional grid-based global disc simulations with local isothermal approximation and compressible-gas model have been done by the GPU version of hydrodynamic code FARGO (GFARGO). To calculate the dust continuum image at submillimetre wavelengths, we combined our hydrodynamical results with a 3D radiative transfer code. By the striking similarities of the calculated and observed submillimetre images, we suggest that the three transition discs can be modelled by a disc possessing a large-scale vortex formed near the disc dead zone edge. Since the larger dust grains (larger than mm in size) are collected in these vortices, the non-axisymmetric submillimetre images of the above transition discs might be interpreted as active planet and planetesimal forming regions situated far (> 50 AU) from the central stars.

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.

(Note: this is a shortened version of the original A&A-style structured abstract). The physical nature of the strong photometric variability of T Tau Sa, the more massive member of the Southern "infrared companion" to T Tau, has long been debated. Intrinsic luminosity variations due to variable accretion were originally proposed but later challenged in favor of apparent fluctuations due to time-variable foreground extinction. In this paper we use the timescale of the variability as a diagnostic for the underlying physical mechanism. Because the IR emission emerging from Sa is dominantly thermal emission from circumstellar dust at <=1500K, we can derive a minimum size of the region responsible for the time-variable emission. In the context of the variable foreground extinction scenario, this region must be (un-) covered within the variability timescale, which implies a minimum velocity for the obscuring foreground material. If this velocity supercedes the local Kepler velocity we can reject foreground extinction as a valid variability mechanism. The variable accretion scenario allows for shorter variability timescales since the variations in luminosity occur on much smaller scales, essentially at the surface of the star, and the disk surface can react almost instantly on the changing irradiation with a higher or lower dust temperature and according brightness. We have detected substantial variations at long wavelengths in T Tau S: +26% within four days at 12.8 micron. We show that this short-term variability cannot be due to variable extinction and instead must be due to variable accretion. Using a radiative transfer model of the Sa disk we show that variable accretion can in principle also account for the much larger (several magnitude) variations observed on timescales of several years. For the long-term variability, however, also variable foreground extinction is a viable mechanism.

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

We present IRS spectra and revised MIPS photometry for the 18 members of the Eta Chamaeleontis cluster. Aged 8 Myr, the Eta Cha cluster is one of the few nearby regions within the 5-10 Myr age range, during which the disk fraction decreases dramatically and giant planet formation must come to an end. For the 15 low-mass members, we measure a disk fraction ~50%, high for their 8 Myr age, and 4 of the 8 disks lack near-IR excesses, consistent with the empirical definition of "transition'' disks. Most of the disks are comparable to geometrically flat disks. The comparison with regions of different ages suggests that at least some of the "transition" disks may represent the normal type of disk around low-mass stars. Therefore, their flattened structure and inner holes may be related to other factors (initial masses of the disk and the star, environment, binarity), rather than to pure time evolution. We analyze the silicate dust in the disk atmosphere, finding moderate crystalline fractions (~10-30%) and typical grain sizes ~1-3 micron, without any characteristic trend in the composition. These results are common to other regions of different ages, suggesting that the initial grain processing occurs very early in the disk lifetime (<1 Myr). Large grain sizes in the disk atmosphere cannot be used as a proxy for age, but are likely related to higher disk turbulence. The dust mineralogy varies between the 8-12micron and the 20-30 micron features, suggesting high temperature dust processing and little radial mixing. Finally, the analysis of IR and optical data on the B9 star Eta Cha reveals that it is probably surrounded by a young debris disk with a large inner hole, instead of being a classical Be star.

The spectral region around 10 micrometer, showing prominent dust emission bands, is commonly used to derive the chemical composition of protoplanetary dust. Different analysis methods have been proposed for this purpose, but so far, no comparative test has been performed to test their validity. We calculated model spectra of disk models with different geometries and central sources, using a 2D radiative transfer code. These spectra were then fitted in a blind test using four different spectral decomposition methods. We studied the effect of disk structure (flared vs. flat), inclination angle, size of the inner disk hole and stellar luminosity on the fitted chemical composition. Our results show that the derived dust compositions by all methods deviate systematically from the real chemical composition. Out of the four tested spectral decomposition methods, our new two-layer temperature distribution method, differs the least from the input dust composition and the results show the weakest systematic effects. The reason for the deviations of the results given by other methods lies in their simplifying assumptions (e.g. single average grain temperature or one component continuum). We also tested the influence of different noise levels on the results of the spectral decomposition methods. We find that, for ground-based observations (8-13 micrometer) the expected uncertainty in the value of the crystallinity is about 11% for a signal-to-noise ratio of 100, while for space-based observations (7-17 micrometer) the uncertainty is about 5%. On the basis of our results, we propose a recipe for the analysis and interpretation of dust spectroscopy data in the mid-infrared which should be especially valuable for analysing Spitzer data and ground-based infrared spectroscopy data in the 10 micrometer window.

(abridged) We present Spitzer infrared observations to constrain disc and dust evolution in young T Tauri stars in MBM 12, a star-forming cloud at high latitude with an age of 2 Myr and a distance of 275 pc. The region contains 12 T Tauri systems, with primary spectral types between K3 and M6; 5 are weak-line and the rest classical T Tauri stars. We first use MIPS and literature photometry to compile spectral energy distributions for each of the 12 members in MBM 12, and derive their IR excesses. The IRS spectra are analysed with the newly developed two-layer temperature distribution (TLTD) spectral decomposition method. For the 7 T Tauri stars with a detected IR excess, we analyse their solid-state features to derive dust properties such as mass-averaged grain size, composition and crystallinity. We find a spatial gradient in the forsterite to enstatite range, with more enstatite present in the warmer regions. The fact that we see a radial dependence of the dust properties indicates that radial mixing is not very efficient in the discs of these young T Tauri stars. The SED analysis shows that the discs in MBM 12, in general, undergo rapid inner disc clearing, while the binary sources have faster discevolution. The dust grains seem to evolve independently from the stellar properties, but are mildly related to disc properties such as flaring and accretion rates.

We present optical, IR and millimeter observations of the solar-type star 13-277, also known as GM Cep, in the 4 Myr-old cluster Tr 37. GM Cep experiences rapid magnitude variations of more than 2 mag at optical wavelengths. We explore the causes of the variability, which seem to be dominated by strong increases in the accretion, being similar to EX-or episodes. The star shows high, variable accretion rates (up to ~10$^{-6}$ Msun/yr), signs of powerful winds, and it is a very fast rotator (Vsini~43 km/s). Its strong mid-IR excesses reveal a very flared disk and/or a remnant envelope, most likely out of hydrostatic equilibrium. The 1.3 millimeter fluxes suggest a relatively massive disk (Mdisk~0.1 Msun). Nevertheless, the millimeter mass is not enough to sustain increased accretion episodes over large timescales, unless the mass is underestimated due to significant grain growth. We finally explore the possibility of GM Cep having a binary companion, which could trigger disk instabilities producing the enhanced accretion episodes.

(abridged) We present spectroscopy data for a sample of 14 FUors and 2 TTauri stars observed with the Spitzer Space Telescope or with the Infrared Space Observatory (ISO). Based on the appearance of the 10 micron silicate feature we define 2 categories of FUors. Objects showing the silicate feature in absorption (Category 1) are still embedded in a dusty and icy envelope. The shape of the 10 micron silicate absorption bands is compared to typical dust compositions of the interstellar medium and found to be in general agreement. Only one object (RNO 1B) appears to be too rich in amorphous pyroxene dust, but a superposed emission feature can explain the observed shape. We derive optical depths and extinction values from the silicate band and additional ice bands at 6.0, 6.8 and 15.2 micron. In particular the analysis of the CO_2 ice band at 15.2 micron allows us to search for evidence for ice processing and constrains whether the absorbing material is physically linked to the central object or in the foreground. For objects showing the silicate feature in emission (Category 2), we argue that the emission comes from the surface layer of accretion disks. Analyzing the dust composition reveals that significant grain growth has already taken place within the accretion disks, but no clear indications for crystallization are present. We discuss how these observational results can be explained in the picture of a young, and highly active accretion disk. Finally, a framework is proposed as to how the two categories of FUors can be understood in a general paradigm of the evolution of young, low-mass stars. Only one object (Parsamian 21) shows PAH emission features. Their shapes, however, are often seen toward evolved stars and we question the object's status as a FUor and discuss other possible classifications.

We investigate the long-term optical-infrared variability of SV Cep, and explain it in the context of an existing UX Ori (UXOR) model. A 25-month monitoring programme was completed with the Infrared Space Observatory in the 3.3-100 um wavelength range. Following a careful data reduction, the infrared light curves were correlated with the variations of SV Cep in the V-band. A remarkable correlation was found between the optical and the far-infrared light curves. In the mid-infrared regime the amplitude of variations is lower, with a hint for a weak anti-correlation with the optical changes. In order to interpret the observations, we modelled the spectral energy distribution of SV Cep assuming a self-shadowed disc with a puffed-up inner rim, using a 2-dimensional radiative transfer code. We found that modifying the height of the inner rim, the wavelength-dependence of the long-term optical-infrared variations is well reproduced, except the mid-infrared domain. The origin of variation of the rim height might be fluctuation in the accretion rate in the outer disc. In order to model the mid-infrared behaviour we tested to add an optically thin envelope to the system, but this model failed to explain the far-infrared variability. Infrared variability is a powerful tool to discriminate between models of the circumstellar environment. The proposed mechanism of variable rim height may not be restricted to UXOR stars, but might be a general characteristic of intermediate-mass young stars.