posted on 2024-03-15, 16:39authored byY Wu, C Baruteau, S Nayakshin
ALMA has spatially resolved over 200 annular structures in protoplanetary discs, many of which are suggestive of the presence of planets. Constraining the mass of these putative planets is quite degenerate for it depends on the disc physical properties, and for simplicity a steady-state is often assumed whereby the planet position is kept fixed and there is a constant source of dust at the outer edge of the disc. Here, we argue against this approach by demonstrating how the planet and dust dynamics can lift degeneracies of such steady-state models. We take main disc parameters from the well-known protoplanetary disc HD 163296 with a suspected planet at R ≈ 86 au as an example. By running gas and dust hydrodynamical simulations post-processed with dust radiative transfer calculations, we first find steady-state disc and planet parameters that reproduce ALMA continuum observations fairly well. For the same disc mass, but now allowing the planet to migrate in the simulation, we find that the planet undergoes runaway migration and reaches the inner disc in ∼0.2 Myr. Further, decreasing the disc mass slows down planet migration, but it then also increases the dust's radial drift, thereby depleting the disc dust faster. We find that the opposing constraints of planet migration and dust drift require the disc mass to be at most 0.025 M⊙, must less massive than previously estimated, and for the dust to be porous rather than compact. We propose that similar analysis should be extended to other sources with suspected planetary companions.
Funding
This research used the ALICE High Performance Computing Facility at the University of Leicester, and DiRAC Data Intensive service at Leicester, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility (www.dirac.ac.uk).
DUSTBUSTERS RISE project (grant agreement number 823823)
Astrophysics Research at the University of Leicester
The FARGO3D code is publicly available from http://fargo.in2p3.fr. The RADMC-3D code is available from https://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d. The python package used for the post-processing of the simulations data and for the analysis of the radiative transfer calculations are available via its GitHub repository: fargo2radmc3d (https://github.com/charango/fargo2radmc3d). The data obtained in our simulations can be made available on reasonable request to the corresponding author.