To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed. In this study we present the first result of the DESTINYS survey (Disk Evolution Study Through Imaging of Nearby Young Stars). DESTINYS is an ESO/SPHERE large program that aims at studying disk evolution in scattered light, mainly focusing on a sample of low-mass stars (<1$M_\odot$) in nearby (~200 pc) star-forming regions. In this particular study we present the observations of the ET Cha (RECX 15) system, a nearby 'old' classical T Tauri star (5-8 Myr, ~100 pc), which is still strongly accreting. We use SPHERE/IRDIS in H-band polarimetric imaging mode to obtain high contrast images of the ET Cha system to search for scattered light from the circumstellar disk as well as thermal emission from close companions. We additionally employ VLT/NACO total intensity archival data taken in 2003. We report here the discovery of a low-mass (sub)stellar companion with SPHERE/IRDIS to ET Cha. We are estimating the mass of this new companion based on photometry. Depending on the system age it is a 5 Myr, 50 $M_{Jup}$ brown dwarf or an 8 Myr, 0.10 $M_\odot$ M-type pre-main-sequence star. We explore possible orbital solutions and discuss the recent dynamic history of the system. Independent of the precise companion mass we find that the presence of the companion likely explains the small size of the disk around ET Cha. The small separation of the binary pair indicates that the disk around the primary component is likely clearing from the outside in, explaining the high accretion rate of the system.

More than a dozen young stars host spiral arms in their surrounding protoplanetary disks. The excitation mechanisms of such arms are under debate. The two leading hypotheses -- companion-disk interaction and gravitational instability (GI) -- predict distinct motion for spirals. By imaging the MWC 758 spiral arm system at two epochs spanning ${\sim}5$ yr using the SPHERE instrument on the Very Large Telescope (VLT), we test the two hypotheses for the first time. We find that the pattern speeds of the spirals are not consistent with the GI origin. Our measurements further evince the existence of a faint "missing planet" driving the disk arms. The average spiral pattern speed is $0.\!^\circ22\pm0.\!^\circ03$ yr$^{-1}$, pointing to a driver at $172_{-14}^{+18}$ au around a $1.9$ $M_\odot$ central star if it is on a circular orbit. In addition, we witness time varying shadowing effects on a global scale that are likely originated from an inner disk.

There is mounting evidence that the composition and structure of planetary systems are intimately linked to their birth environments. During the past decade, several spectral surveys probed the chemistry of the earliest stages of star formation and of late planet-forming disks. However, very little is known about the chemistry of intermediate protostellar stages, i.e. Class I Young Stellar Objects (YSOs), where planet formation may have already begun. We present here the first results of a 3mm spectral survey performed with the IRAM-30m telescope to investigate the chemistry of a sample of seven Class I YSOs located in the Taurus star-forming region. These sources were selected to embrace the wide diversity identified for low-mass protostellar envelope and disk systems. We present detections and upper limits of thirteen small ($N_{\rm atoms}\leq3$) C, N, O, and S carriers - namely CO, HCO$^+$, HCN, HNC, CN, N$_2$H$^+$, C$_2$H, CS, SO, HCS$^+$, C$_2$S, SO$_2$, OCS - and some of their D, $^{13}$C, $^{15}$N, $^{18}$O, $^{17}$O, and $^{34}$S isotopologues. Together, these species provide constraints on gas-phase C/N/O ratios, D- and $^{15}$N-fractionation, source temperature and UV exposure, as well as the overall S-chemistry. We find substantial evidence of chemical differentiation among our source sample, some of which can be traced back to Class I physical parameters, such as the disk-to-envelope mass ratio (proxy for Class I evolutionary stage), the source luminosity, and the UV-field strength. Overall, these first results allow us to start investigating the astrochemistry of Class I objects, however, interferometric observations are needed to differentiate envelope versus disk chemistry.

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.

RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc.

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

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.

To characterize circumstellar systems in high contrast imaging, the fundamental step is to construct a best point spread function (PSF) template for the non-circumstellar signals (i.e., star light and speckles) and separate it from the observation. With existing PSF construction methods, the circumstellar signals (e.g., planets, circumstellar disks) are unavoidably altered by over-fitting and/or self-subtraction, making forward modeling a necessity to recover these signals. We present a forward modeling--free solution to these problems with data imputation using sequential non-negative matrix factorization (DI-sNMF). DI-sNMF first converts this signal separation problem to a "missing data" problem in statistics by flagging the regions which host circumstellar signals as missing data, then attributes PSF signals to these regions. We mathematically prove it to have negligible alteration to circumstellar signals when the imputation region is relatively small, which thus enables precise measurement for these circumstellar objects. We apply it to simulated point source and circumstellar disk observations to demonstrate its proper recovery of them. We apply it to Gemini Planet Imager (GPI) K1-band observations of the debris disk surrounding HR 4796A, finding a tentative trend that the dust is more forward scattering as the wavelength increases. We expect DI-sNMF to be applicable to other general scenarios where the separation of signals is needed.

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.

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.

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa).

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

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.

Debris disks are the natural by-products of the planet formation process. Scattered or polarized light observations are mostly sensitive to small dust grains that are released from the grinding down of bigger planetesimals. High angular resolution observations at optical wavelengths can provide key constraints on the radial and azimuthal distribution of the small dust grains. These constraints can help us better understand where most of the dust grains are released upon collisions. We present SPHERE/ZIMPOL observations of the debris disk around HR 4796 A, and model the radial profiles along several azimuthal angles of the disk with a code that accounts for the effect of stellar radiation pressure. This enables us to derive an appropriate description for the radial and azimuthal distribution of the small dust grains. Even though we only model the radial profiles along (or close to) the semi-major axis of the disk, our best-fit model is not only in good agreement with our observations but also with previously published datasets (from near-IR to sub-mm wavelengths). We find that the reference radius is located at $76.4\pm0.4$ au, and the disk has an eccentricity of $0.076_{-0.010}^{+0.016}$, with the pericenter located on the front side of the disk (north of the star). We find that small dust grains must be preferentially released near the pericenter to explain the observed brightness asymmetry. Even though parent bodies spend more time near the apocenter, the brightness asymmetry implies that collisions happen more frequently near the pericenter of the disk. Our model can successfully reproduce the shape of the outer edge of the disk, without having to invoke an outer planet shepherding the debris disk. With a simple treatment of the effect of the radiation pressure, we conclude that the parent planetesimals are located in a narrow ring of about $3.6$ au in width.

Context. Direct imaging of debris discs gives important information about their nature, their global morphology, and allows us to identify specific structures possibly in connection with the presence of gravitational perturbers. It is the most straightforward technique to observe planetary systems as a whole. Aims. We present the first resolved images of the debris disc around the young F-type star HD 160305, detected in scattered light using the VLT/SPHERE instrument in the near infrared. Methods. We used a post-processing method based on angular differential imaging and synthetic images of debris discs produced with a disc modelling code (GRaTer) to constrain the main characteristics of the disc around HD 160305. All of the point sources in the field of the IRDIS camera were analysed with an astrometric tool to determine whether they are bound objects or background stars. Results. We detect a very inclined (~ 82{\deg}) ring-like debris disc located at a stellocentric distance of about 86au (deprojected width ~27 au). The disc displays a brightness asymmetry between the two sides of the major axis, as can be expected from scattering properties of dust grains. We derive an anisotropic scattering factor g>0.5. A second right-left asymmetry is also observed with respect to the minor axis. We measure a surface brightness ratio of 0.73 $\pm$ 0.18 between the bright and the faint sides. Because of the low signal-to-noise ratio (S/N) of the images we cannot easily discriminate between several possible explanations for this left-right asymmetry, such as perturbations by an unseen planet, the aftermath of the breakup of a massive planetesimal, or the pericenter glow effect due to an eccentric ring. Two epochs of observations allow us to reject the companionship hypothesis for the 15 point sources present in the field.

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed.

In a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and our new Y-band images from the Very Large Telescope SPHERE instrument. We use the knowledge of the disc's radial extent inferred from ALMA observations and the grain size distribution found by SED fitting to generate semi-dynamical dust models of the disc. We compare the models to scattered light and thermal emission data and find that a disc with a maximum of the surface density at 110~au and shallow edges can describe both thermal emission and scattered light observations. This suggests that grains close to the blow-out limit and large grains stem from the same planetesimal population and are mainly influenced by radiation pressure. The influence of inwards transport processes could not be analysed in this study.

Planetary systems hold the imprint of the formation and of the evolution of planets especially at young ages, and in particular at the stage when the gas has dissipated leaving mostly secondary dust grains. The dynamical perturbation of planets in the dust distribution can be revealed with high-contrast imaging in a variety of structures. SPHERE, the high-contrast imaging device installed at the VLT, was designed to search for young giant planets in long period, but is also able to resolve fine details of planetary systems at the scale of astronomical units in the scattered-light regime. As a young and nearby star, NZ Lup was observed in the course of the SPHERE survey. A debris disk had been formerly identified with HST/NICMOS. We observed this system in the near-infrared with the camera in narrow and broad band filters and with the integral field spectrograph. High contrasts are achieved by the mean of pupil tracking combined with angular differential imaging algorithms. The high angular resolution provided by SPHERE allows us to reveal a new feature in the disk which is interpreted as a superimposition of two belts of planetesimals located at stellocentric distances of $\sim$85 and $\sim$115\,au, and with a mutual inclination of about 5$\degb$. Despite the very high inclination of the disk with respect to the line of sight, we conclude that the presence of a gap, that is, a void in the dust distribution between the belts, is likely. We discuss the implication of the existence of two belts and their relative inclination with respect to the presence of planets.

Context. The 51 Eridani system harbors a complex architecture with its primary star forming a hierarchical system with the binary GJ 3305AB at a projected separation of 2000 au, a giant planet orbiting the primary star at 13 au, and a low-mass debris disk around the primary star with possibly a cold component and a warm component inferred from the spectral energy distribution. Aims. We aim to better constrain the orbital parameters of the known giant planet. Methods. We monitored the system over three years from 2015 to 2018 with the VLT/SPHERE exoplanet imaging instrument. Results. We measure an orbital motion for the planet of ~130 mas with a slightly decreasing separation (~10 mas) and find a hint of curvature. This potential curvature is further supported at 3$\sigma$ significance when including literature GPI astrometry corrected for calibration systematics. Fits of the SPHERE and GPI data using three complementary approaches provide broadly similar results. The data suggest an orbital period of 32$^{+17}_{-9}$ yr (i.e. 12$^{+4}_{-2}$ au in semi-major axis), an inclination of 133$^{+14}_{-7}$ deg, an eccentricity of 0.45$^{+0.10}_{-0.15}$, and an argument of periastron passage of 87$^{+34}_{-30}$ deg [mod 180 deg]. The time at periastron passage and the longitude of node exhibit bimodal distributions because we do not detect yet if the planet is accelerating or decelerating along its orbit. Given the inclinations of the planet's orbit and of the stellar rotation axis (134-144 deg), we infer alignment or misalignment within 18 deg for the star-planet spin-orbit. Further astrometric monitoring in the next 3-4 years is required to confirm at a higher significance the curvature in the planet's motion, determine if the planet is accelerating or decelerating on its orbit, and further constrain its orbital parameters and the star-planet spin-orbit.

Context: The complex system HD 100453 AB with a ring-like circumprimary disk and two spiral arms, one of which is pointing to the secondary, is a good laboratory to test spiral formation theories. Aims: To quantify the interaction of HD 100453 B with the circumprimary disk. Methods: Using ALMA band 6 dust continuum and CO isotopologue observations we study the HD 100453 AB system with a spatial resolution of 0.09" x 0.17" at 234 GHz. We use SPH simulations and orbital fitting to investigate the tidal influence of the companion on the disk. Results: We resolve the continuum emission around HD 100453 A into a disk between 0.22" and 0.40" with an inclination of 29.5 deg. and a position angle of 151.0 deg., an unresolved inner disk, and excess mm emission cospatial with the northern spiral arm which was previously detected using scattered light observations. We also detect CO emission from 7 au (well within the disk cavity) out to 1.10", i.e., overlapping with HD 100453 B at least in projection. The outer CO disk PA and inclination differ by up to 10 deg. from the values found for the inner CO disk and the dust continuum emission, which we interpret as due to gravitational interaction with HD 100453 B. Both the spatial extent of the CO disk and the detection of mm emission at the same location as the northern spiral arm are in disagreement with the previously proposed near co-planar orbit of HD 100453 B. Conclusions: We conclude that HD 100453 B has an orbit that is significantly misaligned with the circumprimary disk. Because it is unclear whether such an orbit can explain the observed system geometry we highlight an alternative scenario that explains all detected disk features where another, (yet) undetected, low mass close companion within the disk cavity, shepherds a misaligned inner disk whose slowly precessing shadows excite the spiral arms.

We present 3D smoothed particle hydrodynamics simulations of protoplanetary discs undergoing a flyby by a stellar perturber on a parabolic orbit lying in a plane inclined relative to the disc mid-plane. We model the disc as a mixture of gas and dust, with grains ranging from 1 {\mu}m to 10 cm in size. Exploring different orbital inclinations, periastron distances and mass ratios, we investigate the disc dynamical response during and after the flyby. We find that flybys induce evolving spiral structure in both gas and dust which can persist for thousands of years after periastron. Gas and dust structures induced by the flyby differ because of drag-induced effects on the dust grains. Variations in the accretion rate by up to an order of magnitude occur over a time-scale of order 10 years or less, inducing FU Orionis-like outbursts. The remnant discs are truncated and warped. The dust disc is left more compact than the gas disc, both because of disc truncation and accelerated radial drift of grains induced by the flyby.

Magnetospheric accretion has been thoroughly studied in young stellar systems with full non-evolved accretion disks, but it is poorly documented for transition disk objects with large inner cavities. We aim at characterizing the star-disk interaction and the accretion process onto the central star of LkCa 15, a transition disk system with an inner dust cavity. We obtained quasi-simultaneous photometric and spectropolarimetric observations of the system over several rotational periods. We analyzed the system light curve, as well as changes in spectral continuum and line profile to derive the properties of the accretion flow from the edge of the inner disk to the central star. We also derived magnetic field measurements at the stellar surface. We find that the system exhibits magnetic, photometric, and spectroscopic variability with a period of about 5.70 days. The light curve reveals a periodic dip, which suggests the presence of an inner disk warp that is located at the corotation radius at about 0.06 au from the star. Line profile variations and veiling variability are consistent with a magnetospheric accretion model where the funnel flows reach the star at high latitudes. This leads to the development of an accretion shock close to the magnetic poles. All diagnostics point to a highly inclined inner disk that interacts with the stellar magnetosphere. The spectroscopic and photometric variability of LkCa 15 is remarkably similar to that of AA Tau, the prototype of periodic dippers. We therefore suggest that the origin of the variability is a rotating disk warp that is located at the inner edge of a highly inclined disk close to the star. This contrasts with the moderate inclination of the outer transition disk seen on the large scale and thus provides evidence for a significant misalignment between the inner and outer disks of this planet-forming transition disk system.

We present a methodology to interpret observations of protoplanetary discs where a flyby, also called a tidal encounter, is suspected. In case of a flyby, protoplanetary discs can be significantly disturbed. The resulting dynamical and kinematical signatures can last for several thousands of years after the flyby and hence deeply affect the evolution of the disc. These effects are stronger for closer encounters and more massive perturbers. For the very same flyby parameters, varying the inclination of the perturber's orbit produces a broad range of disc structures: spirals, bridges, warps and cavities. We study this kind of features both in the gas and in the dust for grains ranging from 1 {\mu}m to 10 cm in size. Interestingly, the dust exhibits a different dynamical behaviour compared to the gas because of gas-drag effects. Finally, flybys can also trigger high accretion events in the disc-hosting star, readily similar to FU Orionis-type outbursts. All this information can be used to infer the flyby parameters from an incomplete set of observations at different wavelengths. Therefore, the main scope of our flyby scene investigation (FSI) methodology is to help to interpret recent "puzzling" disc observations.

We present ALMA images of the sub-mm continuum polarisation and spectral index of the protoplanetary ringed disk HD163296. The polarisation fraction at 870{\mu}m is measured to be ~0.9% in the central core and generally increases with radius along the disk major axis. It peaks in the gaps between the dust rings, and the largest value (~4%) is found between rings 1 and 2. The polarisation vectors are aligned with the disk minor axis in the central core, but become more azimuthal in the gaps, twisting by up to +/-9degrees in the gap between rings 1 and 2. These general characteristics are consistent with a model of self-scattered radiation in the ringed structure, without requiring an additional dust alignment mechanism. The 870/1300{\mu}m dust spectral index exhibits minima in the centre and the inner rings, suggesting these regions have high optical depths. However, further refinement of the dust or the disk model at higher resolution is needed to reproduce simultaneously the observed degree of polarisation and the low spectral index.

The disk-outflow connection is thought to play a key role in extracting excess angular momentum from a forming protostar. We present ALMA 0.25" angular resolution observations - in the dust continuum at 1.33 mm and of the molecular line transitions of $^{12}$CO(2-1) and $^{13}$CO(2-1) - of the circumstellar disk and outflow around the T Tauri star HH30, a rare and beautiful example of a pre-main sequence star exhibiting a flared edge-on disk, an optical jet, and a CO molecular outflow. The 1.3 mm continuum emission shows a remarkable elongated morphology along PA=31.2{\deg}+/-0.1{\deg} that has a constant brightness out to a radius of r=75 au. The emission is marginally resolved in the transverse direction, implying an intrinsic vertical width $\leq$24 au and an inclination to the line-of-sight i$\geq$84.8{\deg}. The monopolar outflow, detected in $^{12}$CO, arises from the north-eastern face of the disk from a disk radius r$\leq$22 au and extends up to 5" (or 700 au) above the disk plane. We derive a lower limit to the total mass of the CO cavity/outflow of 1.7E-5 M$_\odot$. The CO cavity morphology is that of a hollow cone with semi-opening angle ~35{\deg}. The derived kinematics are consistent with gas flowing along the conical surface with a constant velocity of 9.3+/-0.7 km/s. We detect small rotation signatures (Vphi x sin(i) $\in$ [0.1;0.5] km/s) in the same sense as the underlying circumstellar disk. From these rotation signatures we infer an average specific angular momentum of the outflow of 38+/-15 au km/s at altitudes z<250 au. We also report the detection of small amplitude wiggling (1.2{\deg}) of the CO axis around an average inclination to the line of sight of i=91{\deg}. The derived morphology and kinematics of the CO cavity are compatible with expectations from a slow disk wind, originating either through photo-evaporation or magneto-centrifugal processes.

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.

In the present study we aim to investigate the circumstellar environment of the spectroscopic binary T Tauri star CS Cha. From unresolved mid- to far-infrared photometry it is predicted that CS Cha hosts a disk with a large cavity. In addition, SED modeling suggests significant dust settling, pointing towards an evolved disk that may show signs of ongoing or completed planet formation. We observed CS Cha with the high contrast imager VLT/SPHERE in polarimetric differential imaging mode to resolve the circumbinary disk in near infrared scattered light. These observations were followed-up by VLT/NACO L-band observations and complemented by archival VLT/NACO K-band and HST/WFPC2 I-band data. We resolve the compact circumbinary disk around CS Cha for the first time in scattered light. We find a smooth, low inclination disk with an outer radius of $\sim$55 au (at 165 pc). We do not detect the inner cavity but find an upper limit for the cavity size of $\sim$15 au. Furthermore, we find a faint co-moving companion with a projected separation of 210 au from the central binary outside of the circumbinary disk. The companion is detected in polarized light and shows an extreme degree of polarization (13.7$\pm$0.4 \% in J-band). The companion's J- and H-band magnitudes are compatible with masses of a few M$_\mathrm{Jup}$. However, K-, L- and I-band data draw this conclusion into question. We explore with radiative transfer modeling whether an unresolved circum-companion disk can be responsible for the high polarization and complex photometry. We find that the set of observations is best explained by a heavily extincted low mass ($\sim 20 \mathrm{M}_\mathrm{Jup}$) brown dwarf or high mass planet with an unresolved disk and dust envelope.

Debris disks are usually detected through the infrared excess over the photospheric level of their host star. The most favorable stars for disk detection are those with spectral types between A and K, while the statistics for debris disks detected around low-mass M-type stars is very low, either because they are rare or because they are more difficult to detect. Terrestrial planets, on the other hand, may be common around M-type stars. Here, we report on the discovery of an extended (likely) debris disk around the M-dwarf GSC 07396-00759. The star is a wide companion of the close accreting binary V4046 Sgr. The system probably is a member of the $\beta$ Pictoris Moving Group. We resolve the disk in scattered light, exploiting high-contrast, high-resolution imagery with the two near-infrared subsystems of the VLT/SPHERE instrument, operating in the YJ bands and the H2H3 doublet. The disk is clearly detected up to 1.5" ($\sim110$ au) from the star and appears as a ring, with an inclination $i\sim83$ degree, and a peak density position at $\sim 70$ au. The spatial extension of the disk suggests that the dust dynamics is affected by a strong stellar wind, showing similarities with the AU Mic system that has also been resolved with SPHERE. The images show faint asymmetric structures at the widest separation in the northwest side. We also set an upper limit for the presence of giant planets to $2 M_J$. Finally, we note that the 2 resolved disks around M-type stars of 30 such stars observed with SPHERE are viewed close to edge-on, suggesting that a significant population of debris disks around M dwarfs could remain undetected because of an unfavorable orientation.

Debris disks are the intrinsic by-products of the star and planet formation processes. Most likely due to instrumental limitations and their natural faintness, little is known about debris disks around low-mass stars, especially when it comes to spatially resolved observations. We present new VLT/SPHERE IRDIS Dual-Polarization Imaging (DPI) observations in which we detect the dust ring around the M2 spectral type star TWA\,7. Combined with additional Angular Differential Imaging observations we aim at a fine characterization of the debris disk and setting constraints on the presence of low-mass planets. We model the SPHERE DPI observations and constrain the location of the small dust grains, as well as the spectral energy distribution of the debris disk, using the results inferred from the observations, and perform simple N-body simulations. We find that the dust density distribution peaks at 25 au, with a very shallow outer power-law slope, and that the disk has an inclination of 13 degrees with a position angle of 90 degrees East of North. We also report low signal-to-noise detections of an outer belt at a distance of ~52 au from the star, of a spiral arm in the Southern side of the star, and of a possible dusty clump at 3.9 au. These findings seem to persist over timescales of at least a year. Using the intensity images, we do not detect any planets in the close vicinity of the star, but the sensitivity reaches Jovian planet mass upper limits. We find that the SED is best reproduced with an inner disk at 7 au and another belt at 25 au. We report the detections of several unexpected features in the disk around TWA\,7. A yet undetected 100 M$_\oplus$ planet with a semi-major axis at 20-30 au could possibly explain the outer belt as well as the spiral arm. We conclude that stellar winds are unlikely to be responsible for the spiral arm.

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 885$\mu$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from $\sim$0.3-20M$_{\oplus}$. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in $M_{dust}$ to $M_{star}$, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early M-stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas:dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M-dwarfs.

In this work we present ALMA continuum observations at 880 $\mu$m of 30 sub-mm cores previously identified with APEX/LABOCA at 870$\mu$m in the Barnard 30 cloud. The main goal is to characterize the youngest and lowest mass population in the cloud. As a result, we report the detection of five (out of 30) spatially unresolved sources with ALMA, with estimated masses between 0.9 and 67 M$_{\rm Jup}$. From these five sources, only two show gas emission. The analysis of multi-wavelength photometry from these two objects, namely B30-LB14 and B30-LB19, is consistent with one Class II- and one Class I low-mass stellar object, respectively. The gas emission is consistent with a rotating disk in the case of B30-LB14, and with an oblate rotating envelope with infall signatures in the case of LB19. The remaining three ALMA detections do not have infrared counterparts and can be classified as either deeply embedded objects or as starless cores if B30 members. In the former case, two of them (LB08 and LB31) show internal luminosity upper limits consistent with Very Low Luminosity objects, while we do not have enough information for LB10. In the starless core scenario, and taking into account the estimated masses from ALMA and the APEX/LABOCA cores, we estimate final masses for the central objects in the substellar domain, so they could be classified as pre-BD core candidates.

We present observations of the Herbig Ae star HD169142 with VLT/SPHERE instruments InfraRed Dual-band Imager and Spectrograph (IRDIS) ($K1K2$ and $H2H3$ bands) and the Integral Field Spectrograph (IFS) ($Y$, $J$ and $H$ bands). We detect several bright blobs at $\sim$180 mas separation from the star, and a faint arc-like structure in the IFS data. Our reference differential imaging (RDI) data analysis also finds a bright ring at the same separation. We show, using a simulation based on polarized light data, that these blobs are actually part of the ring at 180 mas. These results demonstrate that the earlier detections of blobs in the $H$ and $K_S$ bands at these separations in Biller et al. as potential planet/substellar companions are actually tracing a bright ring with a Keplerian motion. Moreover, we detect in the images an additional bright structure at $\sim$93 mas separation and position angle of 355$^{\circ}$, at a location very close to previous detections. It appears point-like in the $YJ$ and $K$ bands but is more extended in the $H$ band. We also marginally detect an inner ring in the RDI data at $\sim$100 mas. Follow-up observations are necessary to confirm the detection and the nature of this source and structure.

DZ Cha is a weak-lined T Tauri star (WTTS) surrounded by a bright protoplanetary disc with evidence of inner disc clearing. Its narrow $\Ha$ line and infrared spectral energy distribution suggest that DZ Cha may be a photoevaporating disc. We aim to analyse the DZ Cha star + disc system to identify the mechanism driving the evolution of this object. We have analysed three epochs of high resolution optical spectroscopy, photometry from the UV up to the sub-mm regime, infrared spectroscopy, and J-band imaging polarimetry observations of DZ Cha. Combining our analysis with previous studies we find no signatures of accretion in the $\Ha$ line profile in nine epochs covering a time baseline of $\sim20$ years. The optical spectra are dominated by chromospheric emission lines, but they also show emission from the forbidden lines [SII] 4068 and [OI] 6300$\,\AA$ that indicate a disc outflow. The polarized images reveal a dust depleted cavity of $\sim7$ au in radius and two spiral-like features, and we derive a disc dust mass limit of $M_\mathrm{dust}<3\MEarth$ from the sub-mm photometry. No stellar ($M_\star > 80 \MJup$) companions are detected down to $0\farcs07$ ($\sim 8$ au, projected). The negligible accretion rate, small cavity, and forbidden line emission strongly suggests that DZ Cha is currently at the initial stages of disc clearing by photoevaporation. At this point the inner disc has drained and the inner wall of the truncated outer disc is directly exposed to the stellar radiation. We argue that other mechanisms like planet formation or binarity cannot explain the observed properties of DZ Cha. The scarcity of objects like this one is in line with the dispersal timescale ($\lesssim 10^5$ yr) predicted by this theory. DZ Cha is therefore an ideal target to study the initial stages of photoevaporation.

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.

A large number of systems harboring a debris disk show evidence for a double belt architecture. One hypothesis for explaining the gap between the belts is the presence of one or more planets dynamically carving it. This work aims to investigate this scenario in systems harboring two components debris disks. All the targets in the sample were observed with the SPHERE instrument which performs high-contrast direct imaging. Positions of the inner and outer belts were estimated by SED fitting of the infrared excesses or, when available, from resolved images of the disk. Very few planets have been observed so far in debris disks gaps and we intended to test if such non-detections depend on the observational limits of the present instruments. This aim is achieved by deriving theoretical predictions of masses, eccentricities and semi-major axes of planets able to open the observed gaps and comparing such parameters with detection limits obtained with SPHERE. The relation between the gap and the planet is due to the chaotic zone around the orbit of the planet. The radial extent of this zone depends on the mass ratio between the planet and the star, on the semi-major axis and on the eccentricity of the planet and it can be estimated analytically. We apply the formalism to the case of one planet on a circular or eccentric orbit. We then consider multi-planetary systems: 2 and 3 equal-mass planets on circular orbits and 2 equal-mass planets on eccentric orbits in a packed configuration. We then compare each couple of values (M,a), derived from the dynamical analysis of single and multiple planetary models, with the detection limits obtained with SPHERE. Our results show that the apparent lack of planets in gaps between double belts could be explained by the presence of a system of two or more planets possibly of low mass and on an eccentric orbits whose sizes are below the present detection limits.

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

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

We report the results of an extended spectropolarimetric and photometric monitoring of the weak-line T Tauri star TAP 26, carried out within the MaTYSSE programme with the ESPaDOnS spectropolarimeter at the 3.6 m Canada-France-Hawaii Telescope. Applying Zeeman-Doppler Imaging to our observations, concentrating in 2015 November and 2016 January and spanning 72 d in total, 16 d in 2015 November and 13 d in 2016 January, we reconstruct surface brightness and magnetic field maps for both epochs and demonstrate that both distributions exhibit temporal evolution not explained by differential rotation alone. We report the detection of a hot Jupiter (hJ) around TAP 26 using three different methods, two using Zeeman-Doppler Imaging (ZDI) and one Gaussian-Process Regression (GPR), with a false-alarm probability smaller than 6.10^-4. However, as a result of the aliasing related to the observing window, the orbital period cannot be uniquely determined; the orbital period with highest likelihood is 10.79 +/- 0.14 d followed by 8.99 +/- 0.09 d. Assuming the most likely period, and that the planet orbits in the stellar equatorial plane, we obtain that the planet has a minimum mass M.sin(i) of 1.66 +/- 0.31 M_Jup and orbits at 0.0968 +/- 0.0032 au from its host star. This new detection suggests that disc type II migration is efficient at generating newborn hJs, and that hJs may be more frequent around young T Tauri stars than around mature stars (or that the MaTYSSE sample is biased towards hJ-hosting stars).

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.

The large cavities observed in the dust and gas distributions of transition disks may be explained by planet-disk interactions. At ~145 pc, 2MASS J16042165-2130284 (J1604) is a 5-12 Myr old transitional disk with different gap sizes in the mm- and $\mu$m-sized dust distributions (outer edges at ~79 and at ~63 au, respectively). Its $^{12}$CO emission shows a ~30 au cavity. This radial structure suggests that giant planets are sculpting this disk. We aim to constrain the masses and locations of plausible giant planets around J1604. We observed J1604 with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) at the Very Large Telescope (VLT), in IRDIFS\_EXT, pupil-stabilized mode, obtaining YJH- band images with the integral field spectrograph (IFS) and K1K2-band images with the Infra-Red Dual-beam Imager and Spectrograph (IRDIS). The dataset was processed exploiting the angular differential imaging (ADI) technique with high-contrast algorithms. Our observations reach a contrast of $\Delta K, YH$ ~12 mag from 0.15" to 0.80" (~22 to 115 au), but no planet candidate is detected. The disk is directly imaged in scattered light at all bands from Y to K, and it shows a red color. This indicates that the dust particles in the disk surface are mainly $\gtrsim0.3\,\mu$m-sized grains. We confirm the sharp dip/decrement in scattered light in agreement with polarized light observations. Comparing our images with a radiative transfer model we argue that the southern side of the disk is most likely the nearest. This work represents the deepest search yet for companions around J1604. We reach a mass sensitivity of $\gtrsim 2-3M_{Jup}$ from ~22 to ~115 au according to a hot start scenario. We propose that a brown dwarf orbiting inside of ~15 au and additional Jovian planets at larger radii could account for the observed properties of J1604 while explaining our lack of detection.

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

We present H-band (1.6 {\mu}m) scattered light observations of the transitional disk RX J1615.3-3255, located in the ~1 Myr old Lupus association. From a polarized intensity image, taken with the HiCIAO instrument of the Subaru Telescope, we deduce the position angle and the inclination angle of the disk. The disk is found to extend out to 68 $\pm$ 12 AU in scattered light and no clear structure is observed. Our inner working angle of 24 AU does not allow us to detect a central decrease in intensity similar to that seen at 30 AU in the 880 {\mu}m continuum observations. We compare the observations with multiple disk models based on the Spectral Energy Distribution (SED) and submm interferometry and find that an inner rim of the outer disk at 30 AU containing small silicate grains produces a polarized intensity signal which is an order of magnitude larger than observed. We show that a model in which the small dust grains extend smoothly into the cavity found for large grains is closer to the actual H-band observations. A comparison of models with different dust size distributions suggests that the dust in the disk might have undergone significant processing compared to the interstellar medium.

In 2015, we initiated a survey of Scorpius-Centaurus A-F stars that are predicted to host warm-inner and cold-outer belts of debris similar to the case of the system HR~8799. The survey aims to resolve the disks and detect planets responsible for the disk morphology. In this paper, we study the F-type star HIP~67497 and present a first-order modelisation of the disk in order to derive its main properties. We used the near-infrared integral field spectrograph (IFS) and dual-band imager IRDIS of VLT/SPHERE to obtain angular-differential imaging observations of the circumstellar environnement of HIP~67497. We removed the stellar halo with PCA and TLOCI algorithms. We modeled the disk emission with the GRaTeR code. We resolve a ring-like structure that extends up to $\sim$450 mas ($\sim$50 au) from the star in the IRDIS and IFS data. It is best reproduced by models of a non-eccentric ring with an inclination of $80\pm1^{\circ}$, a position angle of $-93\pm1^{\circ}$, and a semi-major axis of $59\pm3$ au. We also detect an additional, but fainter, arc-like structure with a larger extension (0.65 arcsec) South of the ring that we model as a second belt of debris at $\sim$130 au. We detect 10 candidate companions at separations $\geq$1''. We estimate the mass of putative perturbers responsible for the disk morphology and compare it to our detection limits. Additional data are needed to find those perturbers, and to relate our images to large-scale structures seen with HST/STIS.

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.

Context: Gaps, cavities and rings in circumstellar disks are signposts of disk evolution and planet-disk interactions. We follow the recent suggestion that Herbig Ae/Be disks with a flared disk harbour a cavity, and investigate the disk around HD~97048. Aims: We aim to resolve the 34$\pm$ 4 au central cavity predicted by Maaskant et al. (2013) and to investigate the structure of the disk. Methods: We image the disk around HD~97048 using ALMA at 0.85~mm and 2.94~mm, and ATCA (multiple frequencies) observations. Our observations also include the 12CO J=1-0, 12CO J=3-2 and HCO+ J=4-3 emission lines. Results: A central cavity in the disk around HD~97048 is resolved with a 40-46 au radius. Additional radial structure present in the surface brightness profile can be accounted for either by an opacity gap at ~90 au or by an extra emitting ring at ~150 au. The continuum emission tracing the dust in the disk is detected out to 355 au. The 12CO J=3-2 disk is detected 2.4 times farther out. The 12CO emission can be traced down to $\approx$ 10 au scales. Non-Keplerian kinematics are detected inside the cavity via the HCO+ J=4-3 velocity map. The mm spectral index measured from ATCA observations suggests that grain growth has occurred in the HD~97048 disk. Finally, we resolve a highly inclined disk out to 150 au around the nearby 0.5~$M_{\odot}$ binary ISO-ChaI 126. Conclusions: The data presented here reveal a cavity in the disk of HD 97048, and prominent radial structure in the surface brightness. The cavity size varies for different continuum frequencies and gas tracers. The gas inside the cavity follows non-Keplerian kinematics seen in HCO+ emission. The variable cavity size along with the kinematical signature suggests the presence of a substellar companion or massive planet inside the cavity.

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.

Models of planet formation and disc evolution predict a variety of observables in the dust structure of protoplanetary discs. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 and Band-7 observations of the transition disc Sz\,91 showing that the continuum emission at 870$\mu$m, which is dominated by emission from large dust grains, is localized in an optically thin narrow ring. We find that most of the emission ($\sim95\%$) is concentrated in a ring located at 110 au from the central star that is only about 44 au wide. In contrast, the $^{12}\mathrm{CO}$ (2-1) emission peaks closer to the star and is detected up to $\sim488$ au from the star. The concentration of large grains in a ring-like structure while the gas disc extends much further in and further out is in qualitative agreement with predictions of hydrodynamical models of planet-disc interactions including radial drift and gas drag.

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

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

Several studies discussing imaging polarimetry observations of protoplanetary disks use the so-called radial Stokes parameters Q_phi and U_phi to discuss the results. This approach has the advantage of providing a direct measure of the noise in the polarized images under the assumption that the polarization is azimuthal only, i.e., perpendicular to the direction towards the illuminating source. However, a detailed study of the validity of this assumption is currently missing. We aim to test whether departures from azimuthal polarization can naturally be produced by scattering processes in optically thick protoplanetary disks at near infrared wavelengths. We use the radiative transfer code MCFOST to create a generic model of a transition disk using different grain size distributions and dust masses. From these models we generate synthetic polarized images at 2.2\mum. We find that even for moderate inclinations (e.g., i = 40degr), multiple scattering alone can produce significant (up to ~4.5% of the Q_phi image) non-azimuthal polarization reflected in the U_phi images. We also find that different grain populations can naturally produce radial polarization (negative values in the Q_phi images). Our results suggest that caution is recommended when interpreting polarized images by only analyzing the Q_phi and U_phi images. We find that there can be astrophysical signal in the U_phi images and negative values in the Q_phi images, which indicate departures from azimuthal polarization. If significant signal is detected in the U_phi images, we recommend to check the standard Q and U images to look for departures from azimuthal polarization. On the positive side, signal in the U_phi images once all instrumental and data-reduction artifacts have been corrected for means that there is more information to be extracted regarding the dust population and particle density.

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri stars (wTTSs) V819 Tau and V830 Tau within the MaTYSSE programme, involving the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope. At ~3 Myr, both stars dissipated their discs recently and are interesting objects for probing star and planet formation. Profile distortions and Zeeman signatures are detected in the unpolarized and circularly-polarized lines, whose rotational modulation we modelled using tomographic imaging, yielding brightness and magnetic maps for both stars. We find that the large-scale magnetic fields of V819 Tau and V830 Tau are mostly poloidal and can be approximated at large radii by 350-400 G dipoles tilted at ~30 degrees to the rotation axis. They are significantly weaker than the field of GQ Lup, an accreting classical T Tauri star (cTTS) with similar mass and age which can be used to compare the magnetic properties of wTTSs and cTTSs. The reconstructed brightness maps of both stars include cool spots and warm plages. Surface differential rotation is small, typically ~4.4x smaller than on the Sun, in agreement with previous results on wTTSs. Using our Doppler images to model the activity jitter and filter it out from the radial velocity (RV) curves, we obtain RV residuals with dispersions of 0.033 and 0.104 km/s for V819 Tau and V830 Tau respectively. RV residuals suggest that a hot Jupiter may be orbiting V830 Tau, though additional data are needed to confirm this preliminary result. We find no evidence for close-in giant planet around V819 Tau.

Protoplanetary discs are now routinely observed and exoplanets, after the numerous indirect discoveries, are starting to be directly imaged. To better understand the planet formation process, the next step is the detection of forming planets or of signposts of young planets still in their disc, such as gaps. A spectacular example is the ALMA science verification image of HL Tau showing numerous gaps and rings in its disc. To study the observability of planet gaps, we ran 3D hydrodynamical simulations of a gas and dust disc containing a 5 M J gap-opening planet and characterised the spatial distribution of migrating, growing and fragmenting dust grains. We then computed the corresponding synthetic images for ALMA. For a value of the dust fragmentation threshold of 15 m s --1 for the collisional velocity, we identify for the first time a self-induced dust pile up in simulations taking fragmentation into account. This feature, in addition to the easily detected planet gap, causes a second apparent gap that could be mistaken for the signature of a second planet. It is therefore essential to be cautious in the interpretation of gap detections.

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.

We model the dust evolution in protoplanetary disks with full 3D, Smoothed Particle Hydrodynamics (SPH), two-phase (gas+dust) hydrodynamical simulations. The gas+dust dynamics, where aerodynamic drag leads to the vertical settling and radial migration of grains, is consistently treated. In a previous work, we characterized the spatial distribution of non-growing dust grains of different sizes in a disk containing a gap-opening planet and investigated the gap's detectability with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we take into account the effects of grain growth and fragmentation and study their impact on the distribution of solids in the disk. We show that rapid grain growth in the two accumulation zones around planet gaps is strongly affected by fragmentation. We discuss the consequences for ALMA observations.

Context. Z\,CMa is a binary composed of an embedded Herbig Be and an FU Ori class star separated by $\sim100$ au. Observational evidence indicate a complex environment in which each star has a circumstellar disk and drives a jet, and the whole system is embedded in a large dusty envelope. Aims. We aim to probe the circumbinary environment of Z\,CMa in the inner 400 au in scattered light. Methods. We use high contrast imaging polarimetry with VLT/NaCo at $H$ and $K_s$ bands. Results. The central binary is resolved in both bands. The polarized images show three bright and complex structures: a common dust envelope, a sharp extended feature previously reported in direct light, and an intriguing bright clump located $0\farcs3$ south of the binary, which appears spatially connected to the sharp extended feature. Conclusions.We detect orbital motion when compared to previous observations, and report a new outburst driven by the Herbig star. Our observations reveal the complex inner environment of Z\,CMa with unprecedented detail and contrast.

We present ALMA (Cycle 0) band-6 and band-3 observations of the transition disk Sz\,91. The disk inclination and position angle are determined to be $i=49.5\degr\pm3.5\degr$ and $\mathrm{PA}=18.2\degr\pm3.5\degr$ and the dusty and gaseous disk are detected up to $\sim220$ au and $\sim400$ au from the star, respectively. Most importantly, our continuum observations indicate that the cavity size in the mm-sized dust distribution must be $\sim97$ au in radius, the largest cavity observed around a T Tauri star. Our data clearly confirms the presence of \co(2-1) well inside the dust cavity. Based on these observational constrains we developed a disk model that simultaneously accounts for the \co and continuum observations (i.e., gaseous and dusty disk). According to our model, most of the millimeter emission comes from a ring located between 97 and 140 au. We also find that the dust cavity is divided into an innermost region largely depleted of dust particles ranging from the dust sublimation radius up to 85 au, and a second, moderately dust-depleted region, extending from 85 to 97 au. The extremely large size of the dust cavity, the presence of gas and small dust particles within the cavity and the accretion rate of Sz\,91 are consistent with the formation of multiple (giant) planets.

Context : The properties of the inner disks of bright Herbig AeBe stars have been studied with near infrared (NIR) interferometry and high resolution spectroscopy. The continuum and a few molecular gas species have been studied close to the central star; however, sensitivity problems limit direct information about the inner disks of the fainter T Tauri stars. Aims : Our aim is to measure some of the properties of the inner regions of disks surrounding southern T Tauri stars. Methods : We performed a survey with the PIONIER recombiner instrument at H-band of 21 T Tauri stars. The baselines used ranged from 11 m to 129 m, corresponding to a maximum resolution of 3mas (0.45 au at 150 pc). Results : Thirteen disks are resolved well and the visibility curves are fully sampled as a function of baseline in the range 45-130 m for these 13 objects. A simple qualitative examination of visibility profiles allows us to identify a rapid drop-off in the visibilities at short baselines in 8 resolved disks. This is indicative of a significant contribution from an extended contribution of light from the disk. We demonstrate that this component is compatible with scattered light, providing strong support to a prediction made by Pinte et al. (2008). The amplitude of the drop-off and the amount of dust thermal emission changes from source to source suggesting that each disk is different. A by-product of the survey is the identification of a new milli-arcsec separation binary: WW Cha. Spectroscopic and interferometric data of AK Sco have also been fitted with a binary and disk model. Conclusions : Visibility data are reproduced well when thermal emission and scattering form dust are fully considered. The inner radii measured are consistent with the expected dust sublimation radii. Modelling of AK Sco suggests a likely coplanarity between the disk and the binary's orbital plane

Inner cavities and annular gaps in circumstellar disks are possible signposts of giant planet formation. The young star HD 142527 hosts a massive protoplanetary disk with a large cavity that extends up to 140 au from the central star, as seen in continuum images at infrared and millimeter wavelengths. Estimates of the survival of gas inside disk cavities are needed to discriminate between clearing scenarios. We present a spatially and spectrally resolved carbon monoxide isotopologue observations of the gas-rich disk HD 142527, in the J=2-1 line of 12CO, 13CO and C18O, obtained with the Atacama Large Millimeter Array (ALMA). We detect emission coming from inside the dust-depleted cavity in all three isotopologues. Based on our analysis of the gas in the dust cavity, the 12CO emission is optically thick, while 13CO and C18O emission are both optically thin. The total mass of residual gas inside the cavity is about 1.5-2 Jupiter masses. We model the gas with an axisymmetric disk model. Our best fit model shows that the cavity radius is much smaller in CO than it is in millimeter continuum and scattered light observations, with a gas cavity that does not extend beyond 105 au (at 3-sigma). The gap wall at its outer edge is diffuse and smooth in the gas distribution, while in dust continuum it is manifestly sharper. The inclination angle, as estimated from the high velocity channel maps, is 28+/-0.5 degrees, higher than in previous estimates, assuming a fix central star mass of 2.2 Solar masses.

Molecular gas constitutes the dominant mass component of protoplanetary discs. To date, these sources have not been studied comprehensively at the longest far-infrared and shortest submillimetre wavelengths. This paper presents Herschel SPIRE FTS spectroscopic observations toward 18 protoplanetary discs, covering the entire 450-1540 GHz (666-195 $\mu$m) range at R~400-1300. The spectra reveal clear detections of the dust continuum and, in six targets, a significant amount of spectral line emission primarily attributable to $^{12}$CO rotational lines. Other targets exhibit little to no detectable spectral lines. Low signal-to-noise detections also include signatures from $^{13}$CO, [CI] and HCN. For completeness, we present upper limits of non-detected lines in all targets, including low-energy transitions of H2O and CH$^+$ molecules. The ten $^{12}$CO lines that fall within the SPIRE FTS bands trace energy levels of ~50-500 K. Combined with lower and higher energy lines from the literature, we compare the CO rotational line energy distribution with detailed physical-chemical models, for sources where these are available and published. Our 13CO line detections in the disc around Herbig Be star HD 100546 exceed, by factors of ~10-30, the values predicted by a model that matches a wealth of other observational constraints, including the SPIRE $^{12}$CO ladder. To explain the observed $^{12}$CO/$^{13}$CO ratio, it may be necessary to consider the combined effects of optical depth and isotope selective (photo)chemical processes. Considering the full sample of 18 objects, we find that the strongest line emission is observed in discs around Herbig Ae/Be stars, although not all show line emission. In addition, two of the six T Tauri objects exhibit detectable $^{12}$CO lines in the SPIRE range.

With Herschel/PACS 134 low mass members of the Taurus star-forming region spanning the M4-L0 spectral type range and covering the transition from low mass stars to brown dwarfs were observed. Combining the new Herschel results with other programs, a total of 150 of the 154 M4-L0 Taurus members members have observations with Herschel. Among the 150 targets, 70um flux densities were measured for 7 of the 7 ClassI objects, 48 of the 67 ClassII members, and 3 of the 76 ClassIII targets. For the detected ClassII objects, the median 70um flux density level declines with spectral type, however, the distribution of excess relative to central object flux density does not change across the stellar/substellar boundary in the M4-L0 range. Connecting the 70um TBOSS values with the results from K0-M3 ClassII members results in the first comprehensive census of far-IR emission across the full mass spectrum of the stellar and substellar population of a star-forming region, and the median flux density declines with spectral type in a trend analogous to the flux density decline expected for the central objects. SEDs were constructed for all TBOSS targets covering the optical to far-IR range and extending to the submm/mm for a subset of sources. Based on an initial exploration of the impact of different physical parameters; inclination, scale height and flaring have the largest influence on the PACS flux densities. From the 24um to 70um spectral index of the SEDs, 5 new candidate transition disks were identified. The steep 24um to 70um slope for a subset of 8 TBOSS targets may be an indication of truncated disks in these systems.Two examples of mixed pair systems that include secondaries with disks were measured. Finally, comparing the TBOSS results with a Herschel study of Ophiuchus brown dwarfs reveals a lower fraction of disks around the Taurus substellar population.

Investigating the evolution of protoplanetary disks is crucial for our understanding of star and planet formation. Several theoretical and observational studies have been performed in the last decades to advance this knowledge. FT Tauri is a young star in the Taurus star forming region that was included in a number of spectroscopic and photometric surveys. We investigate the properties of the star, the circumstellar disk, and the accretion and ejection processes and propose a consistent gas and dust model also as a reference for future observational studies. We performed a multi-wavelength data analysis to derive the basic stellar and disk properties, as well as mass accretion/outflow rate from TNG-Dolores, WHT-Liris, NOT-Notcam, Keck-Nirspec, and Herschel-Pacs spectra. From the literature, we compiled a complete Spectral Energy Distribution. We then performed detailed disk modeling using the MCFOST and ProDiMo codes. Multi-wavelengths spectroscopic and photometric measurements were compared with the reddened predictions of the codes in order to constrain the disk properties. This object can serve as a benchmark for primordial disks with significant mass accretion rate, high gas content and typical size.

We carried out a 12CO(3-2) survey of 52 southern stars with a wide range of IR excesses (LIR/L*) using the single dish telescopes APEX and ASTE. The main aims were (1) to characterize the evolution of molecular gas in circumstellar disks using LIR/L* values as a proxy of disk dust evolution, and (2) to identify new gas-rich disk systems suitable for detailed study with ALMA. About 60% of the sample (31 systems) have LIR/L* > 0.01 typical of T-Tauri or Herbig AeBe stars, and the rest (21 systems) have LIR/L* < 0.01 typical of debris disks. We detect CO(3-2) emission from 20 systems, and 18 (90%) of these have LIR/L* > 0.01. However, the spectra of only four of the newly detected systems appear free of contamination from background or foreground emission from molecular clouds. These include the early-type stars HD 104237 (A4/5V, 116 pc) and HD 98922 (A2 III, 507 pc, as determined in this work), where our observations reveal the presence of CO-rich circumstellar disks for the first time. Of the other detected sources, many could harbor gaseous circumstellar disks, but our data are inconclusive. For these two newly discovered gas-rich disks, we present radiative transfer models that simultaneously reproduce their spectral energy distributions and the 12CO(3-2) line profiles. For both of these systems, the data are fit well by geometrically flat disks, placing them in the small class of non-flaring disks with significant molecular gas reservoirs.

HD 135344B is an accreting (pre-) transition disk that displays the emission of warm CO extending tens of AU inside its 30 AU dust cavity. We used the dust radiative transfer code MCFOST and the thermochemical code ProDiMo to derive the disk structure from the simultaneous modeling of the spectral energy distribution (SED), VLT/CRIRES CO P(10) 4.75 micron, Herschel/PACS [O I] 63 micron, Spitzer-IRS, and JCMT 12CO J=3-2 spectra, VLTI/PIONIER H-band visibilities, and constraints from (sub-)mm continuum interferometry and near-IR imaging. We found a disk model able to describe the current observations simultaneously. This disk has the following structure. (1) To reproduce the SED, the near-IR interferometry data, and the CO ro-vibrational emission, refractory grains (we suggest carbon) are present inside the silicate sublimation radius (0.08<R<0.2 AU). (2) The dust cavity (R<30 AU) is filled with gas, the surface density of this gas must increase with radius to fit the CO P(10) line profile, a small gap of a few AU in the gas is compatible with current data, and a large gap in the gas is not likely. (4) The gas/dust ratio inside the cavity is > 100 to account for the 870 micron continuum upper limit and the CO P(10) line flux. (5) The gas/dust ratio at 30<R<200 AU is < 10 to simultaneously describe the [O I] 63 micron line flux and the CO P(10) line profile. (6) In the outer disk, most of the mass should be located in the midplane, and a significant fraction of the dust is in large grains. Conclusions: Simultaneous modeling of the gas and dust is required to break the model degeneracies and constrain the disk structure. An increasing gas surface density with radius in the inner dust cavity echoes the effect of a migrating Jovian planet. The low gas mass (a few MJupiter) in the HD 135344B's disk suggests that it is an evolved disk that has already lost a large portion of its mass.

In view of both the size of its gap and the previously reported asymmetries and near-infrared spiral arms, the transition disk of the Herbig Fe star HD 142527 constitutes a remarkable case study. This paper focuses on the morphology of the outer disk through ALMA observations of $^{12}$CO J=2-1, $^{12}$CO J=3-2 and $^{13}$CO J=2-1. Both $^{12}$CO J=2-1 and $^{12}$CO J=3-2 show spiral features of different sizes. The innermost spiral arm (S1) is a radio counterpart of the first near-infrared spiral observed by Fukagawa et al. (2006), but it is shifted radially outward. However, the most conspicuous CO spiral arm (S2) lies at the outskirts of the disk and had not been detected before. It corresponds to a cold density structure, with both brightness and excitation temperatures of order 13$\pm$2 K and conspicuous in the $^{12}$CO J=2-1 peak-intensity map, but faint in $^{12}$CO J=3-2. There is also a faint counterarm (S3), point-symmetrical of S2 with respect to the star. These three spirals are modelled separately with two different formulae that approximate the loci of density maxima in acoustic waves due to embedded planets. S1 could be fit relatively well with these formulae, compared to S2 and S3. Alternative scenarios such as gravitational instability or external tidal interaction are discussed. The impact of channelization on spectrally and spatially resolved peak intensity maps is also briefly addressed.

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

HD 166191 has been identified by several studies as hosting a rare and extremely bright warm debris disc with an additional outer cool disc component. However, an alternative interpretation is that the star hosts a disc that is currently in transition between a full gas disc and a largely gas-free debris disc. With the help of new optical to mid-IR spectra and Herschel imaging, we argue that the latter interpretation is supported in several ways: i) we show that HD 166191 is co-moving with the ~4 Myr-old Herbig Ae star HD 163296, suggesting that the two have the same age, ii) the disc spectrum of HD 166191 is well matched by a standard radiative transfer model of a gaseous protoplanetary disc with an inner hole, and iii) the HD 166191 mid-IR silicate feature is more consistent with similarly primordial objects. We note some potential issues with the debris disc interpretation that should be considered for such extreme objects, whose lifetime at the current brightness is mush shorter than the stellar age, or in the case of the outer component requires a mass comparable to the solid component of the Solar nebula. These aspects individually and collectively argue that HD 166191 is a 4-5 Myr old star that hosts a gaseous transition disc. Though it does not argue in favour of either scenario, we find strong evidence for 3-5 um disc variability. We place HD 166191 in context with discs at different evolutionary stages, showing that it is a potentially important object for understanding the protoplanetary to debris disc transition.

The Herschel Space Observatory was used to observe ~ 120 pre-main-sequence stars in Taurus as part of the GASPS Open Time Key project. PACS was used to measure the continuum as well as several gas tracers such as [OI] 63 \mu m, [OI] 145 \mu m, [CII] 158 \mu m, OH, H2O and CO. The strongest line seen is [OI] at 63 \mu m. We find a clear correlation between the strength of the [OI] 63 \mu m line and the 63 \mu m continuum for disk sources. In outflow sources, the line emission can be up to 20 times stronger than in disk sources, suggesting that the line emission is dominated by the outflow. The tight correlation seen for disk sources suggests that the emission arises from the inner disk ($<$ 50 AU) and lower surface layers of the disk where the gas and dust are coupled. The [OI] 63 \mu m is fainter in transitional stars than in normal Class II disks. Simple SED models indicate that the dust responsible for the continuum emission is colder in these disks, leading to weaker line emission. [CII] 158 \mu m emission is only detected in strong outflow sources. The observed line ratios of [OI] 63 \mu m to [OI] 145 \mu m are in the regime where we are insensitive to the gas-to-dust ratio, neither can we discriminate between shock or PDR emission. We detect no Class III object in [OI] 63 \mu m and only three in continuum, at least one of which is a candidate debris disk.

Among the numerous known extrasolar planets, only a handful have been imaged directly so far, at large orbital radii and in rather evolved systems. The Atacama Large Millimeter/submillimeter Array (ALMA) will have the capacity to observe these wide planetary systems at a younger age, thus bringing a better understanding of the planet formation process. Here we explore the ability of ALMA to detect the gaps carved by planets on wide orbits.

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.

Gaseous giant planet formation is thought to occur in the first few million years following stellar birth. Models predict that giant planet formation carves a deep gap in the dust component (shallower in the gas). Infrared observations of the disk around the young star HD142527, at ~140pc, found an inner disk ~10AU in radius, surrounded by a particularly large gap, with a disrupted outer disk beyond 140AU, indicative of a perturbing planetary-mass body at ~90 AU. From radio observations, the bulk mass is molecular and lies in the outer disk, whose continuum emission has a horseshoe morphology. The vigorous stellar accretion rate would deplete the inner disk in less than a year, so in order to sustain the observed accretion, matter must flow from the outer-disk into the cavity and cross the gap. In dynamical models, the putative protoplanets channel outer-disk material into gap-crossing bridges that feed stellar accretion through the inner disk. Here we report observations with the Atacama Large Millimetre Array (ALMA) that reveal diffuse CO gas inside the gap, with denser HCO+ gas along gap-crossing filaments, and that confirm the horseshoe morphology of the outer disk. The estimated flow rate of the gas is in the range 7E-9 to 2E-7 Msun/yr, which is sufficient to maintain accretion onto the star at the present rate.

Circumstellar disc evolution is paramount for the understanding of planet formation. The GASPS program aims at determining the circumstellar gas and solid mass around ~250 pre-main-sequence Herbig Ae and TTauri stars. We aim to understand the origin and nature of the circumstellar matter orbiting 51 Oph, a young (<1 Myr) luminous B9.5 star. We obtained continuum and line observations with the PACS instrument on board the Herschel Space Observatory and continuum data at 1.2 mm with the IRAM 30m telescope. The SED and line fluxes were modelled using the physico-chemo radiative transfer code ProDiMo. We detected a strong emission by OI at 63 microns using the Herschel Space Observatory. The [OI] emission at 145 microns, the [CII] emission at 158 microns, the high-J CO emissions, and the warm water emissions were not detected. Continuum emission was detected at 1.2 mm. The continuum from the near- to the far-infrared and the [OI] emission are well explained by the emission from a compact hydrostatic disc model with a gas mass of 5E-6 MSun, 100 times that of the solid mass. However, this model fails to match the continuum millimeter flux, which hints at a cold outer disc with a mass in solids of 1E-6 MSun or free-free emission from a photoevaporative disc wind. This outer disc can either be devoid of gas and/or is to cold to emit in the [OI] line. A very flat extended disc model (Rout=400 AU) with a fixed vertical structure and dust settling matches all photometric points and most of the [O I] flux. The observations can be explained by an extended flat disc where dust grains have settled. However, a flat gas disc cannot be reproduced by hydrostatic disc models. The low mass of the 51 Oph inner disc in gas and dust may be explained either by the fast dissipation of an initial massive disc or by a very small initial disc mass.

(Abridged) Circumstellar disks are believed to be the birthplace of planets and are expected to dissipate on a timescale of a few Myr. The processes responsible for the removal of the dust and gas will strongly modify the radial distribution of the dust and consequently the SED. In particular, a young planet will open a gap, resulting in an inner disk dominating the near-IR emission and an outer disk emitting mostly in the far-IR. We analyze a full set of data (including VLTI/Pionier, VLTI/Midi, and VLT/NaCo/Sam) to constrain the structure of the transition disk around TCha. We used the Mcfost radiative transfer code to simultaneously model the SED and the interferometric observations. We find that the dust responsible for the emission in excess in the near-IR must have a narrow temperature distribution with a maximum close to the silicate sublimation temperature. This translates into a narrow inner dusty disk (0.07-0.11 AU). We find that the outer disk starts at about 12 AU and is partially resolved by the Pionier, Sam, and Midi instruments. We show that the Sam closure phases, interpreted as the signature of a candidate companion, may actually trace the asymmetry generated by forward scattering by dust grains in the upper layers of the outer disk. These observations help constrain the inclination and position angle of the outer disk. The presence of matter inside the gap is difficult to assess with present-day observations. Our model suggests the outer disk contaminates the interferometric signature of any potential companion that could be responsible for the gap opening, and such a companion still has to be unambiguously detected. We stress the difficulty to observe point sources in bright massive disks, and the consequent need to account for disk asymmetries (e.g. anisotropic scattering) in model-dependent search for companions.

[Abridged] Aims: We provide predictions for ALMA observations of planet gaps that account for the specific spatial distribution of dust that results from consistent gas+dust dynamics. Methods: In a previous work, we ran full 3D, two-fluid Smoothed Particle Hydrodynamics (SPH) simulations of a planet embedded in a gas+dust T Tauri disk for different planet masses and grain sizes. In this work, the resulting dust distributions are passed to the Monte Carlo radiative transfer code MCFOST to construct synthetic images in the ALMA wavebands. We then use the ALMA simulator to produce images that include thermal and phase noise for a range of angular resolutions, wavelengths, and integration times, as well as for different inclinations, declinations and distances. We also produce images which assume that gas and dust are well mixed with a gas-to-dust ratio of 100 to compare with previous ALMA predictions, all made under this hypothesis. Results: Our findings clearly demonstrate the importance of correctly incorporating the dust dynamics. We show that the gap carved by a 1 M_J planet orbiting at 40 AU is visible with a much higher contrast than the well-mixed assumption would predict. In the case of a 5 M_J planet, we clearly see a deficit in dust emission in the inner disk, and point out the risk of interpreting the resulting image as that of a transition disk with an inner hole if observed in unfavorable conditions. Planet signatures are fainter in more distant disks but declination or inclination to the line-of-sight have little effect on ALMA's ability to resolve the gaps. Conclusions: ALMA has the potential to see signposts of planets in disks of nearby star-forming regions. We present optimized observing parameters to detect them in the case of 1 and 5 M_J planets on 40 AU orbits.

The vestiges of planet formation have been observed in debris disks harboring young and massive gaseous giants. The process of giant planet formation is terminated by the dissipation of gas in the protoplanetary disk. The gas-rich disk around HD142527 features a small inner disk, a large gap from \sim10 to \sim140AU, and a massive outer disk extending out to \sim300AU. The gap could have been carved-out by a giant planet. We have imaged the outer regions of this gap using the adaptive-optics camera NICI on Gemini South. Our images reveal that the disk is dynamically perturbed. The outer boundary of the roughly elliptical gap appears to be composed of several segments of spiral arms. The stellar position is offset by 0.17+-0.02" from the centroid of the cavity, consistent with earlier imaging at coarser resolutions. These transient morphological features are expected in the context of disk evolution in the presence of a perturbing body located inside the cavity. We perform hydro-dynamical simulations of the dynamical clearing of a gap in a disk. A 10Mjup body in a circular orbit at r = 90AU, perturbs the whole disks, even after thousands of orbits. By then the model disk has an eccentric and irregular cavity, flanked by tightly wound spiral arms, but it is still evolving far from steady state. A particular transient configuration that is a qualitative match to HD142527 is seen at 1.7Myr.

We present Herschel PACS photometry of seventeen B- to M-type stars in the 30 Myr-old Tucana-Horologium Association. This work is part of the Herschel Open Time Key Programme "Gas in Protoplanetary Systems" (GASPS). Six of the seventeen targets were found to have infrared excesses significantly greater than the expected stellar IR fluxes, including a previously unknown disk around HD30051. These six debris disks were fitted with single-temperature blackbody models to estimate the temperatures and abundances of the dust in the systems. For the five stars that show excess emission in the Herschel PACS photometry and also have Spitzer IRS spectra, we fit the data with models of optically thin debris disks with realistic grain properties in order to better estimate the disk parameters. The model is determined by a set of six parameters: surface density index, grain size distribution index, minimum and maximum grain sizes, and the inner and outer radii of the disk. The best fitting parameters give us constraints on the geometry of the dust in these systems, as well as lower limits to the total dust masses. The HD105 disk was further constrained by fitting marginally resolved PACS 70 micron imaging.

Line spectra of 68 Taurus T Tauri stars were obtained with the Herschel-PACS (Photodetector Array Camera & Spectrometer) instrument as part of the GASPS (Gas Evolution in Protoplanetary Systems) survey of protoplanetary discs. A careful examination of the line scans centred on the [OI] 63.18 microns fine-structure line unveiled a line at 63.32 micron in some of these spectra. We identify this line with a transition of ortho-water. It is detected confidently (i.e., >3 sigma) in eight sources, i.e., 24% of the sub-sample with gas-rich discs. Several statistical tests were used to search for correlations with other disc and stellar parameters such as line fluxes of [OI] 6300 Armstrong and 63.18 microns; X-ray luminosity and continuum levels at 63 microns and 850 microns. Correlations are found between the water line fluxes and the [OI] 63.18 microns line luminosity, the dust continuum, and possibly with the stellar X-ray luminosity. This is the first time that this line of warm water vapour has been detected in protoplanetary discs. We discuss its origins, in particular whether it comes from the inner disc and/or disc surface or from shocks in outflows and jets. Our analysis favours a disc origin, with the observed water vapour line produced within 2-3AU from the central stars, where the gas temperature is of the order of 500-600 K.

HD 181327 is a young F5/F6V star belonging to the Beta Pictoris moving group (12 Myr). It harbors an optically thin belt of circumstellar material at 90 AU. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. We obtained far-IR observations with the Herschel/PACS instrument, and 3.2 mm observations with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS images that break the degeneracy between the disk geometry and the dust properties. We use the radiative transfer code GRaTer to compute a large grid of models, and we identify the grain models that best reproduce the Spectral Energy Distribution through a Bayesian analysis. We attempt to detect the [OI] and [CII] lines with PACS spectroscopy, providing observables to our photochemical code ProDiMo. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected in the far-IR with PACS and the disk is resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of 0.05 M\oplus (up to 1 mm). We discuss evidences that the grains are fluffy aggregates. The upper limits on the gas lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than 17 M\oplus. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry.

We present deep 1.2 millimeter photometry of 37 stars in the young (5 Myr) Upper Scorpius OB association, sensitive to ~4 x 10^-3 Mjup of cool millimeter dust. Disks around four low- and solar-mass stars are detected, as well as one debris disk around an intermediate mass star, with dust masses ranging from 3.6 x 10^-3 -- 1.0 x 10^-1 Mjup. The source with the most massive disk exhibits a transition-disk spectral energy distribution. Combining our results with previous studies, we find the millimeter-detection fraction of Class II sources has significantly decreased from younger ages, and comparison with near-infrared and Halpha measurements indicates the present disks have undergone significant evolution in composition or structure at all radii. The disks of Upper Scorpius represent the tail-end of the depletion of primordial disks; while a few near-solar mass stars may still sustain giant planet formation, this process has finished around higher mass stars

Despite its importance in the thermal-balance of the gas and in the determination of primeval planetary atmospheres, the chemistry in protoplanetary discs remains poorly constrained with only a handful of detected species. We observed the emission from disc around the Herbig Be star HD 100546 with the PACS instrument in the spectroscopic mode on board the Herschel Space Telescope as part of the Gas in Protoplanetary Systems (GASPS) programme and used archival data from the DIGIT programme to search for the rotational emission of CH+. We detected in both datasets an emission line centred at 72.16 micron that most likely corresponds to the J=5-4 rotational emission of CH+. The J=3-2 and 6-5 transitions are also detected albeit with lower confidence. Other CH+ rotational lines in the PACS observations are blended with water lines. A rotational diagram analysis shows that the CH+ gas is warm at 323 (+2320/-151) K with a mass of 3e-14-5e-12 M_Sun. We modelled the CH+ chemistry with the chemo-physical code ProDiMo using a disc density structure and grain parameters that match continuum observations and near- and mid-infrared interferometric data. The model suggests that CH+ is most abundant at the location of the disc rim at 10-13 AU from the star where the gas is warm, consistent with previous observations of hot CO gas emission.

Studying the physical conditions in circumstellar disks is a crucial step toward understanding planet formation. Of particular interest is the case of HD 100546, a Herbig Be star that presents a gap within the first 13 AU of its protoplanetary disk, that may originate in the dynamical interactions of a forming planet. We gathered a large amount of new interferometric data using the AMBER/VLTI instrument in the H- and K-bands to spatially resolve the warm inner disk and constrain its structure. Then, combining these measurements with photometric observations, we analyze the circumstellar environment of HD 100546 in the light of a passive disk model based on 3D Monte-Carlo radiative transfer. Finally, we use hydrodynamical simulations of gap formation by planets to predict the radial surface density profile of the disk and test the hypothesis of ongoing planet formation. The SED and the NIR interferometric data are adequately reproduced by our model. We show that the H- and K-band emissions are coming mostly from the inner edge of the internal dust disk, located near 0.24 AU from the star, i.e., at the dust sublimation radius in our model. We directly measure an inclination of $33^{\circ} \pm 11^{\circ}$ and a position angle of $140^{\circ} \pm 16^{\circ}$ for the inner disk. This is similar to the values found for the outer disk ($i \simeq 42^{\circ}$, $PA \simeq 145^{\circ}$), suggesting that both disks may be coplanar. We finally show that 1 to 8 Jupiter mass planets located at $\sim 8$ AU from the star would have enough time to create the gap and the required surface density jump of three orders of magnitude between the inner and outer disk. However, no information on the amount of matter left in the gap is available, which precludes us from setting precise limits on the planet mass, for now.

The transition between massive Class II circumstellar disks and Class III debris disks, with dust residuals, has not yet been clearly understood. Disks are expected to dissipate with time, and dust clearing in the inner regions can be the consequence of several mechanisms. Planetary formation is one of them that will possibly open a gap inside the disk. According to recent models based on photometric observations, T Cha is expected to present a large gap within its disk, meaning that an inner dusty disk is supposed to have survived close to the star. We investigate this scenario with new near-infrared interferometric observations. We observed T Cha in the H and K bands using the AMBER instrument at VLTI and used the MCFOST radiative transfer code to model the SED of T Cha and the interferometric observations simultaneously and to test the scenario of an inner dusty structure. We also used a toy model of a binary to check that a companion close to the star can reproduce our observations. The scenario of a close (few mas) companion cannot satisfactorily reproduce the visibilities and SED, while a disk model with a large gap and an inner ring producing the bulk of the emission (in H and K-bands) close to 0.1 AU is able to account for all the observations. With this study, the presence of an optically thick inner dusty disk close to the star and dominating the H and K- bands emission is confirmed. According to our model, the large gap extends up to ~ 7.5 AU. This points toward a companion (located at several AU) gap-opening scenario to explain the morphology of T Cha.

(Abridged) Dust particles evolve in size and lattice structure in protoplanetary disks, due to coagulation, fragmentation and crystallization, and are radially and vertically mixed in disks. This paper aims at determining the mineralogical composition and size distribution of the dust grains in disks around 58 T Tauri stars observed with Spitzer/IRS. We present a spectral decomposition model that reproduces the IRS spectra over the full spectral range. The model assumes two dust populations: a warm component responsible for the 10\mu m emission arising from the disk inner regions and a colder component responsible for the 20-30\mu m emission, arising from more distant regions. We show evidence for a significant size distribution flattening compared to the typical MRN distribution, providing an explanation for the usual boxy 10\mu m feature profile generally observed. We reexamine the crystallinity paradox, observationally identified by Olofsson et al. (2009), and we find a simultaneous enrichment of the crystallinity in both the warm and cold regions, while grain sizes in both components are uncorrelated. Our modeling results do not show evidence for any correlations between the crystallinity and either the star spectral type, or the X-ray luminosity (for a subset of the sample). The size distribution flattening may suggests that grain coagulation is a slightly more effective process than fragmentation in disk atmospheres, and that this imbalance may last over most of the T Tauri phase. This result may also point toward small grain depletion via strong stellar winds or radiation pressure in the upper layers of disk. The non negligible cold crystallinity fractions suggests efficient radial mixing processes in order to distribute crystalline grains at large distances from the central object, along with possible nebular shocks in outer regions of disks that can thermally anneal amorphous grains.

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

Aims. The brown dwarf (BD) formation process has not yet been completely understood. To shed more light on the differences and similarities between star and BD formation processes, we study and compare the disk fraction among both kinds of objects over a large angular region in the Taurus cloud. In addition, we examine the spatial distribution of stars and BD relative to the underlying molecular gas Methods. In this paper, we present new and updated photometry data from the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope on 43 BDs in the Taurus cloud, and recalculate of the BD disk fraction in this region. We also useed recently available CO mm data to study the spatial distribution of stars and BDs relative to the cloud's molecular gas. Results. We find that the disk fraction among BDs in the Taurus cloud is 41 \pm 12%, a value statistically consistent with the one among TTS (58 \pm 9%). We find that BDs in transition from a state where they have a disk to a diskless state are rare, and we study one isolated example of a transitional disk with an inner radius of \approx 0.1 AU (CFHT BD Tau 12, found via its relatively small mid-IR excess compared to most members of Taurus that have disks. We find that BDs are statistically found in regions of similar molecular gas surface density to those associated with stars. Furthermore, we find that the gas column density distribution is almost identical for stellar and substellar objects with and without disks.

We present high-resolution spectroscopic mid-infrared observations of the circumstellar disk around the Herbig Ae star HD97048 obtained with the VLT Imager and Spectrometer for the mid-InfraRed (VISIR). We conducted observations of mid-infrared pure rotational lines of molecular hydrogen (H2) as a tracer of warm gas in the disk surface layers. In a previous paper, we reported the detection of the S(1) pure rotational line of H2 at 17.035 microns and argued it is arising from the inner regions of the disk around the star. We used VISIR on the VLT for a more comprehensive study based on complementary observations of the other mid-infrared molecular transitions, namely S(2) and S(4) at 12.278 microns and 8.025 microns respectively, to investigate the physical properties of the molecular gas in the circumstellar disk around HD97048. We do not detect neither the S(2) line nor the S(4) H2 line from the disk of HD97048, but we derive upper limits on the integrated line fluxes which allows us to estimate an upper limit on the gas excitation temperature, T_ex < 570 K. This limit on the temperature is consistent with the assumptions previously used in the analysis of the S(1) line, and allows us to set stronger contraints on the mass of warm gas in the inner regions of the disk. Indeed, we estimate the mass of warm gas to be lower than 0.1 M_Jup. We also discuss the probable physical mechanisms which could be responsible of the excitation of H2 in the disk of HD97048.

We present a study of the circumstellar environment of IRAS 04158+2805 based on multi-wavelength observations and models. Images in the optical and near-infrared, a polarisation map in the optical, and mid-infrared spectra were obtained with VLT-FORS1, CFHT-IR, and Spitzer-IRS. Additionally we used an X-ray spectrum observed with Chandra. We interpret the observations in terms of a central star surrounded by an axisymmetric circumstellar disc, but without an envelope, to test the validity of this simple geometry. We estimate the structural properties of the disc and its gas and dust content. We modelled the dust disc with a 3D continuum radiative transfer code, MCFOST, based on a Monte-Carlo method that provides synthetic scattered light images and polarisation maps, as well as spectral energy distributions. We find that the disc images and spectral energy distribution narrowly constrain many of the disc model parameters, such as a total dust mass of 1.0-1.75x10^-4 sollar masses and an inclination of 62-63 degrees. The maximum grain size required to fit all available data is of the order of 1.6-2.8 microns although the upper end of this range is loosely constrained. The observed optical polarisation map is reproduced well by the same disc model, suggesting that the geometry we find is adequate and the optical properties are representative of the visible dust content. We compare the inferred dust column density to the gas column density derived from the X-ray spectrum and find a gas-to-dust ratio along the line of sight that is consistent with the ISM value. To our knowledge, this measurement is the first to directly compare dust and gas column densities in a protoplanetary disc.

We review recent progress in high-resolution imaging of scattered light from disks around young stellar objects. Many new disks have been discovered or imaged in scattered light, and improved instrumentation and observing techniques have led to better disk images at optical, near-infrared, and thermal-infrared wavelengths. Multi-wavelength datasets are particularly valuable, as dust particle properties have wavelength dependencies. Modeling the changes in scattered-light images with wavelength gives direct information on the dust properties. This has now been done for several different disks. The results indicate that modest grain growth has taken place in some of these systems. Scattered-light images also provide useful constraints on the disk structure, especially when combined with long-wavelength SEDs. There are tentative suggestions in some disks that the dust may have begun to settle. The next few years should see this work extended to many more disks; this will clarify our understanding of the evolution of protoplanetary dust and disks.

We report adaptive optics (AO) observations of the young and nearby association MBM 12 obtained with the Canada-France-Hawaii Telescope. Our main observational result is the discovery of six new binary systems, LkHa 264, E 0255+2018, RX J0255.4+2005, S18, MBM 12-10, RX J0255.3+1915, and the confirmation of HD 17332, already known as a binary. We also detected a possible quadruple system. It is composed of the close binary LkHa 263 AB (separation of 0.41 ''), of LkH\alpha 262 located 15.25 '' from LkHa 263 A, and of LkHa 263 C, located 4.1 '' from LkH\alpha 263 A. A preliminary study of the binary fraction suggests a binary excess in the MBM 12 association as compared to the field and IC 348. Because of the high binarity rate, previous estimations of spectral types and measurements of IR excesses for several candidate members of MBM 12 have to be revised. LkH\alpha 263 C is a nebulous object that we interpret as a disk oriented almost perfectly edge-on and seen in scattered light. This object has already been reported by Jayawardhana et al. (2002). Scattered light models allow us to estimate some of the structural parameters (i.e. inclination, diameter and to a lesser extent dust mass) of the circumstellar disk. We find an inclination of 89^o and a outer radius for the disk, 165 AU if the distance to MBM 12 is 275 pc. With the present data set, we do not attempt to re-assess the distance to MBM 12. We estimate however that the distance to the candidate member RX J0255.3+1915 is d > 175 pc.