{"id":"24a08a3f-f2fc-4b12-9190-76dedda02ccb","arxiv_id":"2512.21668","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Color correlations of quark pairs in static 2Q, 3Q, and 4Q systems all decay along the same universal curve when plotted against the assumed flux-tube path length, approaching random color at large separation.","lead":"This lattice QCD paper measures how the color of two quarks inside static three- and four-quark systems becomes shared with the gluon field as the quarks separate. It reports that the color correlation between any quark pair is controlled by the path length along an assumed flux tube, with a common curve for two-, three-, and four-quark systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flux-tube length L is defined from the same Y/double-Y templates used in the interpolating operators; no test rules out an operator-geometry artifact behind the universal F(L) collapse.","rationale":"The reader identified the definition of L and excited-state contamination as the weakest assumptions. I agree that the unvalidated definition of L is the most load-bearing issue, but I sharpen it as a potential operator-geometry circularity: because L is taken along the same junction paths used in the interpolating operators, a universal curve in L could be baked into the construction. This is not an accusation of error; it is a concrete gap in the evidence. The paper offers some independent support (cross-system comparisons, reconstruction with 2Q data), but it never varies the definition of L or the operator. The deviations in Fig. 14, attributed to X-type geometry, show that the L template is already imperfect. A simple same-L consistency test using the existing data could falsify or support the universality; an alternative-operator test would cleanly separate physical flux-tube dependence from operator dependence. Since these tests are needed before the universality claim is accepted, the reader's CONDITIONAL verdict is appropriate; my concern does not move it.","tokens_in":21297,"tokens_out":9010,"duration_ms":108859,"concrete_test":"Re-analyze the 3Q and twisted-4Q data: group all configurations by integer L; for each L bin, test whether F values from different d,h (and different pairs) agree within jackknife errors (χ²/dof). Additionally, for two representative geometries (one 3Q, one twisted 4Q), recompute the density matrix with an alternative junctionless (Δ-link) baryon/tetraquark operator at several T and check that the extracted F(L) is unchanged and T-independent. If same-L points scatter, or the Δ-link result shifts, the claimed universality is an artifact of the operator/coordinate choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central universality claim is controlled by the coordinate L, which is defined from the same Y/double-Y junction templates used to build the creating operators in Eqs. (36) and (42). The paper does not test whether the collapse of F onto a single curve survives a change of this coordinate. For 3Q, L is the sum of the two arms from each measured quark to the assumed junction; for 4Q twisted, it is the double-Y/X path. Because the epsilon-tensor operators distribute color along exactly those paths, the reduced density matrix of the pair is kinematically tied to the lengths of those same arms; a monotone dependence of F on L may be in part a property of the interpolating operator rather than of the ground-state flux tube. Evidence that this is not merely notational is visible in Fig. 14, where d=4, small-h points deviate from the universal F1 curve; the authors attribute this to the X-type profile, i.e., to a mismatch between the assumed L template and the actual tube. Thus the 'universality' is only demonstrated for a particular, unvalidated choice of L, and the paper never compares alternative definitions (e.g., straight-line distance, actual energy-density path, or a different junction point). The absence of a T-plateau analysis further weakens the extraction of the ground-state density matrix via Eq. (44).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":21652,"tokens_out":3266,"duration_ms":38156,"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":[{"comment":"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.","section":"Sec. III.B, IV.C; Eqs. (36), (42)"},{"comment":"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.","section":"Sec. II.C, Eq. (44)"},{"comment":"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).","section":"Sec. II.B, Eqs. (13)–(18)"},{"comment":"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.","section":"Sec. II.D"}],"minor_comments":[{"comment":"Typos: 'universarity' should be 'universality'; 'protucts' in Sec. II.C; 'st atic' in the header; 'Renyi' should be 'Rényi'.","section":"Abstract and Secs. IV, VI"},{"comment":"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.","section":"Figs. 15–17"},{"comment":"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.","section":"Sec. II.B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a natural continuation of the authors' prior work, and the cross-system comparison is the main new contribution. The largest risk is that the universal F(L) collapse is an artifact of defining L from the same Y/double-Y templates used in the operators; this needs a dedicated robustness test. The absence of configuration counts, T-plateau checks, and error bars on many figures is also a barrier to assessing the strength of the claim. I recommend major revision rather than rejection, because the central idea is interesting and the planar 4Q reconstruction using the independently measured F_1 provides a non-tautological anchor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a real lattice calculation of a new observable, and the cross-system comparison is genuinely new. The universal flux-tube-length claim is plausible but not nailed down; the main worry is that L is defined from the same Y/double-Y junction templates that build the interpolating operators.\n\nThe paper does several things well. The formalism is laid out clearly, the planar 4Q analysis — reconstructing rho-bar-3 and rho-1 using the independently measured 2Q F1 — is a non-circular check that strengthens the ansatz, and the entanglement-entropy section gives a useful diagnostic for flip-flop. Credit where due: new lattice data, new cross-system comparison, and an honest statement of the systematic limitations they do mention.\n\nThe soft spots are real but not fatal. First, L is not coordinate-independent in a way the paper tests: the same Y/double-Y paths appear in the epsilon-tensor operators and in the definition of L. The authors never try alternative definitions (straight-line distance, energy-density path, different junction point). The deviations in Fig. 14 for large d and small h, attributed to the X-type profile, look exactly like a symptom of that coordinate mismatch. Second, ground-state dominance is assumed via Eq. (44), but no T-plateau analysis is shown; that is straightforward to check. Third, the calculation is quenched, at a single lattice spacing, with no stated number of gauge configurations, and many key figures have no visible error bars. That makes the \"universality\" look more precise than the data support. Fourth, the ansatz is imported from the authors' own 2Q work, so the one-parameter mixture is not an independent prediction — but the cross-system comparison mitigates the circularity.\n\nBottom line: a serious referee should engage with this. The central physical idea — that color screening is controlled by flux-tube path length rather than direct distance — is worth testing. Right now it is demonstrated for one choice of L at one quenched spacing. I would send it to peer review and ask for error bars, a T-plateau check, and either a test of alternative L definitions or a softened universality claim.\n\nWho is it for: lattice QCD and hadron-spectroscopy people working on multiquark states and flux tubes. I would cite it if I worked in that area.","headline":"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.","tokens_in":22126,"tokens_out":2533,"would_cite":true,"duration_ms":28425,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Gc"],"model":"deepseek-v4-flash","headline":"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.","keywords":["color correlation","reduced density matrix","static quarks","multiquark systems","flux tube","lattice QCD","color screening","entanglement entropy"],"falsifier":"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.","tokens_in":21152,"feed_emoji":"⚛️","tokens_out":3735,"duration_ms":35709,"temperature":0.7,"pith_summary":"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.","feed_headline":"Quark color screening follows one universal flux-tube curve","feed_subtitle":"For 2Q, 3Q, and 4Q ground states, the same length along the gluon flux tube, not quark distance, sets the color-loss curve.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Color loss in quarks scales with flux-tube length, not distance","Universal color-screening law emerges for multiquark systems","Quark pairs in 2Q, 3Q, 4Q share one color-screening curve","Flux-tube path, not distance, sets quark color screening","Universal color screening: flux-tube length is the key"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Color loss in quarks scales with flux-tube length, not distance","Universal color-screening law emerges for multiquark systems","Quark pairs in 2Q, 3Q, 4Q share one color-screening curve","Flux-tube path, not distance, sets quark color screening","Universal color screening: flux-tube length is the key"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000537,"raw_usage":{"total_tokens":2476,"prompt_tokens":867,"completion_tokens":1609,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":1527}},"tokens_in":611,"tokens_out":1609,"duration_ms":9920,"temperature":1.0,"reasoning_tokens":1527,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T14:01:08.476112+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}