REVIEW 4 major objections 3 minor 21 references
The color correlation between two static quarks is governed by the flux-tube path length between them, and the screening curve is universal across 2Q, 3Q, and 4Q systems.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 14:01 UTC pith:RJSIQRNM
load-bearing objection First lattice computation of two-body color density matrices in static 3Q/4Q systems, with an interesting but not fully proven universal flux-tube-length dependence. the 4 major comments →
Lattice QCD study of color correlations between quarks in static multiquark systems
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central discovery is that the two-body color density matrix of any static quark pair in 2Q, 3Q, and 4Q ground states is well described by the ansatz rho(L) = F(L) rho_MC + (1 - F(L)) rho_random, where L is the minimal flux-tube path length between the pair (Y-type for 3Q, double-Y or X-type for 4Q), rho_MC is the maximally correlated color configuration expected at zero size, and rho_random is the uniform color mixture. The residual rate F(L) extracted from lattice data - F_1 for 2Q, F_3bar for QQ pairs in 3Q/4Q, and F_c for Q-Qbar pairs in 4Q - all collapse onto the same curve as a function of flux-tube length, establishing a universal color-screening law along the gluonic flux tube. Th
What carries the argument
The key object is the reduced two-body color density matrix rho, extracted from lattice correlators with color measurements inserted at the midpoint of the temporal Wilson lines, and projected onto the ground state through the standard Euclidean-time construction. The analysis is carried by an ansatz decomposing rho into a maximally correlated (MC) part and a random part, with F(L) as the residual rate. The flux-tube path length L is defined from assumed Y-type (3Q) and double-Y/X-type (4Q) geometries, and this coordinate converts a complicated color-entanglement problem into a single monotone function of L.
Load-bearing premise
The entire scheme rests on the assumption that the ground-state flux tube has the assumed Y-type or double-Y/X-type geometry with sharp junctions, so that a well-defined path length L exists; if the actual flux-tube profile differs from these templates, the observed universal collapse of F onto a single L-curve could be an artifact of the coordinate definition rather than a physical regularity.
What would settle it
A concrete test: vary the quark geometry while keeping L (the length along the assumed flux tube) fixed, e.g., compare 3Q configurations with different (d,h) pairs that yield the same L; a significant spread in the extracted F values at equal L would falsify the geometric claim. Alternatively, extract the decay rate of F(L) and compare it with the independently measured string tension - a mismatch would indicate that the screening length is not set by the confinement scale.
If this is right
- Color correlations in any static multiquark ground state can be predicted from the flux-tube path length alone, using the universal curve, without a separate dynamical calculation.
- The random-color limit at large L provides a quantitative criterion: a quark pair with no color correlation signals that the pair belongs to different singlet clusters (mesons), enabling lattice identification of genuine tetraquarks versus two-meson states.
- The ansatz rho = F rho_MC + (1-F) rho_random appears to be generally valid for ground-state color density matrices in the confined phase, so the same F(L) analysis can extend to other NQ configurations.
- Entanglement entropy constructed from the color density matrix acts as a sharp indicator of flip-flop: twisted 4Q systems show near-maximal entropy even at small size, while planar systems show entropy consistent with two-meson structure.
Where Pith is reading between the lines
- If the universal L-curve is genuine, a natural testable extension is that the same curve controls quark-antiquark correlations in excited or hybrid multiquark states, with the MC part replaced by the excited-state color structure, as already seen in 2Q hybrid systems.
- The universality suggests that the flux-tube path length, not the direct distance, is the correct coordinate for color screening, which could be tied to the Wilson-loop area law; comparing the fitted decay rate of F(L) with the independently measured string tension would test this geometric interpretation directly.
- The random-vs-connected signature could be inverted as a diagnostic for experiments: measuring color correlations in a candidate tetraquark state and comparing them with the lattice curves may reveal whether it is a genuine multiquark or a meson molecule.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper computes two-body reduced color density matrices for static 3Q (QQQ) and 4Q (Q Qbar Q Qbar) systems using quenched lattice QCD at beta = 5.8 on a 32^3 x 32 lattice. The authors analyze these matrices with the ansatz rho = F rho_MC + (1-F) rho_rand, where rho_MC is a maximally correlated color state and rho_rand is the random color state. They find that the residual rate F extracted from the data depends only on the flux-tube path length L, defined from assumed Y-type (3Q) and double-Y/X-type (4Q) geometries, and that F(L) follows a common curve for 2Q, 3Q, and 4Q ground-state systems. The paper also uses the second Rényi entropy of the reduced density matrix as a diagnostic of flip-flop in planar 4Q configurations.
Significance. If robust, the central claim is conceptually interesting: color correlations between static quarks are controlled by the length of the gluonic flux-tube path between them rather than by their direct spatial separation. This would provide a practical way to diagnose the internal color structure of multiquark systems, e.g., distinguishing connected tetraquark configurations from two-meson states. The paper also contains a useful cross-system comparison and a non-tautological reconstruction check in planar 4Q systems using the independently measured 2Q F_1. However, the universality claim rests on an assumed flux-tube template and a one-parameter ansatz, and the current evidence does not yet rule out coordinate-choice or systematic artifacts.
major comments (4)
- [Sec. III.B, IV.C; Eqs. (36), (42)] The flux-tube path length L is defined from the same Y-type and double-Y/X-type junction templates used to construct the interpolating operators in Eqs. (36) and (42). The paper does not test whether the observed universal collapse of F onto a single L-curve survives a change of this coordinate, e.g., using a different junction position, a straight-line distance, or a path extracted from the energy-density profile. The evidence in Fig. 14 that d=4, small-h points deviate from the universal F_1 curve and are attributed to the X-type profile shows that the L definition is not innocent. This is load-bearing for the universality claim.
- [Sec. II.C, Eq. (44)] The extraction of the ground-state reduced density matrix via Eq. (44) assumes that \(\tilde C_{ab,cd}^{NQ}\) couples only to the ground state. No T-plateau analysis is shown for F or for the density-matrix elements, and the temporal separation T used is not stated. Without plateaus or an explicit demonstration that excited-state contamination is negligible, the extracted rho and hence F may not represent the ground state. The paper also does not state the number of gauge configurations, and several key figures (e.g., Figs. 4, 5, 11–17) have no visible error bars, making the statistical significance of the claimed collapse difficult to assess.
- [Sec. II.B, Eqs. (13)–(18)] The ansatz predicts a specific matrix structure: zero off-diagonal elements, equal diagonal components within the antitriplet block, and equal components within the sextet block. The paper only shows the traces rho_3bar and rho_6, and does not validate the full matrix structure. Consequently, the extracted F is a projection onto the ansatz rather than a test that the ansatz actually describes the data. To support the claim that the color configuration is represented by the ansatz, the off-diagonal suppression and the degeneracies should be shown (or an explicit goodness-of-fit test provided).
- [Sec. II.D] All results are obtained at a single lattice spacing (beta = 5.8, a = 0.14 fm) in the quenched approximation. The paper acknowledges finite-volume effects but does not address discretization effects or sea-quark effects. For a quantitative claim of universality of F(L) across 2Q, 3Q, and 4Q systems, at least a second lattice spacing or an estimate of systematic uncertainties is needed. As written, the claim is based on one ensemble at one cutoff.
minor comments (3)
- [Abstract and Secs. IV, VI] Typos: 'universarity' should be 'universality'; 'protucts' in Sec. II.C; 'st atic' in the header; 'Renyi' should be 'Rényi'.
- [Figs. 15–17] The red dotted line labeled 'MC[2Q+2Q]' takes the same value (3/9 or 1/9) as the random limit in Figs. 15 and 17. The text explains the distinction, but the figure captions should state explicitly that the MC value coincides with the random value in these cases to avoid confusion.
- [Sec. II.B] The ansatz in Eq. (6) is imported from the authors' previous work. This is not a defect, but the paper should state explicitly that the cross-system comparison of F(L) is the novel test, not the ansatz form itself.
Circularity Check
No significant circularity: the universal F(L) collapse is an empirical observation with cross-system checks, and the self-cited ansatz is tested against direct lattice data rather than assumed as the result.
full rationale
The central claim—that the residual rates F1^{2Q}, Fbar3^{3Q}, Fbar3^{4Q}, and Fc^{4Q} collapse onto a common curve as a function of the flux-tube path length L—is an observed, not constructed, regularity. L is defined geometrically from assumed Y-type or double-Y/X-type flux-tube templates (Sec. III B and Sec. IV C), not extracted or fitted from the F data. The paper even reports deviations from the universal curve (Fig. 14, attributed to X-type flux-tube profiles), which would not occur if L were being tuned to force agreement. The reduced density matrix is computed directly from lattice correlators via Eq. (44), and the ansatz of Eq. (6) is used as an interpretive decomposition after the fact; its validity is tested by the single-valuedness of the color components and by agreement with the independently measured 2Q F1^{2Q}. The planar 4Q reconstruction using F1^{2Q} from Ref. [10] (Figs. 15–17) is a genuine non-circular cross-check: the 2Q input is independent of the 4Q correlators, and the comparison succeeds only in the connected-4Q region. The self-citations to Refs. [10–12] supply the ansatz and prior 2Q data, but they are not load-bearing in the sense of forcing the multiquark result; the multiquark universality is demonstrated by the lattice data themselves. Remaining concerns—the lack of an explicit T-plateau test, finite-volume effects, and the sensitivity of L to the assumed flux-tube geometry—are systematic-uncertainty or model-choice issues, not circular reductions of the derivation to its inputs. No equation in the paper reduces to another by construction, and no fitted parameter is renamed as a prediction.
Axiom & Free-Parameter Ledger
free parameters (2)
- Residual rate F_MC(R) (F_3bar or F_c) =
extracted per configuration from lattice rho via Eq. (21) or (34)
- Flux-tube path length L =
computed per (d,h) from assumed Y/double-Y/X junction geometry
axioms (6)
- ad hoc to paper Reduced color density matrix has the ansatz form rho = F rho_MC + (1-F) rho_rand for any quark pair.
- domain assumption The static-light correlator at time T is dominated by the ground state, so <C~> proportional to e^{-E0 T} rho.
- ad hoc to paper Flux tubes in ground-state 3Q/4Q systems have Y-type/double-Y (X-type) geometry with minimal-length junctions.
- domain assumption Coulomb gauge fixing gives a meaningful color basis for the density matrix.
- domain assumption Quenched Wilson lattice QCD at beta=5.8 (a=0.14 fm) with L=32^3 approximates confined QCD enough for the conclusions.
- standard math SU(3) color decomposition into 3bar/6 (QQ) and 1/8 (Q Qbar) with maximally mixed random state.
read the original abstract
We study the color correlation between two static quarks in 3Q ($QQQ$) and 4Q ($QQ\bar Q\bar Q$) multiquark systems at $T=0$ based on the reduced two-body density matrices $\rho$ in color space. We perform quenched lattice QCD calculations with the Coulomb gauge adopting the standard Wilson gauge action, and the spatial volume is $L^3 = 32^3$ at $\beta = 5.8$, which corresponds to the lattice spacing $a=0.14$ fm and the system volume $L^3=4.5^3$ fm$^3$. We evaluate the two-body color density matrix $\rho$ of static quarks, and investigate the dependence of color correlations on the quarks' spatial configuration. As a result, we find that the color correlations depend on the minimal path length along a flux tube which connects two quarks under consideration. The color correlation between quarks quenches because of color leak into the gluon field (flux tube) and finally approaches the random color configuration in the large distance limit. We also find a ``universality'' in the flux-tube path length dependence of the color correlation for 2Q, 3Q, and 4Q ground-state systems. Our results show that the color correlations of end-point quarks can be a clue to clarify the internal structures of hadrons, including exotic (multiquark) hadrons.
Figures
Reference graph
Works this paper leans on
-
[1]
A color configuration of any QQ pair is 3 decomposed into color antitriplet and sextet states
QQ color correlation in 3Q and 4Q systems First, we consider the case when we measure the color correlation between two quarks ( QQ) in multiquark (3Q and 4Q) systems. A color configuration of any QQ pair is 3 decomposed into color antitriplet and sextet states. Let ˆρ¯3i and ˆρ6i be the density operators for color antitriplet states |¯3i⟩ and sextet state...
-
[2]
initial” MC color configuration of a Q ¯Q pair would be neither |1⟩ nor |8i⟩, but a combination of them. Especially when four quarks form a “genuine
Q ¯Q color correlation in 4Q systems The color configuration of a Q ¯Q pair in 4Q ( QQ ¯Q ¯Q) systems is classified into singlet and octet color configura- tions, |1⟩ and |8i⟩. Then we define the density operators for a pair forming a color singlet state |1⟩ and an octet state |8i⟩ in the Coulomb gauge as ˆρ1 = |1⟩⟨1|, ˆρ8i = |8i⟩⟨8i| (i = 1 ∼ N 2 c − 1). (22...
-
[3]
An Introduction to the Confinement Problem,
J. Greensite, “An Introduction to the Confinement Problem,” Lecture Notes in Physics (Springer, 2011) doi:10.1007/978-3-642-14382-3, and references therein
-
[4]
G. S. Bali, K. Schilling and C. Schlichter, Phys. Rev. D 51, 5165 (1995) doi:10.1103/PhysRevD.51.5165 [hep- lat/9409005]
arXiv 1995
-
[5]
V. G. Bornyakov et al. [DIK Collaboration], Phys. Rev. D 70, 054506 (2004) doi:10.1103/PhysRevD.70.054506 [hep-lat/0401026]
Pith/arXiv arXiv 2004
-
[6]
N. Cardoso, M. Cardoso and P. Bicudo, Phys. Rev. D 84, 054508 (2011) doi:10.1103/PhysRevD.84.054508 [arXiv:1107.1355 [hep-lat]]
Pith/arXiv arXiv 2011
-
[7]
G. Tiktopoulos, Phys. Lett. 66B, 271 (1977). doi:10.1016/0370-2693(77)90878-4 15
-
[8]
J. Greensite and C. B. Thorn, JHEP 0202, 014 (2002) doi:10.1088/1126-6708/2002/02/014 [hep-ph/0112326]
Pith/arXiv arXiv 2002
-
[9]
T. T. Takahashi, H. Suganuma, Y. Nemoto and H. Matsufuru, Phys. Rev. D 65, 114509 (2002) doi:10.1103/PhysRevD.65.114509 [arXiv:hep- lat/0204011 [hep-lat]]
arXiv 2002
-
[10]
F. Okiharu, H. Suganuma and T. T. Taka- hashi, Phys. Rev. D 72, 014505 (2005) doi:10.1103/PhysRevD.72.014505 [arXiv:hep- lat/0412012 [hep-lat]]
arXiv 2005
-
[11]
F. Okiharu, H. Suganuma and T. T. Taka- hashi, Phys. Rev. Lett. 94, 192001 (2005) doi:10.1103/PhysRevLett.94.192001 [arXiv:hep- lat/0407001 [hep-lat]]
arXiv 2005
-
[12]
T. T. Takahashi and Y. Kanada-En’yo, Phys. Rev. D 100, no.11, 114502 (2019) doi:10.1103/PhysRevD.100.114502 [arXiv:1910.00859 [hep-lat]]
Pith/arXiv arXiv 2019
-
[13]
T. T. Takahashi and Y. Kanada-En’yo, Phys. Rev. D 103, no.3, 034504 (2021) doi:10.1103/PhysRevD.103.034504 [arXiv:2011.10950 [hep-lat]]
Pith/arXiv arXiv 2021
-
[14]
T. T. Takahashi and Y. Kanada-En’yo, Phys. Rev. D 111, no.1, 014505 (2025) doi:10.1103/PhysRevD.111.014505 [arXiv:2411.13833 [hep-lat]]
Pith/arXiv arXiv 2025
-
[15]
P. Calabrese and J. L. Cardy, J. Stat. Mech. 0406, P06002 (2004) doi:10.1088/1742-5468/2004/06/P06002 [arXiv:hep-th/0405152 [hep-th]]
Pith/arXiv arXiv 2004
- [16]
-
[17]
A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006) doi:10.1103/PhysRevLett.96.110404 [arXiv:hep- th/0510092 [hep-th]]
arXiv 2006
- [18]
-
[19]
S. N. Solodukhin, Living Rev. Rel. 14, 8 (2011) doi:10.12942/lrr-2011-8 [arXiv:1104.3712 [hep-th]]
Pith/arXiv arXiv 2011
-
[20]
W. Donnelly, Class. Quant. Grav. 31, no.21, 214003 (2014) doi:10.1088/0264-9381/31/21/214003 [arXiv:1406.7304 [hep-th]]
Pith/arXiv arXiv 2014
-
[21]
R. Amorosso, S. Syritsyn and R. Venugopalan, JHEP 12, 177 (2024) doi:10.1007/JHEP12(2024)177 [arXiv:2410.00112 [hep-lat]]
Pith/arXiv arXiv 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.