{"id":"1cbf4ef3-c797-4f09-8cc5-16caed426595","arxiv_id":"2504.13352","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"MHD equations are derived from non-equilibrium thermodynamics, yielding generalized Ohm, Fourier, and Newton constitutive laws that reduce to Chapman-Enskog form in the ideal gas limit.","lead":"This tutorial derives magnetohrodynamics (MHD) from non-equilibrium thermodynamics, extending the standard equations beyond the ideal gas limit. It shows that Ohm's and Fourier's laws gain electron chemical potential terms in dense plasmas, with the ideal-gas Chapman-Enskog result as a limit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (35) reduction is not justified: for a two-species plasma d_e + sum d_i = 0, so 'd_e >> sum d_i' cannot hold; the generalized Ohm/Fourier laws may fail when ion chemical potential gradients are not negligible.","rationale":"The paper is a clear tutorial that correctly assembles conservation laws, entropy production, Onsager-Casimir symmetries, and Green-Kubo/Irving-Kirkwood connections. The algebraic structure of Eqs. (49), (53), (68) and the ideal-gas reduction to Chapman-Enskog are internally coherent. However, the single-fluid reduction at Eq. (35) is the hinge of the derivation, and its stated justification is flawed: for a two-species plasma, sum_alpha d_alpha = 0 makes d_e >> sum_i d_i impossible; the actual small parameter is the charge-to-mass ratio and the specific-chemical-potential ratio. A correct reduction introduces (mu_e - mu_i)/(z_e - z_i), so Eq. (49) is an approximation that requires |grad mu_i| << |grad mu_e|. This is likely true for a pure electron-proton plasma in the small mass-ratio limit, but it is not guaranteed for the broad class of 'non-ideal' or molecular plasmas the tutorial advertises. The reader's weakest assumption correctly targets this step; I refine it by noting the first inequality is mathematically impossible as stated. The priority claim ('first complete derivation') is also somewhat assertive given de Groot-Mazur's prior treatment, but that is secondary to the physical closure issue. Overall, a CONDITIONAL verdict remains appropriate: the derivation is sound within a restricted electron-ion single-fluid regime, but the scope and the inequality assumptions need to be restated and quantified.","tokens_in":20777,"tokens_out":9524,"duration_ms":88645,"concrete_test":"Analytically re-derive Eq. (35) for a two-species electron-ion plasma using only the exact relation d_i = -d_e and j = (z_e - z_i)d_e. Show that the entropy production becomes sigma_S = ... + j dot [E'/T - (1/(z_e - z_i)) grad((mu_e - mu_i)/T)] - ... . Then compute the ratio R = |grad(mu_i/T)| / |grad(mu_e/T)| for a model partially ionized hydrogen plasma (e, p, H) or a strongly coupled ion mixture using a simple equation of state. If R is not much smaller than 1 (or if z_i/z_e is not much smaller than 1), Eqs. (49) and (53) must be modified, and the manuscript's claim of a complete single-fluid closure for non-ideal plasmas is overstated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central step that turns the multi-species entropy production (34) into the single-fluid form (35) is the assertion in Sec. II.C that 'd_e >> sum_i d_i, so j ≈ z_e d_e' and 'µ_e >> sum µ_i'. The first inequality is internally inconsistent with the definition d_alpha = rho_alpha(V_alpha - V): summing over all species gives sum_alpha d_alpha = 0, so for an electron-ion plasma with one ion species, d_i = -d_e exactly, and d_e >> |sum_i d_i| is false. The current is dominated by electrons not because d_e >> d_i, but because |z_e| >> |z_i| (charge-to-mass ratio). More importantly, the entropy production contains sum_i d_i · grad(mu_i/T), which cannot be discarded merely because the current is electron-dominated. A correct two-species reduction gives the force conjugate to j as E'/T - [1/(z_e - z_i)] grad[(mu_e - mu_i)/T], not E'/T - (1/z_e) grad(mu_e/T). The paper's Eq. (49) therefore implicitly assumes |grad mu_i| << |grad mu_e| and |z_i| << |z_e|; the latter is true for protons/electrons, but the former is an additional, unquantified assumption that can fail in partially ionized or molecular plasmas with comparable species masses, or in strongly coupled ion-mixture regimes. This directly affects the paper's advertised generality ('valid beyond the usual ideal gas approximation') and the claimed first complete derivation of closed MHD equations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a tutorial derivation of single-fluid magnetohydrodynamic (MHD) equations from non-equilibrium thermodynamics, aiming to be valid beyond the ideal-gas approximation. It derives conservation laws, constructs the entropy production rate, obtains linear constitutive relations for the current density, heat flux, and viscous stress tensor (generalized Ohm, Fourier, and Newton laws), reduces the number of independent transport coefficients using symmetry arguments and Onsager-Casimir relations, and connects the transport coefficients and equation-of-state quantities to microscopic dynamics through the Irving-Kirkwood procedure and Green-Kubo relations. The authors show that the Chapman-Enskog form of MHD is recovered in the ideal-gas limit, and they discuss connections to extended and resistive MHD models.","tokens_in":21184,"tokens_out":13975,"duration_ms":128880,"significance":"If correct, the derivation would provide a useful structural framework for dense-plasma MHD, showing that non-ideal equations of state modify Ohm's and Fourier's laws through chemical-potential gradients and introduce bulk and cross viscosity terms. The symmetry analysis is a genuine strength: the reduction from 225 phenomenological coefficients to 16 independent ones via inversion symmetry, rotational symmetry about the magnetic field, and Onsager-Casimir relations is explicit and checkable. The Green-Kubo and Irving-Kirkwood connections also give a principled route to compute the coefficients from particle trajectories. However, the central single-fluid reduction of the multi-species entropy production contains an internally inconsistent inequality, and the resulting generalized Ohm's and Fourier's laws are derived under unquantified assumptions about ion chemical-potential gradients. This affects the main advertised claim that the equations are valid beyond the usual ideal-gas approximation.","major_comments":[{"comment":"The reduction from Eq. (34) to Eq. (35) rests on the statement that 'd_e ≫ Σ_i d_i'. This is inconsistent with the definition d_α = ρ_α(V_α − V) in Eq. (7), because Σ_α d_α = 0 identically, so for a single ion species Σ_i d_i = −d_e and the magnitudes are exactly equal. The current is electron-dominated because |z_e| ≫ |z_i|, not because the electron diffusion flux is larger than the ion diffusion flux. The reduction should use Σ_α d_α = 0 to eliminate the ion fluxes; the exact two-species entropy production contains the force (1/(z_e − z_i))∇[(μ_e − μ_i)/T] conjugate to j. The simplified form in Eq. (35) follows only after additionally assuming |z_i| ≪ |z_e| and |∇(μ_i/T)| ≪ |∇(μ_e/T)|, neither of which is stated or quantified.","section":"Section II.C, Eq. (35)"},{"comment":"Because the reduction in Eq. (35) is not generally valid, the generalized Ohm's law and Fourier's law as written depend only on ∇(μ_e/T), whereas the correct single-fluid reduction of the species sum involves ∇[(μ_e − μ_i)/T]/(z_e − z_i). For an electron–proton plasma this reduces to the paper's form, but in plasmas with heavier multiply charged ions, multiple ion species of comparable mass, or strong ion–ion coupling, the ion chemical-potential gradient need not be negligible. The paper should either adopt the exact two-species form or give explicit, quantitative conditions under which the simplified form applies. This qualification bears directly on the abstract's claim that the derivation is 'valid beyond the usual ideal gas approximation'.","section":"Eqs. (49) and (53)"}],"minor_comments":[{"comment":"There is a typographical error: 'Champman-Enskog' should be 'Chapman-Enskog'.","section":"Section IV.B"},{"comment":"The text following Eq. (65) lists 'η1, η2, η3, η4 and η5', but only η0 through η4 are defined; the list should refer to η0, η1, η2, η3, η4. Also, in Eq. (68c) the symbol 'ηB_o' should be 'ηB_0'.","section":"Section II.D.5"},{"comment":"The paper itself states that Eq. (97) for the excess electron chemical potential is only approximate because derivatives of pex depend on higher-order correlation functions, and that classical descriptions of dense electron-ion systems suffer Coulomb collapse. These are important limitations of the 'closed and self-consistent' description advertised in the introduction, and they should be reflected more prominently in the abstract or conclusions.","section":"Section III.C"},{"comment":"The sentence 'the average may be dropped so that instantaneous representations of the fluxes are obtained' is imprecise: what is needed is the instantaneous phase-space expression from which the ensemble average is then taken, not the dropping of the average from an averaged quantity. Rephrasing would avoid confusion.","section":"Section III.A"}],"recommendation":"major_revision","confidential_remarks":"The paper is a tutorial with no new data, but its structural derivation is pedagogically useful. The main derivation is fixable: replace the invalid inequality d_e ≫ Σ_i d_i with the exact constraint Σ_α d_α = 0 and restate the generalized Ohm and Fourier laws with the appropriate two-species chemical-potential combination, then quantify the electron-dominance conditions. The authors may also wish to soften the 'first complete derivation' claim, given the prior literature on non-equilibrium thermodynamics of electromagnetic systems that they themselves cite."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this as a tutorial, not a new result. Its value is organizational: it maps the non-equilibrium-thermodynamics route to single-fluid MHD, with symmetry reductions, Green-Kubo/Irving-Kirkwood links, and a clean recovery of Chapman-Enskog in the ideal gas limit. The generalized Ohm/Fourier laws with electron chemical potential gradients are presented clearly, and the viscosity tensor with bulk/cross terms is a nice addition. The internal algebra checks out and assumptions are mostly stated. But the reduction from the multi-species entropy production to the single-fluid form at Eq. (35) is not justified as written: d_e >> sum_i d_i cannot hold because sum_alpha d_alpha = 0, so d_i = -d_e for a two-species plasma. The correct conjugate force is E'/T - [1/(z_e - z_i)] grad[(mu_e - mu_i)/T], and the paper's form follows only under additional conditions (|z_e| >> |z_i| and |grad mu_i| << |grad mu_e|) that are plausible for proton-electron plasmas but unstated. This matters for the paper's advertised generality, since those conditions fail in molecular or partially ionized plasmas and ion mixtures. I'd ask the authors to fix this and quantify the regime of validity. The 'first complete derivation' claim should also be softened; most of the framework is in de Groot and Mazur, and the novel pieces are incremental. No circularity: transport coefficients and EOS are external inputs. No new numbers, but that's fine for a tutorial. For dense-plasma researchers, this is a useful reference and worth serious refereeing, provided the reduction is corrected. Recommend: engage, with revision required on Eq. (35).","headline":"A useful tutorial synthesis, but the single-fluid reduction at Eq. (35) is flawed as written and needs correction before this can be relied on.","tokens_in":21665,"tokens_out":4629,"would_cite":false,"duration_ms":39620,"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":"This paper derives closed magnetohrodynamic equations from non-equilibrium thermodynamics, valid beyond the ideal gas limit, with the standard dilute-gas kinetic-theory closure recovered as the ideal gas case.","keywords":["magnetohydrodynamics","non-equilibrium thermodynamics","generalized Ohm's law","Fourier's law","Newton's law","electron chemical potential","Green-Kubo relations","transport coefficients"],"falsifier":"Compute, from a molecular dynamics simulation of a dense plasma in a steady state with a fixed density gradient, uniform temperature, and no applied electric field in the flowing frame, the current density that appears. The generalized Ohm's law predicts a current proportional to $\\nabla(\\mu_e/T)$; the standard ideal-gas MHD prediction is zero. If the measured current vanishes or lacks the predicted chemical-potential dependence, the central claim fails.","tokens_in":20619,"feed_emoji":"🧲","tokens_out":10071,"duration_ms":85844,"temperature":0.7,"pith_summary":"The paper sets out to show that a closed single-fluid set of magnetohydrodynamic equations can be derived from non-equilibrium thermodynamics alone, without assuming an ideal gas equation of state. It claims that for non-ideal plasmas the constitutive relations take structurally different forms: Ohm's law and Fourier's law contain a gradient of the electron chemical potential over temperature, and Newton's law contains bulk and cross viscosity terms that vanish for a monatomic ideal gas. If correct, this means dense, strongly coupled, or degenerate plasmas cannot be modeled by inserting non-ideal transport coefficients into the standard ideal-gas MHD structure; the equations themselves must change. The derivation also supplies a microscopic closure, connecting every transport coefficient to equilibrium time correlations and the excess equation of state to pair correlations.","feed_headline":"MHD equations now hold for non-ideal plasmas","feed_subtitle":"A derivation from entropy production puts chemical-potential gradients into Ohm's and Fourier's laws, recovering standard kinetic-theory…","key_machinery":"The load-bearing object is the single-fluid entropy production rate after the electron-mass-ratio expansion, Eq. (35): $\\sigma_S = -(1/T^2)\\, q\\cdot\\nabla T + j\\cdot[E'/T - (1/z_e)\\nabla(\\mu_e/T)] - (1/T)\\hat{P}:\\nabla V$. This bilinear form identifies conjugate force-flux pairs, fixes the diffusive entropy flux, and dictates which thermodynamic forces enter the linear constitutive relations. Inversion and rotational symmetry about the magnetic field eliminate vector-tensor coupling and reduce each second-rank transport tensor to three coefficients, while time-reversal symmetry relates the electrothermal and thermoelectric tensors. The viscosity tensor is reduced to seven independent coefficients by symmetry and written in components parallel, perpendicular, and cross to the magnetic field. The macroscopic fluxes are expressed in terms of particle phase-space averages by the Irving-Kirkwood procedure, and their equilibrium autocorrelation functions yield the transport coefficients through Green-Kubo relations, closing the system.","core_discovery":"The central discovery is that the entropy production rate of a quasi-neutral electron-ion plasma, after expanding in the small electron-to-ion mass ratio, can be written in single-fluid form with heat flux, electrical current, and viscous stress as the only dissipative fluxes. From this bilinear form, linear constitutive relations follow, and spatial symmetry decouples vector and tensor forces while time-reversal symmetry reduces the coupling between current and heat flux. The resulting generalized Ohm's law is $j = \\sigma \\cdot [E' - (T/z_e)\\nabla(\\mu_e/T)] + \\phi \\cdot \\nabla T$, Fourier's law is $q = -\\lambda \\cdot \\nabla T + \\varphi \\cdot [E' - (T/z_e)\\nabla(\\mu_e/T)]$, and the viscous stress tensor has seven independent coefficients including bulk viscosity and a cross term. The paper shows explicitly that substituting the ideal-gas expression for $\\nabla(\\mu_e/T)$ reproduces the standard kinetic-theory forms of Ohm's and Fourier's laws, and that dropping bulk and cross viscosity reproduces the conventional magnetized-plasma form of Newton's law. Thus the traditional MHD equations appear as the ideal gas limit of a more general non-ideal closure.","pith_inferences":["The single-fluid closure rests on electron-dominated diffusion, so in partially ionized or molecular plasmas with comparable species masses, the reduction to Eq. (35) would need explicit multi-species diffusion terms and the generalized Ohm's law would likely require modification.","If the chemical-potential gradient term is real, then a dense plasma in a steady state with a density gradient, uniform temperature, and no applied electric field should exhibit a diffusive current; a molecular dynamics simulation could isolate this term from the usual resistive response.","The appearance of bulk and cross viscosity suggests that MHD closures for plasmas with electronic excitation or molecular degrees of freedom need at least these additional coefficients, and Green-Kubo evaluation from molecular dynamics would give their magnitude and dynamical importance.","The derivation implies that dense-plasma MHD codes using the ideal-gas form of Ohm's law are missing a term of the same order as the non-ideal part of the electron chemical potential; estimates from equation-of-state data could show where this correction matters."],"forward_implications":["In dense, strongly coupled, or degenerate plasmas, Ohm's law and Fourier's law contain a term proportional to $\\nabla(\\mu_e/T)$; using the pressure-gradient form of ideal-gas MHD would misrepresent the driving force for current and heat flow.","Newton's law for a non-ideal plasma includes bulk viscosity and a cross-viscosity coefficient that vanish for a monatomic ideal gas, so these terms should appear whenever internal degrees of freedom such as electronic excitation or molecular states are active.","Because every transport coefficient is expressed as a Green-Kubo autocorrelation of phase-space-averaged fluxes, the full MHD closure for a given plasma can in principle be obtained from equilibrium molecular dynamics trajectories without a kinetic-theory solution.","The standard dilute-gas MHD equations are exactly recovered when the ideal-gas equation of state is inserted, including an explicit mapping between the generalized electrothermal and thermoelectric coefficients and the conventional ones.","Phenomena often labeled two-fluid effects, including the Hall effect, appear in the single-fluid description through the off-diagonal components of the resistivity tensor; only the inertial correction is genuinely second order and omitted."],"supporting_citations":[{"why":"It supplies the standard dilute-gas kinetic-theory derivation that the general equations must reproduce in the ideal gas limit.","marker":"[7]"},{"why":"It supplies the conventional magnetized-plasma form of viscous stress to which the non-ideal Newton's law is compared.","marker":"[8]"},{"why":"It provides the non-equilibrium thermodynamics framework, including entropy production and linear constitutive relations, on which the derivation is built.","marker":"[9]"},{"why":"It supplies the reciprocal relations and regression hypothesis used to relate coupled transport coefficients.","marker":"[10]"},{"why":"It extends the reciprocal relations to magnetized systems through time-reversal symmetry, which relates the electrothermal and thermoelectric tensors.","marker":"[11]"},{"why":"It provides the phase-space averaging procedure that expresses macroscopic fluxes as averages of particle positions, velocities, and forces.","marker":"[12]"},{"why":"It supplies one of the Green-Kubo relations connecting transport coefficients to equilibrium autocorrelation functions.","marker":"[13]"},{"why":"It supplies the complementary Green-Kubo relation used for the same microscopic closure.","marker":"[14]"},{"why":"It provides the standard kinetic-theory transport equations used as the explicit ideal-gas comparison and for mapping the generalized coefficients.","marker":"[27]"}],"fun_headline_variants":["MHD upgrades: Ohm, Fourier get chemical-potential terms","Non-ideal plasmas get their own MHD equations","Entropy production yields generalized Ohm and Fourier","MHD beyond ideal gas: laws from thermodynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The single-fluid closure assumes that the electron diffusion flux and the electron specific chemical potential per unit mass are large enough to dominate the ion contributions, so that the multi-species entropy production can be rewritten with electron quantities alone.","fun_headline_variants_meta":{"raw":{"variants":["MHD upgrades: Ohm, Fourier get chemical-potential terms","Non-ideal plasmas get their own MHD equations","Entropy production yields generalized Ohm and Fourier","MHD beyond ideal gas: laws from thermodynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000273,"raw_usage":{"total_tokens":1641,"prompt_tokens":953,"completion_tokens":688,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":625}},"tokens_in":569,"tokens_out":688,"duration_ms":7146,"temperature":1.0,"reasoning_tokens":625,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:10:48.085923+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, from a molecular dynamics simulation of a dense plasma in a steady state with a fixed density gradient, uniform temperature, and no applied electric field in the flowing frame, the current density that appears. The generalized Ohm's law predicts a current proportional to $\\nabla(\\mu_e/T)$; the standard ideal-gas MHD prediction is zero. If the measured current vanishes or lacks the predicted chemical-potential dependence, the central claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the conventional magnetized-plasma form of viscous stress to which the non-ideal Newton's law is compared."},{"cited_title":"Alfv \\'e n ,\\ @noop journal journal Nature \\ volume 150 ,\\ pages 405 ( year 1942 ) NoStop","cited_arxiv_id":null,"evidence_quote":"It provides the non-equilibrium thermodynamics framework, including entropy production and linear constitutive relations, on which the derivation is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the reciprocal relations and regression hypothesis used to relate coupled transport coefficients."},{"cited_title":"Beresnyak ,\\ @noop journal journal Living Reviews in Computational Astrophysics \\ volume 5 ,\\ pages 2 ( year 2019 ) NoStop","cited_arxiv_id":null,"evidence_quote":"It extends the reciprocal relations to magnetized systems through time-reversal symmetry, which relates the electrothermal and thermoelectric tensors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the phase-space averaging procedure that expresses macroscopic fluxes as averages of particle positions, velocities, and forces."},{"cited_title":"Burlaga ,\\ @noop journal journal Space science reviews \\ volume 39 ,\\ pages 255 ( year 1984 ) NoStop","cited_arxiv_id":null,"evidence_quote":"It supplies one of the Green-Kubo relations connecting transport coefficients to equilibrium autocorrelation functions."},{"cited_title":"Wolff \\ and\\ author A","cited_arxiv_id":null,"evidence_quote":"It provides the standard kinetic-theory transport equations used as the explicit ideal-gas comparison and for mapping the generalized coefficients."}],"review_version":1}