Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains.

Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD142527 and HD100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr44. With a fairly axially symmetric ring in the submm radio continuum, DoAr44 is likely also a warped system. We constrain the warp geometry by comparing radiative transfer predictions with the DPI data in H band (Q_\phi(H)) and with a re-processing of archival 336GHz ALMA observations. The observed DPI shadows have coincident radio counterparts, but the intensity drops are much deeper in Q_\phi(H) (~88%), compared to the shallow drops at 336GHz (~24%). Radiative transfer predictions with an inner disk tilt of ~30+-5deg approximately account for the observations. ALMA long-baseline observations should allow the observation of the warped gas kinematics inside the cavity of DoAr44.