{"id":"36c7fbbd-85da-4fad-ab94-6fbdb5aff423","arxiv_id":"2412.07742","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A generalized effective potential for non-axisymmetric MHD perturbations in rotating plasmas is derived, and its negative regions between Alfvénic resonances are shown to confine both high-frequency local MRI and global low-frequency magneto-curvature modes.","lead":"This paper derives a generalized effective potential for perturbations in differentially rotating magnetized plasma disks, and uses it to show that two families of non-axisymmetric instabilities live in a potential well between Alfvénic resonance points. The work offers a unified diagnostic for disk stability and could inform models of accretion and turbulent transport in astrophysical and laboratory plasmas.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The non-axisymmetric 'potential well' is evaluated at the unstable eigenvalue, so Re(U)<0 is automatically implied by any localized solution; it cannot serve as an independent onset criterion.","rationale":"The derivation of Eq. 7 is exact and the m=0 reduction to the MRI bound is genuine independent support, so I am not objecting to the algebraic core. The reader's weakest_assumption already identifies the post-hoc use of omega; I would sharpen it: for a complex solution, the negative Re(U) at the mode location is not merely post-hoc but mathematically forced. At a maximum of R=|Psi|, Re(U)=R''/R-(theta')^2<0, so any localized solution, whether growing, decaying, or neutral, shows a negative well. Consequently the non-axisymmetric 'onset' claim in the abstract and summary is currently unfalsifiable, while the axisymmetric criterion U(x,0)<0 remains substantive. The paper should therefore be accepted only conditionally: either derive a mode-independent condition for m!=0, for example by evaluating U at real frequencies consistent with the resonance condition and comparing the predicted unstable window with the shooting eigenvalues, or explicitly present the well as a diagnostic of the already-computed mode rather than as an onset criterion. The conjugate-mode check in concrete_test would settle whether the well is instability-specific. I keep the reader's CONDITIONAL verdict and note that no code or data are provided for the numerical results.","tokens_in":7354,"tokens_out":9277,"duration_ms":91656,"concrete_test":"For the Tanh2, m=1, k=k1, VA/V0=0.2 case, rerun the shooting solver with the complex-conjugate frequency omega* = omega_r - i gamma, i.e. the damped counterpart of the reported growing mode. Since all coefficients entering Eq. 7 are real analytic functions of omega, this damped eigenmode gives exactly the same Re(U) as the growing mode, including the negative well. If the same negative well appears for a decaying mode, then Re(U)<0 is not an instability criterion, and an independent, eigenvalue-free onset condition for m!=0 must be supplied.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is transferring the U(x,omega)<0 instability criterion from the m=0 case to non-axisymmetric modes. For m=0 the paper evaluates U at the marginal frequency omega=0 (Eq. 9) and obtains a genuine, mode-independent necessary condition, omega_A^2 < -omega_s^2. For m!=0, however, the U displayed in Figs. 5-6 is constructed from the complex eigenvalue omega=omega_r+i gamma returned by the shooting solver. This makes the negative Re(U) well a mathematical consequence of the solution rather than an independent predictor: for any localized mode satisfying Psi''=U Psi, writing Psi=R exp(i theta) and taking a point where R=|Psi| has a maximum gives R'=0 and R''<0, so Re(U)=Re(Psi''/Psi)=R''/R-(theta')^2<0. Thus Re(U)<0 somewhere is guaranteed by localization alone; stable and unstable modes alike would show it. The abstract, Section IIIB, and the summary treat the negative well as the onset mechanism for non-axisymmetric instability, but as it stands this is post-hoc and unfalsifiable. What is missing is an eigenvalue-independent condition for m!=0, analogous to the axisymmetric U(x,0)<0 calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7579,"tokens_out":10498,"duration_ms":95072,"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":[{"comment":"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":"Section IIIB, Figs. 5-6"},{"comment":"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":"Section II, paragraph after Eq. (9)"},{"comment":"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).","section":"Section IIIB, resonance definitions"}],"minor_comments":[{"comment":"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":"Eqs. (4) and (8)"},{"comment":"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.","section":"Section IIIB, discussion of Fig. 6(a)"},{"comment":"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.","section":"Fig. 2 caption"},{"comment":"The references to 'left' and 'right' are ambiguous; use 'inner' and 'outer' boundary instead.","section":"Figs. 2 and 4"}],"recommendation":"major_revision","confidential_remarks":"The central derivation and the m=0 criterion are solid, and the numerical results appear credible, but the non-axisymmetric interpretation as currently worded is circular in the sense that the potential is built from the mode it explains. I would support publication after the authors reframe the non-axisymmetric potential as a diagnostic or supply an independent criterion; I do not see this as grounds for rejection, but the major comments should be addressed in full."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is the transformation of the linearized ideal MHD equations into a Schrödinger-like form with a generalized effective potential U(x,ω,m), Eqs. 6–8. That algebra is clean, and I checked the m=0 limit: it reproduces the standard axisymmetric MRI condition ω_A² < −ω_s². The authors also show that the second derivative of the rotation profile enters the potential, and they benchmark their shooting solver against NIMROD, which is honest work. The two non-axisymmetric branches were already found in their earlier paper, so the genuinely new pieces are the potential form and the curvature free-energy term.\n\nBut the non-axisymmetric part has a load-bearing problem, and the stress-test note is right. For m≠0, the U shown in Figs. 5 and 6 is evaluated at the complex eigenvalue ω returned by the solver. That makes the negative potential well a mathematical consequence of having a localized solution, not an independent predictor. At any point where |Ψ| has a maximum, Re(U)=R''/R−(θ')² is automatically negative for any localized mode, stable or unstable. So Re(U)<0 somewhere is guaranteed by localization alone; it cannot tell you which modes go unstable. The paper's abstract, Section IIIB, and summary treat the well as the onset mechanism for non-axisymmetric MRI and magneto-curvature modes, but that claim is unfalsifiable as presented. What is missing is an eigenvalue-independent condition for m≠0, analogous to the U(x,0)<0 calculation that works so well for m=0.\n\nThere are also smaller issues. No code or data are provided, so the NIMROD comparison is hard to verify independently. The claim that the curvature term is the dominant free energy is plausible but again is only established after solving for the mode, so it is more diagnostic than explanatory. The paper would be strengthened if the authors either explicitly framed the non-axisymmetric potential as a diagnostic tool or found a frequency-independent stability criterion.\n\nWho gets value from this? People working on global stability of rotating MHD plasmas, especially those interested in the mathematical structure of the perturbation equations. The m=0 part is solid and citable. The non-axisymmetric interpretive claims need reworking before I would rely on them.\n\nRecommendation: yes, it deserves peer review. The algebraic core is sound and the paper should be published after the authors clarify the status of the non-axisymmetric potential. A competent referee can sort this out, but only if they push on the post-hoc issue.","headline":"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.","tokens_in":8110,"tokens_out":1655,"would_cite":true,"duration_ms":17479,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.30.Cv","97.10.Gz"],"model":"deepseek-v4-flash","headline":"A generalized effective potential, built from rotation gradients and curvature, predicts the onset of both axisymmetric and non-axisymmetric instabilities in differentially rotating plasma.","keywords":["magnetorotational instability","effective potential","non-axisymmetric modes","Alfvénic resonance","differential rotation","global stability","accretion disks","ideal MHD"],"falsifier":"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.","tokens_in":7127,"feed_emoji":"🌀","tokens_out":13318,"duration_ms":99690,"temperature":0.7,"pith_summary":"This paper derives a single effective potential, $U(x,\\omega,m) = A''/A - C/A^2$, that governs the global linear stability of a current-free, differentially rotating plasma. Cast as a Schrödinger-type equation $\\Psi'' - U\\Psi = 0$, the stability problem becomes a question of where the potential turns negative: instabilities live in potential wells confined between two Alfvénic resonances, or between a resonance and a boundary. In the axisymmetric limit the potential collapses to the familiar magnetorotational instability (MRI) criterion $\\omega_A^2 < -\\omega_s^2$, and for non-axisymmetric modes it separates two distinct branches, a high-frequency localized MRI and a low-frequency global magneto-curvature mode. The payoff is a single diagnostic function of radius, frequency, and azimuthal mode number that organizes the onset of very different instabilities in accretion disks and rotating laboratory plasmas.","feed_headline":"One potential predicts where plasma disks go unstable","feed_subtitle":"It recovers the standard MRI rule and maps the wells that confine two non-axisymmetric modes.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"The companion paper that supplies the linearized radial equation, the complex-frequency shooting solver, and the original classification of the two non-axisymmetric mode branches.","marker":"[25]"},{"why":"Source of the criterion that a negative effective potential somewhere in the domain guarantees an unstable global mode, the interpretive backbone of the stability analysis.","marker":"[27]"},{"why":"The local axisymmetric MRI dispersion relation whose instability condition the potential must reproduce in the $m=0$ limit.","marker":"[3]"},{"why":"Showed localized non-axisymmetric MRI structures confined between Alfvénic resonant points in Cartesian geometry, grounding the resonance-confinement picture.","marker":"[22]"},{"why":"Extended Alfvénic-resonance confinement of non-axisymmetric modes to cylindrical shear flows, the direct predecessor of the potential-well picture.","marker":"[23]"},{"why":"The initial-value code against which the shooting solver's eigenvalues are benchmarked.","marker":"[28]"},{"why":"The hydrodynamical analogue in which waves are amplified in a Rossby region and damped outside it, motivating the confinement interpretation.","marker":"[21]"}],"fun_headline_variants":["A single potential maps disk stability","Unified potential predicts MRI and new modes","Potential wells confine non-axisymmetric instabilities","Rotation curvature drives novel plasma modes","Generalized potential links MRI to curvature modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["A single potential maps disk stability","Unified potential predicts MRI and new modes","Potential wells confine non-axisymmetric instabilities","Rotation curvature drives novel plasma modes","Generalized potential links MRI to curvature modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000263,"raw_usage":{"total_tokens":1620,"prompt_tokens":986,"completion_tokens":634,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":572}},"tokens_in":602,"tokens_out":634,"duration_ms":6606,"temperature":1.0,"reasoning_tokens":572,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:33:21.247404+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Ebrahimi and author M","cited_arxiv_id":null,"evidence_quote":"Source of the criterion that a negative effective potential somewhere in the domain guarantees an unstable global mode, the interpretive backbone of the stability analysis."},{"cited_title":"Velikhov , journal Sov","cited_arxiv_id":null,"evidence_quote":"The local axisymmetric MRI dispersion relation whose instability condition the potential must reproduce in the $m=0$ limit."},{"cited_title":"Glatzel , journal Monthly Notices of the Royal Astronomical Society volume 228 , pages 77 ( year 1987 )","cited_arxiv_id":null,"evidence_quote":"Showed localized non-axisymmetric MRI structures confined between Alfvénic resonant points in Cartesian geometry, grounding the resonance-confinement picture."},{"cited_title":"Ono , author T","cited_arxiv_id":null,"evidence_quote":"Extended Alfvénic-resonance confinement of non-axisymmetric modes to cylindrical shear flows, the direct predecessor of the potential-well picture."},{"cited_title":"Chandrasekhar , title Ch IX: The Stability of Couette Flow in Hydromagnetics ( publisher Dover Publ","cited_arxiv_id":null,"evidence_quote":"The initial-value code against which the shooting solver's eigenvalues are benchmarked."},{"cited_title":"Lovelace , author H","cited_arxiv_id":null,"evidence_quote":"The hydrodynamical analogue in which waves are amplified in a Rossby region and damped outside it, motivating the confinement interpretation."}],"review_version":1}