# Beyond static torques: Self-consistent migration and dynamical trapping of low-mass planets in dusty disks **Authors:** Jose L. Gomez, Pablo Benitez-Llambay, Ximena S. Ramos, Octavio M. Guilera **arXiv:** [2610.27071](https://arxiv.org/abs/2610.27071) **Date:** 2026-10-19 **Summary page:** https://www.summarizepaper.com/fr/arxiv-id/2610.27071v1/ ## Key points - Objectif : étudier comment les solides des disques protoplanétaires modifient la migration des planètes de faible masse et peuvent empêcher leur chute vers l’étoile. - Méthode : 144 simulations hydrodynamiques bidimensionnelles avec FARGO3D suivent la rétroaction entre migration planétaire, distribution des solides et couples gravitationnels. - Trois régimes apparaissent : migration vers l’intérieur, migration vers l’extérieur parfois jusqu’à l’emballement, et confinement radial oscillatoire. - Résultat central : dans les disques riches en solides, notamment avec un rapport poussières/gaz de 0,05 et un nombre de Stokes de 0,55, des planètes de 2,67 à 4 masses terrestres peuvent être piégées dans des attracteurs dynamiques. - Deux simulations de planètes de 4 masses terrestres parties de positions différentes convergent vers une même distance moyenne d’environ 15,5 unités astronomiques ; d’autres confinements apparaissent plus loin dans les disques moins denses. - Mécanisme : un filament de solides devant la planète et une cavité derrière elle produisent un couple positif. Leur destruction puis leur reconstruction entraînent une alternance de migrations vers l’intérieur et vers l’extérieur. - Les particules de nombre de Stokes 0,10 favorisent presque toujours la migration vers l’intérieur ; les planètes de 10 masses terrestres présentent davantage d’instabilités et d’emballements, sans seuil de masse universel. - Les couples calculés sur orbite fixe ne suffisent pas à prédire l’évolution des planètes ; un confinement peut émerger sans maximum de pression préexistant, notamment au-delà de la ligne des glaces. - Limites : modèle bidimensionnel, population unique de solides et thermodynamique simplifiée. Certaines simulations deviennent peu fiables au-delà d’une excentricité d’environ 0,2, et la stabilité à très long terme reste à vérifier. ## Summary Cet article étudie la manière dont la composante solide des disques protoplanétaires modifie la migration des planètes de faible masse et peut empêcher leur chute rapide vers l’étoile. Les modèles classiques prédisent généralement une migration de type I vers l’intérieur, sous l’effet du couple gravitationnel négatif exercé par le gaz. Des travaux antérieurs ont montré que les poussières peuvent produire un couple positif compensateur, mais ils portaient principalement sur des planètes maintenues sur une orbite fixe. Les auteurs examinent ici la rétroaction complète : le déplacement de la planète transforme la distribution des solides, laquelle modifie à son tour les couples et l’évolution orbitale. L’étude repose sur 144 simulations hydrodynamiques bidimensionnelles réalisées avec FARGO3D. Le gaz et les solides sont modélisés comme deux fluides : le premier est localement isotherme, tandis que le second est dépourvu de pression. Le modèle intègre les échanges de quantité de mouvement entre ces fluides ainsi que la diffusion turbulente des poussières. Six paramètres sont explorés : la position initiale de la planète, fixée à 1 ou 10 rayons de référence ; sa masse, égale à 2,67, 4 ou 10 masses terrestres ; le nombre de Stokes des particules, égal à 0,10, 0,26 ou 0,55 ; le rapport poussières/gaz, fixé à 0,01 ou 0,05 ; la viscosité turbulente, fixée à 0,0001 ou 0,001 ; et deux normalisations de la densité surfacique du gaz, différant d’un facteur proche de huit. Un maillage de 1 024 cellules radiales et de 2 048 cellules azimutales accompagne la planète afin de limiter la diffusion numérique et de préserver les structures coorbitales. Trois comportements principaux se dégagent : une migration monotone vers l’intérieur, une migration vers l’extérieur pouvant s’accélérer jusqu’à l’emballement et un confinement radial oscillatoire. Ce dernier constitue le résultat central de l’étude. Dans les environnements riches en solides, notamment lorsque le rapport poussières/gaz atteint 0,05 et que le nombre de Stokes des particules vaut 0,55, les planètes de 2,67 à 4 masses terrestres peuvent être capturées dans des attracteurs dynamiques. Deux simulations de référence, portant sur une planète de 4 masses terrestres initialement située à 1 ou 10 rayons de référence, convergent ainsi vers une même position moyenne de 5,16 rayons. Avec l’échelle physique adoptée, soit 3 unités astronomiques par rayon de référence, cette position correspond à environ 15,5 unités astronomiques. Dans les disques moins denses, des confinements apparaissent également entre 14 et 27 rayons de référence, notamment autour de 60 unités astronomiques. Le mécanisme de piégeage repose sur la régénération cyclique d’un filament dense de solides devant la planète et d’une cavité derrière elle. Cette asymétrie produit un couple positif. Lorsque la migration vers l’extérieur devient trop rapide, le filament ne se reconstitue plus assez vite : le couple gazeux reprend alors le dessus et entraîne la planète vers l’intérieur. La reconstruction de l’asymétrie relance ensuite le mouvement vers l’extérieur. La rétroaction des poussières sur le gaz favorise ce cycle malgré la diffusion turbulente. À l’inverse, les particules dont le nombre de Stokes vaut 0,10 conduisent presque toujours à une migration vers l’intérieur, avec seulement 7 exceptions sur 48 simulations. Les cœurs de 10 masses terrestres présentent davantage d’instabilités et d’emballements, sans que cette masse constitue pour autant un seuil universel. Les auteurs concluent que les couples calculés pour des planètes sur orbite fixe ne suffisent pas à prédire le devenir des embryons planétaires. Un confinement peut émerger sans maximum de pression préexistant, en particulier au-delà de la ligne des glaces. La généralisation de ce résultat reste toutefois incertaine : le modèle est bidimensionnel, ne considère qu’une seule population de solides à nombre de Stokes constant et repose sur une thermodynamique simplifiée. La prise en compte de la structure verticale pourrait affaiblir les asymétries et modifier les seuils de masse. Enfin, certaines simulations deviennent peu fiables lorsque l’excentricité dépasse environ 0,2, car les perturbations débordent alors du domaine numérique. La stabilité de ces mécanismes à très long terme nécessite donc des études supplémentaires.