{"id":"69659714-94f7-4aac-a405-9cf6e9c74da5","arxiv_id":"2607.21879","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A purely electric field sources momentum anisotropy in a relativistic plasma through a diffusion-shear coupling absent from standard Israel-Stewart hydrodynamics.","lead":"This paper reports that a strong electric field alone can create momentum anisotropy in a relativistic two-component plasma, even with no flow gradient. The effect comes from a kinetic-theory framework imported from the authors' prior work, and it becomes subleading under Bjorken expansion.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (2)–(3) contain dimensional inconsistencies that make the numerical central claim irreproducible from the printed manuscript.","rationale":"The reader's weakest assumption is that the transport equations from [3] are correct. I sharpen this: the manuscript as printed is internally inconsistent, so the equations used for the numerics are not identifiable. This does not prove the physics is wrong, but it means the central claim cannot be independently checked from the text. The concern is load-bearing because Eq. (3)'s V E source is the sole driver of the homogeneous anisotropy; any typo in that equation or in the relaxation terms could materially change π_xx/ε. A re-derivation and recomputation with corrected equations would settle whether the O(0.1) estimate survives. This supports the reader's CONDITIONAL verdict rather than strengthening or reversing it.","tokens_in":4488,"tokens_out":21026,"duration_ms":199546,"concrete_test":"Recompute the homogeneous π_xx/ε evolution using the dimensionally consistent forms of Eqs. (2)–(3) (τ_V \\dot V + V = σ_E E + nonlinear terms; τπ \\dot π + π = (8/5)τπ V^{<μ}E^{ν>} + (8/15)ετπ σ^{μν} + ...) with ε0=1000 fm^-4, η/s=1, σ_T/σ_T^+-=0.1, and E0=20 fm^-2, explicitly specifying |q| and n/T^3 (as implied by a classical massless gas). Check whether π_xx/ε reaches O(0.1) and matches Fig. 2; if it does not, the central quantitative claim is unsupported. Also compare the coefficient of V^{<μ>}E^{ν>} obtained by an independent 14-moment re-derivation with (8/5)τπ.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends entirely on Eq. (3)'s source term (8/5)τπ V^{<μ}E^{ν>}, but the printed equations are not self-consistent enough to verify that term or the numerical results. In Eq. (2), τ_V \\dot V has dimension fm^-3, while Γ_V V and G_E E have dimension fm^-4; moreover, the linear Ohm's law quoted in Sec. 2a is τ_V \\dot V + V = σ_E E, which cannot be recovered from Eq. (2) with Γ_V ≈ 1/τ_V. In Eq. (3), the term (8/15) ε σ^{μν} has dimension fm^-5, whereas all other terms in that equation are fm^-4; its Bjorken reduction in Eq. (4), (8/45) ε/τ, implies the intended term is actually (8/15) ε τπ σ^{μν}, i.e. a missing factor τπ in the displayed Eq. (3). Thus the exact equations solved to produce Figs. 1–3 are not unambiguously specified. Since the quantitative claim (E0~20 fm^-2 ⇒ π_xx/ε~0.1) is sensitive to the coefficient of V^{<μ}E^{ν>} and to the relaxation rates, a mistyped factor could change the result. Delegation to [3] does not resolve this: the manuscript must stand alone.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":4878,"tokens_out":16909,"duration_ms":157576,"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":[{"comment":"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.","section":"Sec. 1, Eqs. (2)–(3), and Sec. 3, Eq. (4)"},{"comment":"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.","section":"Sec. 2a"},{"comment":"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.","section":"Sec. 1, transport-coefficient list"},{"comment":"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.","section":"Abstract and Sec. 4"}],"minor_comments":[{"comment":"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.","section":"Sec. 2, after Eq. (5)"},{"comment":"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.","section":"Sec. 3, Eq. (4) and surrounding text"},{"comment":"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.","section":"Sec. 1, Eq. (1)"},{"comment":"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.","section":"Abstract and Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper leans heavily on the self-cited companion paper [3] for the derivation and for c_σ. Given the dimensional inconsistencies in Eqs. (2)–(3), I recommend that the editor ask the authors to confirm that the printed equations match exactly those solved numerically and those in [3]. The central homogeneous claim may survive after the equations are corrected, but the manuscript as it stands is not reproducible from the printed text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know about this paper. First, the physical effect it reports—a purely electric field sourcing momentum anisotropy through the V_q^{<mu}E^{nu>} term in the shear equation—is real and interesting. The numerical demonstration in homogeneous and Bjorken settings is new as far as I can tell. Second, the printed equations are not self-consistent enough to reproduce the numbers. Eq. (2) has a dimensional mismatch: tau_{V_q} \\dot V is fm^-3, while Gamma_{V_q} V is fm^-4. Eq. (3) has (8/15) epsilon sigma^{mu nu} at fm^-5 while every other term is fm^-4; the Bjorken reduction in Eq. (4) implies the intended term is (8/15) epsilon tau_pi sigma^{mu nu}. These are not cosmetic typos—they change the evolution and the claimed E0 ~ 20 fm^-2 giving pi_xx/epsilon ~ 0.1. The paper also omits c_sigma and asserts an oscillatory regime without derivation. No code is provided, so the figures can't be independently checked from the text.\n\nWhat the paper does well: it isolates a clean mechanism and shows with a simple ODE setup that the coupling can be sizable. The Bjorken section is a nice sanity check showing the effect is subleading under expansion. The transport coefficients are closed-form and come from the authors' earlier PRD derivation [3], so this is a solid extension of an established program, not a speculative new framework.\n\nWhere it's soft: the manuscript is not self-contained. All the heavy lifting is delegated to [3], which is fine as a companion letter, but then the equations printed here must be exactly those solved. They aren't, at least not transparently. Also, the \"linearized Ohm's law\" quoted in Sec. 2a doesn't follow from Eq. (2) unless Gamma_{V_q} is 1/tau_{V_q}, which conflicts with the listed relation tau_{V_q} = (1-alpha_{V_q})/Gamma_{V_q}. These issues are fixable—likely just typos or missing factors—but they need to be corrected before the numerical results can be trusted.\n\nMy take: the underlying physics is probably right, and the paper has a worthwhile message for the heavy-ion/MHD community. It deserves a serious referee, but the referee should be asked to verify the equations and the reproducibility of the figures. I would not cite it in its current form.\n\nRecommendation: send to peer review with a request for major revision focusing on the equations and the numerical implementation.","headline":"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.","tokens_in":5277,"tokens_out":4682,"would_cite":false,"duration_ms":37841,"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":"A purely electric field can anisotropize a relativistic plasma even when no velocity gradient exists.","keywords":["resistive magnetohydrodynamics","momentum anisotropy","charge diffusion current","shear-stress tensor","Boltzmann-Vlasov equation","14-moment approximation","Bjorken flow","heavy-ion collisions"],"falsifier":"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.","tokens_in":4423,"feed_emoji":"⚡","tokens_out":6360,"duration_ms":59961,"temperature":0.7,"pith_summary":"The paper derives a closed second-order resistive magnetohydrodynamics for a two-component ultrarelativistic plasma from the Boltzmann–Vlasov equation, and argues that the shear-stress tensor carries a source term proportional to the product of the net-charge diffusion current and the electric field. That source means a homogeneous, locally neutral plasma with no flow gradients can develop a transient momentum anisotropy from the electric field alone. The authors estimate that E0 around 20 fm^-2 (about 10^19 Gauss) produces pi_xx/epsilon of order 0.1, comparable to ordinary flow-induced shear in heavy-ion collisions. Under Bjorken expansion the coupling persists but is subleading to the expansion source.","feed_headline":"Electric fields alone can create momentum anisotropy","feed_subtitle":"Kinetic theory shows ~10^19 Gauss fields produce anisotropy as large as flow-driven shear in heavy-ion collisions.","key_machinery":"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.","core_discovery":"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,","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Electric fields alone drive momentum anisotropy","No flow needed: electric field creates anisotropy","Field-induced momentum anisotropy from kinetic theory","Electric fields source anisotropy without shear","Momentum anisotropy from electric fields, not flow"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Electric fields alone drive momentum anisotropy","No flow needed: electric field creates anisotropy","Field-induced momentum anisotropy from kinetic theory","Electric fields source anisotropy without shear","Momentum anisotropy from electric fields, not flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000156,"raw_usage":{"total_tokens":1015,"prompt_tokens":668,"completion_tokens":347,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":412,"completion_tokens_details":{"reasoning_tokens":284}},"tokens_in":412,"tokens_out":347,"duration_ms":3740,"temperature":1.0,"reasoning_tokens":284,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T06:25:01.657389+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}