REVIEW 2 major objections 4 minor 29 references
Local three-dimensional simulations of the convective overstability in protoplanetary discs
T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Abstract and body disagree: the promised disc simulations do not appear, and the full text is a spiral vortex-sheet proof.
desk verdict The abstract promises 3D COS simulations, but the body is an unrelated math paper on Prandtl spirals—desk reject as submitted, though the math itself looks serious. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
For the body's actual argument, the load-bearing object is an explicit conformal map that sends a vertical strip onto the exterior of the logarithmic spiral, turning the vortex-sheet uniqueness problem into a non-standard boundary-value problem for holomorphic functions on the strip; the proof then uses the Schwarz reflection principle and the identity theorem to show that two candidate velocity fields must coincide. The abstract's advertised results would instead have rested on three-dimensional shearing-box simulation machinery with a cooling parameter, but no such machinery appears anywhere in the full text.
What would settle it
Searching the full text for the abstract's key terms — 'convective overstability', 'shearing box', 'zonal flow', 'Péclet' — finds none of them, which immediately settles that the abstract's claims are not present in this manuscript. For the body's theorem itself, constructing a divergence-free velocity field that satisfies the decay condition at the origin but violates the velocity matching condition would exhibit exactly the non-uniqueness the theorem is designed to rule out.
Extended reading notes
Core claim
The actual content in this manuscript is a pure-mathematics theorem: given a divergence-free velocity field whose distributional vorticity is supported on a single logarithmic (Prandtl) spiral, and assuming the velocity matching condition and a decay condition at the spiral's origin, the velocity field is unique; the same holds for vorticity supported on a union of concentric logarithmic spirals. The proof reduces the problem via an explicit conformal map from the spiral exterior to a strip, where uniqueness is settled by the Schwarz reflection principle under non-standard boundary conditions, and as a by-product the paper re-derives the known velocity and complex-potential formulas. This is a self-contained result in vortex-sheet theory, but it is entirely disconnected from the abstract's claims about convective overstability in protoplanetary discs.
Load-bearing premise
The load-bearing premise of the abstract's claims is that the body text is the same paper as the abstract; it is actually a different mathematics paper, so the simulation results are unsupported by anything in this manuscript.
Editorial extensions
If this is right
- If the uniqueness theorem holds, the explicit velocity formulas for the Prandtl spiral are the only divergence-free fields compatible with the spiral vorticity once the matching and decay conditions are imposed.
- The pressure matching condition is necessary: without it, the candidate velocity field does not exist, so not every logarithmic spiral is a weak solution of the two-dimensional Euler equations.
- The conformal-map construction yields an alternative derivation of the velocity and complex-potential formulas, independent of earlier contour-integration methods.
- Adding any holomorphic function that vanishes at the origin produces infinitely many divergence-free velocity fields with the same vorticity, which is why the velocity matching condition is indispensable.
- For families of concentric spirals, uniqueness extends, and the pressure matching condition becomes a discrete system of equations in the spiral parameters.
Reading between the lines
- The abstract and the full text are two different documents, so any evaluation that relies on the advertised simulation results should read the document that actually reports them; this manuscript supplies no such report.
- If the convective-overstability results were provided in a corrected submission, the described diversity of saturated states would imply that dust-concentration predictions depend sharply on Reynolds and Péclet numbers, so claims about solid accumulation would need to be regime-specific.
- The conformal-map approach to uniqueness could plausibly extend to other families of vortex sheets, such as spirals with different growth rates or smooth vorticity distributions, giving a broader uniqueness theory beyond logarithmic spirals.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The submission is presented with an astro-ph.EP title and an abstract reporting local three-dimensional SNOOPY simulations of the convective overstability (COS) in protoplanetary discs, with claimed transitions from weakly nonlinear and wave-turbulent states, through bursty cycles involving zonal flows, planar vortices, and elliptical instability, to persistent zonal flows as Reynolds number and Péclet number are varied. The full text, however, is a mathematics paper, arXiv:2508.03554v1 [math.AP], by Cieślak, Kokocki, and Kosewski, on the uniqueness of divergence-free velocity fields associated with Prandtl-spiral vorticity. The body contains no disc model, no simulation setup, no SNOOPY code, no cooling prescription, no numerical diagnostics, and no mention of convective overstability; it develops a conformal-map proof of a uniqueness theorem for vortex-sheet velocity fields. The abstract's claims therefore have no support anywhere in the manuscript.
Significance. If the abstract's results were present and correct, they would be a valuable contribution: local 3D COS simulations showing Re- and Pe-dependent saturated states would bear directly on the viability of COS-induced vortices as dust-concentration sites and on the relation between COS and the subcritical baroclinic instability. None of that content is in the manuscript. On its own terms, the body appears to be a serious mathematical contribution: it constructs an explicit conformal map from a strip to the exterior of a logarithmic spiral, reduces uniqueness to a non-standard boundary value problem, and identifies in Remark 1.2 the precise role of the velocity matching condition. These mathematical strengths do not rescue the submission, because the body addresses a different question and the claimed simulation results remain completely unverifiable.
major comments (2)
- [Full text (header and body)] The central claims of the abstract—3D SNOOPY simulations of convective overstability, bursty cycles involving zonal flows, planar vortices, and elliptical instability, and the Re/Pe dependence of saturated states—are not stated, derived, or supported anywhere in the body. The body opens with the header 'arXiv:2508.03554v1 [math.AP]' and is a mathematics paper on uniqueness of divergence-free velocity fields whose vorticity is supported on Prandtl spirals. It contains no model equations, no grid or box parameters, no cooling scheme, no diagnostics, no numerical data, and no reference to protoplanetary discs or convective overstability. The abstract's results therefore rest entirely on content that is absent from the manuscript, so the submission does not present the paper it describes.
- [Abstract vs. Sections 1–8] The claimed parameter study is not merely under-reported; it is absent. The abstract states that 'at higher Re, but low Peclet number (Pe), we obtain bursty cycles...' and that 'for larger Pe... zonal flows can persist,' yet no definition of Re or Pe, no diagnostic for zonal flows, no elliptical-instability diagnostic, and no figures or tables of simulation outcomes appear in Sections 1–8. The body's Theorem 1.1 and Propositions 3.2, 4.1, 6.1, 7.1, and 8.1 concern velocity potentials for spiral vortex sheets, not convective overstability. This is a load-bearing gap in the submitted artifact as a whole.
minor comments (4)
- [Abstract] The abstract contains the typo 'axisymmetic' for 'axisymmetric'.
- [Full text header] The body is labeled arXiv:2508.03554v1 [math.AP], while the submission is arXiv:2508.03557 (astro-ph.EP); the manuscript must be associated with a single consistent identifier, title, and author list.
- [Section 5, Eq. (5.3) area] The alternative derivation of the velocity formula is explicitly conditional on the pressure matching condition (4.1) imported from [3]; this should be stated as a limitation if the mathematical body is to be read as a standalone paper.
- [Throughout] The full text contains encoding artifacts (e.g., the author string 'PRZEMYS/suppress LAW KOSEWSKI' and scattered garbled characters) that make the text difficult to read.
Circularity Check
No circular derivation in the supplied body; the artifact-level abstract/body mismatch is an absence-of-content issue, not circularity.
full rationale
The supplied manuscript body is arXiv:2508.03554 [math.AP], a mathematics paper by Cieślak, Kokocki, and Kosewski on the uniqueness of divergence-free velocity fields supported on Prandtl spiral vortex sheets. Its claimed result is a conditional uniqueness theorem, and the proof is self-contained: the uniqueness argument in Sections 3 and 6 uses only the stated jump, matching, and decay conditions together with an explicit conformal map to a strip; it does not assume the target velocity formula. Section 5 re-derives the known formula (1.3) under the pressure matching condition (4.1) imported from the authors' earlier work [3], but the paper explicitly presents this as an alternative derivation and verification, not as a new prediction, so no fitted input is being renamed as an output. The heavy citation of [3] supplies the Prandtl spiral construction and parameter conditions; it does not supply the uniqueness conclusion, which is proven here. Thus there is no circular step within the mathematical body. However, the abstract provided for this submission is for arXiv:2508.03557 [astro-ph.EP], a convective-overstability simulation paper by Teed and Latter, and none of the abstract's claims about SNOOPY runs, zonal flows, or elliptical instability appear in the body. That is a serious missing-support problem at the artifact level, but it is an absence of content rather than a derivation that reduces to its own inputs, so it does not raise the circularity score.
Assumptions & free parameters
free parameters (3)
- Reynolds number Re =
varied across runs (not fitted)
- Peclet number Pe =
varied across runs (not fitted)
- Numerical grid, box size, and resolution of the 3D runs =
not stated in the provided text
assumptions (5)
- ad hoc to paper The velocity matching condition (1.14) and the decay condition at the origin (1.8) define the admissible class of velocity fields
- domain assumption The pressure matching condition (4.1), taken from the authors' prior work [3, eq. (1.18)], holds
- domain assumption Self-similar velocity ansatz u(z,t) = (1/t)U(z/t) is inherited from [3]
- standard math Standard results of complex analysis: Schwarz reflection principle and the identity theorem
- domain assumption The local shearing-box model with an imposed adverse entropy gradient and Peclet-number cooling represents COS-relevant disc regions
Cite this review
Pith. "Pith review of Local three-dimensional simulations of the convective overstability in protoplanetary discs." pith.science (2026). https://pith.science/paper/4SKREDDL
@misc{pith2026250803557,
author = {Pith},
title = {Pith review of: Local three-dimensional simulations of the convective overstability in protoplanetary discs},
year = {2026},
howpublished = {\url{https://pith.science/paper/4SKREDDL}},
note = {Machine review of arXiv:2508.03557}
}
read the original abstract
At certain radii protoplanetary discs may sustain a form of oscillatory convection (`convective overstability'; COS) due to localised adverse entropy gradients. The resulting hydrodynamical activity can produce coherent structures, such as zonal flows and vortices, that may concentrate solid material and aid their further coagulation. In this paper we extend previous axisymmetric runs by performing local three-dimensional simulations of the COS, using the code SNOOPY. As parameters are varied, we characterise how the various axisymmetric COS saturated states are transformed in 3D, while also tracking their interrelationship with the subcritical baroclinic instability. In particular, at low Reynolds number (Re) our 3D simulations exhibit similar weakly nonlinear and wave turbulent states to our earlier axisymmetic runs. At higher Re, but low Peclet number (Pe), we obtain bursty cycles involving the creation of zonal flows, the subsequent development of planar vortices, and their destruction by elliptical instability. For larger Pe, however, zonal flows can persist, alongside weaker more elongated vortices. These results further reveal the diversity of the COS's behaviour, and show that solid accumulation via COS-induced vortices may not be straightforward.
Reference graph
Works this paper leans on
-
[1]
Chang E., Youdin A. N., 2024, @doi [ ] 10.3847/1538-4357/ad81f6 , https://ui.adsabs.harvard.edu/abs/2024ApJ...976..100C 976, 100
-
[2]
534, Protostars and Planets VII
Drazkowska J., et al., 2023, in Inutsuka S., Aikawa Y., Muto T., Tomida K., Tamura M., eds, Astronomical Society of the Pacific Conference Series Vol. 534, Protostars and Planets VII. p. 717 ( @eprint arXiv 2203.09759 ), @doi 10.48550/arXiv.2203.09759
-
[3]
pp 391--413, @doi 10.1051/eas/1982035
Fromang S., Lesur G., 2019, in EAS Publications Series. pp 391--413, @doi 10.1051/eas/1982035
-
[4]
Goldreich P., Lynden-Bell D., 1965, @doi [ ] 10.1093/mnras/130.2.125 , https://ui.adsabs.harvard.edu/abs/1965MNRAS.130..125G 130, 125
-
[5]
Hawley J. F., Gammie C. F., Balbus S. A., 1995, @doi [ ] 10.1086/175311 , https://ui.adsabs.harvard.edu/abs/1995ApJ...440..742H 440, 742
doi:10.1086/175311 1995
-
[6]
Klahr H. H., Bodenheimer P., 2003, @doi [ ] 10.1086/344743 , https://ui.adsabs.harvard.edu/abs/2003ApJ...582..869K 582, 869
doi:10.1086/344743 2003
-
[7]
Klahr H., Hubbard A., 2014, @doi [ ] 10.1088/0004-637X/788/1/21 , https://ui.adsabs.harvard.edu/abs/2014ApJ...788...21K 788, 21
-
[8]
Klahr H., Baehr H., Melon Fuksman J. D., 2023, @doi [arXiv e-prints] 10.48550/arXiv.2305.08165 , https://ui.adsabs.harvard.edu/abs/2023arXiv230508165K p. arXiv:2305.08165
work page Pith review arXiv doi:10.48550/arxiv.2305.08165 2023
Show all 29 references
-
[9]
N., 2016, @doi [ ] 10.1093/mnras/stv2449 , https://ui.adsabs.harvard.edu/abs/2016MNRAS.455.2608L 455, 2608
Latter H. N., 2016, @doi [ ] 10.1093/mnras/stv2449 , https://ui.adsabs.harvard.edu/abs/2016MNRAS.455.2608L 455, 2608
2016 doi
-
[10]
N., Papaloizou J., 2017, @doi [ ] 10.1093/mnras/stx2038 , https://ui.adsabs.harvard.edu/abs/2017MNRAS.472.1432L 472, 1432
Latter H. N., Papaloizou J., 2017, @doi [ ] 10.1093/mnras/stx2038 , https://ui.adsabs.harvard.edu/abs/2017MNRAS.472.1432L 472, 1432
2017 doi
-
[11]
Lehmann M., Lin M.-K., 2024, @doi [ ] 10.3847/1538-4357/ad496c , https://ui.adsabs.harvard.edu/abs/2024ApJ...970...15L 970, 15
2024 doi
-
[12]
Lehmann M., Lin M.-K., 2025, @doi [ ] 10.3847/1538-4357/adbf06 , https://ui.adsabs.harvard.edu/abs/2025ApJ...983..135L 983, 135
2025 doi
-
[13]
Lesur G., 2021, @doi [Journal of Plasma Physics] 10.1017/S0022377820001002 , https://ui.adsabs.harvard.edu/abs/2021JPlPh..87a2001P 87, 205870101
2021 doi
-
[14]
Lesur G., Longaretti P.-Y., 2005, Astronomy & Astrophysics, 444, 25
2005
-
[15]
Lesur G., Longaretti P.-Y., 2007, Monthly Notices of the Royal Astronomical Society, 378, 1471
2007
-
[16]
Lesur G., Papaloizou J. C. B., 2009, @doi [ ] 10.1051/0004-6361/200811577 , https://ui.adsabs.harvard.edu/abs/2009A&A...498....1L 498, 1
2009 doi
-
[17]
Lesur G., Papaloizou J. C. B., 2010, @doi [ ] 10.1051/0004-6361/200913594 , https://ui.adsabs.harvard.edu/abs/2010A&A...513A..60L 513, A60
2010 doi
- [18]
-
[19]
Lithwick Y., 2007, @doi [ ] 10.1086/522074 , https://ui.adsabs.harvard.edu/abs/2007ApJ...670..789L 670, 789
2007 doi
-
[20]
Lyra W., 2014, The Astrophysical Journal, 789, 77
2014
-
[21]
M., 2019, @doi [ ] 10.1088/1538-3873/aaf5ff , https://ui.adsabs.harvard.edu/abs/2019PASP..131g2001L 131, 072001
Lyra W., Umurhan O. M., 2019, @doi [ ] 10.1088/1538-3873/aaf5ff , https://ui.adsabs.harvard.edu/abs/2019PASP..131g2001L 131, 072001
2019 doi
-
[22]
B., Umurhan O
Lyra W., Yang C.-C., Simon J. B., Umurhan O. M., Youdin A. N., 2024, @doi [ ] 10.3847/2041-8213/ad5af6 , https://ui.adsabs.harvard.edu/abs/2024ApJ...970L..19L 970, L19
2024 doi
-
[23]
R., Julien K., Stewart G
Petersen M. R., Julien K., Stewart G. R., 2007, @doi [ ] 10.1086/511513 , https://ui.adsabs.harvard.edu/abs/2007ApJ...658.1236P 658, 1236
2007 doi
-
[24]
Raettig N., Lyra W., Klahr H., 2021, @doi [ ] 10.3847/1538-4357/abf739 , https://ui.adsabs.harvard.edu/abs/2021ApJ...913...92R 913, 92
2021 doi
-
[25]
D., Papaloizou J
Railton A. D., Papaloizou J. C. B., 2014, @doi [ ] 10.1093/mnras/stu2060 , https://ui.adsabs.harvard.edu/abs/2014MNRAS.445.4409R 445, 4409
2014 doi
-
[26]
J., Latter H
Teed R. J., Latter H. N., 2021, @doi [ ] 10.1093/mnras/stab2311 , https://ui.adsabs.harvard.edu/abs/2021MNRAS.507.5523T 507, 5523
2021 doi
-
[27]
J., Fromang S., Gammie C., Klahr H., Lesur G., Wardle M., Bai X
Turner N. J., Fromang S., Gammie C., Klahr H., Lesur G., Wardle M., Bai X. N., 2014, in Beuther H., Klessen R. S., Dullemond C. P., Henning T., eds, Protostars and Planets VI. p. 411 ( @eprint arXiv 1401.7306 ), @doi 10.2458/azu\_uapress\_9780816531240-ch018
2014 arXiv
-
[28]
I., 2016, @doi [ ] 10.1093/mnras/stw2238 , https://ui.adsabs.harvard.edu/abs/2016MNRAS.463.3725V 463, 3725
Vanon R., Ogilvie G. I., 2016, @doi [ ] 10.1093/mnras/stw2238 , https://ui.adsabs.harvard.edu/abs/2016MNRAS.463.3725V 463, 3725
2016 doi
-
[29]
write newline
" write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.stat...
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.