REVIEW 3 major objections 4 minor 32 references
A generalized effective potential for differentially rotating plasmas
T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read A generalized effective potential, built from rotation gradients and curvature, predicts the onset of both axisymmetric and non-axisymmetric instabilities in differentially rotating plasma.
desk verdict The generalized effective potential is a correct and useful piece of algebra for m=0, but for non-axisymmetric modes the potential well is built from the very eigenvalue it claims to explain, so the onset story is post-hoc rather than predictive. 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
The central object is the generalized effective potential $U(x,\omega,m) = A''/A - C/A^2$, obtained by rewriting the second-order radial differential equation in the Schrödinger-compatible form $\Psi'' - U\Psi = 0$ through $\Psi = A u$ with $A = \sqrt{\Gamma(x)x/Q(x)}$, $\Gamma(x) = \bar{\omega}^2 - \omega_A^2$, and $Q(x) = k^2 x + m^2$. The machinery packs every destabilizing ingredient, the Doppler-shifted frequency, the Alfvén frequency, rotation shear $\omega_s^2$, curvature $\omega_c^2$, and the radial derivatives of the rotation profile, into one frequency- and mode-dependent function whose sign at each radius decides stability. The stability principle borrowed from potential-well reasoning is that an unstable global mode exists when the potential is negative somewhere between its turning points; here the turning points are the Alfvénic resonances or the domain boundaries.
What would settle it
Run the paper's shooting solver on a rotation profile engineered to have a deep negative region in $\mathrm{Re}(U)$ between the Alfvénic resonances at several trial frequencies and show that no growing non-axisymmetric eigenmode exists; that would break the claimed equivalence between a negative potential well and instability. A growing mode found where $\mathrm{Re}(U)$ is positive across the whole domain would break the converse.
Extended reading notes
Core claim
The central claim is that Equation (7), $U(x,\omega,m) = A''/A - C/A^2$ with $\Gamma = \bar{\omega}^2 - \omega_A^2$ and $Q = k^2 x + m^2$, is a generalized effective potential for non-axisymmetric disturbances: the change of variable $\Psi = A u$ with $A = \sqrt{\Gamma x/Q}$ turns the radial stability equation into $\Psi'' - U\Psi = 0$. Instability sets in where $\mathrm{Re}(U) < 0$, and the mode is confined to the region between two Alfvénic resonances (where the Doppler-shifted frequency equals the Alfvén frequency) or between one resonance and a boundary. In the axisymmetric limit with a purely axial field, the potential reduces to Equation (9), whose negativity at marginal stability reproduces the MRI condition $\omega_A^2 < -\omega_s^2$. For $m=1$ the paper finds two distinct branches, a high-frequency localized MRI and a lower-frequency global magneto-curvature (MCI) mode, each traced back to negative regions of the potential built from the rotation gradients and rotation curvature. The free energy for the non-axisymmetric modes comes from the first and second derivatives of the global rotation profile, so even Rayleigh-stable profiles with $1 < q < 2$ can be non-axisymmetrically unstable when $q(r)$ is localized.
Load-bearing premise
The load-bearing premise is that a negative dip in the potential, evaluated at the mode's own frequency, is what makes the mode unstable; for non-axisymmetric modes that frequency is supplied by the very eigenvalue the potential is meant to explain, so the well describes the found mode rather than independently predicting instability.
Editorial extensions
If this is right
- The axisymmetric limit of the potential recovers the standard MRI condition $\omega_A^2 < -\omega_s^2$, so one diagnostic connects the classical instability to its non-axisymmetric relatives.
- Non-axisymmetric instabilities appear as modes confined in a negative potential well between two Alfvénic resonances, or between a resonance and a boundary, with the mode evanescent outside the well.
- Two distinct $m=1$ branches are separated by the potential picture: a high-frequency localized MRI and a low-frequency global magneto-curvature (MCI) mode.
- Rotation profiles with localized curvature in $q(r)$ can be non-axisymmetrically unstable even at zero magnetic field, despite being Rayleigh-stable.
- Well depth tracks mode localization: deeper wells give more localized modes with smaller growth rates, shallower wells give more global modes with larger growth rates.
Reading between the lines
- A parameter-free, $\omega$-independent version of the criterion would make the potential a true predictive tool, since the well is currently drawn using the eigenfrequency it is meant to explain; the paper provides no such bound.
- The same construction should extend to compressible, resistive, or toroidal settings, where the Alfvénic-resonance confinement may shift; a direct test is whether resistive modes in liquid-metal Couette experiments sit where the ideal potential is negative.
- Because the magneto-curvature mode is global and persists at stronger fields, it may be a more robust source of turbulence and dynamo action in accretion disks than the localized MRI branch, a consequence the paper does not explore.
- The zero-field instability found for localized $q(r)$ profiles suggests a purely hydrodynamic, curvature-driven non-axisymmetric route to instability that could seed magnetic fields without MRI; this possibility is not tested in the paper.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives a normal-form effective potential U(x,ω,m) for linear ideal MHD perturbations of current-free differentially rotating cylindrical plasmas by transforming the eigenvalue ODE (Eq. 4) into Ψ''−UΨ=0 (Eqs. 6–7). It shows that in the axisymmetric limit the instability condition U(x,0)<0 reproduces the standard MRI criterion ω_A²<−ω_s². It then solves the eigenvalue problem numerically for m=0 and m=1 modes with several rotation profiles, identifies two non-axisymmetric branches (a high-frequency local MRI and a low-frequency global magneto-curvature mode), and interprets their confinement in terms of a negative Re(U) well between Alfvénic resonances or between a resonance and a boundary.
Significance. If the predictive interpretation were fully supported, the paper would provide a useful global diagnostic for non-axisymmetric instabilities in accretion disks and laboratory experiments. The algebraic reduction and the m=0 marginal-stability condition are sound, and the numerical shooting solver is benchmarked against NIMROD; the identification of two distinct non-axisymmetric branches with different frequency characters is interesting. However, the non-axisymmetric potential is evaluated at the eigenvalue it is used to explain, so the claimed onset mechanism is not an independent predictor; this limits the significance unless the authors either supply an eigenvalue-independent criterion or explicitly reframe the potential as a post-hoc diagnostic.
major comments (3)
- [Section IIIB, Figs. 5-6] The non-axisymmetric potential U(x,ω,m) in Eq. (7) is evaluated at the complex eigenfrequency returned by the shooting solver, so the negative Re(U) well is not independent evidence for the onset of the mode. For any solution of Ψ''=UΨ, writing Ψ=R exp(iθ) and choosing a point where R=|Ψ| has a local maximum gives R'=0, R''<0, and hence Re(U)=R''/R−(θ')^2<0. Thus any localized mode, stable or unstable, produces such a well. The paper should either construct an ω-independent necessary condition for m≠0 analogous to the U(x,0)<0 calculation for m=0, or explicitly label the non-axisymmetric well as a diagnostic consequence of the mode rather than an onset mechanism.
- [Section II, paragraph after Eq. (9)] The statement that U(x,ω)<0 somewhere guarantees an unstable global mode is transferred from the axisymmetric marginal analysis to the frequency-dependent complex potential without proof. For m=0 the evaluation at ω²=0 yields a genuine, mode-independent necessary condition, but for m≠0 U is complex and depends on ω, so the cited Schrödinger-operator intuition (Ref. [27]) does not directly apply. The authors need to state the conditions under which the criterion holds for non-self-adjoint eigenvalue-dependent potentials, or restrict the claim to the axisymmetric case.
- [Section IIIB, resonance definitions] For an unstable mode with complex frequency ω=ω_r+iγ, Γ=(ω_r−mΩ+iγ)^2−ω_A^2 cannot vanish on the real x-axis: the imaginary part vanishes only at co-rotation, where the real part is −γ²−ω_A²<0. The paper nevertheless defines Alfvénic resonances as the real points satisfying Re(bar-ω)^2−ω_A^2=0 and uses them as boundaries of the potential well. This is an approximation whose accuracy should be stated (for example, a condition on γ/ω_A), otherwise the confinement picture is not rigorously tied to singularities of Eq. (7).
minor comments (4)
- [Eqs. (4) and (8)] The last term in Eqs. (4) and (8) is typeset as (Ωbar-ω+... )^2 Q(x) k^2 / Γ(x); to reproduce Eq. (9) it should be (Ωbar-ω+... )^2 k^2 / [Q(x) Γ(x)]. Please fix the typography.
- [Section IIIB, discussion of Fig. 6(a)] With bar-ω=ω−mΩ, the displayed identity U1=Γ''/(2Γ)=((mΩ')^2+m bar-ω Ω'')/Γ has the wrong sign on the second term; the correct expression is ((mΩ')^2−m bar-ω Ω'')/Γ. Please check.
- [Fig. 2 caption] The caption should state whether U is evaluated at the computed eigenvalue or at marginal ω²=0; the text distinguishes these cases for m=0, and the figure would be clearer if the choice were explicit.
- [Figs. 2 and 4] The references to 'left' and 'right' are ambiguous; use 'inner' and 'outer' boundary instead.
Circularity Check
Non-axisymmetric potential well is evaluated at the shooting-solver eigenvalue, so Re(U)<0 is automatic for any localized solution; the m=0 criterion is independent.
-
self definitional
[Section III B, Eqs. 6-7 and Figs. 5-6; cf. Eq. 9 and the m=0 criterion after Eq. 9.]
"We find that in the Alfv´enic propagating region where the mode is confined between the Alfv´enic resonance points, the potential (ℜ(U )) is negative and changes sign as it is evanescent outside this region."
The non-axisymmetric potential is not evaluated at a marginal or otherwise mode-independent frequency: the paper first solves Eq. 4 with a complex eigenvalue shooting solver, obtains the eigenvalue omega=omega_r+i gamma of the mode it wants to explain, and then constructs U(x,omega,m) from Eq. 7. Because Eq. 6 defines U=Psi''/Psi, any localized solution automatically yields Re(U)=R''/R-(theta')^2<0 at a maximum of R=|Psi| regardless of stability. Thus the negative well in Figs. 5-6 is guaranteed by the construction U=Psi''/Psi and cannot serve as an independent onset condition. The axisymmetric case is different: Eq. 9 is evaluated at omega=0 and gives the mode-independent condition omega_A^2<-omega_s^2, so that part is self-contained.
full rationale
The axisymmetric analysis is self-contained: Eq. 9 evaluated at the marginal frequency omega=0 yields omega_A^2<-omega_s^2 without using any eigenvalue, so that result is independent support. The non-axisymmetric interpretation, which is a central claim of the paper, is circular in the narrower sense: U(x,omega,m) is computed from the very complex eigenvalue found by the shooting solver, and Eq. 6 then makes Re(U)<0 automatic at a maximum of |Psi| for any localized solution. The negative potential well displayed in Figs. 5-6 is therefore a post-hoc diagnostic restatement of the mode, not an independent predictor of onset. Reliance on Ref. 25 for the ODE in Eq. 4 is a self-citation, but that equation is an algebraic consequence of Eqs. 1-3 and is not itself a circularity; I do not count it. Score 6 reflects the partial circularity of the central non-axisymmetric claim, while the m=0 MRI condition and the numerical mode searches retain independent content.
Assumptions & free parameters
free parameters (2)
- Rotation profile parameters (a1, a2, r1, r2) =
a1=0.84 (Tanh), a1=0.9, a2=0.6 (Tanh2), r1=1, r2=1.5 (mKep)
- Wave numbers m and k =
m=1, k=pi/4
assumptions (4)
- domain assumption The linearized ODE (Eq. 4) from Ref. 25 is valid.
- domain assumption Incompressibility and ideal MHD with uniform vertical magnetic field.
- standard math The normal mode ansatz, where perturbations scale as exp(i(m phi + k z - omega t)).
- domain assumption The criterion that an unstable global mode requires U(x,omega) < 0 somewhere in the domain, from Ref. 27.
Cite this review
Pith. "Pith review of A generalized effective potential for differentially rotating plasmas." pith.science (2026). https://pith.science/paper/XNSVWT3B
@misc{pith2026241207742,
author = {Pith},
title = {Pith review of: A generalized effective potential for differentially rotating plasmas},
year = {2026},
howpublished = {\url{https://pith.science/paper/XNSVWT3B}},
note = {Machine review of arXiv:2412.07742}
}
read the original abstract
Global stability of differentially rotating plasma is investigated using a generalized effective potential. We first, for a current-free system, obtain a general form of an effective potential in terms of the free energies of global curvature and gradients of rotation for non-axisymmetric disturbances. We then examine the stability of differentially rotating disks for several rotation profiles and present the associated effective potential for the onset of these instabilities in the MHD regime. In particular, results for global axisymmetric magnetorotational instability (MRI) as well as local and global non-axisymmetric modes are presented. The latter constitute two distinct non-axisymmetric modes, a high frequency local MRI and a global low-frequency non-axisymmetric mode (the magneto-curvature mode, introduced in Ebrahimi&Pharr, ApJ 2022), confined either between two Alfv\'enic resonances or an Alfv\'enic resonance and a boundary.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[27]
author F. Ebrahimi and author M. Pharr , journal The Astrophysical Journal volume 936 , pages 145 ( year 2022 )
work page 2022
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archive author booktitle chapter collaboration edition editor eid eprint howpublished institution isbn issn journal key month note number numpages organization pages publisher school series title type url volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
-
[3]
author E. Velikhov , journal Sov. Phys. JETP volume 36 , pages 995 ( year 1959 )
work page 1959
-
[4]
author S. Chandrasekhar , journal Proc. Natl. Acad. Sci. volume 46 , pages 253 ( year 1960 )
work page 1960
-
[5]
author S. A. Balbus and author J. F. Hawley , journal Astrophysical Journal, Part 1 (ISSN 0004-637X), vol. 376, July 20, 1991, p. 214-233. volume 376 , pages 214 ( year 1991 )
work page 1991
-
[6]
author F. Rincon , author G. Ogilvie , and author M. Proctor , journal Physical review letters volume 98 , pages 254502 ( year 2007 )
work page 2007
-
[7]
author F. Ebrahimi , author S. C. Prager , and author D. D. Schnack , journal The Astrophysical Journal volume 698 , pages 233 ( year 2009 ), ://doi.org/10.1088/0004-637x/698/1/233
Show all 32 references
-
[8]
Ebrahimi and author E
author F. Ebrahimi and author E. G. Blackman , journal Monthly Notices of the Royal Astronomical Society volume 459 , pages 1422 ( year 2016 ), ISSN issn 0035-8711 , https://academic.oup.com/mnras/article-pdf/459/2/1422/8040856/stw724.pdf , ://doi.org/10.1093/mnras/stw724
2016 doi
-
[9]
Bhat , author F
author P. Bhat , author F. Ebrahimi , and author E. G. Blackman , journal Monthly Notices of the Royal Astronomical Society volume 462 , pages 818 ( year 2016 ), ISSN issn 0035-8711 , https://academic.oup.com/mnras/article-pdf/462/1/818/18470481/stw1619.pdf , ://doi.org/10.109...
2016 doi
-
[10]
author D. R. Sisan , author N. Mujica , author W. A. Tillotson , author Y.-M. Huang , author W. Dorland , author A. B. Hassam , author T. M. Antonsen , and author D. P. Lathrop , journal Physical Review Letters volume 93 , pages 114502 ( year 2004 )
2004
-
[11]
Spence , author A
author E. Spence , author A. Roach , author E. Edlund , author P. Sloboda , and author H. Ji , journal Physics of Plasmas volume 19 ( year 2012 )
2012
-
[12]
author K. J. Caspary , author D. Choi , author F. Ebrahimi , author E. P. Gilson , author J. Goodman , and author H. Ji , journal Physical Review E volume 97 , pages 063110 ( year 2018 )
2018
-
[13]
Choi , author F
author D. Choi , author F. Ebrahimi , author K. J. Caspary , author E. P. Gilson , author J. Goodman , and author H. Ji , journal Physical Review E volume 100 , pages 033116 ( year 2019 )
2019
-
[14]
Mishra , author G
author A. Mishra , author G. Mamatsashvili , and author F. Stefani , journal Physical Review Fluids volume 7 , pages 064802 ( year 2022 )
2022
-
[15]
Wang , author E
author Y. Wang , author E. P. Gilson , author F. Ebrahimi , author J. Goodman , author K. J. Caspary , author H. W. Winarto , and author H. Ji , journal Nature communications volume 13 , pages 4679 ( year 2022 a )
2022
-
[16]
Wang , author E
author Y. Wang , author E. P. Gilson , author F. Ebrahimi , author J. Goodman , and author H. Ji , journal Physical review letters volume 129 , pages 115001 ( year 2022 b )
2022
-
[17]
Cairns , journal Journal of Fluid Mechanics volume 92 , pages 1 ( year 1979 )
author R. Cairns , journal Journal of Fluid Mechanics volume 92 , pages 1 ( year 1979 )
1979
-
[18]
Papaloizou and author J
author J. Papaloizou and author J. Pringle , journal Monthly Notices of the Royal Astronomical Society volume 208 , pages 721 ( year 1984 )
1984
-
[19]
Goldreich , author J
author P. Goldreich , author J. Goodman , and author R. Narayan , journal Monthly Notices of the Royal Astronomical Society volume 221 , pages 339 ( year 1986 )
1986
-
[20]
Goodman , author R
author J. Goodman , author R. Narayan , and author P. Goldreich , journal Monthly Notices of the Royal Astronomical Society volume 225 , pages 695 ( year 1987 )
1987
-
[21]
Lovelace , author H
author R. Lovelace , author H. Li , author S. Colgate , and author A. Nelson , journal The Astrophysical Journal volume 513 , pages 805 ( year 1999 )
1999
-
[22]
Glatzel , journal Monthly Notices of the Royal Astronomical Society volume 228 , pages 77 ( year 1987 )
author W. Glatzel , journal Monthly Notices of the Royal Astronomical Society volume 228 , pages 77 ( year 1987 )
1987
-
[23]
Ono , author T
author T. Ono , author T. Muto , author T. Takeuchi , and author H. Nomura , journal The Astrophysical Journal volume 823 , pages 84 ( year 2016 )
2016
-
[24]
Matsumoto and author T
author R. Matsumoto and author T. Tajima , type Tech. Rep. , institution Univ. of Texas, Austin, TX (United States). Institute for Fusion Studies ( year 1995 )
1995
-
[25]
author G. I. Ogilvie and author J. E. Pringle , journal Monthly Notices of the Royal Astronomical Society volume 279 , pages 152 ( year 1996 ), ISSN issn 0035-8711 , https://academic.oup.com/mnras/article-pdf/279/1/152/4170650/279-1-152.pdf , ://doi.org/10.1093/mnras/279.1.152
1996 doi
-
[26]
Goedbloed and author R
author H. Goedbloed and author R. Keppens , journal arXiv preprint arXiv:2201.11551 ( year 2022 )
2022 arXiv
-
[28]
Chandrasekhar , title Ch IX: The Stability of Couette Flow in Hydromagnetics ( publisher Dover Publ
author S. Chandrasekhar , title Ch IX: The Stability of Couette Flow in Hydromagnetics ( publisher Dover Publ. , year 2006 ), p. pages 382–427
2006
-
[29]
Pino and author S
author J. Pino and author S. Mahajan , journal The Astrophysical Journal volume 678 , pages 1223 ( year 2008 )
2008
-
[30]
Sovinec , author A
author C. Sovinec , author A. Glasser , author T. Gianakon , author D. Barnes , author R. Nebel , author S. Kruger , author D. Schnack , author S. Plimpton , author A. Tarditi , and author M. Chu , journal Journal of Computational Physics volume 195 , pages 355 ( year 2004 ), ...
2004
-
[31]
Curry and author R
author C. Curry and author R. E. Pudritz , journal Monthly Notices of the Royal Astronomical Society volume 281 , pages 119 ( year 1996 )
1996
-
[32]
Wang , author F
author Y. Wang , author F. Ebrahimi , author H. Lu , author J. Goodman , author E. P. Gilson , and author H. Ji , journal arXiv preprint arXiv:2411.02361 ( year 2024 )
2024 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.