{"id":"f2481c5d-bc46-4c58-8446-6ee536749fac","arxiv_id":"2607.18581","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a quasiparticle kinetic model, expansion-induced momentum anisotropy reduces shear and bulk viscosities of baryon-asymmetric QCD matter, while baryon asymmetry increases them.","lead":"This paper calculates how shear and bulk viscosities of the quark-gluon plasma change when the plasma is both baryon-asymmetric and momentum-anisotropic from fast expansion. It predicts anisotropy lowers the viscosities while baryon asymmetry raises them, shifting the expected sound attenuation, flow character, and conformal-symmetry measures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Scalar η extracted via Eq. (27) is not uniquely defined for an anisotropic plasma; Eq. (28) uses an implicit angular average, so the central 'anisotropy decreases η' claim is ambiguous.","rationale":"The reader's conditional verdict is reasonable, but the weakest assumption is not primarily the relaxation-time approximation or the O(ξ) truncation. The more fundamental problem is the definition of the quantity being computed: a single scalar shear viscosity cannot describe a uniaxial anisotropic fluid. The paper's Eq. (27) is the isotropic Navier-Stokes form, while the background distribution has a preferred direction. The angular integrals in the derivation implicitly average c over the sphere, which is equivalent to choosing the isotropic part of the response. That choice should be stated and the results reinterpreted. This does not necessarily overturn the qualitative conclusion, but it is a necessary clarification before the central claim can be accepted. The reader's concern about τ is a modeling uncertainty; the O(ξ) concern is quantitative; the current concern is definitional and affects the meaning of every plotted η. I therefore keep the conditional verdict, i.e., no change to the reader's verdict.","tokens_in":20202,"tokens_out":11461,"duration_ms":142087,"concrete_test":"Set n=ẑ, use Eq. (26), and compute the effective shear viscosity for two pure-shear velocity gradients: W^A with only W_{xy}=W_{yx}=1 and W^B with only W_{xz}=W_{zx}=1 (both traceless). Form η_eff = -∆T^{ij}W_{ij}/(W_{kl}W_{kl}) from the full angular integrals in Eq. (26) without replacing c by 1/3. If η_eff(A) ≠ η_eff(B) at O(ξ), the single η in Eq. (28) is not a well-defined tensorial viscosity. Then rerun Figures 5, 6, and 8 using a stated projection (e.g., transverse projector) and check whether the qualitative result survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result is that expansion-induced anisotropy decreases η and ζ. But the paper never defines a unique η for an anisotropic medium. Eq. (27) assumes the isotropic constitutive relation ∆T^{ij} = -η W^{ij} - ζδ^{ij}∂_l u_l. Substituting the anisotropic distributions (1)-(3) into Eq. (26) produces angular integrands containing c(α,θ,φ)=(p·n)^2/p^2 with n fixed. For fixed n, ∫dΩ p^i p^j p^k p^l c is not proportional to the isotropic rank-4 tensor; it contains n^i n^j, δ^{ij}n^k n^l, and n^i n^j n^k n^l terms. Consequently ∆T^{ij} has different shear-response coefficients for, e.g., W_{xy} and W_{xz}. The closed forms (28)-(29) are obtained by replacing c with its full-solid-angle average 1/3, an implicit isotropic projection. The paper does not state this projection or justify it. Thus 'η' in Figures 5, 6, 8 is an angle-averaged coefficient, not the shear viscosity in the usual tensorial sense. The central claim 'anisotropy decreases η' is therefore ambiguous until a projection is specified. This is a structural issue, independent of the relaxation-time approximation and the O(ξ) truncation.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the viscous transport properties of an expanding, baryon-asymmetric quark-gluon plasma with a weak momentum-space anisotropy. Working in the relaxation-time approximation (RTA) to the relativistic Boltzmann equation with Romatschke–Strickland deformed distributions expanded to first order in the anisotropy parameter ξ, and with temperature-, chemical-potential-, and ξ-dependent quasiparticle masses, the authors derive closed-form expressions for the shear viscosity η and bulk viscosity ζ (Eqs. (29) and (38)) and, using these, compute the Prandtl number, Reynolds number, η/s, and ζ/s. The central claims are that expansion-induced anisotropy reduces both η and ζ, that finite baryon asymmetry increases them, and that the combined effects leave distinct signatures on the derived observables—with η/s approaching the conjectured 1/4π bound in the baryonless anisotropic case and ζ/s exhibiting a temperature dip whose depth and position shift with μ and ξ. The paper is essentially a calculational survey, with numerical results displayed in eight figures for T ∈ [0.16, 0.64] GeV, μ = 0 or 0.1 GeV, and ξ = 0 or 0.4.","tokens_in":20594,"tokens_out":30661,"duration_ms":323961,"significance":"If the results hold, the calculation provides a self-contained RTA/quasiparticle description of momentum transport in the early-stage QGP, quantifying two competing effects: anisotropy (from longitudinal expansion) suppresses η and ζ, while baryon asymmetry enhances them. The closed-form expressions (29) and (38) are easy to reuse, the calculation is parameter-driven (ξ, μ, T are inputs; no output fitting), and the qualitative predictions—most notably the reduction of η/s toward 1/4π with anisotropy and the opposite shift from μ—are falsifiable against more complete kinetic-theory or hydrodynamics treatments. A strength is the self-consistency of the setup: the ξ-dependent quasiparticle masses are computed from the same deformed distributions that enter the transport integrals, so the mass and transport inputs are not independently fitted. The significance is moderated by the issues raised in the major comments: the extraction of scalar viscous coefficients from an anisotropic medium requires an explicit projection prescription, and the O(ξ) results are used at values of ξ where the truncation error is untested.","major_comments":[{"comment":"For a fixed anisotropy axis n, the angular integrands in (26) contain c(α,θ,φ)=(n·p̂)². The tensor ∫dΩ p̂^i p̂^j p̂^k p̂^l (n·p̂)² is not proportional to the isotropic rank-4 tensor; contracted with traceless W_kl it yields different coefficients for different shear components (e.g., W_xy vs W_xz) plus a trace contribution, so no unique scalar η satisfying the isotropic relation (27) exists. The closed forms (28)–(29) are obtained by replacing c with its angular average 1/3; the same applies to ζ through (30)–(33) → (37)–(38). This implicit isotropic projection is never stated or justified. Please specify the projection used to define η and ζ (e.g., a definite contraction with W^ij, or report the full set of anisotropic coefficients η∥, η⊥ and the corresponding bulk response), and show that the closed forms follow. As written, the quantities labeled η and ζ in Figs. 5–8 are angle-average","section":"§2, Eqs. (26)–(29), (37)–(38)"},{"comment":"All analytic results are first order in ξ, but the paper plots and interprets results at ξ = 0.4–0.6, where ξ² = 0.16–0.36 is not negligible compared with unity. No truncation-error estimate is given, and the isotropic baryonless limit ξ = 0, μ = 0 is not benchmarked against existing RTA/quasiparticle calculations, despite being the baseline for all four scenarios. Since the quasiparticle masses (56)–(57) are also linear in ξ, the trends reported in Figs. 1–8 may receive significant O(ξ²) corrections. Please (i) compute η and ζ at ξ = 0.4 and 0.6 using the exact Romatschke–Strickland distributions (1)–(3) without the ξ-expansion and compare with (29) and (38); (ii) report the ξ = 0, μ = 0 values against a published limit.","section":"§2, Eqs. (4)–(6), (29), (38); §5, Figs. 1–8"},{"comment":"The relaxation times (16)–(17) are momentum-independent isotropic-equilibrium expressions; their use in an anisotropic, baryon-asymmetric plasma is assumed without discussion of possible ξ- or μ-dependence, and no sensitivity test is provided. In addition, the derivatives (22)–(24) omit terms proportional to ∂ω/∂T and ∂ω/∂μ arising from the ξ-dependent quasiparticle masses (56)–(57) that enter the distributions; such terms are known to contribute to bulk viscosity in quasiparticle models. Because the absolute values of η/s, Pr, Re, and the quantitative size of the baryon-asymmetry enhancement depend on these choices, I ask that the assumptions be stated explicitly and their impact be assessed—e.g., by including the mass-derivative terms, or by testing an anisotropic and/or momentum-dependent τ at one representative temperature.","section":"§2, Eqs. (13)–(17), (22)–(24)"}],"minor_comments":[{"comment":"The numerical evaluation uses ε, P, and the speed-of-sound squared (∂P/∂ε) that appear in (31)–(38), but the explicit expressions used for these quantities are not given. Please provide them so the results are reproducible.","section":"§4–5, Eqs. (31)–(38)"},{"comment":"The Landau–Lifshitz conditions are imposed species-by-species; only the total ΔT^00 = 0 is required by energy-momentum conservation. Please justify the per-species choice or show that it is equivalent to the total condition.","section":"§2, Eqs. (34)–(36)"},{"comment":"The definition of C_P in (42) is not the standard constant-pressure specific heat T(∂s/∂T)_P, and it is not stated whether the derivative is taken at fixed μ or fixed P. Since the Prandtl number (39) also uses κ from Ref. [10], please clarify the consistency of the two calculations.","section":"§3.1, Eq. (42)"},{"comment":"Please verify the momentum powers and coefficients in (54)–(55). Direct substitution of (4)–(6) into (52) gives a gluon ξβ term of the form −ξβg²/(9π²)∫dp p³/ω_g f_g(1+f_g), whereas Eq. (54) shows p² f_g(1+f_g); several terms in (55) appear to have the same issue. If an ultrarelativistic approximation ω ≈ p is intended, state it explicitly.","section":"§4, Eqs. (54)–(55)"},{"comment":"The notation f_f for both flavor and function makes Eqs. (4)–(6) and (22)–(24) hard to parse; a cleaner notation (e.g., f_q, f_qbar) would help. The inset panels in Figs. 5 and 7 are too small to read, and the curves in Fig. 8(b) are nearly coincident at low T, making the claimed ordering at the dip hard to verify.","section":"Throughout; Figs. 5–8"}],"recommendation":"major_revision","confidential_remarks":"The central ambiguity in the viscosity definition (major comment 1) should be resolved before acceptance. The manuscript also leans heavily on the authors' own prior work (Refs. [8]–[11], [16]–[19], [30]) for the framework and for the thermal conductivity used in Pr, so independent validation of the baseline ξ = 0, μ = 0 values would substantially strengthen the paper. If the authors choose to report angle-averaged coefficients, they should be explicit that standard viscous hydrodynamics with a single scalar η assumes isotropy, and that their η/s comparisons to the 1/4π bound have that interpretation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives the first explicit relaxation-time-approximation calculation of shear and bulk viscosities for a Romatschke-Strickland anisotropic QGP with finite baryon chemical potential, using xi-dependent quasiparticle masses. The derivation is algebraically detailed and the results are genuine model predictions, not fits. That part is solid and useful.\n\nWhat is actually new: the combination of expansion-induced anisotropy and finite baryon asymmetry for eta and zeta, plus the derived Prandtl number, Reynolds number, eta/s, and zeta/s. It extends the same group's charge/heat transport calculation (ref. [10]) to momentum transport, and it goes a step further by letting the quasiparticle masses depend on the same anisotropic distributions used for the transport integrals. No data fitting, no fabricated validation. Credit where due.\n\nThe main soft spot is the one the stress-test note flags, and I think it holds up. Equation (27) assumes the isotropic relation Delta T^ij = -eta W^ij - zeta delta^ij div u. But inserting a distribution with a fixed anisotropy direction n leaves angular integrands containing (p dot n)^2. After integrating over angles, the coefficient of W^ij is not a scalar; it has tensor pieces built from n^i n^j. The closed forms in (28)-(29) effectively replace (p dot n)^2 with its average p^2/3, an implicit isotropic projection that is never stated or justified. So the scalar eta in the figures is an angle-averaged coefficient, not the shear viscosity in the usual tensorial sense. The qualitative direction—anisotropy lowers eta and zeta, baryon asymmetry raises them—probably survives, but the quantitative values are conditional on the projection, and the paper should say which projection.\n\nOther soft spots, in proportion: the relaxation times in (16)-(17) are taken from isotropic equilibrium and applied to an anisotropic medium with xi-dependent masses; that sets the overall scale and could shift numbers, but it is a standard approximation and not fatal. The O(xi) expansion is used at xi = 0.4-0.6, where xi^2 is not tiny; an error estimate or a check at smaller xi would have been easy and reassuring. And there is no baseline comparison to known xi=0, mu=0 results, again easy to add.\n\nNone of this is load-bearing enough to desk-reject. It is a competent, honest extension that needs a definitional cleanup and a few robustness checks. Send it to peer review; a referee can ask for the projection to be made explicit and for an estimate of higher-order terms. I would cite it once those are in place.","headline":"A careful RTA calculation of shear and bulk viscosities in an anisotropic, baryon-asymmetric QGP, but the scalar eta is an implicit angular average that needs to be stated and justified before the central numbers hold up.","tokens_in":21027,"tokens_out":4284,"would_cite":true,"duration_ms":51188,"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":"Expansion-induced anisotropy reduces shear and bulk viscosity of hot QCD matter; baryon asymmetry increases both.","keywords":["quark-gluon plasma","shear viscosity","bulk viscosity","momentum anisotropy","baryon chemical potential","relaxation time approximation","Prandtl number","Reynolds number"],"falsifier":"A lattice QCD calculation of shear and bulk viscosities for a quark-gluon plasma with finite quark chemical potential and an anisotropic distribution—or a kinetic-theory computation without the small-ξ truncation—would settle whether η and ζ indeed decrease with ξ and increase with μ at the quantitative level claimed here.","tokens_in":20111,"feed_emoji":"🌀","tokens_out":5272,"duration_ms":55920,"temperature":0.7,"pith_summary":"Solving the relativistic transport equation in the relaxation-time approximation for a quark-gluon plasma with both a small momentum-space anisotropy and a finite baryon chemical potential, the authors find that anisotropy suppresses the shear viscosity η and the bulk viscosity ζ, while baryon asymmetry enhances them. The mechanism is a deformation of the parton distribution functions and a corresponding shift in the quasiparticle masses, both of which feed into the dissipative part of the energy-momentum tensor. Because η and ζ set the scale of dissipative corrections in heavy-ion phenomenology, the predicted changes in the Prandtl number, Reynolds number, and specific viscosities give concrete, temperature-dependent signatures of early-time anisotropy and net baryon density. The paper extends known isotropic and baryonless results to a more realistic early-stage fireball that is simultaneously expanding and baryon-rich.","feed_headline":"Anisotropy lowers QGP viscosity; baryon density raises it","feed_subtitle":"Solving the kinetic equation in a deformed plasma shows expansion cuts shear and bulk viscosity while baryon asymmetry compensates.","key_machinery":"The central objects are the deformed phase-space distribution functions of quarks, antiquarks, and gluons, obtained from an isotropic distribution by the substitution p → √(p² + ξ(p·n)²), normalized by √(1+ξ), and expanded to first order in the anisotropy parameter ξ. These distributions determine both the transport coefficients—through the dissipative part of the energy-momentum tensor in the relaxation-time approximation—and the quasiparticle masses of the medium, which become functions of temperature, chemical potential, and ξ. The shear and bulk viscosities are then read off as the coefficients of the traceless and trace parts of the momentum-space gradient terms.","core_discovery":"The central claim is that in a baryon-asymmetric QCD plasma, a weak expansion-induced anisotropy (parameter ξ > 0) lowers both the shear viscosity η and the bulk viscosity ζ relative to the isotropic plasma, while increasing the quark chemical potential at fixed ξ raises both coefficients. The paper derives closed-form integral expressions for η and ζ from the deformed quark, antiquark, and gluon distributions and the ξ-dependent quasiparticle masses, and shows numerically that the trends persist across temperatures from 0.16 to 0.64 GeV. The same computation yields the Prandtl number, Reynolds number, and specific viscosities, all of which show a stronger response to anisotropy than to bary","pith_inferences":["Because the relaxation times are taken from isotropic equilibrium and momentum-independent, the quantitative magnitudes of the viscosity shifts should be tested against a calculation with momentum-dependent relaxation; the qualitative sign pattern is likely robust to that change.","The same deformed-distribution machinery should also modify electrical and thermal conductivities in a baryon-asymmetric anisotropic plasma, so the pattern of transport-coefficient suppression with ξ is a testable family of predictions.","Bayesian extractions of η/s and ζ/s from heavy-ion data at different collision energies—where baryon chemical potential varies—could distinguish the anisotropic model from an isotropic baryonless baseline.","Extending the first-order-in-ξ expansion to second order at ξ ≈ 0.5 could quantify the error of the truncation and check whether the decreasing trends survive."],"forward_implications":["Hydrodynamic models of the pre-equilibrium stage should use lower η and ζ when momentum anisotropy is present, reducing dissipative corrections during the early expansion.","At finite baryon density, viscosities rise relative to baryonless matter, partially counteracting the anisotropy effect and shifting the balance with temperature.","The Prandtl number exceeding unity in all scenarios means momentum diffusion dominates sound attenuation, with anisotropy weakening that dominance and baryon asymmetry barely changing it.","The specific shear viscosity approaches the conjectured 1/(4π) bound in the baryonless anisotropic case, so early-time anisotropy makes the plasma appear more nearly perfect.","The temperature at which ζ/s dips is modified by both ξ and μ, giving a concrete signature of where the plasma is closest to conformal symmetry."],"fun_headline_variants":["Expansion anisotropy trims QGP viscosity, baryons boost it","Deformed QGP: shear and bulk viscosity drop with anisotropy","Baryon asymmetry compensates viscosity loss from anisotropy","Hot QCD: anisotropy lowers viscosity, baryon density raises it","QGP viscosity: expansion cuts, baryons restore"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The relaxation times for quarks, antiquarks, and gluons are taken from isotropic equilibrium and assumed momentum-independent, which fixes the overall scale of both viscosities; if the true relaxation times depend on anisotropy, chemical potential, or momentum, the magnitude—and possibly the ordering—of the effects would change.","fun_headline_variants_meta":{"raw":{"variants":["Expansion anisotropy trims QGP viscosity, baryons boost it","Deformed QGP: shear and bulk viscosity drop with anisotropy","Baryon asymmetry compensates viscosity loss from anisotropy","Hot QCD: anisotropy lowers viscosity, baryon density raises it","QGP viscosity: expansion cuts, baryons restore"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000151,"raw_usage":{"total_tokens":1084,"prompt_tokens":836,"completion_tokens":248,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":164}},"tokens_in":580,"tokens_out":248,"duration_ms":4006,"temperature":1.0,"reasoning_tokens":164,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T14:57:16.269377+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice QCD calculation of shear and bulk viscosities for a quark-gluon plasma with finite quark chemical potential and an anisotropic distribution—or a kinetic-theory computation without the small-ξ truncation—would settle whether η and ζ indeed decrease with ξ and increase with μ at the quantitative level claimed here.","supporting_citations":[],"review_version":1}