{"id":"9ea9231b-d925-4162-ae5e-5601bbbb51c1","arxiv_id":"2501.15658","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A three-component Chapman-Enskog model with running coupling and screening mass predicts higher shear viscosity when quarks are included and a decrease of eta/s as the QGP cools.","lead":"This paper calculates how the stickiness of quark-gluon plasma depends on temperature and time using a kinetic theory that treats gluons, quarks, and antiquarks as separate species. It finds that adding quark species raises the viscosity, and that as the plasma cools the viscosity falls toward the quantum lower bound.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (16) keeps only the ∾t-channel pole for every 2→2 amplitude, but the full gg and qq amplitudes also contain ∾u-channel singularities; dropping these and all finite terms leaves the reported η(T) quantitatively unsupported without a numerical check.","rationale":"The reader’s weakest assumption is essentially correct: the numerical results hinge on Eq. (16), and the paper never checks the size of the neglected terms. My stress-test sharpens the objection by noting that the omitted terms are not only finite terms: for identical particles, the full amplitudes contain ∾u-channel poles that are also divergent and are suppressed nowhere in the integration range. This makes the t-channel-only reduction even less safe than the reader’s phrasing suggests. The proposed numerical test would settle the issue directly, because all integrals in the Chapman–Enskog collision kernel are one-dimensional once the Mandelstam variables are fixed, so the comparison is straightforward and unambiguous. I do not recommend moving to REJECT: the qualitative physics, more species raising the phase-space for scattering and a falling T lowering both η and η/s, is plausible and may survive the test. However, the paper’s central quantitative claims, including the detailed comparison with AMPT in Fig. 3, cannot be assessed until the full cross sections are integrated. The reader’s CONDITIONAL verdict is therefore the right one, and the condition should include the numerical check described above.","tokens_in":9601,"tokens_out":6062,"duration_ms":57224,"concrete_test":"Numerically evaluate T_lk and J_lk in Eq. (5) on a temperature grid from 150 to 550 MeV using the full amplitudes of Eqs. (7)–(10), with every singular channel Debye-screened: 1/∾t^2 → 1/(∾t − m_D^2)^2 and 1/∾u^2 → 1/(∾u − m_D^2)^2, while retaining the finite terms. Then recompute η(T) and η/s(T) from Eq. (12) with the same running g(T), m_D(T), and molar fractions. If the resulting curves differ from Figs. 1–3 by more than roughly 20% anywhere in the deconfined window, or if the ∾u-channel contribution to T_lk/J_lk is comparable to the ∾t-channel contribution, Eq. (16) is not a controlled approximation. A minimal first step is to tabulate the four C′′_lk integrals at one benchmark temperature, e.g., T = 350 MeV.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative results in Figs. 1–3 and Eq. (12) are determined entirely by the screened-Coulomb cross sections of Eq. (16). That replacement is not the “divergent channels only” approximation it claims to be. For gg scattering, Eq. (10) contains both ∾t- and ∾u-channel singular terms, symmetric under ∾t ↔ ∾u; for identical qq scattering, Eq. (8) likewise has ∾t- and ∾u-channel poles. Eq. (16) keeps only the ∾t pole for every channel. Consequently, the ∾u → 0 endpoint ∾t → −∾s of the collision integrals in Eq. (5) is left unscreened, and it contributes to T_lk with weight 8/3 + 4∾t/∾s → −4/3, which is of the same logarithmic order as the kept forward term. The finite terms omitted from Eq. (9), such as the −1/(∾s∾u) piece in gq scattering, are also integrated over the full ∾t range and need not be negligible at m_D/T ≈ 2 and α_s ≈ 0.2–0.34. Because C′′_lk is linear in dσ/d∾t, every curve in Figs. 1–3 is controlled by this single unchecked approximation. The claimed quark-induced enhancement and the claimed alignment with pQCD are therefore not yet quantitatively established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript (arXiv:2501.15658, nucl-th) proposes a multi-component Chapman-Enskog (CE) framework for the shear viscosity of a fully thermalized quark-gluon plasma (QGP), treating the system as a three-component mixture of gluons, quarks, and antiquarks. The authors use species-specific elastic 2→2 cross sections from the Standard Model, combine them with a running strong coupling g(T) and running Debye mass m_D(T), and then impose a time-dependent cooling law T(t) = 550 MeV exp(-1.204 t/fm). The central claims are that the inclusion of (anti-)quarks enhances both η and η/s relative to a gluon-only description, that η and η/s decrease over time as the plasma cools, and that these results align with perturbative QCD predictions. The CE results are compared with AMPT transport-model calculations, and the discrepancy is attributed to AMPT's constant cross-section assumption.","tokens_in":9900,"tokens_out":8371,"duration_ms":66020,"significance":"If the central claims were quantitatively established, this would be a useful analytical tool for estimating QGP transport coefficients with species-resolved degrees of freedom and temperature-dependent running parameters. The paper has notable strengths: it makes the multi-component CE approach explicit, includes a straightforward comparison with AMPT, and uses standard pQCD inputs rather than fitting to the target η values. However, the quantitative results rest on an unchecked approximation of the scattering amplitudes, and several internal inconsistencies make the plotted curves non-reproducible from the text. For these reasons, the significance of the claimed 'alignment with pQCD' and the quark-induced enhancement cannot yet be assessed with confidence.","major_comments":[{"comment":"The central numerical results in Figs. 1–3 are produced with the screened t-channel-only cross sections in Eq. (16), but the manuscript gives no quantitative justification for discarding the remaining terms of the full 2→2 amplitudes in Eqs. (7)–(10). Because the collision kernel C''_lk in Eqs. (4)–(5) is linear in dσ/dt, every plotted curve is controlled by this approximation. In particular, Eq. (9) for gq scattering contains a 1/(s u) term, and the gg and identical-fermion amplitudes in Eqs. (10) and (8) have u-channel poles at u = 0 (t = -s), which are not regularized because the Debye mass in Eq. (13) screens only the t-channel. The integrals in Eq. (5) run over the full t ∈ [-s, 0] range, so the u-pole region is unscreened and is of the same logarithmic order as the kept forward pole at m_D/T ≈ 2 and α_s ≈ 0.2–0.34. Please provide a numerical comparison of η(T) computed from the complete amplitudes versus the approximation in Eq. (16), or a bound on the omitted contributions. Without such a test, the claimed quark-induced enhancement and the claimed quantitative agreement with pQCD are not established.","section":"§3, Eq. (16); §2.2, Eqs. (7)–(10)"},{"comment":"There is an internal inconsistency in the overall factor of γ0. Equation (3) reads η = T γ0/(10σ)(x1 C0,1 + x2 C0,2 + x3 C0,3), whereas the C0,i given in Eqs. (2) and (11) already contain one power of γ0; substituting them into Eq. (3) yields a factor γ0^2, as indeed appears in Eq. (12). Please clarify which formula is used to produce the figures. In addition, the stated small-z limit γ0 = 1600/z^2 is not the standard Bessel asymptotics for z→0, where K3(z)/K2(z) ~ 4/z and hence γ0 ~ -40/z; since z = m/T for massless partons, the finiteness and numerical value of the plotted η require an explicit explanation.","section":"§2.1, Eqs. (3), (12); text after Eq. (1)"},{"comment":"The paper never reports the values of the integrals T_lk and J_lk in Eq. (5) or any other numerical intermediate, so the curves in Figs. 1–3 are not independently verifiable from the text. Please provide a table of η(T) and η/s(T) for representative temperatures, or a reproducible specification of the integration procedure, the cutoff or regularization details, and the parameter values used, so that the results can be checked by a reader.","section":"§2.1–§3, Eqs. (5), Figs. 1–3"}],"minor_comments":[{"comment":"The manuscript contains numerous typos, including 'incorpoartes' in the abstract, 'viscsoity' in the introduction, 'significally' in §3, 'adjustemtns' in §4, and 'T able 1' before Table 1; a careful proofread is needed.","section":"Abstract, §1, §3, §4"},{"comment":"The sentence 'where g ∼ m_D/T is defined as (1/g^2(T))^{-1/2}' is circular and confusing; please define g(T) unambiguously, e.g., g(T) = sqrt(4π α_s(T)).","section":"§2.3, Eq. (15)"},{"comment":"The hybrid construction 'N = 1 AMPT viscosity and N = 3 entropy density' is unusual and should be explicitly justified, since the comparison in Fig. 2 mixes different model assumptions in a way that is not explained in the main text.","section":"Fig. 2 caption, §4"},{"comment":"The cooling law in Eq. (20) is introduced as an exponential fit without derivation, uncertainty, or comparison to hydrodynamic expansion models; the manuscript should clearly label this as a phenomenological input and discuss its sensitivity.","section":"§4.1, Eq. (20)"},{"comment":"Because the N = 3 CE expansion in Eq. (2) is taken from Ref. [27], which is the author's own earlier preprint, the present paper should either reproduce the derivation in an appendix or explicitly state the dependence; the notation in Eq. (4) (δN,1, ΔN,1, and the ≃ symbol) is also under-specified.","section":"§2.1, Eq. (2), Ref. [27]"}],"recommendation":"major_revision","confidential_remarks":"The paper is a self-contained follow-up to the author's previous preprint, but the numerical results currently cannot be reproduced from the text because the collision-kernel integrals are not evaluated and the γ0 normalization is inconsistent between equations. The central approximation of keeping only screened t-channel poles is unchecked, so the claim of quantitative agreement with pQCD is stronger than the evidence supports. The topic is suitable for a nuclear theory journal, but the manuscript needs a substantial revision with numerical validation before it can be considered for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a follow-up application of the author's multi-component Chapman-Enskog work (Ref. 27) to a thermal QGP with running Debye mass and coupling, compared against AMPT. The qualitative finding—that quarks raise eta and eta/s and that cooling lowers both—is plausible and consistent with having more degrees of freedom and stronger screening. The paper is clearly written, transparent about the model setup, and the comparison with AMPT is a useful sanity check.\n\nThe soft spots are real and load-bearing. The central quantitative results in Eq. (12) and Figs. 1–3 are computed from the screened-Coulomb cross sections of Eq. (16), which keep only the t-channel pole for every 2→2 amplitude. The full amplitudes in Eqs. (7)–(10) also contain u-channel singularities and finite terms; for gg and identical-fermion scattering, the u-channel pole is left unscreened and contributes at the same logarithmic order as the kept t-channel term. The paper gives no numerical check that these neglected terms are small at the relevant m_D/T ~ 2 and alpha_s ~ 0.2–0.34. Since the collision kernels are linear in dsigma/dt, every curve in Figs. 1–3 is controlled by this single unchecked approximation. That is the main flaw.\n\nTwo other issues: the integrals in Eq. (5) are never evaluated or tabulated, and no code or data are provided, so the plotted curves are not verifiable from the text. And there is an internal inconsistency between Eq. (3), which has a factor gamma0, and Eq. (12), which has gamma0^2. That is not a trivial typo; it changes the magnitude of eta by roughly an order of magnitude and must be fixed.\n\nThe claim that the results \"align with perturbative QCD\" is made without a quantitative benchmark. No comparison to known pQCD results for eta/s(T) is given, so the claim is not yet supported.\n\nThe framework itself is a natural extension of the author's prior work, and the novelty is in the application rather than the formalism. That is fine, but it means the paper stands or falls on whether the approximations are checked. A serious referee should ask for: (1) a derivation of the screened cross sections from the full amplitudes, including the u-channel and finite terms, justified numerically; (2) an explicit evaluation of the collision kernel integrals, even as a small code release; (3) a correction of the gamma0 factor; and (4) a pQCD baseline curve. With those, the paper would be a useful contribution. As written, the numbers are unsupported.\n\nI would send this to peer review—the question is real and the framework is legitimate—but the revision request should be heavy. Not something I'd cite in its current form.","headline":"A plausible extension of the author's own CE framework to QGP with running couplings, but the numbers rest on an unchecked screened-pole approximation and an internal gamma0 inconsistency; needs a careful referee.","tokens_in":10525,"tokens_out":2453,"would_cite":false,"duration_ms":21700,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Mh","25.75.-q"],"model":"deepseek-v4-flash","headline":"This paper claims that a three-species Chapman-Enskog treatment of the quark-gluon plasma, with gluons plus (anti-)quarks, running QCD parameters, and time-dependent cooling, yields shear viscosity and eta/s values that are higher than…","keywords":["quark-gluon plasma","shear viscosity","eta/s","Chapman-Enskog","Debye screening","running coupling","AMPT","heavy-ion collisions"],"falsifier":"Recompute $\\eta$ and $\\eta/s$ with the same Chapman-Enskog kernels but using the full Standard Model amplitudes of Eqs. (7)-(10) with the same running $m_D(T)$ and $\\alpha_s(T)$, at temperatures from 200 to 500 MeV; if the resulting values differ significantly from the screened-Coulomb-only results, or if the quark-on versus quark-off ordering reverses, the central claim fails. Alternatively, run a Green-Kubo or transport simulation with the full amplitudes under the same cooling law and compare the time-dependent $\\eta/s$ curves.","tokens_in":2086,"feed_emoji":"⚛️","tokens_out":3923,"duration_ms":68083,"temperature":0.7,"pith_summary":"The paper tries to establish that the shear viscosity $\\eta$ and the ratio $\\eta/s$ of a fully thermalized quark-gluon plasma are larger when quarks and antiquarks are treated as separate species alongside gluons, rather than folded into an effective gluon gas. It further claims that with a temperature-dependent running Debye mass and running strong coupling, plus an exponential cooling law, both quantities decrease over time as the plasma expands and cools, approaching the KSS lower bound $\\eta/s = 1/(4\\pi)$ from above. These results, obtained from an $N=3$ Chapman-Enskog formula, are presented as more consistent with perturbative QCD expectations than a constant-cross-section transport model whose one-component approximation shows trends that disagree with Bayesian analyses at later times. A sympathetic reader would care because transport coefficients like $\\eta/s$ are the key link between QGP microscopic interactions and observable collective flow in heavy-ion collisions.","feed_headline":"Quark flavors raise QGP shear viscosity above gluon-only estimates","feed_subtitle":"Multi-species Chapman-Enskog with running QCD parameters lowers eta/s as the plasma cools toward the KSS bound.","key_machinery":"The central object is the $N=3$ Chapman-Enskog viscosity expression Eq. (12), built from iterative $C_{0,k}$ coefficients and the linearized collision kernels $C''_{lk}$ of Eqs. (4)-(5). The kernels are linear in the differential cross sections $d\\sigma_{lk}/d\\hat t$, and the paper evaluates them using the screened-Coulomb approximation of Eq. (16): only the divergent $t$-channel survives, so $d\\sigma_{lk}/d\\hat t \\simeq c_{lk}\\,\\alpha_s^2/(\\hat t - m_D^2)^2$, with $m_D(T)$ from Eq. (15), $\\alpha_s(T)$ from Eq. (14), and $\\Lambda_T = 30\\,\\mathrm{MeV}$. The time evolution enters through the exponential cooling law $T(t)$, which carries the viscosity curves from early hot high-viscosity conditions toward the KSS bound as the plasma expands.","core_discovery":"Using a multi-component Chapman-Enskog method in full equilibrium, the paper derives an $N=3$ shear viscosity formula, Eq. (12), for a plasma of gluons, quarks, and antiquarks, built from iterative partial viscosities and linearized collision kernels $C''_{lk}$. Restricting all elastic $2\\to 2$ parton scatterings to their divergent $t$-channel amplitudes, screened by a temperature-dependent Debye mass $m_D(T)$ with the QCD coupling $\\alpha_s(T)$ running, the paper finds that switching on three quark flavors raises both $\\eta$ and $\\eta/s$ relative to the quark-off gluon-gas case. Substituting the cooling law $T(t) = 550\\,\\mathrm{MeV}\\exp(-1.204\\,t/\\mathrm{fm})$ makes both quantities fall monotonically with time, with $\\eta/s$ drifting toward the KSS bound. These trends are contrasted with calculations that use a constant cross section and a one-component gluon-gas viscosity, which the paper argues are oversimplified and disagree with pQCD expectations at lower temperatures.","pith_inferences":["A direct numerical check using the full Standard Model amplitudes of Eqs. (7)-(10), rather than only the screened $t$-channel terms, would show whether the quark-induced enhancement survives when non-divergent momentum-dependent terms are included; the paper does not perform this check.","The same three-species Chapman-Enskog machinery could be extended to partial chemical equilibrium with quark fugacities, where the quark-to-gluon ratio changes over time and would likely alter the cooling trajectory of $\\eta/s$.","Because the $t$-channel-only approximation emphasizes small-momentum-transfer forward scattering, the true $\\eta$ could be lower than presented; a Green-Kubo or transport simulation with the full amplitudes at $T \\sim 200$-$500$ MeV would quantify the difference.","The exponential cooling law with boundary temperatures 550 MeV at $t=0$ and 150 MeV at $t=10$ fm/$c$ is a strong modeling assumption; replacing it with a hydrodynamic cooling profile could shift the time-dependent curves even if the equilibrium transport coefficients are unchanged."],"forward_implications":["Including (anti-)quark partial viscosities raises $\\eta$ and $\\eta/s$ at every temperature considered, so single-species gluon-gas estimates likely underestimate the plasma's resistance to shear.","With running $m_D(T)$ and $\\alpha_s(T)$, both $\\eta$ and $\\eta/s$ decrease as the QGP cools, moving toward the KSS lower bound and implying the plasma appears more nearly perfect at later times.","The $N=3$ framework with running parameters produces trends that are closer to perturbative QCD expectations than constant-cross-section one-component estimates, particularly at lower temperatures and later times.","At $N_f = 0$ the formalism reduces to the $N=1$ gluon-gas limit, providing a smooth interpolation between the multi-component and one-component descriptions.","If correct, the enhanced early-time $\\eta/s$ implies that hydrodynamic modeling of heavy-ion collisions should incorporate flavor-resolved transport coefficients rather than a single effective gluon-gas viscosity."],"supporting_citations":[{"why":"Supplies the $N=3$ iterative Chapman-Enskog expansions for $C_{0,k}$ and the anisotropic collision kernels used in Eq. (12).","marker":"[27]"},{"why":"Provides the AMPT full-equilibrium comparison data and the statement that constant-cross-section results show trends that disagree with Bayesian and pQCD expectations.","marker":"[15]"},{"why":"Gives the Debye mass formula and the screened $t$-channel regularization for pQCD parton scattering amplitudes.","marker":"[22]"},{"why":"Provides the Standard Model cross-section formulas, Eqs. (7)-(10), for quark-antiquark, quark-quark, gluon-quark, and gluon-gluon elastic channels.","marker":"[29]"},{"why":"Supplies the running gauge coupling $g(T)$ expression, Eq. (14), used to define $\\alpha_s(T)$ and the temperature-dependent Debye mass.","marker":"[32]"},{"why":"Supports the consistency of the Chapman-Enskog method with Green-Kubo results for isotropic and anisotropic scatterings, justifying the method's use for QGP viscosity.","marker":"[14]"},{"why":"Defines the ZPC model whose two-body forward-angle scattering cross section underlies the comparison baseline.","marker":"[19]"}],"fun_headline_variants":["Quark flavors boost QGP shear viscosity predictions","Multi-species model raises quark-gluon plasma viscosity","Quark and antiquark contributions lift QGP viscosity","Full equilibrium CE framework yields higher eta/s","Chapman-Enskog with quarks increases viscosity estimates"],"cache_read_input_tokens":12416,"weakest_assumption_plain":"The numerical results rest on the assumption that in every elastic $2\\to 2$ parton channel only the divergent $t$-channel survives, so all differential cross sections collapse to the screened-Coulomb form $d\\sigma_{lk}/d\\hat t \\simeq c_{lk}\\,\\alpha_s^2/(\\hat t - m_D^2)^2$; if the neglected non-divergent terms in the full amplitudes contribute substantially at these temperatures, every viscosity curve changes.","fun_headline_variants_meta":{"raw":{"variants":["Quark flavors boost QGP shear viscosity predictions","Multi-species model raises quark-gluon plasma viscosity","Quark and antiquark contributions lift QGP viscosity","Full equilibrium CE framework yields higher eta/s","Chapman-Enskog with quarks increases viscosity estimates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000233,"raw_usage":{"total_tokens":1510,"prompt_tokens":981,"completion_tokens":529,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":453}},"tokens_in":597,"tokens_out":529,"duration_ms":5875,"temperature":1.0,"reasoning_tokens":453,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:04:33.575352+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute $\\eta$ and $\\eta/s$ with the same Chapman-Enskog kernels but using the full Standard Model amplitudes of Eqs. (7)-(10) with the same running $m_D(T)$ and $\\alpha_s(T)$, at temperatures from 200 to 500 MeV; if the resulting values differ significantly from the screened-Coulomb-only results, or if the quark-on versus quark-off ordering reverses, the central claim fails. Alternatively, run a Green-Kubo or transport simulation with the full amplitudes under the same cooling law and compare the time-dependent $\\eta/s$ curves.","supporting_citations":[{"cited_title":"Shear Viscosity of an $N$-Component Gas Mixture using the Chapman--Enskog Method under Anisotropic Scatterings","cited_arxiv_id":"2406.07764","evidence_quote":"Supplies the $N=3$ iterative Chapman-Enskog expansions for $C_{0,k}$ and the anisotropic collision kernels used in Eq. (12)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Standard Model cross-section formulas, Eqs. (7)-(10), for quark-antiquark, quark-quark, gluon-quark, and gluon-gluon elastic channels."}],"review_version":1}