{"id":"87662e08-22cd-4ba4-895c-b024f058e5b8","arxiv_id":"2411.13833","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Excited hybrid quark-antiquark channels are color-octet at short distance and exponentially relax to random color, with screening masses ordered B(Σ_g^+) ~ B(Π_u) < B(Π_u') ~ B(Δ_g) < B(Δ_g') < B(Σ_g^+').","lead":"This lattice QCD study measures how the color charge of a static quark and antiquark leaks into the gluon field as the pair is separated, in the ground state and in five excited hybrid channels. It finds that excited pairs start color-octet at short distance and all channels relax to a random color mixture at large separation, with channel-dependent screening speeds.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The hierarchy's largest screening mass B(Σ_g^+') is extracted from an ad hoc constant offset δ, so the headline ordering is not yet supported.","rationale":"The Reader's verdict is CONDITIONAL, and I agree with that verdict, but for a somewhat different reason. The weakest-assumption cited by the Reader — that the density matrix is diagonal with equal octet entries — is actually not a free assumption for any physical state with total color charge zero. A color-singlet full state gives a reduced density matrix invariant under global color rotations; by Schur's lemma on 3⊗3bar = 1⊕8, the matrix must be αP_s+βP_o, i.e., diagonal in the {singlet, octet} basis with all eight octet probabilities equal. So lack of a direct off-diagonal test for excited channels is not, by itself, load-bearing. The real weak point is more specific: the top of the claimed B hierarchy, B(Σ_g^+')=0.489(21), is extracted after manually adding an R-independent offset δ=0.2043 to F(R;Σ_g^+'). This offset is not present in the ansatz and has no quoted uncertainty; it has the signature of excited-state contamination from the Σ_g^+ ground state. A fit with δ free, or a two-state fit, could shift B(Σ_g^') substantially, potentially erasing the separation from B(Δ_g^')=0.418(1) on which the ordering B(Δ_g')<B(Σ_g^') depends. Additionally, all B values come from one lattice spacing and one volume, with no continuum or infinite-volume extrapolation; the paper explicitly defers this. The qualitative finding — singlet at R→0 for Σ_g^+, octet at R→0 for hybrid channels, all approaching 1:8 random mixture — is well supported by the plotted diagonal data and consistent with the constituent-gluon picture. The quantitative hierarchy, however, is not yet solid. I would keep the verdict CONDITIONAL. The single check that would most directly settle my concern is a refit of Σ_g^+' with δ as a free parameter and full error propagation.","tokens_in":13744,"tokens_out":16564,"duration_ms":174261,"concrete_test":"Refit the Σ_g^+' channel two ways with full covariance: (i) F(R)=A exp(-BR) over 4≤R≤7 without the offset; (ii) F(R)=A exp(-BR)-δ with δ a free parameter, propagating δ's uncertainty into B. Then compare B(Σ_g^+') with B(Δ_g^')=0.418(1). If the two values are no longer separated by more than 2σ, or if fit (i) gives a substantially different slope, the claimed hierarchy should be downgraded. As a complementary check, compare δ with the ground-state overlap from the GEVP; if they match, a proper two-state fit should replace the ad hoc shift.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the screening-mass hierarchy B(Σ_g^+)∼B(Π_u)<B(Π_u')∼B(Δ_g)<B(Δ_g')<B(Σ_g^+') (Eq. (29), Table I). The top of this hierarchy is the only part that rests on a quantity outside the model of Eqs. (8)-(17): for Σ_g^+' the authors find that the measured F(R) is contaminated by an R-independent offset and define \\tilde F(R)=F(R)+δ with δ=0.2043, selected by a numerical fit with no quoted uncertainty and no physical origin in the ansatz. The presence of such an offset is exactly what one expects from imperfect isolation of the excited Σ_g^+' state, which shares quantum numbers with the much lighter Σ_g^+ ground state. If δ is even moderately uncertain, the slope B(Σ_g^+')=0.489(21) changes, and the separation from B(Δ_g^')=0.418(1) may disappear. Since the ansatz's diagonal/equal-octet form is largely protected by global color symmetry, the most load-bearing weakness is not the ansatz shape but this unmodeled offset in the channel that fixes the top of the hierarchy.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the reduced color density matrix ρ(R) of a static quark-antiquark pair in quenched lattice QCD on a 32^3×32 lattice at β=5.8, for the ground-state channel Σ_g^+ and five channels with gluonic excitations (Σ_g^+', Π_u, Π_u', Δ_g, Δ_g'). From the singlet and octet components of ρ, the authors define F(R;Γ), the fraction of the initial color configuration remaining at separation R, and fit F(R;Γ)=A(Γ)exp(−B(Γ)R) to extract channel-dependent screening masses B(Γ). They report that all six channels approach the random 1:8 singlet-to-octet mixture at large R, that all excited channels start from an octet configuration at small R, and that the screening masses satisfy B(Σ_g^+)≈B(Π_u)<B(Π_u')≈B(Δ_g)<B(Δ_g')<B(Σ_g^+'). A comparison with finite-temperature results from a previous work is also presented.","tokens_in":14022,"tokens_out":6225,"duration_ms":63408,"significance":"The qualitative finding that color correlation leaks into the gluon field and flows to the random color configuration in every channel, including excited channels, is clearly supported by the lattice data and is a useful, physically appealing result. The small-R octet dominance in all excited channels is consistent with the constituent-gluon picture and with earlier hybrid-potential calculations. The construction of a reduced density matrix for excited static q-qbar systems is a natural and valuable extension of the authors' earlier work. However, the quantitative screening-mass hierarchy, which is the main new quantitative claim, rests on one channel whose extraction requires an ad hoc offset, and the quoted B values are not accompanied by fit-range or systematic uncertainties. With those caveats addressed, the paper would be a solid contribution; in its present form the quantitative claim is not yet fully supported.","major_comments":[{"comment":"The top of the claimed hierarchy in Eq. (29), namely B(Σ_g^+')=0.489(21), is obtained by fitting F(R;Σ_g^+')+δ with δ=0.2043. This offset is introduced ad hoc: no uncertainty is quoted for δ, no physical origin is given, and no stability test of B(Σ_g^+') with respect to δ is reported. Since Σ_g^+' shares quantum numbers with the much lighter Σ_g^+ ground state, a constant contamination from imperfect isolation of the excited state is a concrete and plausible mechanism for such an offset. A moderate uncertainty in δ could change B(Σ_g^+') enough to erase the separation from B(Δ_g')=0.418(1). Please provide a correlated fit that includes δ with an uncertainty, or a variational check of the operator overlap, or a systematic scan over δ, and state explicitly how the ordering in Eq. (29) changes as δ is varied within its uncertainty.","section":"Sec. III.C, Eq. (28), Table I"},{"comment":"The screening masses B(Γ) are extracted from the range 4≤R≤7 lattice units, chosen because the log plot appears linear there. No fit-range systematics are given for any channel, and the Σ_g^+' data turn upward for R>1.2 fm after the offset correction. Because the quoted B values cover only a factor of about 2.5 and several entries are close, e.g., B(Δ_g)=0.322(2) and B(Π_u')=0.338(5), a table varying R_min and R_max for every channel (and δ for Σ_g^+') is needed to establish that the hierarchy in Eq. (29) is not an artifact of the chosen fitting window.","section":"Sec. III.C, Fig. 10"}],"minor_comments":[{"comment":"The one-parameter form of the reduced density matrix in Eqs. (10)-(17) can be justified more directly than by citing the earlier fits: for a total color-singlet q-qbar-gluon state, the reduced density matrix is invariant under simultaneous color rotations of the quark indices, and hence, by Schur's lemma, it must be a linear combination of the singlet and octet projectors. This makes the diagonal form with equal octet entries exact at every R, and stating it would preempt concerns about the ansatz.","section":"Sec. II.B"},{"comment":"There are typographical errors: 'quar ks' appears in the title, 'matirices' appears in Sec. IV, and 'Fig. 11 in Ref. [6]' should be 'Fig. 11 of Ref. [6]'.","section":"Sec. IV"},{"comment":"The construction of O(n)_Γ is only shown explicitly for Γ=Σ_g^+; for the excited channels the reader is referred to Refs. [9,10]. A brief sentence or a small appendix listing the operator shapes for Π_u, Δ_g, and the primed channels would improve reproducibility.","section":"Sec. II.C"},{"comment":"The temperature in Eq. (30) is set to T=250 MeV, but the lattice data shown are at T=220, 244, and 275 MeV. The choice of 250 MeV should be justified, or the comparison should be reframed as a qualitative consistency check between two interpolated values.","section":"Sec. III.D, Fig. 12"},{"comment":"Table I lists only B(Γ); the fitted amplitudes A(Γ) and the offset δ for Σ_g^+' are not quoted. Reporting all fitted parameters with uncertainties would allow readers to reproduce the curves in Fig. 10 and to assess the quality of the exponential description.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The paper is a natural continuation of the authors' previous work, and the self-citations to Refs. [6,7] are appropriate in that context. The qualitative results are convincing and suitable for the journal; the main issue is that the headline quantitative hierarchy depends on an unmodeled offset in the Σ_g^+' channel. I believe this is fixable with a robustness analysis rather than requiring a fundamentally different approach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is the extension of your reduced-density-matrix method to five excited gluonic channels, and that part is genuinely worth having. The data clearly show the excited q-qbar pairs sitting in a color-octet configuration at short distance and relaxing toward the random 1:8 singlet-to-octet mixture at large R. That qualitative picture is solid and matches the constituent-gluon intuition. The comparison with your earlier finite-T results is a nice consistency check, and the paper is honest about the single lattice spacing and finite volume.\n\nThe soft spot is exactly the one the stress-test flags. The screening-mass hierarchy in Eq. (29) puts B(Σ_g^+') at the top, but that value comes from F(R) corrected by an additive offset δ=0.2043 that has no quoted uncertainty and no physical origin in the ansatz. You can see in Fig. 8 that the raw Σ_g^+' data do not follow the exponential form, so you modify the quantity before fitting. That makes the largest screening mass the least reliable number in the table, and if δ is even moderately uncertain the separation from B(Δ_g') probably disappears. The fix is straightforward: quote a correlated fit of F(R) = A exp(-BR) − δ with all parameters and errors, or show that the ordering survives different treatments of the offset.\n\nThe ansatz itself is a separate, lesser concern. You validated it for the ground state in Ref. [6], but for the excited channels you do not directly test whether off-diagonal color matrix elements vanish. I think the diagonal form is plausible and probably protected by the color-symmetry arguments, but because the screening masses are extracted from this single function, a direct check on the excited-state density matrix would strengthen the quantitative claims.\n\nBottom line: the paper deserves a serious referee. The qualitative conclusions are supported by the plots, and the new quantitative information is useful for hybrid meson studies. But the specific hierarchy, especially the top end, should be treated as provisional until the offset is handled more carefully. I would send this to review and ask for a more systematic treatment of the fitting systematics before publication.","headline":"A useful new lattice measurement of color screening in hybrid q-qbar channels, with a credible qualitative picture but a screening-mass hierarchy that rests on a fitted offset in the highest channel.","tokens_in":14499,"tokens_out":1578,"would_cite":true,"duration_ms":17168,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V05","81T25"],"pacs":["12.38.Gc","11.15.Ha"],"model":"deepseek-v4-flash","headline":"Hybrid quark–antiquark pairs begin as color octets and lose color correlation exponentially as they separate.","keywords":["lattice QCD","static quark-antiquark pair","color density matrix","gluonic excitation","hybrid hadron","color screening","flux tube","quenched approximation"],"falsifier":"On the same $32^3$ ensemble at $\\beta=5.8$, compute the off-diagonal color matrix element $\\rho_{s,a_i}$ for the $\\Pi_u$ channel from the un-normalized $L^{(n)}_{ij,kl}$ correlators; if any such element is nonzero beyond noise, the one-parameter diagonal ansatz fails for the excited channels and the quoted $B(\\Gamma)$ values are not the screening masses of the paper's model.","tokens_in":13551,"feed_emoji":"🎨","tokens_out":5810,"duration_ms":53144,"temperature":0.7,"pith_summary":"This paper asks how the color of a static quark–antiquark pair changes when the gluon field between them is excited, and it computes the reduced density matrix in color space from quenched lattice QCD in six channels: the ground-state $\\Sigma_g^+$ and five hybrid channels ($\\Sigma_g^+{}'$, $\\Pi_u$, $\\Pi_u'$, $\\Delta_g$, $\\Delta_g'$). The central finding is that in every hybrid channel the pair is a pure color octet at short separation, and that as the separation grows the color correlation is exponentially screened toward the completely random mixture in which singlet and octet components appear in the ratio $1:8$. The paper extracts a screening mass $B(\\Gamma)$ for each channel and finds the pattern $B(\\Sigma_g^+)\\sim B(\\Pi_u) < B(\\Pi_u')\\sim B(\\Delta_g) < B(\\Delta_g') < B(\\Sigma_g^+{}')$. A sympathetic reader should care because this gives a quantitative, gauge-invariant picture of how color leaks from quarks into the flux tube, and it supports the constituent-gluon picture of hybrid hadrons.","feed_headline":"Hybrid quark pairs: octet at short range, random at long range","feed_subtitle":"Lattice QCD shows color leaking exponentially into the gluon field, with screening masses ordering all six channels.","key_machinery":"The load-bearing object is the reduced two-body color density matrix $\\rho_{ij,kl}(R)$ for a static $q\\bar q$ pair, obtained by integrating out gluon degrees of freedom, together with the one-parameter ansatz $\\rho^{\\rm ansatz}(R)=F(R)\\hat\\rho_{\\rm initial}+(1-F(R))\\hat\\rho_{\\rm rand}$, where $\\hat\\rho_{\\rm rand}=\\frac{1}{N_c^2}\\hat I$ is the maximally mixed color state. A single fraction $F(R;\\Gamma)$ carries all the separation dependence, and the screening mass $B(\\Gamma)$ is extracted from the exponential slope of $F$. Numerically, the machinery is quenched lattice QCD in Coulomb gauge at $\\beta=5.8$ on a $32^3\\times 32$ lattice, with operators $O_\\Gamma^{(n)}(R,T)$ carrying the channel quantum numbers and a generalized eigenvalue problem used to isolate excited-state signals.","core_discovery":"The paper extends the color-density-matrix ansatz of its earlier work to excited $q\\bar q$ channels. In the ansatz, $\\rho(R)$ is a convex combination of an initial color configuration and the random configuration $\\frac{1}{N_c^2}\\hat I$; for the ground state the initial configuration is color singlet, while for each hybrid channel it is color octet. The lattice data show that at $R\\to0$ the excited channels are indeed pure octet, consistent with a constituent gluon bound to an octet $q\\bar q$ pair. As $R$ grows, all channels lose correlation: $F(R;\\Gamma)$, the fraction of the initial state, decays as $A(\\Gamma)\\exp(-B(\\Gamma)R)$, with fitted values in lattice units such as $B(\\Sigma_g^+)=0.200(5)$ and $B(\\Pi_u)=0.211(2)$, and with the ordering $B(\\Sigma_g^+)\\sim B(\\Pi_u) < B(\\Pi_u')\\sim B(\\Delta_g) < B(\\Delta_g') < B(\\Sigma_g^+{}')$. The $\\Sigma_g^+{}'$ channel requires an additional constant offset $\\delta=0.2043$ before the exponential falloff appears. In addition, a thermal density matrix built from the six zero-temperature channels reproduces the finite-temperature lattice data for $R<0.8$ fm.","pith_inferences":["If the ansatz survives a direct test, the screening mass $B(\\Gamma)$ could be connected quantitatively to gluelump masses, giving a spectroscopic interpretation of why higher channels screen faster; the paper only raises this as a possibility.","A direct lattice computation of off-diagonal color matrix elements would determine where the one-parameter ansatz breaks; the current paper validates the diagonal form only for the ground-state channel.","The constant offset $\\delta$ in the $\\Sigma_g^+{}'$ channel could be either a finite-volume artifact or a sign of a second, nearly degenerate state; a volume scan would distinguish these.","Continuum and dynamical-fermion extrapolations would show whether the hierarchy of screening masses persists beyond the quenched, single-lattice-spacing setting studied here."],"forward_implications":["Hybrid $q\\bar q$ systems are color octets at short range, so their constituent-gluon interpretation is supported by direct color-space measurement.","Color screening by the flux tube is not special to the singlet ground state: every excited channel also flows to the random $1:8$ singlet-to-octet mixture.","The close screening masses of $\\Sigma_g^+$ and $\\Pi_u$ suggest the $\\Pi_u$ mode is a simple fundamental gluonic excitation that does not accelerate color leak, while higher excitations randomize quark color faster.","The zero-temperature six-channel density matrix predicts finite-temperature color correlation up to $R\\approx0.8$ fm, so low-lying $T=0$ physics accounts for the sub-$T_c$ behavior.","The same reduced-density-matrix construction transfers to multiquark systems, which the authors identify as the next target."],"supporting_citations":[{"why":"Supplies the reduced-density-matrix ansatz and the ground-state color-screening analysis that this paper extends to excited channels.","marker":"[6]"},{"why":"Provides the finite-temperature color-density-matrix data that the thermal reconstruction in Sec. III D must match.","marker":"[7]"},{"why":"Supplies the hybrid $q\\bar q$ potential calculation and the operator construction whose quantum numbers label the six channels.","marker":"[9]"},{"why":"Provides excited-state operators used to construct $O_\\Gamma^{(n)}$ in the lattice formalism.","marker":"[10]"},{"why":"Supplies the generalized eigenvalue method used to isolate excited-state signals.","marker":"[19]"},{"why":"Provides the companion lattice technique for extracting excited-state overlaps used in the operator optimization.","marker":"[20]"}],"fun_headline_variants":["Color leaks from excited quark pairs into gluon field","Screening mass orders all six gluonic-excitation channels","Color density matrices reveal exponential screening in hybrids","Excited quark pairs: color leaks into gluon field","Screening mass quantifies color leak in six hybrid channels"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the quark and antiquark colors are described at every separation by a single number, the fraction of color singlet or color octet, so that all mixing between different color states and all differences among the eight octet possibilities are assumed absent; the paper tests this for the ground state but not directly for the excited channels.","fun_headline_variants_meta":{"raw":{"variants":["Color leaks from excited quark pairs into gluon field","Screening mass orders all six gluonic-excitation channels","Color density matrices reveal exponential screening in hybrids","Excited quark pairs: color leaks into gluon field","Screening mass quantifies color leak in six hybrid channels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000848,"raw_usage":{"total_tokens":3762,"prompt_tokens":1093,"completion_tokens":2669,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":2592}},"tokens_in":709,"tokens_out":2669,"duration_ms":17935,"temperature":1.0,"reasoning_tokens":2592,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:49:04.867843+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On the same $32^3$ ensemble at $\\beta=5.8$, compute the off-diagonal color matrix element $\\rho_{s,a_i}$ for the $\\Pi_u$ channel from the un-normalized $L^{(n)}_{ij,kl}$ correlators; if any such element is nonzero beyond noise, the one-parameter diagonal ansatz fails for the excited channels and the quoted $B(\\Gamma)$ values are not the screening masses of the paper's model.","supporting_citations":[{"cited_title":"Lattice QCD study of static quark and antiquark correlations via entanglement entropies","cited_arxiv_id":"1910.00859","evidence_quote":"Supplies the reduced-density-matrix ansatz and the ground-state color-screening analysis that this paper extends to excited channels."},{"cited_title":"Lattice QCD study of static quark and antiquark correlations at finite $T$ via entanglement entropies","cited_arxiv_id":"2011.10950","evidence_quote":"Provides the finite-temperature color-density-matrix data that the thermal reconstruction in Sec. III D must match."},{"cited_title":"Perantonis and C","cited_arxiv_id":null,"evidence_quote":"Provides the companion lattice technique for extracting excited-state overlaps used in the operator optimization."}],"review_version":1}