REVIEW 4 major objections 4 minor 6 references
Momentum anisotropy from Resistive Magnetohydrodynamics
T0 review · 4 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read A purely electric field can anisotropize a relativistic plasma even when no velocity gradient exists.
desk verdict Interesting physical effect—an electric field alone can source momentum anisotropy—but the printed equations are dimensionally inconsistent and the numerical claims cannot be verified as the text stands. 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 load-bearing object is the source term V_q^{<mu}E^{nu>} — the symmetric traceless product of the net-charge diffusion current and the electric four-field — in the shear-stress evolution equation. It is the only term that can generate shear stress in the homogeneous, locally neutral, B=0 limit, and it closes a feedback loop because shear stress feeds back on the diffusion current through the Omega_Epi term. All transport coefficients are explicit closed functions of the particle density and the intra- and inter-species cross sections; the 14-moment ansatz for the distribution functions and the gradient expansion eliminating the relative heat flow are the enabling approximations.
What would settle it
A moment-truncation-free numerical solution of the Boltzmann–Vlasov equation for a homogeneous, locally neutral, B=0 plasma of massless charged particles with no velocity gradient, initialized at E0≈20 fm^-2, sigma_T/sigma_+-_T=0.1, eta/s=1: if pi_xx/epsilon does not rise to ~0.1 and relax as in Figure 2, the 14-moment closure overstates the field-shear coupling. A second check is an independent derivation of the coefficient of V_q^{<mu}E^{nu>} in the shear equation; any value different from (8/5)tau_pi would falsify the central claim.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is a kinetic-theory derivation in which the 14-moment closure of the Boltzmann–Vlasov equation for massless, oppositely charged particles produces a coupled system for the diffusion current and the shear stress. The structural result is the term V_q^{<mu}E^{nu>} in the shear relaxation equation: it lets an electric field source momentum anisotropy directly, independently of the velocity-gradient terms sigma^{mu nu}. In a homogeneous neutral plasma with B=0, this produces a transient pi_xx/epsilon ~ O(0.1) for E0 ~ 20 fm^-2, relaxing once Maxwell depletion removes the field. The same equations give a relaxation-type Ohm's law at moderate fields,
Load-bearing premise
The equations inherit all transport coefficients and the elimination of relative heat flow from an earlier derivation without re-derivation; if that closure is inaccurate at field strengths up to 100 fm^-2, the field-induced anisotropy could be an artifact of the truncation.
Editorial extensions
If this is right
- Momentum anisotropy in heavy-ion collisions can be sourced electromagnetically, not only by velocity gradients; E0≈20 fm^-2 yields pi_xx/epsilon≈0.1, the same order as typical flow-driven shear corrections.
- Relaxation-type Ohm's law with a field-independent late-time conductivity remains accurate up to E0≈30 fm^-2; beyond that, nonlinear corrections to the diffusion current become visible.
- For large viscosity (eta/s≈5), the charge-current response becomes underdamped and oscillatory, a regime absent from standard second-order viscous hydrodynamics.
- Under Bjorken expansion, the field-induced contribution to pi_xx/epsilon is subleading to the 8epsilon/(45tau) expansion source, bounding how much of observed anisotropic flow could come from electric fields.
Reading between the lines
- If the coupling survives at finite magnetic field and net charge, early-time flow harmonics in heavy-ion collisions could carry an electromagnetic imprint; extending the derivation to B≠0 would quantify this.
- The predicted underdamped oscillations in the charge current at high viscosity could leave an imprint on electromagnetic emission (photons or dileptons), giving an observable beyond the anisotropy itself.
- The homogeneous setup is a clean target for a moment-truncation-free Boltzmann–Vlasov simulation; matching Figure 2 would validate the closure, while a miss would identify the truncation as the source of the effect.
- Because all coefficients are closed-form, the E0≈20 fm^-2 threshold translates into a physical field strength for any plasma once the cross sections are specified, for instance from lattice or kinetic input.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript reports a compact second-order resistive magnetohrodynamic framework for an ultrarelativistic, two-component, locally neutral plasma, obtained from the Boltzmann–Vlasov equation with the 14-moment approximation (with the detailed derivation delegated to the authors' companion paper [3]). The central claims are: (i) the charge diffusion and shear stress equations are coupled by electric-field terms, in particular the source V^{<μ}E^{ν>} in the shear equation; (ii) in a homogeneous system, an electric field alone can generate a transient momentum anisotropy π_xx/ε ~ O(0.1) for E_0 ~ 20 fm^{-2}; (iii) under Bjorken expansion this field-induced anisotropy persists but is subleading to the expansion source; and (iv) large viscosity can produce an underdamped oscillatory regime not present in standard Israel–Stewart theory. The paper gives numerical solutions for the homogeneous and Bjorken cases, using scipy's solve_ivp and explicitly stated initial conditions and parameter choices.
Significance. If the central claim is correct, the paper identifies a genuinely new microscopic channel for generating momentum anisotropy without any flow gradient, with a concrete quantitative estimate (E0~20 fm^{-2} -> π_xx/ε~O(0.1)) and a falsifiable Bjorken-flow prediction. The paper's numerical setup is transparent, and the traceless projection coefficient (8/5 -> 16/15) is internally consistent between Eq. (3) and Eq. (4), which suggests the underlying derivation in [3] may be sound. However, as printed, the equations are not dimensionally self-consistent, one transport coefficient (c_σ) is left unspecified, and the oscillatory regime is asserted but never demonstrated. These issues prevent an independent check of the quantitative central claim, so the manuscript requires substantive revision.
major comments (4)
- [Sec. 1, Eqs. (2)–(3), and Sec. 3, Eq. (4)] The printed equations are dimensionally inconsistent. From Maxwell's equation \.dot E^μ = -V_q^μ in Sec. 2, V_q has dimension fm^{-3}; then in Eq. (2), τ_{Vq}\dot V is fm^{-3} while Γ_{Vq}V and G_E E are fm^{-4}. In Eq. (3), the LHS τ_π\dot π + π has dimension fm^{-4}, but the terms (8/15)εσ, (4/3)θπ, and (10/7)σπ have dimension fm^{-5}; Eq. (4) corresponds instead to the divided form \.dot π + Σ_ππ = (16/15)V E + (8/45)ε/τ - ..., which is not the equation obtained by dividing the printed Eq. (3) by τ_π. The exact ODE system integrated for Figs. 1–3 is therefore not unambiguously specified. Please provide a single, dimensionally consistent form (multiplied or divided) for both equations, state which was numerically solved, and verify that the printed equations reproduce the displayed reductions.
- [Sec. 2a] The claim that switching off Ω_Eπ and Γ_NL in Eq. (2) reduces it to the linear Ohm's law τ_{Vq}\dot V + V = σ_E E is not obtainable from the printed equation. With those coefficients set to zero, Eq. (2) is τ_{Vq}\dot V + Γ_{Vq}V = G_E E; dividing by Γ_{Vq} gives (τ_{Vq}/Γ_{Vq})\dot V + V = σ_E E, and since τ_{Vq}= (1-α_{Vq})/Γ_{Vq}, the relaxation coefficient is τ_{Vq}/Γ_{Vq} = τ_{Vq}^2/(1-α_{Vq}), not τ_{Vq}. This affects the baseline against which Fig. 1 is compared and also the homogeneous shear source through V_q(t). Please clarify whether the intended equation is τ_{Vq}\dot V + V = ..., or \dot V + Γ_{Vq}V = ..., or define Γ_{Vq} as dimensionless.
- [Sec. 1, transport-coefficient list] The text states that every transport coefficient is given in closed microscopic form, but c_σ in Eq. (2) is never defined. Since c_σ multiplies σ^μ_ν V^ν and therefore enters the Bjorken calculation used for Fig. 3, its absence prevents independent verification of the expanding-plasma results. Please provide the explicit expression or, if it is not needed in the limits studied, state so and remove the blanket claim.
- [Abstract and Sec. 4] The underdamped oscillatory regime 'absent from the standard Israel–Stewart formulation' is stated in the abstract and repeated in Sec. 4, but no derivation, threshold condition, or numerical illustration is given anywhere in the paper. The only comment in Sec. 3 is that the oscillatory regime is not reached in Bjorken flow. If this claim is to remain part of the abstract, it must be supported by an explicit analysis (e.g., a linear-stability condition or a figure showing oscillations) or be removed as an overclaim.
minor comments (4)
- [Sec. 2, after Eq. (5)] For the homogeneous case, the energy density ε(t) used as the denominator in π_xx/ε is not specified. If ε evolves due to Joule heating (E·V_q) or remains constant, this should be stated, and the corresponding equation should be given.
- [Sec. 3, Eq. (4) and surrounding text] The Bjorken system is not fully specified: Eq. (4) gives only π_xx and E_x, while Eq. (5) involves π_yy and π_ηη. The assumptions π_xx=π_yy and tracelessness (τ^2π^{ηη} = -2π_xx) should be stated explicitly so the system is closed.
- [Sec. 1, Eq. (1)] The collision terms C[f±, f±'] and C[f±, f∓'] are written in a shorthand that leaves the momentum arguments of the two distribution functions implicit. A concrete notation would improve clarity.
- [Abstract and Introduction] The abstract says 'We derive ...', while Sec. 1 says 'we summarize such a derivation [3]'. Please align the wording to avoid implying that the full derivation is contained here, and consider citing [3] in the abstract.
Circularity Check
No significant circularity; central anisotropy result is a forward solution of previously derived equations, not a fit.
full rationale
The paper's central claim is that an electric field alone sources shear stress through Eq. (3)'s V^{<mu} E^{nu>} term. This term is not fitted to the anisotropy it is used to predict; its coefficient and all transport coefficients are stated as closed-form functions of microscopic cross sections, and the numerical inputs E0, eta/s, and sigma_T/sigma_T^{+-} are chosen, not tuned to match pi_xx/epsilon. The homogeneous and Bjorken results are forward integrations of the ODE system; no output value is fed back into the equations to force the quoted O(0.1) anisotropy. The late-time Ohmic convergence in Fig. 1 is a built-in consistency property (sigma_E = G_E/Gamma_Vq), not a fitted prediction. The main caveat is that Eqs. (2)-(3) are imported from the authors' prior work [3] without derivation, and the printed equations contain apparent typographical/dimensional inconsistencies (e.g., the displayed (8/15) epsilon sigma^{mu nu} in Eq. (3) has the wrong dimension unless a tau_pi factor is missing, and the displayed Ohm's law is not the Gamma_Vq=1 reduction of Eq. (2)). These are correctness/reproducibility concerns, not circular reductions: they do not make the prediction equal to the input. No step in the paper reduces a predicted quantity to a fitted parameter or to a self-citation that is itself the claimed result.
Assumptions & free parameters
free parameters (4)
- sigma_T / sigma_T^+- =
0.1
- eta/s =
1 (also 0.5, 5)
- E0 =
10-100 fm^-2 (20 for central estimate)
- epsilon0, tau0 =
1000 fm^-4, 0.1 fm
assumptions (5)
- domain assumption Boltzmann-Vlasov equation with binary elastic collisions and constant transport cross sections
- domain assumption 14-moment truncation and Landau frame matching
- domain assumption Elimination of the relative energy-diffusion current via gradient expansion
- domain assumption Local charge neutrality (n_q=0) and zero magnetic field
- domain assumption Bjorken boost invariance for the expanding case
Cite this review
Pith. "Pith review of Momentum anisotropy from Resistive Magnetohydrodynamics." pith.science (2026). https://pith.science/paper/4TMSP2DI
@misc{pith2026260721879,
author = {Pith},
title = {Pith review of: Momentum anisotropy from Resistive Magnetohydrodynamics},
year = {2026},
howpublished = {\url{https://pith.science/paper/4TMSP2DI}},
note = {Machine review of arXiv:2607.21879}
}
read the original abstract
We derive relativistic resistive magnetohydrodynamics framework for a two-component ultrarelativistic plasma of massless, oppositely charged particles directly from the Boltzmann-Vlasov equation using the 14-moment approximation. The resulting second-order equations couple the net-charge diffusion current to the shear-stress tensor through the electric field, with all transport coefficients given in closed microscopic form. In the homogeneous limit, the charge-current dynamics is well described by relaxation-type Ohm's law for moderate field strengths, while large viscosity drives the system into an underdamped oscillatory regime absent from the standard Israel-Stewart formulation. Most strikingly, a purely electric field generates sizable momentum anisotropy even without any underlying flow gradient. Under Bjorken expansion this field-induced anisotropy persists but becomes subleading to the hydrodynamic expansion source.
Figures
Reference graph
Works this paper leans on
-
[3]
K. Kushwah, G. S. Denicol, C. V . P. de Brito, Relativistic resistive magnetohydrodynamics for a two-component plasma, Phys. Rev. D 113 (3) (2026) 036021. arXiv:2511.14787, doi:10.1103/PhysRevD.113.036021
arXiv 2026
-
[1]
V . Skokov, A. Y . Illarionov, V . Toneev, Estimate of the magnetic field strength in heavy-ion collisions, Int. J. Mod. Phys. A 24 (2009) 5925–5932. doi:10.1142/S0217751X09047570
-
[2]
J. D. Bjorken, Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Region, Phys. Rev. D 27 (1983) 140–151. doi:10.1103/PhysRevD.27.140
-
[4]
K. Kushwah, G. S. Denicol, Relativistic dissipative magnetohydrodynamics from the Boltzmann equation for a two-component gas, Phys. Rev. D 109 (9) (2024) 096021. arXiv:2402.01597, doi:10.1103/PhysRevD.109.096021
arXiv 2024
-
[5]
Israel, J
W. Israel, J. M. Stewart, Transient relativistic thermodynamics and kinetic theory, Annals Phys. 118 (1979) 341–
1979
-
[372]
doi:10.1016/0003-4916(79)90130-1. 4
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.