The recent ALMA observations of the disk surrounding HL Tau reveal a very complex dust spatial distribution. We present a radiative transfer model accounting for the observed gaps and bright rings as well as radial changes of the emissivity index. We find that the dust density is depleted by at least a factor of 10 in the main gaps compared to the surrounding rings. Ring masses range from 10-100 M ⊕ in dust, and we find that each of the deepest gaps is consistent with the removal of up to 40 M ⊕ of dust. If this material has accumulated into rocky bodies, these would be close to the point of runaway gas accretion. Our model indicates that the outermost ring is depleted in millimeter grains compared to the central rings. This suggests faster grain growth in the central regions and/or radial migration of the larger grains. The morphology of the gaps observed by ALMA -well separated and showing a high degree of contrast with the bright rings over all azimuths -indicates that the millimeter dust disk is geometrically thin (scale height ≈ 1 AU at 100 AU) and that a large amount of settling of large grains has already occurred. Assuming a standard dust settling model, we find that the observations are consistent with a turbulent viscosity coefficient of a few 10 −4 . We estimate the gas/dust ratio in this thin layer to be of the order of 5 if the initial ratio is 100. The HCO + and CO emission is consistent with gas in Keplerian motion around a 1.7 M star at radii from ≤ 10 − 120 AU.

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

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. The first paper of this series focuses on the assumptions about the shape of the disk, the dust opacities, dust settling, and polycyclic aromatic hydrocarbons (PAHs). In particular, we propose new standard dust opacities for disk models, we present a simplified treatment of PAHs in radiative equilibrium which is sufficient to reproduce the PAH emission features, and we suggest using a simple yet physically justified treatment of dust settling. We roughly adjust parameters to obtain a model that predicts continuum and line observations that resemble typical multi-wavelength continuum and line observations of Class II T Tauri stars. 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 mainstream continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63 μm, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties, i.e. large grains, often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger near-to far-IR 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 efficiently shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities in comparison to evolved dust. 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 H 2 , as additional constraints to determine a number of key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

Discs of gas and dust surrounding young stars are the birthplace of planets. However, the direct detection of protoplanets forming within discs has proved elusive to date. We present the detection of a large, localized deviation from Keplerian velocity in the protoplanetary disc surrounding the young star HD 163296. The observed velocity pattern is consistent with the dynamical effect of a two-Jupiter-mass planet orbiting at a radius ≈ 260 au from the star.

Aims. We present a new continuum 3D radiative transfer code, MCFOST, based on a Monte-Carlo method. MCFOST can be used to calculate (i) monochromatic images in scattered light and/or thermal emission; (ii) polarisation maps; (iii) interferometric visibilities; (iv) spectral energy distributions; and (v) dust temperature distributions of protoplanetary disks. Methods. Several improvements to the standard Monte Carlo method are implemented in MCFOST to increase efficiency and reduce convergence time, including wavelength distribution adjustments, mean intensity calculations, and an adaptive sampling of the radiation field. The reliability and efficiency of the code are tested against a previously-defined benchmark, using a 2D disk configuration. No significant difference (no more than 10% and usually much less) is found between the temperatures and SEDs calculated by MCFOST and by other codes included in the benchmark. Results. A study of the lowest disk mass detectable by Spitzer, around young stars, is presented and the colours of "representative" parametric disks compared to recent IRAC and MIPS Spitzer colours of solar-like young stars located in nearby star-forming regions.

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

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

Context. Accretion and ejection are complex and related processes that vary on various timescales in young stars. Aims. We intend to investigate the accretion and outflow dynamics and their interaction from observations of the classical T Tauri star AA Tau. Methods. From a long time series of high resolution (R = 115 000) HARPS spectra and simultaneous broad-band photometry, we report new evidence for magnetospheric accretion as well as ejection processes in the nearly edge-on classical T Tauri star AA Tau. Results. AA Tau's light curve is modulated with a period of 8.22 d. The recurrent luminosity dips are due to the periodic occultation of the central star by the magnetically-warped inner disk edge located at about 9 R . Balmer line profiles exhibit a clear rotational modulation of high-velocity redshifted absorption components with a period of 8.22 days as well, with a maximum strength when the main accretion funnel flow passes through the line of sight. At the same time, the luminosity of the system decreases by about 1 mag, indicative of circumstellar absorption of the stellar photosphere by the magnetically-warped, corotating inner disk edge. The photospheric and He I radial velocities also exhibit periodic variations, and the veiling is modulated by the appearance of the accretion shock at the bottom of the accretion funnel. Diagnostics of hot winds and their temporal behaviour are also presented. Conclusions. The peculiar geometry of the young AA Tau system (nearly edge-on) allows us to uniquely probe the acretion-ejection region close to the star. We find that most spectral and photometric diagnostics vary as expected from models of magneticallychannelled accretion in young stars, with a large scale magnetosphere tilted by 20• onto the star's spin axis. We also find evidence for time variability of the magnetospheric accretion flow on a timescale of a few rotational periods.