{"id":"68f7f3df-fd19-4dd1-88d2-fb435289625f","arxiv_id":"2511.14790","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For cold plasmas, the frequency-dependent bulk viscosity follows the Mandelstam-Leontovich/Drude form with high accuracy, and can exceed shear viscosity by orders of magnitude.","lead":"This paper derives bulk viscosity for cold plasmas by solving a kinetic equation for ionization-recombination, and finds that the old Mandelstam-Leontovich formula is nearly exact for such plasmas. It applies the result to the solar chromosphere, arguing that bulk viscosity is essential for acoustic heating models.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ML exactness is proven only for a single ionizable species; for multi-species cold plasmas Eq. (81) generically has multiple relaxation modes, so the abstract's 'exact for cold plasmas' is unsupported.","rationale":"The reader's weakest_assumption is the empirical uncertainty in the ionization cross-section parameters (w, C_W) from McGowan-Clarke. That concern affects the absolute value of τ and the solar-heating predictions, but it does not threaten the ML functional form: for pure hydrogen, the ML shape is exact for any rate coefficient, and for the solar cocktail hydrogen dominates. The more load-bearing issue is the paper's generalization of the exactness claim from a single ionizable species to 'cold plasmas' generally. The matrix structure of Eq. (81) shows that multi-species cold plasmas with comparable ionization contributions have multiple relaxation times, so the single-pole ML formula cannot be exact in general. The solar example works only because of hydrogen dominance, and the paper does not provide a general proof or a quantitative criterion for when the single-mode reduction is valid. This is an internal-scope inconsistency rather than a mere disagreement with external consensus: the abstract claims more than the derivation supports. The reader's conditional verdict remains appropriate, but the required condition should be that the plasma is dominated by one ionizable species, not merely that T<<I_a. The proposed test with a two-alkali cocktail would settle whether the multi-species deviation is practically significant or negligible in the intended regime.","tokens_in":18645,"tokens_out":7897,"duration_ms":81001,"concrete_test":"Solve Eq. (81) for a synthetic cold two-alkali cocktail, e.g., Na and K with comparable abundances and temperature/density chosen via the Saha equation so that both α_a are non-negligible (while still T<<I_a). Compute the exact γ(ω) over a wide frequency range and fit the ML formula (108) with two free parameters (Γ0,∞ and τ). If the maximum deviation between the exact and ML Argand plots exceeds a few percent, the abstract's general claim fails. Additionally, derive the spectral dominance condition (e.g., max_a A_a ι_a / sum_a A_a ι_a ≈ 1) under which the matrix equation reduces to a single effective relaxation time.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, as stated in the abstract, is that the Mandelstam-Leontovich (ML) single-relaxation-time formula is 'actually an exact solution in the case of cold plasmas.' The derivation proves this exactly for pure hydrogen (Section VI, Eqs. (151)-(164)), where the response function reduces to a single pole. For a multi-species cocktail, however, the linearized kinetic equation (81) is a matrix equation with species-dependent relaxation times τ_a. Its solution is generically a sum of relaxation modes, not a single-pole rational function. The numerical agreement shown for the solar chromosphere is not evidence for the general claim; it follows because hydrogen dominates the ionization balance at the chosen height: α_He is exponentially small at T~0.5 eV, and metal abundances are trace. The paper itself attributes the agreement to 'the number of a single type ions dominates' (Section V). Thus the abstract overstates the scope of the result. The claim 'ML approximation has an excellent accuracy for cold plasmas with T<<I_a' is only established for plasmas where one ionizable species dominates, not for arbitrary cold plasmas. This is a correctness/scope risk: a reader applying the formula to a two-alkali cold plasma with comparable ionization potentials would get a misleading single-relaxation-time description.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the complex polytropic index and frequency-dependent bulk viscosity for a cold plasma (T much smaller than the first ionization potentials) from a linearized kinetic equation for electron-impact ionization/recombination. The ionization rates are based on the Wannier near-threshold cross-section, and the resulting matrix system (Eq. (81)) is solved for a six-species H/He/C/Mg/Si/Fe cocktail representing the solar chromosphere at a chosen height. The result is compared with the Mandelstam-Leontovich (ML) single-relaxation-time formula. For pure hydrogen and, more generally, for a single ionizable species with neutral noble-gas diluent, the ML form is shown to be exact (Sec. VI). For the chromospheric cocktail, the Argand plot of the complex polytropic index agrees with the ML fit to about 2% in the extracted relaxation time. The paper also presents thermodynamic formulas for alkali-noble-gas cocktails and discusses implications for acoustic heating of the solar chromosphere. Source code and data are deposited on Zenodo.","tokens_in":18906,"tokens_out":6971,"duration_ms":66931,"significance":"If the derivation is correct, the paper gives a rare first-principles calculation of a kinetic coefficient—bulk viscosity of cold plasmas—expressed in terms of Saha equilibrium, the Wannier cross-section, and a matrix inversion. The microscopic derivation of the phenomenological Mandelstam-Leontovich/Drude form is conceptually valuable, and the explicit hydrogen/noble-gas solution (Eq. (154)) is a clean, testable result. The availability of the code on Zenodo is a further strength. The main limitation is that the exact ML statement is proven only for a single dominant ionizable species; for a general multi-species cold plasma, the matrix equation (81) generically yields multiple relaxation modes, so the unqualified abstract claim overstates the result.","major_comments":[{"comment":"The abstract claims that the Mandelstam-Leontovich approximation is 'actually an exact solution in the case of cold plasmas.' This is too broad. Exact single-relaxation-time behavior is established only for a single ionizable species (Sec. VI, Eqs. (151)-(163)), and for the chromospheric cocktail the conclusion is numerical, with the text itself attributing the agreement to the fact that 'the number of a single type ions dominates' (Sec. V). For arbitrary cold plasmas with several comparable ionizable species, Eq. (81) is a matrix equation with species-dependent relaxation times tau_a, and its solution generically contains multiple relaxation poles, not a single ML pole. The abstract and the summarizing sentence in Sec. V.B should be qualified, e.g., by stating that the ML form is exact when one ionizable species dominates and is an excellent approximation for hydrogen-dominated mixtures","section":"Abstract; Sec. V.B; Sec. VI"},{"comment":"The central linearized kinetic equation (70) is introduced after the sentence 'Omitting the complete derivation.' This is a load-bearing step: Eqs. (69)-(71) and (81) all follow from it, and the accuracy of all later results depends on it. As written, the reader cannot verify the signs, the treatment of the Saha-density derivative (65), or the temperature-oscillation relation (69). The full derivation should be supplied in an appendix (or a referenced companion paper) before the manuscript is published.","section":"Sec. IV.B, Eq. (70)"},{"comment":"The ionization rate (39) is obtained by integrating the Wannier near-threshold cross-section (37) with parameters w≈1.18 and C_W≈2.7 taken from a single early experimental study of H(1s) (McGowan & Clarke 1968). The same Wannier parameters are then applied to all elements in the chromospheric cocktail (He, C, Mg, Si, Fe). The sensitivity of the extracted relaxation time and bulk viscosity to these parameters is not quantified, and the neglect of radiative recombination and multi-step ionization is not justified for the chromospheric conditions. Since the solar application is a headline result, a sensitivity analysis, or at least a careful discussion of why these processes are negligible at the chosen height in the AL08 profile, is needed.","section":"Sec. III, Eq. (39)"}],"minor_comments":[{"comment":"In the definition of the matrix \\hat M_B, the last displayed row appears as 'D_{a_max} L_B', missing the prefactor \\hat\\Gamma_{a_max} that appears in all other rows. Please correct the typo.","section":"Eq. (88)"},{"comment":"The quantity c_e is used in Eq. (81) (in the term (\\iota_a + c_v + 1/c_e)) before it is defined in Eq. (95). Please move the definition earlier or add a note at first use.","section":"Sec. IV.C, Eq. (81)"},{"comment":"The caption says 'ML fit with fitting parameter \\Gamma_{0,\\infty}=0.654', but \\Gamma_{0,\\infty} is computed thermodynamically (Eq. (107)), not fitted. In the Argand plot, the fitting parameter is actually \\tau, and \\tau cancels in the dimensionless plot. Please clarify what is fitted in the figure.","section":"Fig. 1 caption"},{"comment":"The abstract states that magnetic diffusivity, heat conductivity, shear viscosity, and bulk viscosity 'all play an almost equally important part' in wave damping. No calculation of magnetic diffusivity or heat conductivity is presented in this paper; the analysis in Sec. VII compares only bulk and shear viscosity. Please either remove the statement or cite the companion work where the other coefficients are computed.","section":"Abstract; Sec. VII"},{"comment":"There are several typos and duplicated words, e.g., 'according according' in Sec. V.B, 'Naval nozzle' (should probably be 'de Laval nozzle') in Sec. VII.A, and the duplicated line p=Re(\\hat p) in Eqs. (52)-(53). A careful proofread is recommended.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially interesting, and the single-species exactness result is a solid contribution. The main issues are fixable: (i) the abstract overstates the scope of the ML exactness, (ii) the central linearized equation is not derived, and (iii) the solar application rests on a single-experiment cross-section with no sensitivity test. I recommend major revision rather than rejection because the core derivation appears sound for the claimed single-dominant-species regime and the code is provided for reproducibility."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper deserves a serious referee. The genuinely new piece is a first-principles kinetic derivation of the frequency-dependent bulk viscosity of cold partially ionized plasmas, together with the demonstration that the Mandelstam-Leontovich (Drude) single-relaxation-time form is exact for pure hydrogen and numerically excellent for the solar chromosphere cocktail. Earlier work by the same group covered hydrogen polytropic indices; the multi-species matrix formulation and the direct bulk-viscosity result are the step forward.\n\nWhat it does well: the single-species case is proven, not assumed. The algebra closes to a one-pole response (Eqs. 151–164), the three independent determinations of the relaxation time agree to about 2%, and the low- and high-frequency limits correctly reproduce the thermodynamic γ0 and the monatomic 5/3. Code and data are deposited on Zenodo. That is real, checkable evidence.\n\nThe main problem is the abstract. It claims ML is “actually an exact solution in the case of cold plasmas.” That is too broad. For multi-species plasmas, Eq. (81) is a matrix equation with species-dependent relaxation times, whose generic solution is a sum of relaxation modes, not a single-pole rational function. Exactness is proven for one ionizable species. The multi-species agreement shown is for a cocktail where hydrogen dominates the ionization balance — the paper itself attributes the agreement to “the number of a single type ions dominates.” Anyone applying the formula to a two-alkali mixture with comparable ionization potentials would get a misleading answer. This is a scoping error in the claims, not a flaw in the math, and it is fixable.\n\nSecond soft spot: the paper explicitly omits the derivation of the linearized kinetic equation (“Omitting the complete derivation...”). Since Eq. (81) is the load-bearing object, that derivation belongs in an appendix or supplement. Moderate.\n\nThird: the solar application rests on one atmospheric model and on near-threshold Wannier cross-section parameters (w≈1.18, C_W≈2.7) from one experimental study, with no sensitivity analysis; radiative recombination and second ionization are neglected. That is acceptable for an illustration, but “studied thoroughly” overstates it, and the heating numbers should be read as indicative. Minor.\n\nThe stress-test concern is correct and should be pressed in review. The math itself holds up. Audience: plasma kinetic theorists, solar-atmosphere modelers, and anyone working with alkali-noble lab plasmas. Recommendation: accept for peer review, conditional on rewriting the abstract, supplying the omitted derivation, and adding a sensitivity check on the cross-section constants. Also add a working Zenodo link.","headline":"Solid kinetic derivation of cold-plasma bulk viscosity with ML exactness proven for hydrogen — but the abstract overclaims the multi-species case.","tokens_in":19455,"tokens_out":4868,"would_cite":true,"duration_ms":44952,"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":"Cold plasmas have a single-relaxation-time bulk viscosity: exact for pure hydrogen, about 2% accurate for six-element chromospheric mixtures.","keywords":["bulk viscosity","cold plasma","complex polytropic index","ionization-recombination kinetics","single relaxation time","near-threshold ionization cross-section","solar chromosphere","acoustic wave damping"],"falsifier":"A laboratory measurement of the frequency dependence of sound absorption in a cold partially ionized gas, for instance a sodium–neon cocktail, should reproduce the predicted semicircular plot of the complex polytropic index and the relation between the low-frequency slope of its imaginary part and the location of the absorption peak; a deviation beyond the few-percent level would rule out the exact single-relaxation-time claim. A direct re-measurement of the hydrogen near-threshold ionization cross-section that disagrees with σ ∝ (ε/I − 1)^1.18 and a coefficient near 2.7 would also undermine t","tokens_in":18470,"feed_emoji":"⚛️","tokens_out":9240,"duration_ms":80985,"temperature":0.7,"pith_summary":"This paper solves the kinetic equation for ionization–recombination in a cold plasma, meaning a plasma whose temperature is far below the first ionization potentials of its constituents. It derives an explicit expression for the frequency-dependent complex polytropic index, and from it the bulk viscosity. The central result is that this frequency dependence is exactly the textbook single-relaxation-time formula: for pure hydrogen the derivation closes exactly, and for a realistic six-element solar chromospheric cocktail the full kinetic calculation and the single-relaxation-time fit agree to about 2% in the relaxation time. The relaxation time is not a free parameter; it follows from the equilibrium ionization state and the near-threshold electron-impact ionization rate. Having bulk viscosity from first principles matters because bulk viscosity is the dominant sound-damping channel in many partially ionized gases, including the solar chromosphere, where ignoring it would miss most of the acoustic heating.","feed_headline":"Cold-plasma bulk viscosity obeys one relaxation-time formula","feed_subtitle":"Relaxation time comes from ionization-recombination kinetics, so sound damping becomes first-principles.","key_machinery":"The engine of the derivation is the complex generalized polytropic index γ̂(ω) ≡ δp̂/δρ̂ / c_N², the frequency-dependent pressure-to-density response of the plasma; its imaginary part is directly proportional to bulk viscosity. The machinery is a set of linearized kinetic equations for the first-ionization fractions α_a of each element, with a near-threshold electron-impact ionization rate σ ∝ (ε/I − 1)^w (w ≈ 1.18) used to compute the relaxation rates, plus the energy-conservation and charge-neutrality constraints that determine the accompanying temperature and electron-density oscillations. Solving this linear system gives γ̂(ω), from which the bulk viscosity and sound damping are obtained","core_discovery":"For a cold plasma (T much less than the first ionization potentials), the paper claims that the frequency-dependent bulk viscosity is exactly described by the single-relaxation-time formula: a complex generalized polytropic index γ̂(ω) whose real part rises from the equilibrium thermodynamic value γ0 to the monatomic value 5/3 and whose imaginary part traces a semicircle in the complex plane. The relaxation time is not fitted but derived from the linearized kinetic equations for the ionization fractions, coupled to temperature and electron-density oscillations through energy conservation and charge neutrality. For pure hydrogen the derivation is closed, making the single-relaxation-time form","pith_inferences":["A natural extension is to regard any reversible internal degree of freedom that equilibrates through a single bottleneck rate—vibrational relaxation in molecular gases, dissociation-recombination, or chemical reactions—as producing a similar Debye-like bulk-viscosity peak; the paper does not make this broader claim.","The computed long relaxation time at chromospheric conditions shifts most of the acoustic spectrum into the dispersionless absorption band, so chromospheric heating models could be simplified from a frequency convolution to a local heating rate, an implication the authors only sketch.","The 2% agreement between the full kinetic calculation and the single-relaxation-time fit could be used as a quantitative probe: deviations should grow as multiple ionization stages or species with very different ionization potentials contribute significantly, and a deliberate test with a two-species laboratory mixture would map that boundary.","One could invert the paper's logic and use observed chromospheric wave damping as a remote constraint on the near-threshold ionization cross-sections of minor elements like magnesium, silicon, and iron, which are less well measured than hydrogen."],"forward_implications":["In the chromospheric model, bulk viscosity dominates sound-wave damping below a critical frequency of about 35 mHz, so realistic acoustic-heating models must include it or they miss the dominant dissipation channel.","The high-frequency damping rate becomes frequency-independent, meaning the acoustic heating power is directly proportional to the wave spectral density in that band.","The closed formula ζ0 = τ p (γ∞ − γ0), with τ derived from the ionization kinetics, gives a parameter-free prediction for sound absorption in alkali–noble-gas mixtures such as sodium–neon, testable in laboratory plasmas.","For pure hydrogen plasma, the single-relaxation-time description is exact, so no multi-relaxation-time corrections are needed in cold hydrogen plasmas.","The same thermodynamic and kinetic framework yields an explicit expression for the convective-instability criterion, connecting bulk viscosity to models of convection in partially ionized atmospheres."],"fun_headline_variants":["Cold-plasma bulk viscosity gets exact single-relaxation-time formula","Exact formula found for cold-plasma bulk viscosity","Bulk viscosity of cold plasmas: single relaxation time proven exact","Ionization-recombination kinetics yields exact bulk viscosity in cold plasmas"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the near-threshold electron-impact ionization rate, with its power-law form and coefficients taken from a single hydrogen measurement, accurately describes every element in the cocktail and that only the first ionization state participates; if this rate is wrong or higher ionization steps or radiative recombination matter, the predicted relaxation time and bulk viscosity shift.","fun_headline_variants_meta":{"raw":{"variants":["Cold-plasma bulk viscosity gets exact single-relaxation-time formula","Exact formula found for cold-plasma bulk viscosity","Bulk viscosity of cold plasmas: single relaxation time proven exact","Ionization-recombination kinetics yields exact bulk viscosity in cold plasmas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000181,"raw_usage":{"total_tokens":1123,"prompt_tokens":701,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":349}},"tokens_in":445,"tokens_out":422,"duration_ms":4447,"temperature":1.0,"reasoning_tokens":349,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:15:10.107445+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A laboratory measurement of the frequency dependence of sound absorption in a cold partially ionized gas, for instance a sodium–neon cocktail, should reproduce the predicted semicircular plot of the complex polytropic index and the relation between the low-frequency slope of its imaginary part and the location of the absorption peak; a deviation beyond the few-percent level would rule out the exact single-relaxation-time claim. A direct re-measurement of the hydrogen near-threshold ionization cross-section that disagrees with σ ∝ (ε/I − 1)^1.18 and a coefficient near 2.7 would also undermine t","supporting_citations":[],"review_version":1}