{"id":"2ef48ca2-e625-4391-8811-7cad435f574e","arxiv_id":"2505.10397","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Linear stability and causality of the second-order magnetohydrodynamics from Ref. [64] are verified for any magnetic field in a locally neutral two-component massless plasma.","lead":"This paper tests whether a recently proposed relativistic magnetohydrodynamics theory for a two-component ultrarelativistic plasma is causal and stable. It concludes the theory is linearly stable for any magnetic field, and shows that a common astrophysical simplification only works for very strong fields.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'always causal and stable' claim is proven only for perturbations strictly parallel or perpendicular to the magnetic field; oblique directions are never analyzed, so the universal conclusion in Sec. VI is not supported by the presented evidence.","rationale":"I read the paper's central claim as unconditional: the second-order MHD from Ref. [64] is linearly causal and stable for any B, any cross-sections, and hence any perturbation direction. The analysis, however, only treats propagation directions along and across the magnetic field. The general equations (42)-(44) are not solved for oblique wave vectors, and no determinant is presented for that case. The reader's verdict flags the intermediate-k numerical gap, which is a related but distinct limitation; I see the unanalyzed oblique sector as the more direct reason the universal statement is not established. This does not mean the claim is false—the special-direction results are credible evidence—but it means the conclusion overreaches. The authors can either close the gap with a Routh-Hurwitz analysis over θ and k^2 or soften the claim. Accordingly, the reader's CONDITIONAL verdict remains appropriate, and I do not recommend changing it.","tokens_in":20775,"tokens_out":21664,"duration_ms":203076,"concrete_test":"Set κb = k cos θ and κ⊥ = k sin θ in the general linearized system (42)-(44), form the characteristic determinant, and apply Routh-Hurwitz with a computer algebra quantifier-elimination pass to test whether every root obeys Im ω ≥ 0 for all k > 0, θ ∈ [0, π/2], B ≥ 0, and Σ, Σ' > 0, with the asymptotic group-velocity bound (45) checked separately. A specific numerical counterexample at θ = π/4 would immediately refute the universal claim; if the symbolic conditions hold for all θ, the claim would be proven.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV explicitly restricts the analysis to κ⊥=0 (Sec. IV A) and κb=0 (Sec. IV B); all dispersion relations in Eqs. (50) and (57) are for these two directions only. The full linear system in Eqs. (42)-(44) contains both κb and κ⊥ and is never solved at arbitrary angle. Because the background magnetic field defines a preferred direction, stability of the parallel and perpendicular sectors does not imply stability of oblique modes: the coefficient matrix depends nontrivially on the angle, and instabilities or superluminal asymptotic velocities can appear only for intermediate angles. The conclusion in Sec. VI states that the formulation 'was shown to be always causal and stable in the linear regime' and that small perturbations 'decrease exponentially and propagate subluminally,' with no direction restriction. This is a proof gap in the central claim, not an internal inconsistency; the paper should either provide the general-angle dispersion relation with a stability and causality proof (e.g., Routh-Hurwitz conditions in k^2) or explicitly qualify the claim to parallel and perpendicular perturbations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes the linear stability and causality of the second-order relativistic magnetohydrodynamic theory derived in Refs. [64,65] for a locally neutral, non-resistive, two-component plasma of massless particles. The authors linearize the coupled conservation, Maxwell, and shear-stress equations around global equilibrium, decompose Fourier-space perturbations in an orthonormal basis adapted to the background magnetic field, and derive dispersion relations for perturbations parallel (κ⊥=0) and perpendicular (κb=0) to the field. They compute small- and large-wavenumber asymptotic expansions, plot intermediate-wavenumber modes, and find stable, subluminal spectra that include oscillatory nonhydrodynamic modes at large magnetic field, in contrast to Israel-Stewart theory. They also study a simplified 'longitudinal shear-stress' truncation used in astrophysics and discuss its range of validity. The paper concludes that the formulation is always linearly causal and stable for any magnetic field.","tokens_in":21011,"tokens_out":7731,"duration_ms":80965,"significance":"If the analysis were complete, this would be a useful and timely result: it would certify the kinetic-theory-based second-order MHD of Ref. [64] as linearly stable and causal without fitting parameters, and it identifies qualitative differences from Israel-Stewart theory that are relevant for heavy-ion and astrophysical plasmas. The paper is explicit about the mode structure, gives analytic dispersion relations and asymptotic expansions, and makes a practically useful comparison with the longitudinal approximation. The main caveat is that the universal conclusion is stronger than the evidence: only two propagation directions are solved, and some transverse modes are analyzed only numerically. For that reason the result, as stated, is not yet fully established, although the gap appears fixable.","major_comments":[{"comment":"The claim in Sec. VI that the formulation 'was shown to be always causal and stable in the linear regime' is not supported for perturbations at arbitrary angle to the magnetic field. Section IV A restricts to κ⊥=0 and Section IV B to κb=0; the full linear system in Eqs. (42) and (44) contains both κb and κ⊥ and is never solved at general angle. Because the background field defines a preferred direction, stability of the parallel and perpendicular sectors does not imply stability of oblique modes. The authors should either analyze the general-angle dispersion relation (for example, using Routh-Hurwitz conditions in k^2 for the determinant of the coupled system) or explicitly qualify the conclusion to perturbations parallel and perpendicular to the background magnetic field.","section":"Sec. IV and Sec. VI"},{"comment":"Three transverse modes are excluded from the analytic large-wavenumber analysis with the statement that they have 'rather intricate analytical form and thus were omitted here,' yet the causality criterion in Eq. (45) is an asymptotic condition and the stability claim is made for all positive Σ, Σ′ and all B. The plotted curves in Fig. 3 sample specific temperatures, magnetic fields, and relaxation coefficients, so they do not prove the asymptotic bound or positivity of imaginary parts in general. The authors should provide analytic large-k expansions (or rigorous bounds) for the omitted roots of Eq. (57b), or state which part of the conclusion relies on numerical exploration.","section":"Sec. IV B, Eq. (57b)"},{"comment":"Even within the two parallel and perpendicular sectors, the stability analysis at intermediate wavenumbers is largely numerical. The paper gives small-k and large-k asymptotics and then plots roots for selected parameter values; it does not give a general algebraic argument (e.g., Routh-Hurwitz applied to the dispersion polynomials as functions of k^2) that all modes have positive imaginary part for every k>0 and every positive Σ, Σ′. Thus the phrasing 'always causal and stable' overstates the proof content, even before the oblique-direction gap is considered.","section":"Sec. IV, Eqs. (50) and (57)"}],"minor_comments":[{"comment":"The last term of the shear-stress decomposition is printed as ∆̃π_{kq}(κ̂⊥^μ q̂^ν + κ̂⊥^μ q̂^ν); the second factor should presumably be κ̂⊥^ν q̂^μ to make the combination symmetric.","section":"Eqs. (28) and (30)"},{"comment":"The concluding paragraph attributes the analyzed theory to Ref. [65], while the rest of the paper consistently refers to Ref. [64] as the derivation used; this citation should be corrected.","section":"Sec. VI"},{"comment":"The text says 'The black stars denote the mode...', but the figure legend shows a black curve labeled 'Longitudinal limit'; the wording should be aligned with the actual plot.","section":"Fig. 6"},{"comment":"Several dispersion relations and asymptotic expansions are quoted without derivation, e.g., Eqs. (50), (52), (53), (57), (59), and (60); an appendix or supplementary file with the algebraic steps would make the analysis easier to verify.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The paper is a legitimate follow-up to the authors' own derivation, and the self-citation in Refs. [64,65] is not by itself a concern. The main issue is the mismatch between the universal claim in the abstract and conclusions and the two-direction analysis actually performed; this is fixable either by a general-angle proof or by careful qualification. I would not reject, since the explicit mode structure and the comparison with the longitudinal approximation are valuable for the target community."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the analysis is genuinely new and mostly sound, but the paper overstates its own conclusion. The theory from Ref. [64] is shown to be causal and stable for perturbations strictly parallel and strictly perpendicular to the background magnetic field. That is the scope of the actual calculation; \"always causal and stable in the linear regime,\" as written in Sec. VI and the abstract, goes beyond it.\n\nWhat the paper does well: no one has checked causality and stability of this second-order MHD formalism before, and the linear mode analysis is the right way to do it. The asymptotic expansions are internally consistent; the large-k velocities are at or below c, and the imaginary parts are positive for all modes the authors compute. The numerical plots at intermediate k for the transverse modes from Eq. (57b) are consistent with stability for the parameters shown. The physics differences from Israel–Stewart are real: the magnetic field introduces oscillatory nonhydrodynamic modes at large B, the Alfvén channel becomes parametrically weakly dissipative at large B, and Sec. V gives a useful quantitative comparison showing when the longitudinal-only approximation used in astrophysics is and is not acceptable. The self-citation to the authors' own derivation is not a problem here; the stability result is derived, not assumed.\n\nSoft spots, in proportion. The main one: the general equations (42)–(44) are written with both κb and κ⊥, but only the two special directions are solved. Because the background B breaks isotropy, stability of the parallel and perpendicular sectors does not automatically cover oblique directions. This is a proof gap, not a demonstrated counterexample; the fix is either a general-angle dispersion analysis (Routh–Hurwitz in k², for instance) or an explicit qualification of the claim. Secondary: the \"any B, any cross-section\" claim at intermediate k rests on asymptotics plus a few parameter-restricted plots, so the universal phrasing in Sec. IV B is stronger than the evidence. Minor: some expansions are stated without derivation, and there are small typos, including the conclusion attributing the theory to Ref. [65] while the rest of the paper uses [64].\n\nWhom it is for: relativists working on second-order MHD for heavy-ion and astro applications. I would send it to a serious referee, with the main request being to close or restate the angle gap in the central claim.","headline":"The stability analysis of the Kushwah–Denicol MHD theory is genuinely new and mostly sound, but the 'always causal and stable' claim goes beyond what is proven: only perturbations parallel and perpendicular to B are analyzed.","tokens_in":21485,"tokens_out":3866,"would_cite":true,"duration_ms":39563,"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 shows that the kinetic-theory-derived magnetohydrodynamics for a locally neutral, non-resistive two-component ultrarelativistic plasma is linearly causal and stable around global equilibrium for any magnetic field strength.","keywords":["relativistic magnetohydrodynamics","linear stability","causality","two-component plasma","shear-stress tensor","Boltzmann-Vlasov equation","Alfvén waves","heavy-ion collisions"],"falsifier":"Solve the full linearized Boltzmann–Vlasov equation for a locally neutral, massless, two-component plasma in a homogeneous magnetic field without the 14-moment truncation and without dropping bulk viscosity and diffusion; if any Fourier mode with $\\mathrm{Im}\\,\\omega<0$ or $|\\partial\\,\\mathrm{Re}\\,\\omega/\\partial k|>1$ appears for any $B$, the claim fails. A cheaper version: add a bulk-viscous relaxation or a diffusion current to the linearized equations (20) and recompute the dispersion relations—an unstable or superluminal root would refute the 'always causal and stable' conclusion.","tokens_in":20579,"feed_emoji":"🧲","tokens_out":14371,"duration_ms":117259,"temperature":0.7,"pith_summary":"This paper asks whether a second-order relativistic magnetohydrodynamics, derived from kinetic theory for a locally neutral, non-resistive two-component plasma of massless charged particles, is a physically usable theory: do small perturbations around global equilibrium decay rather than grow, and do they propagate at or below the speed of light? The authors linearize the equations, decompose the perturbations in Fourier space on an orthonormal basis aligned with the magnetic field, and extract the dispersion relations for perturbations parallel and perpendicular to the field. They find that every hydrodynamic and nonhydrodynamic mode has a positive imaginary part and a subluminal asymptotic group velocity, for any value of the magnetic field and for the positive transport coefficients supplied by the kinetic derivation. This matters because the theory is a candidate for the strongly magnetized early stage of ultrarelativistic heavy-ion collisions and for astrophysical plasmas, where an unstable or acausal formulation would make simulations meaningless.","feed_headline":"Magnetized plasma modes stay causal and stable at any B","feed_subtitle":"A kinetic-theory magnetohydrodynamics for neutral plasmas passes causality and stability at any field strength.","key_machinery":"The load-bearing object is the pair of coupled relaxation equations for the total shear-stress tensor $\\pi^{\\mu\\nu}$ and the relative shear-stress tensor $\\delta\\pi^{\\mu\\nu}$, obtained by a 14-moment truncation of the Boltzmann–Vlasov equation, i.e., a closure that keeps only a limited, fixed set of momentum moments. The magnetic field enters through the frequency $\\omega_0 = 2|q|B/(5T)$ and through the antisymmetric tensor $b^{\\mu\\nu}$ built from the field direction; these couple the two shear tensors and are what turns otherwise purely damped nonhydrodynamic modes into oscillatory ones at large $B$. The analytic machinery is the Fourier-space decomposition of all perturbations in the orthonormal basis $\\{u_0^\\mu, b_0^\\mu, \\hat{\\kappa}_\\perp^\\mu, \\hat{q}^\\mu\\}$, which decouples the linear system into smaller blocks. From those blocks the paper derives the dispersion relations and analyzes their roots in the limits $k\\to 0$ and $k\\to\\infty$; the positive definiteness of the relaxation coefficients $\\Sigma$ and $\\Sigma'$ then forces stability, and the bounded root velocities enforce causality.","core_discovery":"The central claim is that the non-resistive magnetohydrodynamics developed in Ref. [64] is linearly causal and stable around global equilibrium for any magnetic field strength. Linearizing about a static equilibrium and decomposing the perturbations in Fourier space, the paper obtains partially decoupled sets of dispersion relations for longitudinal and transverse perturbations; in the small- and large-wavenumber limits all roots have positive imaginary part, and all asymptotic group velocities obey $\\lim_{k\\to\\infty}|\\partial\\,\\mathrm{Re}\\,\\omega/\\partial k|\\le 1$. The paper also displays the intermediate-wavenumber behavior numerically. A distinctive consequence is that at large magnetic field the otherwise purely damped nonhydrodynamic shear modes acquire oscillatory real parts, a behavior the traditional Israel–Stewart formalism cannot produce, while the Alfvén mode's damping is suppressed as the field grows. When the theory is truncated to the longitudinal component of the shear-stress tensor, as is common in astrophysical applications, the spectrum reduces to a subset of these modes and the Alfvén modes become non-dissipative; the paper argues this approximation is justified only for sufficiently large $B$, and not in the early heavy-ion collision regime where $B_0/T_0^2\\sim 0.2$–$2$.","pith_inferences":["Beyond the paper, an analytic Routh–Hurwitz check of the full dispersion polynomials would settle whether causality and stability hold at every intermediate wavenumber, not just in the asymptotic limits plus numerical plots.","The linear result strongly suggests the linearized equations form a mathematically well-posed initial-value problem, but nonlinear well-posedness and shock stability are separate questions the paper leaves open.","The predicted oscillatory shear relaxation at large $B$ could be looked for in kinetic-theory simulations of quark-gluon plasma or in laboratory plasma analogues, though the paper itself does not propose such a test.","The comparison with the longitudinal approximation implies that accretion-disk simulations that keep only longitudinal shear may underestimate dissipation from Alfvén modes whenever $B$ is not asymptotically large; a perturbative inclusion of the transverse modes would quantify the error."],"forward_implications":["The theory derived in Ref. [64] can be used in numerical simulations of strongly magnetized relativistic plasmas without introducing unstable linear modes or superluminal signal speeds, for any value of the field.","At large magnetic field, the nonhydrodynamic shear modes become oscillatory, so a measurement or simulation of shear-stress relaxation in a strongly magnetized plasma would distinguish this theory from Israel–Stewart-type equations.","The Alfvén mode and the transverse-coupled modes become increasingly dissipationless as $B$ grows, making the plasma approach ideal magnetohydrodynamic behavior at large field strength.","The longitudinal-only shear approximation, used in some astrophysical plasma codes, reproduces only a subset of the full spectrum and misses the damping of Alfvén waves; it is reliable only when $B$ is large compared with $T^2$, not in the early heavy-ion collision regime."],"supporting_citations":[{"why":"Supplies the second-order magnetohydrodynamic equations for a locally neutral two-component plasma whose linear modes are analyzed here.","marker":"[64]"},{"why":"Companion derivation of the same two-component magnetohydrodynamics, cited in the conclusions as the theory whose causality and stability were verified.","marker":"[65]"},{"why":"Originates the Boltzmann–Vlasov moment method, including the 14-moment truncation, that produces the shear-stress relaxation equations.","marker":"[48]"},{"why":"Supplies the relativistic 14-moment approximation used to close the kinetic equations.","marker":"[74]"},{"why":"Introduces the orthonormal Fourier basis and linearization procedure the paper uses to derive dispersion relations.","marker":"[75]"},{"why":"Establishes that linear causality and stability in a static background imply the same in moving backgrounds, justifying the rest-frame analysis.","marker":"[76]"},{"why":"Defines the longitudinal-only shear-stress approximation for strong magnetic fields that the paper compares with the full theory.","marker":"[69]"},{"why":"Earlier Israel–Stewart magnetohydrodynamics causality and stability analysis serving as the baseline contrast for the new theory's modes.","marker":"[66]"}],"fun_headline_variants":["Causal and stable MHD for neutral plasma at any B","No B limit: causal MHD for ultrarelativistic two-component gas","Magnetized neutral gas MHD: stable and causal for all B","Causality and stability hold for this MHD irrespective of B","Any B field is fine: this MHD stays causal and stable"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the untested premise that the truncated kinetic equations it inherits—the 14-moment Boltzmann–Vlasov closure with bulk viscosity and diffusion neglected—faithfully describe the plasma's linear response; if the omitted couplings change the mode structure, the all-field stability and causality results do not transfer to the physical system.","fun_headline_variants_meta":{"raw":{"variants":["Causal and stable MHD for neutral plasma at any B","No B limit: causal MHD for ultrarelativistic two-component gas","Magnetized neutral gas MHD: stable and causal for all B","Causality and stability hold for this MHD irrespective of B","Any B field is fine: this MHD stays causal and stable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001695,"raw_usage":{"total_tokens":6705,"prompt_tokens":927,"completion_tokens":5778,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":5685}},"tokens_in":543,"tokens_out":5778,"duration_ms":37133,"temperature":1.0,"reasoning_tokens":5685,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:09:40.824341+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full linearized Boltzmann–Vlasov equation for a locally neutral, massless, two-component plasma in a homogeneous magnetic field without the 14-moment truncation and without dropping bulk viscosity and diffusion; if any Fourier mode with $\\mathrm{Im}\\,\\omega<0$ or $|\\partial\\,\\mathrm{Re}\\,\\omega/\\partial k|>1$ appears for any $B$, the claim fails. A cheaper version: add a bulk-viscous relaxation or a diffusion current to the linearized equations (20) and recompute the dispersion relations—an unstable or superluminal root would refute the 'always causal and stable' conclusion.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the second-order magnetohydrodynamic equations for a locally neutral two-component plasma whose linear modes are analyzed here."},{"cited_title":"Shear Viscosity of hadronic matter at finite temperature and magnetic field","cited_arxiv_id":"2204.01639","evidence_quote":"Companion derivation of the same two-component magnetohydrodynamics, cited in the conclusions as the theory whose causality and stability were verified."},{"cited_title":"Kushwah and G","cited_arxiv_id":null,"evidence_quote":"Supplies the relativistic 14-moment approximation used to close the kinetic equations."},{"cited_title":"Biswas, A","cited_arxiv_id":null,"evidence_quote":"Establishes that linear causality and stability in a static background imply the same in moving backgrounds, justifying the rest-frame analysis."}],"review_version":1}