Pith. sign in

REVIEW 3 major objections 5 minor 3 references

Formation of multiple dust rings and gaps in protoplanetary discs by a single migrating planet II: radiative discs and observational signatures

T0 review · 3 major / 5 minor · reviewed 2026-07-14 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A single migrating planet in a low-viscosity disc can open multiple dust rings and gaps that stay visible for at least 400,000 years, even with radiative cooling.

desk verdict Solid radiative extension of intermittent-migration multi-ring work; the abstract over-claims “all EOS” because the adaptive case only jumps after a post-hoc mass cut. read the letter →

arxiv 2606.02734 v2 pith:LUDTG3H4 submitted 2026-06-01 astro-ph.EP

Kim M. Weiskopf , Sören C. Meiners , Alexandros Ziampras , Cornelis P. Dullemond This is my paper · ORCID
classification astro-ph.EP
keywords protoplanetarydiscsplanet–discinteractionsplanetarymigrationdustringsandgapsradiativecoolingvortex-assistedintermittent
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks whether one migrating planet, not a chain of planets, can leave behind the multiple dust rings and gaps that millimetre observations routinely show. In high-resolution two-dimensional simulations that include dust and three different treatments of cooling, the planet repeatedly jumps inward, each jump leaving pressure bumps that trap dust. Those rings and gaps remain bright in synthetic continuum images for at least 400 kyr. Cooling decides how many jumps occur and how the planet migrates, yet it does not control how long the dust structures survive; that lifetime is set by the slow-migration intervals between jumps and by the disc’s low viscosity. Large vortices appear only briefly and are rare once cooling is realistic, so most of the observable structure is nearly axisymmetric. The work therefore argues that multi-ring discs can record the migration history of a single planet, provided models use a cooling prescription that is not merely isothermal.

What carries the argument

Intermittent migration jumps: brief runaway (type-III-like) episodes triggered when the planet catches its own inverse-vortensity maximum and the dynamical corotation torque becomes strongly negative, assisted by Rossby-wave vortices and, when cooling is present, by baroclinic forcing that further modifies vortensity around the horseshoe region.

What would settle it

A fully radiative (flux-limited diffusion) simulation of the same planet and disc that produces either no migration jumps or dust rings that dissipate well before 400 kyr would falsify the central claim that the structures are robust under realistic cooling.

Watch

Extended reading notes

Core claim

For every equation of state tested—isothermal, constant β-cooling, and adaptive β—the planet executes one or more intermittent migration jumps, each of which creates dust rings and gaps that remain visible in radiative-transfer images for at least 400 kyr. Cooling strongly shapes the migration track and the number of jumps, but has no measurable effect on the lifetime of the resulting dust structures.

Load-bearing premise

The local adaptive cooling timescale is assumed to capture the baroclinic forcing and vortensity changes that drive the jumps, even though it is known to overestimate vortensity growth compared with full radiative diffusion.

Editorial extensions

If this is right

  • Multi-ring continuum images can be produced by one planet’s migration history rather than by multiple planets or dead-zone edges.
  • The number of rings and the spacing between them encode the number and timing of migration jumps, which themselves depend on the local cooling timescale.
  • Dust asymmetries from large vortices should be rare and short-lived (~90 kyr), so most observed rings are expected to look nearly circular.
  • Structure lifetime is controlled by the duration of slow vortex-assisted migration, not by cooling, so low-viscosity discs keep fossil rings long after the planet has stalled.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If cooling sets the jump count, then discs of different optical depth or stellar irradiation should systematically show different numbers of rings for planets of similar mass.
  • Kinematic signatures of a single super-thermal planet that has stalled after jumps—deep primary gap plus weaker secondary gap—could be searched for in molecular-line data of multi-ring systems.
  • Because the adaptive-β model overestimates vortensity growth, the true window of planet masses that produce jumps may be narrower than reported here.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper uses high-resolution 2-D FARGOCPT hydrodynamical simulations of a single 100 M⊕ (or 62 M⊕) planet migrating in a low-viscosity (α=10^{-4}) disc from 5–500 au, comparing three equations of state (locally isothermal, constant β=1, and adaptive β_fld). Dust is evolved as Lagrangian super-particles and post-processed into synthetic 1.25 mm continuum images. The central claim is that the planet undergoes one or more intermittent migration jumps for all EOS considered, each jump producing pressure maxima/minima that trap dust into rings and gaps remaining visible for at least ~400 kyr; cooling controls the number and timing of jumps (via baroclinic forcing and vortensity evolution) but not structure lifetimes, while large-scale vortices are short-lived (~90 kyr average) and rare under realistic cooling.

Significance. If the results hold, the work supplies a concrete, observationally testable pathway by which a single migrating planet can generate the multi-ring/gap continuum morphologies commonly seen by ALMA/DSHARP without invoking multiple planets or other mechanisms. The systematic comparison across isothermal, constant-β and adaptive-β cooling, the inclusion of multi-size dust, and the production of beam-convolved intensity maps are genuine strengths that go beyond earlier isothermal studies (McNally et al. 2019; Wafflard-Fernandez & Baruteau 2020; Paper I). The explicit demonstration that structure lifetime tracks the duration of slow (vortex-assisted) migration rather than the cooling timescale itself is a useful, falsifiable prediction. The caveats about 2-D geometry, local β_fld versus full FLD, and the absence of self-gravity are stated, so the paper remains a solid contribution provided the claims are accurately scoped.

major comments (3)
  1. [Abstract, §3.2.3–3.2.4, Conclusions] Abstract and §3.2.3–3.2.4: The unqualified statement that “for all equations of state considered, the planet undergoes one or several migration jumps” is not supported by the fiducial adaptive-β run (model adapt, Mp=100 M⊕). That run stalls in the type-II regime without a jump; a jump appears only after the planet mass is lowered post-hoc to 62 M⊕ (adapt-low-mass) “since the focus of this paper lays on planets which experience a migration jump.” Because the number of jumps is itself mass- and cooling-dependent near the thermal-mass boundary, the abstract, strongest claim and conclusions must be revised to state clearly that realistic cooling produces multi-ring structures only for a restricted mass range, and that the adaptive suite required a mass reduction outside the original parameter set.
  2. [§2.1, §4.1] §2.1 (Eqs. 5–7) and §4.1: The adaptive cooling model uses the local β_fld approximation, which the authors themselves note (citing Ziampras et al. 2024) overestimates vortensity growth relative to full flux-limited diffusion. Because baroclinic forcing and the resulting small-scale vortices control whether the planet catches its own IV structures and jumps, this approximation is load-bearing for the claim that realistic cooling still permits intermittent migration. A quantitative estimate of the bias (or at least a stronger caveat that the adaptive jumps may be artificially facilitated) is required before the multi-structure narrative can be applied to real discs.
  3. [§3.2, Figs. 2 and 11] §3.2 and Fig. 2 / Fig. 11: The transition between intermittent and type-II regimes is stated to depend on both Mp/Mth and the local cooling timescale, yet only two planet masses and a single α are explored for the adaptive case. Given that the entire multi-ring claim for the most realistic EOS rests on the lower-mass run, the paper needs either an additional adaptive run at the original 100 M⊕ with a modestly different initial surface-density normalisation (or opacity) that still produces a jump, or an explicit statement that the adaptive multi-gap morphology is not robust across the thermal-mass boundary.
minor comments (5)
  1. [Title page] Title page and affiliations contain duplicated “Germany” and the incomplete phrase “centre for Astronomie (ZAH), Heidelberg University o Germany Germany”.
  2. [Fig. 1] Fig. 1 caption and text: β_fld is plotted at t=0; a second panel or curve showing β_fld after the gap has opened would help the reader judge how strongly the cooling timescale evolves during the vortex-assisted phase.
  3. [§3.3, Fig. 14] §3.3 and Fig. 14: The “noise” in the synthetic images is correctly attributed to particle binning, but a short quantitative note on the effective surface-density resolution (particles per cell) would reassure readers that the reported ring contrasts are not resolution artefacts.
  4. [§3.1, Eq. (15)] Eq. (15) and surrounding text: the dynamical corotation torque is written with IVh/IVp; a brief reminder that IV is inverse vortensity (already defined in Eq. 4) would improve readability for non-specialists.
  5. [Appendix A] Appendix A: the additional constant-β runs (0.01, 3, 7.5) are useful; stating the corresponding Mp/Mth values at the jump radius would make the comparison with the adaptive suite more direct.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: multi-ring and lifetime claims emerge from forward hydro simulations, not from definitions, fits, or load-bearing self-citation chains.

full rationale

The paper reports results of a suite of 2-D FARGOCPT hydrodynamical simulations (isothermal, constant-β, adaptive β_fld) that include Lagrangian dust and post-processed continuum images. Migration jumps, pressure maxima, dust rings/gaps and their ~400 kyr visibility are measured outcomes of those runs, not quantities obtained by construction from fitted parameters or from definitions that already encode the target. The adaptive-β cooling formula (Eqs. 5–7) and the vortex-assisted / intermittent-migration taxonomy are taken from prior literature (including overlapping-author papers) and used as modelling tools; the simulations themselves decide whether a given mass and EOS produce jumps. The post-hoc lowering of planet mass from 100 M⊕ to 62 M⊕ for the adaptive case is an explicit parameter choice to study the phenomenon of interest, not a circular derivation that forces the abstract claim. No uniqueness theorem is invoked, no observational data are fitted and then re-predicted, and no known empirical pattern is merely renamed. The derivation chain is therefore self-contained against external benchmarks and free of the circular patterns listed in the analyser specification.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The central multi-ring claim rests on standard disc hydrodynamics plus a set of modelling choices (viscosity, planet mass, cooling approximations, 2-D geometry) that are conventional in the field but not derived from first principles inside the paper. No new particles or forces are invented; free parameters are the usual simulation knobs.

free parameters (5)
  • α viscosity = 10^{-4}
    Fixed at 10^{-4}; controls gap depth, migration regime and structure lifetime. Chosen to match recent observational estimates rather than derived.
  • planet mass Mp = 100 M⊕ (62 M⊕ variant)
    Set to 100 M⊕ (or 62 M⊕ for adapt-low-mass) so that Mp ~ 1–2 M_th at 50 au; determines whether intermittent jumps occur.
  • constant β = 1
    β=1 chosen because it maximises radiative damping of secondary gaps; other values (0.01, 3, 7.5) explored only in appendix.
  • initial surface density normalisation = 10 g cm^{-2}
    Σ=10 g cm^{-2} at 50 au sets disc mass ~0.12 M⊙ and Toomre Q>1; conventional outer-disc value.
  • dust sizes / Stokes numbers = 3 µm, 30 µm, 0.3 mm
    Three fixed sizes (3 µm, 30 µm, 0.3 mm) corresponding to St=10^{-4}–10^{-2}; chosen to sample ALMA-relevant grains.
assumptions (4)
  • domain assumption Vertically integrated 2-D Navier–Stokes equations with ideal-gas EOS and Shakura–Sunyaev viscosity adequately capture gap opening and intermittent migration.
    Standard in the planet–disc literature; invoked throughout Sect. 2.
  • domain assumption Local β_fld cooling (surface + mid-plane diffusion) is a sufficient proxy for full flux-limited diffusion for the purpose of baroclinic forcing and migration jumps.
    Adopted from Ziampras et al. (2023); authors note it overestimates vortensity growth (Sect. 2.1, 4.1).
  • domain assumption Gas self-gravity can be neglected because Toomre Q>1 everywhere.
    Stated in Sect. 2.2; disables disc–disc indirect terms.
  • domain assumption Dust is passive (no back-reaction) and follows an MRN size distribution with fixed Stokes numbers.
    Sect. 2.2; simplifies particle tracking.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Formation of multiple dust rings and gaps in protoplanetary discs by a single migrating planet II: radiative discs and observational signatures." pith.science (2026). https://pith.science/paper/LUDTG3H4

@misc{pith2026260602734,
  author       = {Pith},
  title        = {Pith review of: Formation of multiple dust rings and gaps in protoplanetary discs by a single migrating planet II: radiative discs and observational signatures},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LUDTG3H4}},
  note         = {Machine review of arXiv:2606.02734}
}
abstract

Dust structures in protoplanetary discs have been widely observed and their creation remains an active field of research. Several possible origins have already been explored, including magneto-hydrodynamics, shadows and planets-disc interactions. The goal of this paper is to investigate whether a single migrating planet in a low-viscosity disc, including radiative processes, is capable of generating observable dust structures. We aim to examine both the lifetime of such structures and potential asymmetries within them. We perform a set of high-resolution, two-dimensional hydrodynamic simulations of migrating planets using three different equations of state: isothermal, constant $\beta$-cooling and an adaptive $\beta$ model. Dust is included in all simulations and the resulting dust density profiles are then post-processed to create radiative transfer images. For all equations of state considered, the planet undergoes one or several migration jumps, each producing dust rings and gaps. The lifetime of these structures depends on the phase of slow migration preceding and occurring between jumps, but in all cases they remain visible for at least 400 kyr. We find that cooling has a deciding effect on the migration behaviour and the number of jumps, but no measurable influence on the lifetime of the dust structures. The structures exhibit relatively few asymmetries, and large-scale vortices persist for an average of only 90 kyr. Our models highlight the capacity of planets to open multiple gaps while migrating, and stress the importance of a realistic cooling model. Care should be taken when interpreting and comparing such models directly to observations.

Figures

Figures reproduced from arXiv: 2606.02734 by the authors.

Figure 1
Figure 1. Dimensionless cooling timescale βfld over the extend of the disc at t = 0 kyr. We additionally plot the approximations for the optically thick and optically thin limit. Name EOS Mp [M⊕] iso isothermal 100 beta β = 1 100 adapt β = βfld 100 adapt-low-mass β = βfld 62 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Semi-major axis (ap) evolution of model iso, beta and adapt. We plot the inverse vortensity in [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Surface density profile of model iso at t = 500 kyr. A sample size of the dust particles is plotted over the profile [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figures from the paper (9 more)
Figure 6
Figure 6. Figure 6: Small scale vortices at the outer horseshoe edge of model [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: Schematic graphic explaining the baroclinic forcing con [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: Normalized, azimuthally averaged gas surface density profile of model iso and model beta at t = T1 and T2. The dot marks the planet’s position [PITH_FULL_IMAGE:figures/full_fig_p007_8.png]
Figure 9
Figure 9. Figure 9: Surface density profile with dust of model [PITH_FULL_IMAGE:figures/full_fig_p007_9.png]
Figure 11
Figure 11. Figure 11: Semi-major axis (ap) evolution of model adapt-low-mass and adapt. At t ∼ 180 kyr, Rp ∼ 41 au the planet enters the type-II regime, but migrates outwards. At this point, Mp = 2.4Mth and βfld has an approximate value of βfld ≈ 0.146, almost the same as the start￾ing val…
Figure 10
Figure 10. Figure 10: Normalized inverse vortensity profile and baroclinic [PITH_FULL_IMAGE:figures/full_fig_p008_10.png]
Figure 13
Figure 13. Figure 13: Surface density profile with dust of model [PITH_FULL_IMAGE:figures/full_fig_p009_13.png]
Figure 14
Figure 14. Figure 14: Intensity Iν of model iso, beta, adapt and adapt-low-mass, convolved with a Gaussian beam at t = 200, 300, 400, and 500 kyr. The white ellipse represents the beam, while the white cross indicates the planet’s position. The planet in model adapt did not undergo a migra…
Figure 15
Figure 15. Figure 15: Comparison of the convolved intensity image of model [PITH_FULL_IMAGE:figures/full_fig_p012_15.png]

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

3 extracted references

  1. [1]

    M., Huang, J., Pérez, L

    Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018, The Astrophysical Journal, 869, L41, aDS Bibcode: 2018ApJ...869L..41A Ansdell, M., Williams, J. P., Trapman, L., et al. 2018, ApJ, 859, 21 Armitage, P. J. 2020, Astrophysics of Planet Formation, 2nd edn. (Cambridge University Press) Armitage, P. J. & Kley, W. 2019, From Protoplanetary Disks to Planet F...

    2018

  2. [2]

    J., Livio, M., & Pringle, J

    Swiss Society for Astrophysics and Astron- omy, aDS Bibcode: 2019SAAS...45.....A Armitage, P. J., Livio, M., & Pringle, J. E. 2001, MNRAS, 324, 705 Bae, J. & Zhu, Z. 2018a, ApJ, 859, 118 Bae, J. & Zhu, Z. 2018b, ApJ, 859, 119 Béthune, W., Lesur, G., & Ferreira, J. 2016, A&A, 589, A87 Crida, A., Baruteau, C., Gonzalez, J.-F., et al. 2025a, The Open Journal...

    2001

  3. [3]

    peaks forβ∼1 (at least among the values tested), in line with the discussion in Sect. 3.2.2. By plotting the individual components ofS, we can identify that this ribbon is indeed driven by termS 1 (top row) and is mostly sensitive to∂T/∂ϕ(see also Fig. 7). Since∂T/∂ϕ=0 in our locally isothermal models, this term is completely absent there, hence the lack ...

    2024

Pith tools

Reviewed July 14, 2026 · model on record in the stance chip above.