{"id":"7db72c30-2144-43b7-ac83-ba38e7a9efc5","arxiv_id":"2504.20138","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a matrix model of the semi-QGP with teen ghost fields, zeta/s is largest at T_d, comparable to eta/s, while eta/s stays well above the AdS/CFT bound.","lead":"This paper computes the shear and bulk viscosities of a pure glue plasma using a matrix model with nontrivial holonomy and new 'teen' ghost fields. The model predicts that the bulk viscosity to entropy ratio peaks near the deconfinement temperature at values comparable to the shear viscosity, then drops sharply.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"At the central parameter point the leading-log quasiparticle expansion is not controlled: α_sN_c(T_d)≈1.5 and ln(κ/g^2N_c)≈1.2, so the claimed ζ/s peak lies outside the calculation's regime of validity.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing premise that fails most directly: the coupling is assumed moderate near T_d, justifying on-shell quasiparticles and leading-log 2→2 kinetic theory. The quantitative check underlines this: at the parameter values used for the central figures, α_sN_c(T_d)≈1.5 and the leading logarithm is barely larger than unity. Thus the claimed ζ/s peak is not a controlled weak-coupling prediction. No independent fatal inconsistency was found beyond the limitations the paper already acknowledges, so the correct response is to keep the CONDITIONAL verdict rather than reject the work outright. The bulk-viscosity result also depends on taking Δv_s^2 from lattice thermodynamics (Appendix B), but that is a modeling choice, not the chief weakness; the coupling hierarchy is the more fundamental condition. My review therefore does not change the reader's verdict but strengthens the reason for conditionality by showing that the central parameter choice numerically violates the stated assumption.","tokens_in":43908,"tokens_out":10412,"duration_ms":119552,"concrete_test":"At T=T_d, w=0.5, κ=64, and with the holonomy solution of Sec. III, compute the imaginary part of the retarded hard-gluon self-energy at on-shell momentum p∼T using the total HTL self-energy of Eq. (58), including the teen contribution of Eq. (61). If the resulting on-shell damping width Γ is an O(1) fraction of T, the on-shell quasiparticle Boltzmann equation used in Secs. V–VI is internally inconsistent; if Γ/T≲0.2, the moderate-coupling premise survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central prediction—ζ/s largest at T_d and comparable to η/s, with η/s≈0.38–0.4 and ζ/s≈1.6–1.7 reported near T_d—is obtained from a leading-logarithmic Boltzmann equation for on-shell gluon and teen quasiparticles. That framework requires weak coupling and a hierarchy of scales. The paper's own running coupling, Eq. (165), with the values used for the headline figures (w=0.5, κ=64), gives α_sN_c = 6π/[11 ln(w2πT/T_d)] ≈ 1.5 at T_d, i.e. 't Hooft coupling λ≈18.8; even for w=1, α_sN_c≈0.93. Consequently ln(κ/g^2N_c)=ln(64/λ)≈1.2 at T_d, so the logarithm in Eq. (138) is not large and the leading-log truncation is not controlled. Equivalently, with N_c=3 the 'soft' scale gT exceeds T at T_d, so the HTL separation P∼gT≪T used in Sec. IV is absent. The statements in Sec. I that the coupling is 'moderate' are inconsistent with this numerical choice. Appendix B's acknowledged discrepancy between kinetic-theory teen pressure and thermodynamic pressure, and Appendix C's zero-mode regularization, are additional internal limitations, but the coupling-size problem alone means the calculated values and the claim that ζ/s peaks at T_d are not established by the calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends a large-N_c effective matrix model of the SU(N_c) holonomy potential by adding adjoint two-dimensional ghost fields, called 'teens,' whose parameters are adjusted to reproduce the lattice Polyakov loop. The authors then compute the shear and bulk viscosities from a linearized Boltzmann equation with 2-to-2 scattering at leading logarithmic order in the presence of nonzero holonomy. The central reported result is that η/s remains above the AdS/CFT bound while ζ/s is largest at T_d, with values near ζ/s ≈ 1.6–1.7, and falls off rapidly by about 2T_d.","tokens_in":44154,"tokens_out":8821,"duration_ms":89346,"significance":"If the calculation were controlled, it would provide a concrete transport framework for the semi-quark-gluon plasma and a falsifiable prediction that the bulk viscosity peaks at deconfinement in pure-glue theory. The manuscript has genuine strengths: the matrix elements are written out explicitly, the sign structure for ghost fields is treated carefully in Appendix A, the cancellation of the δΓ terms through the equations of motion is a nontrivial consistency check around Eq. (59), and the zero-mode problem in the bulk-viscosity collision kernel is addressed numerically in Appendix C. However, the significance is conditional on several heuristic ingredients, and one of them—the validity of the weak-coupling leading-log expansion at T_d—is load-bearing for the main quantitative claim.","major_comments":[{"comment":"The leading-logarithmic expansion is not controlled at the temperatures where the headline prediction is made. With the values used for the figures, w = 0.5 and κ = 64, Eq. (165) gives g^2 N_c = 24π^2/[11 ln(π)] ≈ 18.8 at T = T_d, i.e., α_s N_c ≈ 1.5. The 'soft' momentum scale is then gT > T, so the HTL separation P ~ gT ≪ T used in Sec. IV is absent. Moreover, ln(κ/g^2 N_c) = ln(64/18.8) ≈ 1.2, so the logarithm that controls the truncation in Eq. (138) is not large and the leading-log approximation is not justified. The quoted value ζ/s ≈ 1.6–1.7 at T_d, and the associated claim that ζ/s peaks at T_d, are therefore not established by the calculation as presented.","section":"Sec. V.C, Eq. (138) and Sec. VII, Eq. (196)"},{"comment":"The bulk viscosity is largely inherited from lattice thermodynamics rather than derived from teen quasiparticle dynamics. The source term in Eq. (168) is proportional to Δv_s^2 = 1/3 − v_s^2, and v_s^2 is computed from the lattice pressure fit in Eq. (27). Consequently, the near-T_d enhancement of ζ/s is controlled by the input lattice equation of state. The manuscript should quantify how much of the peak survives if Δv_s^2 is instead computed self-consistently from the matrix model, and should avoid presenting the ζ/s peak as an independent prediction of the teen kinetic theory.","section":"Sec. VI, Eqs. (168), (27), and (196)"},{"comment":"The kinetic-theory and thermodynamic descriptions of the teen fields are not mutually consistent. Appendix B shows that kinetic theory gives E_kt = −T^2 T_d^2/24 and P_kt = −T^2 T_d^2/72, hence E_kt = 3P_kt, while direct thermodynamics gives E_th = P_th = −T^2 T_d^2/24. The bulk-viscosity calculation uses the thermodynamic Δv_s^2 while using teen quasiparticles with kinetic-theory dispersion in the collision kernel. This acknowledged mismatch means the computation is not a closed derivation from a single kinetic model; its quantitative output should be labeled accordingly.","section":"Appendix B and Sec. VI"},{"comment":"The teen-field action in Eq. (33) contains a transverse momentum cutoff at T_d that explicitly violates gauge invariance. The paper argues that the teen self-energy does not enter at leading logarithmic order, but the teen-gluon and teen-teen scattering amplitudes used in Secs. V and VI do rely on the teen propagator with this gauge-violating cutoff. A more explicit demonstration is needed that the gauge-violating terms cannot affect the leading-log results for η/s and ζ/s; otherwise the transport coefficients also rest on an unquantified model assumption.","section":"Sec. IV.A and Sec. V.A"}],"minor_comments":[{"comment":"The text near Fig. 15 states that the parameters are chosen with w = 1 and κ = 64, while the figure caption and Eq. (196) use w = 0.5 and κ = 64; these must be reconciled.","section":"Sec. VI.B, Fig. 15"},{"comment":"The running coupling is written in an ambiguous form; please display g^2 N_c explicitly and state the renormalization scheme and scale used.","section":"Eq. (165)"},{"comment":"The abstract emphasizes the large-N_c limit, but several comparisons are made to lattice data for N_c = 3 and N_c = 5; the large-N_c status of the final results should be stated more precisely.","section":"Abstract and Sec. III"},{"comment":"The use of the Bengali character for the teen field is likely to be confusing to readers; a table of symbols or an alternative notation would improve readability.","section":"Sec. IV.A"},{"comment":"The matrix model does not reproduce the sharp peak in the interaction measure near 1.1T_d, and this limitation should be explicitly connected to the reliability of ζ/s in that temperature region.","section":"Sec. VI.B and Fig. 3"}],"recommendation":"reject","confidential_remarks":"The manuscript is transparent about its assumptions and the computation is careful, but the central quantitative claim—ζ/s ≈ 1.6–1.7 and comparable to η/s at T_d—lies outside the regime where the leading-logarithmic weak-coupling expansion is controlled. This is not a presentation issue; it is a load-bearing validity problem that would require either a next-to-leading-order calculation or a fundamental reframing of the results as uncontrolled estimates. I therefore do not see a path to acceptance within the manuscript's current scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a long, careful leading-log transport calculation in the semi-QGP matrix model, and the first to include the \"teen\" ghost fields in both shear and bulk viscosity. The main qualitative claim — ζ/s rises toward T_d and becomes comparable to η/s there — is new and, within the model, stable across the two eigen-density ansatze. The paper is also honest: it flags the gauge-violating transverse cutoff for teens, the mismatch between kinetic-theory teen pressure and thermodynamics, and the slow approach of the model Polyakov loop to unity.\n\nThe calculation itself is a genuine feat. The matrix elements, color sums, and the cancellation of the δΓ terms via the equations of motion are worked out explicitly, and the paper gives enough detail to be checked. That is worth referee time.\n\nBut the stress-test note is right, and it lands on the load-bearing part. With the headline values w=0.5 and κ=64, Eq. (165) gives g^2 N_c ≈ 18.8 at T_d, i.e. α_s N_c ≈ 1.5 and ln(κ/g^2 N_c) ≈ 1.2. There is no large log, and the soft scale gT is not small compared to T. The introduction calls the coupling \"moderate\", but the numbers used in the figures are not moderate; they are strongly coupled. So the leading-log Boltzmann equation is operating outside its regime of validity, and the reported ζ/s ≈ 1.6–1.7 at T_d is not established by the calculation. This is a criticism of the numeric parameter choice, not of the algebra.\n\nThe other soft spots are softer. The teen fields are ad hoc, but the paper acknowledges it and shows that the problematic self-energy drops out at leading log. The kinetic-theory pressure issue is bypassed by taking Δv_s^2 from lattice thermodynamics; that is a phenomenological patch, but an explicit one. The Polyakov-loop fit is heuristic, and the results depend on κ and w, which the paper discloses.\n\nWho gets value: readers in the semi-QGP program and those comparing effective models of transport near T_d. As a quantitative prediction for ζ/s, I would not rely on it. As a worked-out, partially controlled extension of the matrix model program, it deserves a serious referee. My recommendation: send to peer review, with a request that the authors either choose parameters where the leading-log expansion is actually controlled or reframe the claim as a model extrapolation and discuss the coupling-size problem explicitly.","headline":"Careful new transport calculation in the semi-QGP matrix model with teen ghosts, but the headline ζ/s peak sits in a regime where the leading-log expansion is not controlled (α_s N_c ≈ 1.5 at T_d).","tokens_in":44796,"tokens_out":3065,"would_cite":false,"duration_ms":31137,"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 claims that in a pure-glue matrix model with teen ghost fields, the bulk viscosity-to-entropy ratio peaks at the deconfinement transition at values comparable to shear viscosity, then falls rapidly to negligible by about twice…","keywords":["shear viscosity","bulk viscosity","semi-quark-gluon plasma","Polyakov loop","matrix model","teen ghost fields","leading logarithmic order","pure glue QCD"],"falsifier":"A next-to-leading-logarithmic computation of the same matrix-model Boltzmann equation would settle the leading-log prediction, but the sharpest test is a lattice extraction of the pure-glue $\\zeta/s$ near $1.1T_d$: if the ratio remains below about $0.3$ instead of rising toward the model's $1.6$ to $1.7$, the central claim would be contradicted.","tokens_in":43568,"feed_emoji":"⚛️","tokens_out":6500,"duration_ms":65580,"temperature":0.7,"pith_summary":"This paper tries to establish that a pure-glue plasma made of gluons and effective two-dimensional adjoint ghost fields, called teens, has transport properties shaped by partial deconfinement. Using Boltzmann kinetic theory at leading logarithmic order in a background with nonzero Polyakov-loop holonomy, it claims that the shear viscosity to entropy ratio $\\eta/s$ decreases toward the deconfinement temperature $T_d$ but stays above the conjectured $1/4\\pi$ bound. In contrast, the bulk viscosity to entropy ratio $\\zeta/s$ is largest at $T_d$, comparable to $\\eta/s$, then falls off rapidly and is negligible by about $2T_d$. A sympathetic reader would care because a large $\\zeta/s$ near the transition would affect hydrodynamic descriptions of heavy-ion collisions and would be a distinctive, testable signature of the semi-quark-gluon-plasma picture.","feed_headline":"Bulk viscosity peaks at deconfinement, dies by 2Td","feed_subtitle":"A matrix model with teen ghost fields finds ζ/s comparable to η/s at Td and negligible soon after.","key_machinery":"The central object is the large-$N_c$ matrix model of deconfinement, whose eigenvalue density has the form $\\rho(q)=1+b(T)\\cos(d(T)q)$ and determines the Polyakov loops $\\ell_k$ as well as the free energy. The load-bearing new ingredient is the teen field, defined as a two-dimensional adjoint ghost whose negative pressure cancels the leading gluon pressure near $T_d$; teen propagators and vertices enter the scattering amplitudes on the same footing as gluons. Transport is computed from the linearized Boltzmann equation at leading logarithmic order, keeping only $2\\to 2$ forward scattering with a soft gluon exchanged at momentum $k\\sim gT$, which produces the logarithm $\\ln(\\kappa/g^2 N_c)$ after angular integration. For the bulk viscosity, the source term is proportional to $2p^2(1/3-v_s^2)$, so the deviation of the speed of sound from the conformal value drives the result.","core_discovery":"The central claim is that, in an SU($N_c$) pure glue theory described by an effective matrix model for the thermal holonomy, the shear and bulk viscosities can be computed to leading logarithmic order in weak coupling using gluons and teen ghost fields as quasiparticles. The paper finds that $\\eta/s$ is moderately suppressed toward $T_d$ but remains well above the conjectured lower bound of $1/4\\pi$, whereas $\\zeta/s$ is largest at $T_d$, around $1.6$ to $1.7$ for the central parameter choices, comparable to $\\eta/s$, and then drops sharply so that bulk dissipation is negligible by roughly $2T_d$. The origin of the large bulk viscosity is the nonconformal equation of state, measured by $\\Delta v_s^2 = 1/3 - v_s^2$, rather than thermal masses, so the result differs strongly from ordinary perturbative estimates where the bulk viscosity is strongly suppressed.","pith_inferences":["If a large $\\zeta/s$ near $T_d$ is real, hydrodynamic simulations of deconfined matter that currently employ small bulk viscosity near the transition may need revision; a testable signature would be enhanced entropy production or modified radial flow in the near-transition region.","The same machinery, applied to QCD with dynamical quarks where the Polyakov loop at the chiral transition is much smaller, should produce a substantially smaller $\\eta/s$; computing $\\zeta/s$ in that chiral semi-QGP is the natural next step.","The quoted values at $T_d$ span roughly $\\zeta/s \\approx 1.1$ to $7$ when the parameters $\\kappa$ and $w$ are varied within the ranges considered, so a next-to-leading-logarithmic computation is the cleanest way to sharpen the prediction.","Because the model's bulk viscosity is driven by the equation of state rather than by scattering dynamics, lattice determinations of the pure-glue speed of sound near $T_d$ directly constrain the magnitude of the predicted $\\zeta/s$."],"forward_implications":["At $T_d$, $\\zeta/s \\approx 1.6$ to $1.7$ for the central parameters, comparable to $\\eta/s$; near $1.5T_d$, $\\eta/s \\approx 0.38$ to $0.4$ and $\\zeta/s \\approx 0.076$ to $0.087$, with the shear value at the high end of current lattice estimates.","$\\zeta/s$ increases strongly as $T\\to T_d$ and becomes negligible by about $2T_d$, so bulk dissipation is a near-transition phenomenon in this model.","The inclusion of gluon-teen and teen-teen scattering raises the shear viscosity near $T_d$ by about six percent relative to gluon-gluon scattering alone, because the ghost fields act as additional scatterers.","The temperature-dependent ratio $\\zeta/\\eta$ divided by $(1/3-v_s^2)^2$ lies around $56$ to $65$ at $T_d$, larger than the relaxation-time-approximation rule of thumb, and rises to a maximum near $2.5T_d$ before falling.","Even with large choices of the logarithmic parameter $\\kappa$, the model keeps $\\eta/s$ above the conjectured $1/4\\pi$ bound throughout the deconfined phase."],"supporting_citations":[{"why":"Defines the teen ghost fields whose negative pressure cancels the gluon pressure near the transition, the key model ingredient being tested.","marker":"[44]"},{"why":"Supplies the large-$N_c$ matrix-model solution for the eigenvalue density and the vanishing of higher Polyakov loops at $T_d$.","marker":"[33]"},{"why":"Provides the leading-logarithmic Boltzmann collision framework used to compute the shear viscosity.","marker":"[45, 46]"},{"why":"Supplies the bulk-viscosity method and the treatment of zero modes in the collision kernel.","marker":"[47]"},{"why":"Gives the nonzero-holonomy hard thermal loops and Debye mass expressions used in the gluon and teen self-energies.","marker":"[28]"},{"why":"Provides the lattice SU(3) pressure used to compute the entropy density and the conformality measure $\\Delta v_s^2$.","marker":"[58]"},{"why":"Supplies lattice renormalized Polyakov loop values used to fit the temperature dependence of the matrix-model solution.","marker":"[54]"}],"fun_headline_variants":["Teen ghosts make bulk viscosity spike at deconfinement","Bulk viscosity rivals shear at Td, then dies by 2Td","Matrix model finds bulk viscosity peak at Td, quick drop","Pure glue: bulk viscosity peaks at deconfinement, then fades"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the gauge coupling being moderate even at $T_d$, so that gluons and teen fields are well-defined quasiparticles whose leading-logarithmic $2\\to 2$ scattering controls transport; if the effective coupling is large near $T_d$, the Boltzmann computation and the predicted large $\\zeta/s$ do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Teen ghosts make bulk viscosity spike at deconfinement","Bulk viscosity rivals shear at Td, then dies by 2Td","Matrix model finds bulk viscosity peak at Td, quick drop","Pure glue: bulk viscosity peaks at deconfinement, then fades"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000536,"raw_usage":{"total_tokens":2589,"prompt_tokens":972,"completion_tokens":1617,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":1539}},"tokens_in":588,"tokens_out":1617,"duration_ms":11869,"temperature":1.0,"reasoning_tokens":1539,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:36:24.144724+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A next-to-leading-logarithmic computation of the same matrix-model Boltzmann equation would settle the leading-log prediction, but the sharpest test is a lattice extraction of the pure-glue $\\zeta/s$ near $1.1T_d$: if the ratio remains below about $0.3$ instead of rising toward the model's $1.6$ to $1.7$, the central claim would be contradicted.","supporting_citations":[{"cited_title":"Matrix model for deconfinement in an SU(2) gauge theory in 2+1 dimensions","cited_arxiv_id":"1306.2943","evidence_quote":"Defines the teen ghost fields whose negative pressure cancels the gluon pressure near the transition, the key model ingredient being tested."},{"cited_title":"Zero interface tension at the deconfining phase transition for a matrix model of a $SU(\\infty)$ gauge theory","cited_arxiv_id":"1301.7432","evidence_quote":"Supplies the bulk-viscosity method and the treatment of zero modes in the collision kernel."}],"review_version":1}