P. Pinilla

Asymmetric features in the protoplanetary disk MWC 758

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 SPHERE/VLT at 1.04 m to resolve scattered light at a smaller inner working angle (0.093 00 ) and a higher angular resolution (0.027 00 ) than previously achieved. Results. We observe polarized scattered light within 0.53 00 (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.

ALMA HINTS AT THE PRESENCE OF TWO COMPANIONS IN THE DISK AROUND HD 100546

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

Gas and dust structures in protoplanetary disks hosting multiple planets

Context. Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of particle trapping and grain growth in pressure bumps originating from the disk interaction with a planetary companion. The multiple features observed in some transition disks, such as SR 21, suggest the presence of more than one planet. Aims. We aim to study the gas and dust distributions of a disk hosting two massive planets as a function of different disk and dust parameters. Observational signatures, such as spectral energy distributions, sub-millimetre, and polarised images, are simulated for various parameters. Methods. Two dimensional hydrodynamical and one dimensional dust evolution numerical simulations are performed for a disk interacting with two massive planets. Adopting the previously determined dust distribution, and assuming an axisymmetric disk model, radiative transfer simulations are used to produce spectral energy distributions and synthetic images in polarised intensity at 1.6 μm and sub-millimetre wavelengths (850 μm). We analyse possible scenarios that can lead to gas azimuthal asymmetries. Results. We confirm that planets can lead to particle trapping, although for a disk with high viscosity (αturb = 10-2), the planet should be more massive than 5 MJup and dust fragmentation should occur with low efficiency (vf ~ 30 m s-1). This will lead to a ring-like feature as observed in transition disks in the millimetre. When trapping occurs, we find that a smooth distribution of micron-sized grains throughout the disk, sometimes observed in scattered light, can only happen if the combination of planet mass and turbulence is such that small grains are not fully filtered out. A high disk viscosity (αturb = 10-2) ensures a replenishment of the cavity in micron-sized dust, while for lower viscosity (αturb = 10-3), the planet mass is constrained to be less than 5 MJup. In these cases, the gas distribution is likely to show low-amplitude azimuthal asymmetries caused by disk eccentricity rather than by long-lived vortices.

Scattered light images of spiral arms in marginally gravitationally unstable discs with an embedded planet

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

A COMPACT CONCENTRATION OF LARGE GRAINS IN THE HD 142527 PROTOPLANETARY DUST TRAP

10^{-2}

Direct detection of scattered light gaps in the transitional disk around HD 97048 with VLT/SPHERE

) can trigger gravitational instability in the outer disc. The spirals created by this instability and the density waves launched by the planet can overlap resulting in large-scale, more open spiral arms in the outer disc. The contrast of these spirals is well above the detection limit of current telescopes.

Direct imaging of the water snow line at the time of planet formation using two Alma continuum bands

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply conned and furthermore the trapped grains should follow a distribution that is markedly dierent from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD 142527, from 34 GHz to 700 GHz, that reveal a compact concentration of cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of mm-sized particles. The emission peaks at wavelengths shorter than 1 mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is ecient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reect the gas distribution. Subject headings: Protoplanetary disks | Planet-disk interactions | Stars: individual: (HD 142527)

Dust trapping by spiral arms in gravitationally unstable protostellar discs

Aims. We studied the well-known circumstellar disk around the Herbig Ae/Be star HD97048 with high angular resolution to reveal undetected structures in the disk which may be indicative of disk evolutionary processes such as planet formation. Methods. 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. Results. 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 similar to 270 au. Using the surface height profile we measured the scattering phase function of the disk and found that it is consistent with theoretical models of compact dust aggregates. We discuss the origin of the detected features and find that low mass (\textless= 1 M-Jup) nascent planets are a possible explanation.

Steepening of the 820 μm continuum surface brightness profile signals dust evolution in TW Hydrae’s disk

Molecular snow lines in protoplanetary disks have been studied theoretically for decades because of their importance in shaping planetary architectures and compositions. The water snow line lies in the planet formation region at < 10 AU, and so far its location has been estimated only indirectly from spatially-unresolved spectroscopy. This work presents a proof-of-concept method to directly image the water snow line in protoplanetary disks through its physical and chemical imprint in the local dust properties. We adopt a physical disk model that includes dust coagulation, fragmentation, drift, and a change in fragmentation velocities of a factor 10 between dry silicates and icy grains as found by laboratory work. We find that the presence of a water snow line leads to a sharp discontinuity in the radial profile of the dust emission spectral index {\alpha}_mm, due to replenishment of small grains through fragmentation. We use the ALMA simulator to demonstrate that this effect can be observed in protoplanetary disks using spatially-resolved ALMA images in two continuum bands. We explore the model dependence on the disk viscosity and find that the spectral index reveals the water snow line for a wide range of conditions, with opposite trends when the emission is optically thin rather than thick. If the disk viscosity is low ({\alpha}_visc < 10^-3) the snow line produces a ring-like structure with a minimum at {\alpha}_mm ~ 2 in the optically thick regime, possibly similar to what has been measured with ALMA in the innermost region of the HL Tau disk.

Testing particle trapping in transition disks with ALMA

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.