{"id":"221a78c9-3af5-46a7-98d7-388739f2a60f","arxiv_id":"1908.08821","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Explicit sufficient stability criteria and a Hall MHD energy principle for XMHD equilibria with flow, derived from energy-Casimir, dynamically accessible, and Lagrangian methods.","lead":"This paper derives new sufficient stability conditions for extended and Hall magnetohydrodynamic plasma equilibria with flow, using Hamiltonian and Lagrangian energy methods. The conditions are conservative checks: failing them does not prove instability, but satisfying them certifies linear stability within the model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The DA second variation (63) is asserted without derivation and is not obviously equivalent to (115); the generic DA stability criteria, including (121), rest on this unchecked algebra.","rationale":"The reader's conditional verdict is reasonable and I would keep it: the paper contains several independent criteria and the EC and Lagrangian sections are plausible enough to merit conditional acceptance. However, the single most load-bearing point is not the inherited bracket/Casimirs identified as the reader's weakest assumption. Those objects are published in the group's earlier papers, and the first-variation equations in the present paper are consistent with them. The decisive step is the transition from (115) to (63) in Section III: without (63) there is no explicit DA stability criterion for generic equilibria, and the appendix's simple toroidal-flow criterion (121) relies on it. This transition is asserted without proof, and it is not pointwise or even superficially obvious; the missing rho^-1 in Eq. (37) shows that algebraic slips of exactly this kind occur elsewhere in the manuscript. The concrete test I propose would settle the issue by checking the integral identity and the claimed MHD/HMHD limits. I marked disagreement with the reader's weakest_assumption not because the Casimir question is irrelevant, but because the DA algebra is more immediately load-bearing for the paper's central claims as stated in the abstract and strongest_claim.","tokens_in":26657,"tokens_out":13540,"duration_ms":134295,"concrete_test":"Independently derive Eq. (63) from Eq. (115) using the definitions (52)-(54), (60)-(62) and the XMHD equilibrium equations, ideally with a computer-algebra system (SymPy or Mathematica), and verify the integral identity exactly. If it holds, set d_e = 0 in both (115) and (63) and compare term by term with the barotropic HMHD expression of Hirota et al. [22], and set d_i = 0 to compare with the Frieman-Rotenberg/MHD limit. If the identities fail or the HMHD/MHD limits disagree, then the DA stability criteria, including Eq. (121), are not established. As a secondary check, re-derive Eqs. (36)-(37) by explicit variation of delta^2 H_C in (15); Eq. (37) currently lacks the rho^-1 factor required by the Euler-Lagrange equation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III presents two expressions for the dynamically accessible second variation of H: the intermediate form (Appendix A, Eq. (115)) and the simplified form (Eq. (63)), obtained, per the text, after 'some tedious but also straightforward manipulations.' No derivation is given. The two forms are not manifestly equivalent: the velocity square in (115) contains -rho^-1 v div(rho zeta), while (63) contains zeta.grad v - v.grad zeta; these differ by density- and divergence-dependent terms, so the claimed integral identity requires nontrivial cancellations using the equilibrium equations. This is the central object of Section III: the only explicit generic stability criterion in the paper, Eq. (121) (c_s^2 - v_phi^2 - d_e^2 |J|^2/rho^2 > 0 for axisymmetric toroidal flow with flux-surface-tangent perturbations), is derived from (63). The asserted d_e -> 0 limit reproducing the HMHD result of Hirota et al. [22] is also not checked term by term. If (63) has a sign error or a missing term, the DA stability claims do not follow and the error would be invisible in the paper's presentation. The surrounding algebra is demonstrably error-prone: Eq. (37), the minimizer for delta v_perp obtained from (15), omits the rho^-1 factor that the Euler-Lagrange equation and the later substitution into (39) both require. This supports treating (63) as unverified rather than as a routine simplification.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops formal stability criteria for extended magnetohrodynamic (XMHD) equilibria using three Hamiltonian-based methods. In the energy-Casimir (EC) approach, sufficient stability conditions are derived for axisymmetric XMHD equilibria with toroidal flow (Eqs. (26)-(27)) and for Hall MHD (HMHD) equilibria (Eqs. (28)-(30)), with a numerical HMHD example. For restricted perturbations, a conditional stability criterion (47) is obtained by partially minimizing the second variation. In the dynamically accessible (DA) approach, a second-order variation of the Hamiltonian is computed for generic equilibria, and a sufficient condition (121) is given for axisymmetric equilibria with toroidal flow under flux-surface-tangent perturbations. In the Lagrangian part, the second-order Lagrangian and Hamiltonian of the linearized two-fluid dynamics are constructed; upon setting the electron mass to zero, a HMHD energy principle (114) is derived including the electron pressure contribution, and the perturbed induction equation is obtained from a consistency condition. The paper emphasizes that all stability criteria are sufficient, not necessary.","tokens_in":26958,"tokens_out":12122,"duration_ms":121498,"significance":"If the derivations are correct, the paper provides explicit Lyapunov-function-based sufficient stability conditions for XMHD and HMHD equilibria with flow, going beyond static MHD. The EC criteria (26)-(30) are explicit and amenable to numerical evaluation, as demonstrated in Figs. 1 and 2. The Lagrangian/Hamiltonian formulation yields a HMHD energy principle with electron pressure and a derivation of the perturbed induction equation from a dynamical constraint, which are valuable new results. The DA section is potentially useful for generic equilibria because it avoids the known Casimir MHD-limit difficulties. However, the strength of the paper rests on several algebraic steps that are asserted rather than shown, and at least one of those steps contains a visible error. The paper’s central claims are plausible, but the load-bearing algebra must be fully verified before the results can be accepted as rigorous.","major_comments":[{"comment":"The minimization of the functional (15) with respect to δv⊥ gives δv⊥ = -ρ^{-1} v⊥ δρ + ρ^{-1}(γ∇δF + μ∇δG) × ∇φ, not the expression printed in Eq. (37), which omits the factor ρ^{-1} multiplying the second term. This is not a harmless typo: after substitution into (38)-(39), the kinetic contribution changes, so the coefficients in (39) and hence the sufficient conditions (44)-(45) and the criterion (47) are not justified as written. Please correct the minimizer and re-derive the subsequent inequalities, or explicitly state a different convention for the variation that makes Eq. (37) correct.","section":"Sec. II.D, Eq. (37)"},{"comment":"The equivalence of the intermediate form (115) and the simplified form (63) is asserted after 'some tedious but also straightforward manipulations' without any derivation. The two forms are not manifestly identical: for example, (115) contains the term -ρ^{-1} v ∇·(ρζ) inside the squared velocity, while (63) contains ζ·∇v - v·∇ζ, so the required cancellations involve the equilibrium equations and are nontrivial. Since Eq. (63) underlies the DA stability discussion and the claimed HMHD and MHD limits, this is a load-bearing algebraic assertion. Please provide a complete derivation, either in the paper or in a supplementary file, and identify which equilibrium equations and integration-by-parts identities are used at each step.","section":"Sec. III, Eq. (63) and Appendix A, Eq. (115)"},{"comment":"The reduction from (115) to (121) for flux-surface-tangent perturbations is described only verbally ('A rigorous proof can be carried out...'), with no term-by-term computation. The criterion c_s^2 - v_φ^2 - d_e^2 |J|^2/ρ^2 > 0 is one of the main DA results, so the fate of each term in (115), including the terms containing ∇g0, ζ×ω, η×B*, and all cross terms involving δB_da, should be displayed explicitly. Without this demonstration the reader cannot exclude hidden indefinite terms that would invalidate the sufficient condition.","section":"Appendix A, Eq. (121)"}],"minor_comments":[{"comment":"The title contains a typographical error: 's tability' should read 'stability'.","section":"Title"},{"comment":"For reproducibility of the numerical examples, the values of the free coefficients g_i, m_i, n_i in Eq. (35) and the computational parameters should be listed, since the stability diagrams in Figs. 1 and 2 depend on them.","section":"Sec. II.C.2"},{"comment":"The condition ⟨(|k_x|^2 - C^{-1})(δx)^2⟩ ≥ 0 is a consequence of the Poincaré inequality and is not an additional restriction on the perturbations; this should be stated explicitly to avoid confusion about its role in the criterion.","section":"Sec. II.D, Eq. (47)"},{"comment":"The notation ζ and η is reused for the Eulerianized displacements after dropping the tildes; a sentence reminding the reader of this change of notation would improve clarity.","section":"Sec. IV, after Eq. (84)"},{"comment":"Reference [57] is a technical report that may not be readily accessible; please indicate whether a journal version or a stable online source exists.","section":"Ref. [57]"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the conceptual framework is sound, but the DA section and the conditional-stability section contain algebraic steps that are not verified, and Eq. (37) has a concrete error. I recommend a major revision during which the authors supply the missing derivations; if the derivations check out, the paper could be acceptable. The EC criteria for the special equilibria and the Lagrangian/HMHD energy principle appear to be the strongest parts of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I gave this a careful read, and here is my take. The paper is a substantial extension of formal stability analysis to extended MHD: it writes down explicit sufficient energy-Casimir criteria for axisymmetric XMHD and Hall MHD equilibria with toroidal flow, derives the dynamically accessible second variation for XMHD, and constructs the Hamiltonian of the linearized two-fluid dynamics in mixed Eulerian-Lagrangian variables, from which a Hall MHD energy principle with electron pressure follows. That last piece is new and looks like the most durable part of the paper. The authors are honest about the conditions: everything is sufficient, often only for symmetric equilibria or constrained perturbations, and they say so. What is new: the explicit inequalities (26)-(27), (28)-(30), the conditional criterion (47), the HMHD energy principle (114), and the form B1 = curl[(zeta - d_i eta) x B]. The EC part is careful and mostly self-contained. The comparison with MHD limits is thoughtful. Now the soft spots. Section III, the dynamically accessible analysis, carries the weight of the paper's most general claim, and that is exactly where the presentation fails. The expression for delta^2 H_da is quoted after 'tedious but also straightforward manipulations' with no derivation. The intermediate form (115) and the final form (63) are not manifestly equal, and the route between them is exactly where a sign error or missing term would be invisible. More concretely, Eq. (37), the minimizer for delta v_perp, is missing a factor of rho^{-1} that both the Euler-Lagrange equation and the subsequent substitution require. It may be a typo, but it sits in a load-bearing derivation, and it makes Eq. (39) and the derived c_phi, c_xi conditions formally unreliable as printed. The claimed d_e -> 0 limit reproducing Hirota et al. is also not checked term by term. This is the one part of the paper I would not trust without checking everything myself. The numerical demonstration is a sanity check, not evidence: no code, no full parameter list, and it only checks sufficient conditions, so it cannot catch an algebraic error in the criteria. The reliance on the noncanonical bracket and axisymmetric Casimirs from earlier papers is a real inheritance, but I do not hold that against them: the bracket is from a published paper and the Casimirs are theirs. That is normal for this literature. Who is this for? Plasma theorists working on formal stability of two-fluid models. It deserves a serious referee, but that referee should be asked to verify Section III line by line. My recommendation: send to peer review, with a request for the derivation of (63) and a fix for Eq. (37), or a withdrawal of the generic DA criterion until it is shown.","headline":"Extends energy-Casimir, dynamically accessible, and Lagrangian stability tools to XMHD and gets a new HMHD energy principle, but the load-bearing DA algebra is asserted, and Eq. (37) has a missing rho^{-1} as printed.","tokens_in":27509,"tokens_out":4047,"would_cite":false,"duration_ms":39344,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives explicit sufficient stability criteria for axisymmetric XMHD and Hall MHD equilibria with toroidal flow, a wavenumber-limited criterion for incompressible perturbations, and a Lagrangian HMHD energy principle that…","keywords":["extended magnetohydrodynamics","Hall MHD","energy-Casimir stability","dynamically accessible variations","Lagrangian stability","noncanonical Hamiltonian structure","axisymmetric plasma equilibria","energy principle"],"falsifier":"Take a numerically constructed XMHD equilibrium satisfying (26)–(27), or an HMHD equilibrium satisfying (28)–(30), and solve the linearized equations with a small axisymmetric perturbation; any growing mode would disprove the claimed Lyapunov sufficiency. A cheaper equivalent is to scan the linear eigenmode spectrum of the operators in (15) and (38) for a positive growth rate in the parameter region where the paper's conditions hold.","tokens_in":26455,"feed_emoji":"🧲","tokens_out":10767,"duration_ms":104148,"temperature":0.7,"pith_summary":"The paper aims to give explicit sufficient stability criteria for stationary plasma equilibria of extended magnetohydrodynamics (XMHD) and Hall MHD (HMHD), including equilibria with macroscopic flow. Three Hamiltonian-based routes are used: the energy-Casimir variational principle, dynamically accessible variations, and a mixed Eulerian–Lagrangian action. The central results are the inequalities (26)–(27) for axisymmetric XMHD equilibria with toroidal flow, (28)–(30) for HMHD, a wavenumber-limited criterion for incompressible perturbations, and the Lagrangian/Hamiltonian stability conditions (93) and (114). If the analysis is right, each criterion is a genuine Lyapunov-functional condition for the linearized dynamics, and the HMHD energy principle including the electron pressure contribution is a new result.","feed_headline":"Criteria now certify stability of flowing two-fluid plasma states","feed_subtitle":"Energy-Casimir, dynamically accessible, and Lagrangian methods yield sufficient conditions including electron pressure.","key_machinery":"The central objects are the noncanonical Poisson bracket (2), the axisymmetric Casimir invariants (9)–(12), and the generalized magnetic field $B^* = B + d_e^2 \\nabla\\times(\\nabla\\times B/\\rho)$ that carries electron inertia. These convert the Hamiltonian into an energy-Casimir functional whose first variation is the equilibrium condition and whose second variation is the candidate Lyapunov functional; the explicit inequalities come from requiring the matrix $\\mathsf A$ in (17) and the surviving squared terms to be positive definite. The dynamically accessible generator $W$ produces variations tangent to the level sets of the Casimirs, giving a second variation for generic equilibria. In the Lagrangian route, the second-order action $L_2$ and its Legendre transform $H_2$ supply the energy principles, and for HMHD the momentum $\\pi_\\eta$ becomes a constraint whose consistency yields the perturbed induction equation.","core_discovery":"Within the noncanonical Hamiltonian description of XMHD, the paper constructs energy-Casimir functionals $H-\\sum_i C_i$ whose first variation vanishes on the equilibria derived in the authors' earlier work and whose second variation, when positive definite, serves as a Lyapunov functional for the linearized dynamics. For purely toroidal flow the indefinite terms in $\\delta^2 H_C$ are eliminated, leaving algebraic conditions on the free functions: (26)–(27) for XMHD and (28)–(30) for HMHD, with numerical satisfaction shown on diverted tokamak-like equilibria. For incompressible perturbations, partial minimization over $\\delta B_\\varphi$ and $\\delta v_\\perp$ produces the condition $c_\\phi>0$, $c_\\xi>0$, together with the wavenumber restriction $|k_x|<k_+$ and a lower bound on $|k_x|$ imposed by the geometry of the domain. The dynamically accessible route, generated by $W=\\int d^3x(g_0\\rho + g_1\\cdot v + g_2\\cdot B^*)$, yields a second variation $\\delta^2 H_{da}$ valid for generic equilibria. The mixed Eulerian–Lagrangian expansion of the quasineutral two-fluid action gives the linearized Hamiltonian $H_2$ with the sufficient condition $-\\int W\\,d^3x \\ge 0$; taking the electron mass to zero produces the HMHD energy principle $-\\int W_{\\rm hmhd}\\,d^3x \\ge 0$, and the perturbed induction equation $B_1=\\nabla\\times[(\\zeta-d_i\\eta)\\times B_0]$ follows from the consistency condition (102) on the canonical momentum $\\pi_\\eta$.","pith_inferences":["A natural extension the paper does not pursue is to test (26)–(30) against direct linear eigenmode solvers; such a test would show how far the sufficient conditions are from necessary ones.","Because the DA second variation avoids the singular Casimir/MHD limit noted in the paper, it may be the most practical of the three tools for numerical stability scans of three-dimensional equilibria.","The mixed Eulerian–Lagrangian derivation keeps the full two-fluid pressure structure before the zero-electron-mass limit, so the resulting energy principle could likely be adapted to non-barotropic closures or to equilibria with electrostatic potentials."],"forward_implications":["Axisymmetric XMHD equilibria with purely toroidal rotation are formally stable whenever the free functions $M,N$ are concave and the electron-inertia-modified toroidal speed stays below the sound speed, conditions (26)–(27).","Hall MHD equilibria with toroidal flow satisfy an analogous explicit criterion (28)–(30); the numerical examples indicate the extra condition (30) is the restrictive one at high $\\beta$ and that increasing the Hall parameter $d_i$ can enlarge the stable region.","For incompressible perturbations, the criteria $c_\\phi>0$, $c_\\xi>0$ with the wavenumber bound $|k_x|<k_+$ provide sufficient stability for equilibria that need not be of the purely toroidal class.","The dynamically accessible second variation supplies a stability test for generic three-dimensional equilibria, and its MHD limit recovers the standard flowing-MHD potential energy without the singularities that plague the Casimir/MHD limit.","The Lagrangian energy principle for the quasineutral two-fluid model, and its HMHD reduction, give a sufficient condition $-\\int W\\,d^3x \\ge 0$ that now includes the electron pressure contribution."],"supporting_citations":[{"why":"Supplies the noncanonical Poisson bracket (2) on which all three stability analyses rest.","marker":"[29]"},{"why":"Supplies the axisymmetric equilibrium equations and the Casimir invariants (9)–(12) used by the energy-Casimir method.","marker":"[39]"},{"why":"Establishes the formal-stability framework: first variation vanishing at equilibrium and second variation serving as a Lyapunov functional.","marker":"[20]"},{"why":"Provides the HMHD dynamically accessible second variation and the perturbed induction equation that this paper generalizes.","marker":"[22]"},{"why":"Supplies the MHD energy-Casimir analogue and the partial-minimization technique used to remove indefinite terms.","marker":"[36]"},{"why":"Supplies the mixed Eulerian–Lagrangian action formalism used to derive $L_2$.","marker":"[35]"},{"why":"Gives the fully Lagrangian MHD stability framework that the mixed Eulerian–Lagrangian derivation extends.","marker":"[56]"},{"why":"The standard MHD energy principle that motivates the sufficient-stability logic and that the HMHD principle extends.","marker":"[1]"}],"fun_headline_variants":["Stability criteria found for flowing two-fluid plasmas","Energy-Casimir method certifies XMHD equilibrium stability","Sufficient conditions for MHD equilibria with flow","Hall MHD energy principle from Lagrangian stability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire analysis inherits, without re-derivation, the noncanonical Poisson bracket (2) and the axisymmetric Casimirs (9)–(12) from Refs. [29] and [39]; if that bracket or those invariants are incomplete, the stated inequalities do not certify stability.","fun_headline_variants_meta":{"raw":{"variants":["Stability criteria found for flowing two-fluid plasmas","Energy-Casimir method certifies XMHD equilibrium stability","Sufficient conditions for MHD equilibria with flow","Hall MHD energy principle from Lagrangian stability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1517,"prompt_tokens":1091,"completion_tokens":426,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":707,"completion_tokens_details":{"reasoning_tokens":363}},"tokens_in":707,"tokens_out":426,"duration_ms":5135,"temperature":1.0,"reasoning_tokens":363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:29:16.436366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a numerically constructed XMHD equilibrium satisfying (26)–(27), or an HMHD equilibrium satisfying (28)–(30), and solve the linearized equations with a small axisymmetric perturbation; any growing mode would disprove the claimed Lyapunov sufficiency. A cheaper equivalent is to scan the linear eigenmode spectrum of the operators in (15) and (38) for a positive growth rate in the parameter region where the paper's conditions hold.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the noncanonical Poisson bracket (2) on which all three stability analyses rest."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the axisymmetric equilibrium equations and the Casimir invariants (9)–(12) used by the energy-Casimir method."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the formal-stability framework: first variation vanishing at equilibrium and second variation serving as a Lyapunov functional."},{"cited_title":"Hirota, Z","cited_arxiv_id":null,"evidence_quote":"Provides the HMHD dynamically accessible second variation and the perturbed induction equation that this paper generalizes."},{"cited_title":"Andreussi, P","cited_arxiv_id":null,"evidence_quote":"Supplies the MHD energy-Casimir analogue and the partial-minimization technique used to remove indefinite terms."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the mixed Eulerian–Lagrangian action formalism used to derive $L_2$."},{"cited_title":"Lingam, G","cited_arxiv_id":null,"evidence_quote":"Gives the fully Lagrangian MHD stability framework that the mixed Eulerian–Lagrangian derivation extends."},{"cited_title":"F ′ = G′ = 0, it is clear that Q > 0 implies δ2HC > 0","cited_arxiv_id":null,"evidence_quote":"The standard MHD energy principle that motivates the sufficient-stability logic and that the HMHD principle extends."}],"review_version":1}