Pith. sign in

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 →

arxiv 2512.21668 v2 pith:RJSIQRNM submitted 2025-12-25 hep-lat

Lattice QCD study of color correlations between quarks in static multiquark systems

classification hep-lat PACS 12.38.Gc
keywords color correlationreduced density matrixstatic quarksmultiquark systemsflux tubelattice QCDcolor screeningentanglement entropy
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper aims to show that in static multiquark systems, the color correlation between any pair of quarks is controlled by the length of the confining flux tube that connects them, rather than the straight-line interquark distance. Using quenched lattice QCD, the authors compute a two-body reduced color density matrix for 3Q and 4Q ground states and find that the residual rate of the maximally correlated color configuration decays with flux-tube path length exactly as it does in the known 2Q case. At large path length, the color configuration approaches a fully random mixture, indicating that color leaks into the gluon field. This gives a simple geometric rule for color screening in complex multiquark states and offers a practical lattice diagnostic for distinguishing genuine tetraquark states from two-meson configurations.

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.

Watch this falsifier — get emailed when new claim-graph text bears on it.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 3 minor

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)
  1. [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.
  2. [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.
  3. [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).
  4. [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)
  1. [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'.
  2. [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.
  3. [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

0 steps flagged

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

2 free parameters · 6 axioms · 0 invented entities

The central quantitative claim rests on a one-parameter ansatz imported from the authors' previous work, a model-dependent definition of the flux-tube path length, and a single quenched lattice ensemble. No new particles or forces are introduced. The main free parameters are the extracted residual rate F and the geometric coordinate L.

free parameters (2)
  • Residual rate F_MC(R) (F_3bar or F_c) = extracted per configuration from lattice rho via Eq. (21) or (34)
    One-parameter interpolation weight in the ansatz rho = F rho_MC + (1-F) rho_rand; it is measured from data, not predicted from first principles.
  • Flux-tube path length L = computed per (d,h) from assumed Y/double-Y/X junction geometry
    Coordinate used for the universality claim; depends on unmeasured junction positions and template flux-tube shapes.
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.
    Eqs. (6), (13), (25); imported from authors' prior 2Q papers [10-12]; restricts rho to a one-parameter family.
  • domain assumption The static-light correlator at time T is dominated by the ground state, so <C~> proportional to e^{-E0 T} rho.
    Eq. (44); no multi-exponential or variational analysis is shown; excited-state contamination could bias extracted F.
  • ad hoc to paper Flux tubes in ground-state 3Q/4Q systems have Y-type/double-Y (X-type) geometry with minimal-length junctions.
    Sec. III B, IV C, Figs. 3 and 10; L is defined along these assumed tubes.
  • domain assumption Coulomb gauge fixing gives a meaningful color basis for the density matrix.
    Sec. II D; rho components are gauge-dependent.
  • 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.
    Sec. II D; one lattice spacing, no sea quarks; authors acknowledge finite-volume/systematic effects.
  • standard math SU(3) color decomposition into 3bar/6 (QQ) and 1/8 (Q Qbar) with maximally mixed random state.
    Eqs. (7)-(12), (22)-(24); group-theoretic input.

pith-pipeline@v1.3.0-alltime-deepseek · 20937 in / 12185 out tokens · 119064 ms · 2026-08-03T14:01:08.476112+00:00 · methodology

0 comments
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

Figures reproduced from arXiv: 2512.21668 by Toru T. Takahashi, Yoshiko Kanada-En'yo.

Figure 1
Figure 1. Figure 1: FIG. 1. The paths Γ [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Three quarks are located at (+ [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Antitriplet and sextet components ( [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Shematic pictures which express the definition of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Quarks and antiquarks are located at (+ [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. In connected 4Q configurations (left), all the quarks [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Schematic pictures which express the definition of [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The data with errorbars labeled “E ( [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: shows the L dependence of the antitriplet components (ρ3¯) in the Q1Q2 pair (the left panel of [PITH_FULL_IMAGE:figures/full_fig_p008_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Singlet components ( [PITH_FULL_IMAGE:figures/full_fig_p009_12.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. The residual rate [PITH_FULL_IMAGE:figures/full_fig_p009_14.png] view at source ↗
Figure 17
Figure 17. Figure 17: shows the singlet components (ρ {1,4} 1 ) mea￾sured in a Q1Q¯ 4 pair in planar 4Q systems plotted as a function of the flux-tube path length L. The black dashed line indicates the random limit value, ρ {1,4} 1 = 1/9. The red dotted line labeled “MC[2Q+2Q]” shows ρ1 = 1/9, the value expected in the MC limit of a disconnected 4Q configuration. In disconnected 4Q configurations, Q1 and Q¯ 4 belong to differe… view at source ↗
Figure 16
Figure 16. Figure 16: shows the singlet components (ρ {1,3} 1 ) mea￾sured in a Q1Q¯ 3 pair in planar 4Q systems plotted as a function of the flux-tube path length L. The black dashed line indicates the random limit value, ρ {1,3} 1 = 1/9. The red dotted line labeled “MC[2Q+2Q]” shows ρ1 = 1, the value expected in the MC limit of a disconnected 4Q con￾figuration. In disconnected 4Q configurations, Q1 and Q¯ 3 belong to the iden… view at source ↗
Figure 18
Figure 18. Figure 18: The quark pair of which we investigate the cor [PITH_FULL_IMAGE:figures/full_fig_p012_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: FIG. 19 [PITH_FULL_IMAGE:figures/full_fig_p012_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: FIG. 20 [PITH_FULL_IMAGE:figures/full_fig_p012_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: FIG. 21 [PITH_FULL_IMAGE:figures/full_fig_p013_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: FIG. 22 [PITH_FULL_IMAGE:figures/full_fig_p013_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: FIG. 23 [PITH_FULL_IMAGE:figures/full_fig_p014_23.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

21 extracted references · 1 canonical work pages

  1. [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. [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. [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. [4]

    G. S. Bali, K. Schilling and C. Schlichter, Phys. Rev. D 51, 5165 (1995) doi:10.1103/PhysRevD.51.5165 [hep- lat/9409005]

  5. [5]

    V. G. Bornyakov et al. [DIK Collaboration], Phys. Rev. D 70, 054506 (2004) doi:10.1103/PhysRevD.70.054506 [hep-lat/0401026]

  6. [6]

    Cardoso, M

    N. Cardoso, M. Cardoso and P. Bicudo, Phys. Rev. D 84, 054508 (2011) doi:10.1103/PhysRevD.84.054508 [arXiv:1107.1355 [hep-lat]]

  7. [7]

    Tiktopoulos, Phys

    G. Tiktopoulos, Phys. Lett. 66B, 271 (1977). doi:10.1016/0370-2693(77)90878-4 15

  8. [8]

    Greensite and C

    J. Greensite and C. B. Thorn, JHEP 0202, 014 (2002) doi:10.1088/1126-6708/2002/02/014 [hep-ph/0112326]

  9. [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]]

  10. [10]

    Okiharu, H

    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]]

  11. [11]

    Okiharu, H

    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]]

  12. [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]]

  13. [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]]

  14. [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]]

  15. [15]

    Calabrese and J

    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]]

  16. [16]

    Ryu and T

    S. Ryu and T. Takayanagi, JHEP 08, 045 (2006) doi:10.1088/1126-6708/2006/08/045 [arXiv:hep- th/0605073 [hep-th]]

  17. [17]

    Kitaev and J

    A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006) doi:10.1103/PhysRevLett.96.110404 [arXiv:hep- th/0510092 [hep-th]]

  18. [18]

    Amico, R

    L. Amico, R. Fazio, A. Osterloh and V. Ve- dral, Rev. Mod. Phys. 80, 517-576 (2008) doi:10.1103/RevModPhys.80.517 [arXiv:quant- ph/0703044 [quant-ph]]

  19. [19]

    S. N. Solodukhin, Living Rev. Rel. 14, 8 (2011) doi:10.12942/lrr-2011-8 [arXiv:1104.3712 [hep-th]]

  20. [20]

    Donnelly, Class

    W. Donnelly, Class. Quant. Grav. 31, no.21, 214003 (2014) doi:10.1088/0264-9381/31/21/214003 [arXiv:1406.7304 [hep-th]]

  21. [21]

    Amorosso, S

    R. Amorosso, S. Syritsyn and R. Venugopalan, JHEP 12, 177 (2024) doi:10.1007/JHEP12(2024)177 [arXiv:2410.00112 [hep-lat]]