pith. sign in

arxiv: 2606.26568 · v1 · pith:JVPVKVSInew · submitted 2026-06-25 · 🌌 astro-ph.EP

Ring Spacing from a Fourth-Order Radial Feedback Green Function in a Keplerian Accretion Disk

Yiwei Bao , Can Cui This is my paper

Pith reviewed 2026-06-26 04:07 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords protoplanetary disksring structuresGreen functionradial feedbackKeplerian accretion diskALMA observationsaccretion diskssurface density
0
0 comments X

The pith

A fourth-order radial feedback closure on disk surface density produces a Green function with damped oscillatory rings spaced about 12 AU.

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

The paper develops an exactly solvable local model for the concentric rings observed by ALMA in protoplanetary disks. Passive advection and diffusion separate into a vertical Sturm-Liouville spectrum and a radial modified-Bessel Green function that remains smooth and non-periodic. Adding a minimal fourth-order radial feedback closure to the ring-averaged surface density alters the Green function, producing a damped oscillatory exterior branch once the decaying spatial roots turn complex. Under an observationally motivated AU scaling this yields an illustrative ring spacing of order 12 AU generated entirely inside the disk. The vertical diffusivity profile changes only the passive kernel and does not source the oscillation.

Core claim

For a localized steady ring source the resulting Green function has a damped oscillatory exterior branch whenever the decaying spatial roots are complex. Under an observationally motivated AU scaling the same dimensionless solution gives an illustrative spacing of order 12 AU, distinct from the fastest-growing temporal wavelength. The model separates the passive transport kernel from the feedback mechanism that selects the observable radial scale.

What carries the argument

The minimal fourth-order radial feedback closure for the ring-averaged surface density, which modifies the passive modified-Bessel kernel to admit complex decaying spatial roots and thereby an oscillatory exterior branch.

If this is right

  • Ring spacing is set by the feedback closure rather than by the passive advection-diffusion kernel alone.
  • The oscillation arises internally to the disk and does not require pre-existing planets.
  • The vertical diffusivity profile affects the background kernel but leaves the source of the Green-function oscillation unchanged.
  • The observable spacing differs from the wavelength of fastest temporal growth.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Varying disk mass or radius while keeping the same closure would predict which systems develop rings and at what characteristic scales.
  • A time-dependent or nonlinear version of the feedback could be used to follow how rings grow, migrate, or merge.
  • Comparing the model spacing to observed rings in disks with independently measured vertical structure would test whether the passive kernel and feedback scale separately as claimed.

Load-bearing premise

A minimal fourth-order radial feedback closure is the appropriate modeling choice for the ring-averaged surface density and produces the desired complex roots.

What would settle it

A derivation from disk microphysics showing that the effective radial feedback is not fourth-order or does not produce complex decaying spatial roots, or ALMA measurements of ring spacings that fail to follow the predicted scaling with orbital distance.

Figures

Figures reproduced from arXiv: 2606.26568 by Can Cui, Yiwei Bao.

Figure 1
Figure 1. Figure 1: FIG. 1. Exact steady Green response of the fourth-order ra [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
read the original abstract

ALMA observations of protoplanetary disks reveal ubiquitous concentric ring structures whose origin remains debated. We present an exactly solvable local model for ring spacing in a Keplerian disk. The passive advection--diffusion problem with a general vertical diffusivity profile separates into a vertical Sturm--Liouville spectrum and a radial modified-Bessel Green function. This passive kernel is smooth and does not by itself generate a periodic ring train. We therefore introduce a minimal fourth-order radial feedback closure for the ring-averaged surface density. For a localized steady ring source, the resulting Green function has a damped oscillatory exterior branch whenever the decaying spatial roots are complex. Under an observationally motivated AU scaling, the same dimensionless solution gives an illustrative spacing of order $12$~AU, distinct from the fastest-growing temporal wavelength. The model separates the passive transport kernel from the feedback mechanism that selects the observable radial scale. Because this mechanism is internal to the disk, it does not require pre-existing planets. The vertical diffusivity profile affects the background kernel but is not the source of the Green-function oscillation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The paper claims to provide an exactly solvable local model for ring spacing in Keplerian accretion disks. The passive advection-diffusion problem with a general vertical diffusivity profile separates into a vertical Sturm-Liouville spectrum and a radial modified-Bessel Green function; this passive kernel is smooth. A minimal fourth-order radial feedback closure is then introduced for the ring-averaged surface density. For a localized steady ring source the resulting Green function possesses a damped oscillatory exterior branch when the decaying spatial roots are complex. Under an observationally motivated AU scaling the same dimensionless solution yields an illustrative spacing of order 12 AU, distinct from the fastest-growing temporal wavelength. The model separates the passive transport kernel from the feedback mechanism that selects the observable radial scale and asserts that this mechanism is internal to the disk, thereby not requiring pre-existing planets.

Significance. If the fourth-order radial feedback closure were derived from vertical structure, turbulence, or a conservation law rather than introduced as a modeling choice, the construction would supply a mathematically clean, internal mechanism capable of generating concentric rings without planets and would cleanly separate the passive kernel from the scale-selection feedback. In its present form the result remains an illustrative mathematical example whose oscillation is imposed by the choice of operator and scaling.

major comments (3)
  1. [Abstract / model-setup paragraph] Abstract (and the paragraph introducing the closure): the fourth-order radial feedback operator and its coefficients are introduced because 'the passive kernel is smooth and does not by itself generate a periodic ring train'; no derivation from the vertical Sturm-Liouville eigenfunctions, diffusivity profile, or any moment closure is supplied. Consequently the damped oscillatory exterior branch follows by construction once the operator is imposed, rendering the claim that the mechanism is 'internal to the disk' dependent on an unmotivated modeling assumption.
  2. [Abstract] Abstract: the reported spacing of order 12 AU is obtained by applying an 'observationally motivated AU scaling' to the dimensionless solution rather than emerging as a parameter-free prediction from the governing equations. The abstract explicitly labels the spacing 'illustrative,' which undercuts any claim that the model predicts a specific radial scale from disk physics alone.
  3. [model formulation (presumed §2-3)] The separation into vertical Sturm-Liouville spectrum and radial Green function is presented as exact for the passive problem, yet the subsequent feedback closure is applied only to the ring-averaged surface density; the manuscript does not demonstrate that this averaging is consistent with the vertical eigenfunction expansion or that the fourth-order operator is the minimal one compatible with any conservation law.
minor comments (1)
  1. Notation for the feedback coefficients and the precise definition of the ring-averaged surface density should be introduced with an equation number at first use to avoid ambiguity when the operator is later applied to the Green function.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the detailed and constructive report. The comments correctly identify that the fourth-order closure is a modeling choice rather than a derived result. We respond to each major comment below and indicate where revisions will be made.

read point-by-point responses

  1. Referee: [Abstract / model-setup paragraph] Abstract (and the paragraph introducing the closure): the fourth-order radial feedback operator and its coefficients are introduced because 'the passive kernel is smooth and does not by itself generate a periodic ring train'; no derivation from the vertical Sturm-Liouville eigenfunctions, diffusivity profile, or any moment closure is supplied. Consequently the damped oscillatory exterior branch follows by construction once the operator is imposed, rendering the claim that the mechanism is 'internal to the disk' dependent on an unmotivated modeling assumption.

    Authors: We agree that the fourth-order operator is posited as a minimal closure rather than derived from the vertical eigenfunctions or a conservation law. The passive kernel is indeed smooth, and the oscillatory exterior solution is a direct consequence of the chosen operator. We have revised the abstract and the paragraph introducing the closure to state explicitly that this is a phenomenological modeling assumption and to qualify the description of the mechanism as internal to the disk. revision: yes

  2. Referee: [Abstract] Abstract: the reported spacing of order 12 AU is obtained by applying an 'observationally motivated AU scaling' to the dimensionless solution rather than emerging as a parameter-free prediction from the governing equations. The abstract explicitly labels the spacing 'illustrative,' which undercuts any claim that the model predicts a specific radial scale from disk physics alone.

    Authors: The abstract already qualifies the result as an 'illustrative spacing' obtained under an 'observationally motivated AU scaling'. The governing equations yield a dimensionless spacing set by the complex roots; the numerical value is provided only for context and is not presented as a parameter-free prediction. The existing wording is therefore appropriate and no revision is required. revision: no

  3. Referee: [model formulation (presumed §2-3)] The separation into vertical Sturm-Liouville spectrum and radial Green function is presented as exact for the passive problem, yet the subsequent feedback closure is applied only to the ring-averaged surface density; the manuscript does not demonstrate that this averaging is consistent with the vertical eigenfunction expansion or that the fourth-order operator is the minimal one compatible with any conservation law.

    Authors: The passive separation is exact. The feedback is applied to the ring-averaged surface density as an effective radial description after vertical integration. We acknowledge that a detailed demonstration of consistency between the averaging and the eigenfunction expansion is not supplied. We will add a short discussion of this approximation and of the status of the fourth-order operator as a minimal closure in the revised manuscript. revision: partial

standing simulated objections not resolved
  • Derivation of the specific fourth-order radial feedback operator from the vertical Sturm-Liouville eigenfunctions, diffusivity profile, or a conservation law (it is introduced as a minimal modeling choice).

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The model rests on the separation of the advection-diffusion operator and on the introduction of an ad-hoc fourth-order feedback closure whose justification is not supplied in the abstract.

free parameters (2)
  • fourth-order feedback coefficients and order
    The specific order and form of the radial feedback closure are introduced to generate complex roots and oscillation; no first-principles derivation is given.
  • AU scaling factor
    An observationally motivated scaling converts the dimensionless solution into a physical spacing of order 12 AU.
axioms (2)
  • domain assumption The passive advection-diffusion problem with general vertical diffusivity separates into a vertical Sturm-Liouville spectrum and a radial modified-Bessel Green function.
    Stated as the starting point that yields a smooth passive kernel.
  • domain assumption The disk is Keplerian.
    Standard assumption for the background rotation profile.
invented entities (1)
  • minimal fourth-order radial feedback closure no independent evidence
    purpose: To produce a damped oscillatory exterior branch in the Green function that selects ring spacing
    Introduced as a modeling device without derivation from microphysical processes; independent evidence is not provided.

pith-pipeline@v0.9.1-grok · 5714 in / 1531 out tokens · 63149 ms · 2026-06-26T04:07:34.449233+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references

  1. [1]

    ALMA Partnership et al., Astrophys. J. Lett.808, L3 (2015), URLhttps://ui.adsabs.harvard.edu/abs/ 2015ApJ...808L...3A/abstract

    2015

  2. [2]

    S. M. Andrews et al., Astrophys. J. Lett.869, L41 (2018), URLhttps://ui.adsabs.harvard.edu/ abs/2018ApJ...869L..41A/abstract

    2018

  3. [3]

    S. M. Andrews et al., Astrophys. J. Lett.820, L40 (2016), URLhttps://ui.adsabs.harvard.edu/ abs/2016ApJ...820L..40A/abstract

    2016

  4. [4]

    Isella et al., Phys

    A. Isella et al., Phys. Rev. Lett.117, 251101 (2016), URLhttps://ui.adsabs.harvard.edu/abs/ 2016PhRvL.117y1101I/abstract

    2016

  5. [5]

    Fedele et al., Astron

    D. Fedele et al., Astron. Astrophys.610, A24 (2018), URLhttps://ui.adsabs.harvard.edu/abs/ 2018A%26A...610A..24F/abstract

    2018

  6. [6]

    Huang et al., Astrophys

    J. Huang et al., Astrophys. J. Lett.869, L42 (2018), URLhttps://ui.adsabs.harvard.edu/abs/ 2018ApJ...869L..42H/abstract

    2018

  7. [7]

    R. Dong, Z. Zhu, R. R. Rafikov, and J. M. Stone, Astro- phys. J. Lett.809, L5 (2015), URLhttps://ui.adsabs. harvard.edu/abs/2015ApJ...809L...5D/abstract

    2015

  8. [8]

    Z. Zhu, R. Dong, J. M. Stone, and R. R. Rafikov, As- trophys. J.813, 88 (2015), URLhttps://ui.adsabs. harvard.edu/abs/2015ApJ...813...88Z/abstract

    2015

  9. [9]

    D. M. Segura-Cox, A. Schmiedeke, J. E. Pineda, I. W. Stephens, M. Fern´ andez-L´ opez, L. W. Looney, P. Caselli, Z.-Y. Li, L. G. Mundy, W. Kwon, et al., Na- ture586, 228 (2020), URLhttps://doi.org/10.1038/ s41586-020-2779-6

    2020

  10. [10]

    Lambrechts and A

    M. Lambrechts and A. Johansen, Astron. Astrophys. 572, A107 (2014), URLhttps://doi.org/10.1051/ 0004-6361/201424343

    2014

  11. [11]

    C. F. Gammie, Astrophys. J.457, 355 (1996), URLhttps://ui.adsabs.harvard.edu/abs/1996ApJ. ..457..355G/abstract

    1996

  12. [12]

    Hasegawa and R

    Y. Hasegawa and R. E. Pudritz, Astrophys. J. Lett. 710, L167 (2010), URLhttps://ui.adsabs.harvard. edu/abs/2010ApJ...710L.167H/abstract

    2010

  13. [13]

    Flock, J

    M. Flock, J. P. Ruge, N. Dzyurkevich, T. Hen- ning, H. Klahr, and S. Wolf, Astron. Astrophys.574, A68 (2015), URLhttps://ui.adsabs.harvard.edu/ abs/2015A%26A...574A..68F/abstract

    2015

  14. [14]

    Zhang, G

    K. Zhang, G. A. Blake, and E. A. Bergin, Astro- phys. J. Lett.806, L7 (2015), URLhttps://ui.adsabs. harvard.edu/abs/2015ApJ...806L...7Z/abstract

    2015

  15. [15]

    Johansen, A

    A. Johansen, A. Youdin, and M.-M. Mac Low, Astrophys. J. Lett.704, L75 (2009), URLhttps://ui.adsabs. harvard.edu/abs/2009ApJ...704L..75J/abstract

    2009

  16. [16]

    Bai, Astrophys

    X.-N. Bai, Astrophys. J.791, 137 (2014), URL https://ui.adsabs.harvard.edu/abs/2014ApJ...791. .137B/abstract

    2014

  17. [17]

    Ercolano and I

    B. Ercolano and I. Pascucci, R. Soc. Open Sci.4, 170114 (2017), URLhttps://ui.adsabs.harvard.edu/ abs/2017RSOS....470114E/abstract

    2017

  18. [18]

    Pinilla and A

    P. Pinilla and A. Youdin, inFormation, Evolution, and Dynamics of Young Solar Systems(Springer, 2017), vol. 445 ofAstrophysics and Space Science Library, p. 91, URLhttps://ui.adsabs.harvard.edu/abs/2017ASSL. .445...91P/abstract

    2017

  19. [19]

    A. N. Youdin and J. Goodman, Astrophys. J.620, 459 (2005), URLhttps://ui.adsabs.harvard.edu/ abs/2005ApJ...620..459Y/abstract

    2005

  20. [20]

    Wu and Y

    Y. Wu and Y. Lithwick, Astrophys. J.923, 123 (2021), URLhttps://ui.adsabs.harvard.edu/abs/ 2021ApJ...923..123W/abstract

    2021

  21. [21]

    T. Ueda, M. Flock, and T. Birnstiel, Astrophys. J. Lett. 914, L38 (2021), URLhttps://ui.adsabs.harvard. edu/abs/2021ApJ...914L..38U/abstract

    2021

  22. [22]

    Ziampras, C

    A. Ziampras, C. P. Dullemond, T. Birnstiel, M. Benisty, and R. P. Nelson, Mon. Not. R. Astron. Soc.540, 1185 (2025), URLhttps://academic.oup.com/mnras/ article/540/1/1185/8138670

    2025

  23. [23]

    H. Salo, J. Schmidt, and F. Spahn, Icarus153, 295 (2001), URLhttps://ui.adsabs.harvard.edu/ abs/2001Icar..153..295S/abstract

    2001

  24. [24]

    Schmidt, H

    J. Schmidt, H. Salo, F. Spahn, and O. Petzschmann, Icarus153, 316 (2001), URLhttps://ui.adsabs. harvard.edu/abs/2001Icar..153..316S/abstract

    2001

  25. [25]

    Rein and H

    H. Rein and H. N. Latter, Mon. Not. R. Astron. Soc. 431, 145 (2013), URLhttps://ui.adsabs.harvard. edu/abs/2013MNRAS.431..145R/abstract

    2013

  26. [26]

    A. M. Hughes, G. Duchˆ ene, and B. C. Matthews, Annu. Rev. Astron. Astrophys.56, 541 (2018), URLhttps://ui.adsabs.harvard.edu/abs/2018ARA% 26A..56..541H/abstract

    2018

  27. [27]

    K. A. Rosenfeld, S. M. Andrews, D. J. Wilner, and H. C. Stempels, Astrophys. J.775, 136 (2013), URLhttps://ui.adsabs.harvard.edu/abs/2013ApJ. ..775..136R/abstract

    2013

  28. [28]

    G. I. Taylor, Proc. R. Soc. London A219, 186 (1953), URLhttps://ui.adsabs.harvard.edu/abs/ 1953RSPSA.219..186T/abstract

    1953

  29. [29]

    RING SP ACING FROM A FOUR TH-ORDER RADIAL FEEDBACK GREEN FUNCTION IN A KEPLERIAN ACCRETION DISK

    A. N. Youdin, inPhysical Processes in Circumstel- lar Disks around Young Stars, edited by P. J. V. Garcia (University of Chicago Press, Chicago, 2011), p. 1, URLhttps://press.uchicago.edu/ucp/books/ book/chicago/P/bo11432060.html. 6 SUPPLEMENT AL MA TERIAL FOR “RING SP ACING FROM A FOUR TH-ORDER RADIAL FEEDBACK GREEN FUNCTION IN A KEPLERIAN ACCRETION DISK...

    2011

  30. [30]

    each reported root satisfiesP(q) = 0

  31. [31]

    the two selected right roots have Req <0 and the two selected left roots have Req >0

  32. [32]

    the solved coefficients satisfy continuity off,f ′, andf ′′ at the source

  33. [33]

    the third-derivative jump satisfiesB[f ′′′(0+)−f ′′′(0−)] = 1

  34. [34]

    These checks verify the algebraic Green function and the extrema reported in the main text

    the analytic maxima satisfyf ′(x) = 0 andf ′′(x)<0. These checks verify the algebraic Green function and the extrema reported in the main text. P ASSIVE VER TICAL TRANSPOR T KERNEL This section gives the passive vertical transport calculation used as the background kernel for the local feedback model. The steady passive tracer equation is v· ∇f= 1 r ∂ ∂r ...

  35. [35]

    a radial feedback closure produces the fourth-order linear operator in Eq. (S6)

  36. [36]

    the homogeneous temporal problem hask 2 ∗ =A/(2B) andγ max =−λ+A 2/(4B)

  37. [37]

    (S13), with continuity off,f ′, and f ′′ and the third-derivative jump in Eq

    the steady ring-source problem is the fourth-order Green equation in Eq. (S13), with continuity off,f ′, and f ′′ and the third-derivative jump in Eq. (S17)

  38. [38]

    for a complex right-decaying root pair, the exterior spacing is ∆r Green = 2π/|Imq +|

  39. [39]

    passive vertical diffusion supplies a background transport kernel and motivates the effective coefficients in the local feedback model; it is not the mathematical source of the feedback spacing