{"id":"fc6447c3-9afd-4503-9714-0326e2568db9","arxiv_id":"2412.10767","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Restricting maximal-center-gauge maximization to a Gaussian-distributed subset of random gauge copies reproduces the SU(2) string tension, with ad hoc corrections needed at beta=2.7.","lead":"This paper studies a gauge-fixing method used to find magnetic vortex structures in lattice QCD, and proposes a modified rule for choosing among many equivalent gauge copies. The proposed rule makes the computed quark-antiquark string tension match established lattice results, supporting the center vortex picture of quark confinement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim is not yet a prediction: the Gaussian restriction, the beta=2.6/2.7 exclusions, and the choice N=2-3 are all calibrated against the literature string tensions, and no statistical errors are reported to show the agreement is significant.","rationale":"The strongest claim is explicitly about prediction. A prediction requires an out-of-sample check. The paper provides none: the same data used to define the selection rules is used to demonstrate agreement. The reader's weakest assumption focused on whether non-Gaussian deviations at beta=2.6/2.7 are physical; my concern is broader and includes that assumption, because the symmetrization and exclusion rules are part of the selection that is tuned. I do not question the existence of the near-linear R_MCG-sigma_cp correlation; the figures are consistent with earlier work (Refs. [15,19]). The issue is evidentiary: 'prediction' is too strong a claim without a fixed protocol and uncertainty quantification. This is consistent with the reader's CONDITIONAL verdict, which asks for error bars and justified selection criteria. I therefore do not change the verdict. It is worth crediting the paper for a clear review and for making the data distributions visible, but the central numerical conclusion remains conditional.","tokens_in":7678,"tokens_out":7302,"duration_ms":64290,"concrete_test":"Hold out one coupling in the claimed range, e.g. beta=2.55 or a new lattice volume at beta=2.6, and prespecify the complete protocol: symmetrize if and only if a predefined skewness/kurtosis threshold is exceeded, exclude copies only if a predefined V(R) reversal criterion is met, take the N=3 highest R_MCG copies in the Gaussian subset, and extract sigma_cp from the large-R slope. Compute sigma_cp with a jackknife over the 200 physical configurations, and compare with the independent literature value. If the prespecified result agrees within 1-2 sigma, the central claim is supported; if it requires further tuning, the claim is a fit, not a prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's conclusion states that maximizing the MCG functional within the Gaussian-distributed subset of R_MCG values provides very good predictions of the string tension for 2.3 <= beta <= 2.6. For this to be a prediction, the selection rule must be fixed before comparing to the target values, and the agreement must be statistically meaningful. Neither condition is met. First, the procedure contains three post hoc choices: (i) for beta=2.6 the ensemble is symmetrized and for beta=2.7 gauge copies are excluded using a criterion derived from the observed reversal of the projected potential (Section 3, Fig. 4); (ii) the maximum is taken over the N highest R_MCG copies within the Gaussian subset, and N is set to 2-3 because Fig. 5 shows agreement with literature at those values, while N=1 underestimates and N>=7 overestimates; (iii) the Gaussian-restricted subset itself is defined after inspecting the distributions in Fig. 2. Second, no error bars are given for sigma_cp, so the deviations from the literature lines in Fig. 5 could be within statistical noise; the 200 gauge copies per physical configuration are not independent, so the effective sample size is much smaller than 20000. Since sigma_cp varies approximately linearly with R_MCG, and N selects the average R_MCG of the chosen copies, there will generally exist some N that reproduces any target string tension lying between the N=1 and large-N limits. The fact that N=2-3 works for all beta is suggestive of a real correlation, but without a prespecified rule and error bars it does not establish the predictive claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reviews the dual-superconductor and center-vortex pictures of confinement and presents new lattice SU(2) simulations in maximal center gauge (MCG) on 24^4 lattices at β = 2.3, 2.4, 2.5, 2.6, and 2.7. For each β the authors generate 200 physical configurations, with 100 random gauge copies each, and analyze the distribution of the MCG functional R_MCG at local maxima. They report that these distributions are approximately Gaussian for β ≤ 2.5, that the center-projected string tension σ_cp varies approximately linearly with R_MCG, and that averaging the N = 2 or N = 3 gauge copies with the highest R_MCG values within the Gaussian subset reproduces literature values of the unprojected string tension. The conclusion states that maximization of Eq. (1) restricted to the Gaussian-distributed copies gives very good predictions of the string tension for 2.3 ≤ β ≤ 2.6, with a modified exclusion also working at β = 2.7.","tokens_in":7950,"tokens_out":4405,"duration_ms":41054,"significance":"If the restricted-maximization rule were fixed before comparison and supported by statistical uncertainties, this would be a substantial step toward resolving the Gribov-copy problem in center projection and would strengthen the center-vortex mechanism as a quantitative description of confinement. The paper's main positive contribution is the documentation of the R_MCG distribution across β and the observed nearly linear relation between σ_cp and R_MCG, which is a useful empirical constraint on MCG-based vortex detection. However, the central predictive claim is not yet established: the prescription is calibrated to literature string tensions after inspecting the data, and the reported agreement is not accompanied by error bars or an independent validation. The manuscript is therefore better characterized as a suggestive empirical study than as a demonstrated prediction.","major_comments":[{"comment":"No statistical uncertainties are reported for σ_cp, and the effective independent sample size is much smaller than the 20,000 gauge fields per β, since these arise from only 200 physical configurations (20 starts × 10 separations) with 100 correlated gauge copies each. Without error bars, the apparent agreement between the center-projected and literature string tensions in Fig. 5 cannot be assessed; deviations such as the variations visible for β = 2.6 could be statistically insignificant. The claim of 'very good predictions' requires a quantitative comparison with uncertainties.","section":"Sec. 3, Figs. 3–5"},{"comment":"The choice N = 2–3 is made after inspecting the comparison with literature values in Fig. 5, not from an a priori rule. Because σ_cp varies approximately linearly with R_MCG and N selects the average R_MCG of the N highest copies, the procedure can reproduce any target string tension lying between the N = 1 and large-N limits by tuning N. The observed agreement at N = 2–3 is therefore calibration, not a prediction, unless the selection rule is fixed before comparison and tested on independent data.","section":"Sec. 3, Fig. 5 and Conclusion"},{"comment":"The Gaussian-restriction prescription is defined post hoc. For β = 2.6 the ensembles are symmetrized following Ref. [19], and for β = 2.7 gauge copies are excluded using a criterion based on the reversal of the projected potential; both steps are introduced only after non-Gaussian deviations are seen in Fig. 2. The assertion that these deviations are artifacts of incorrect vortex detection is unsupported, and no independent diagnostic distinguishes 'incorrect vortex detection' from physical finite-lattice or weak-coupling effects. Thus the claimed range 2.3 ≤ β ≤ 2.6 does not uniformly satisfy the Gaussian condition, since β = 2.6 already shows systematic deviations before symmetrization.","section":"Sec. 3, Figs. 2 and 4"},{"comment":"The literature string tensions used as reference lines [8,16–18] are taken from different lattice actions, volumes, and extrapolation procedures. The manuscript does not discuss how these differences affect the comparison, nor does it compute unprojected string tensions on the same 24^4 ensembles used for the center-projected values. A direct same-ensemble comparison would remove an uncontrolled systematic and make the claimed agreement meaningful.","section":"Sec. 3, Fig. 3"}],"minor_comments":[{"comment":"The caption contains a typo: 'β =∈ {2.3, 2.4, 2.5}' should read 'β ∈ {2.3, 2.4, 2.5}'.","section":"Fig. 3 caption"},{"comment":"The text says 'colorfull' where 'colorful' is intended, and the phrase 'precautious linearity' is unclear; a more standard wording such as 'approximate linearity' would be preferable.","section":"Sec. 2"},{"comment":"The caption contains a stray period and an incomplete sentence after 'gauge copies.'; it should be cleaned up.","section":"Fig. 5 caption"},{"comment":"Ref. [19] lists the arXiv number 1908.09711, which is the same as Ref. [2] (GlueX); this is likely a citation error and should be corrected.","section":"References"},{"comment":"The skewness panel would be more informative with error bars or at least a statement of how skewness uncertainties were estimated, since the claim of 'nearly symmetric' for β = 2.3–2.5 is central to the analysis.","section":"Sec. 3, Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's title and framing suggest a broad review, but the substantive content is a new numerical study with a specific protocol. Given the calibrational issues and missing error bars, I cannot recommend acceptance in the present form. The path to acceptance would be a substantial revision that supplies uncertainties, fixes the N-selection rule a priori, and validates the Gaussian-restriction prescription on independent configurations or at least with a clear out-of-sample test."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new increment here is small. The Gaussian-distribution analysis of the R_MCG ensemble and the restricted-maximization recipe already appear in the authors' Refs. [15] and [19], and this preprint largely follows those papers. What it adds is the beta=2.7 exclusion heuristic and a review framing of the confinement mechanisms. That is not nothing, but it is not a new result on its own. The paper's own abstract and conclusion oversell the novelty by presenting the restricted-Gaussian recipe as the main finding without clearly flagging how much of it is prior work. The citation pattern is otherwise fine; self-citation is justified when the earlier results are the actual basis.  The numerical core does have real merit. The distribution of 20000 local maxima per beta is a solid, reproducible observation, and the nearly linear dependence of sigma_cp on R_MCG is visible in the figures. The symmetrization for beta=2.6 is a reasonable response to the first signs of skewness, and the potential curves show the claimed kink. So the paper is not sloppy or incoherent; the basic lattice work appears competently done.  The soft spots are exactly where the stress-test note lands. The central claim is that maximizing within the Gaussian subset gives 'very good predictions' of the string tension. But the selection rules are post hoc: the Gaussian subset is defined after inspecting the distributions, the beta=2.6 symmetrization and the beta=2.7 exclusion are motivated by the very deviation from Gaussianity that the procedure is meant to treat, and the N=2-3 choice is justified by the agreement with literature at those values. With sigma_cp varying approximately linearly with R_MCG, there will generally exist some N that hits a target value between the N=1 and large-N limits. The fact that N=2-3 works across all beta is suggestive of a real correlation, but without error bars on sigma_cp it is impossible to tell whether the agreement is significant or just within noise. The 200 configurations are not independent, so quoting 20000 copies overstates the sample size.  Who is this for? Lattice practitioners working on center vortices and Gribov copies. They will get a useful recipe and a clear statement of one empirically motivated prescription, but they should not treat it as a confirmed prediction until the selection criteria are fixed in advance and uncertainties are provided.  Recommendation: this deserves a serious referee, but it should come back with major revision. The authors need to (1) report statistical errors on all extracted string tensions, (2) specify the Gaussian-restriction and N-selection rules before comparing to literature, or clearly separate calibration from validation, and (3) state plainly what is new relative to Refs. [15] and [19]. I would not cite it as a standalone result until those issues are addressed.","headline":"A real but incremental lattice observation wrapped in a review, whose central 'prediction' is weakened by post hoc calibration and missing error bars.","tokens_in":8561,"tokens_out":1096,"would_cite":false,"duration_ms":12016,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.-t","12.38.Aw","12.38.Gc"],"model":"deepseek-v4-flash","headline":"Random gauge copies, restricted to their Gaussian distribution, reproduce the QCD string tension on the lattice.","keywords":["color confinement","center vortices","maximal center gauge","Gribov copies","string tension","lattice QCD","monopoles","SU(2) gauge theory"],"falsifier":"Repeat the analysis on a larger lattice or with a different lattice action and check whether the Gaussian-restricted prescription reproduces the unprojected string tension without any symmetrization or copy exclusion; if agreement at $\\beta=2.7$ requires removing the non-Gaussian tails, the claim that the random-copy ensemble contains the full string-tension information would be selection-dependent rather than measured.","tokens_in":7367,"feed_emoji":"🌀","tokens_out":9396,"duration_ms":76002,"temperature":0.7,"pith_summary":"This paper argues that the confining string tension can be read from center-projected Wilson loops after a modified maximal center gauge (MCG) prescription: instead of searching for the global maximum of the gauge functional, one should restrict maximization to the Gaussian-distributed subset of random gauge copies. On a $24^4$ lattice at coupling strengths $\\beta = 2.3$ to $2.6$, the local maxima of the MCG functional over random copies are nearly symmetric Gaussians, and the string tensions extracted from the high-$R_{MCG}$ tail of these distributions agree with the unprojected string tensions in the literature. The same procedure works at $\\beta = 2.7$ after excluding copies whose potentials indicate mis-detected vortices. If right, this gives a practical answer to the Gribov-copy problem in vortex detection and strengthens the center vortex model of confinement.","feed_headline":"Random gauge copies in a Gaussian window reproduce the string tension","feed_subtitle":"Maximizing only over the Gaussian spread of gauge copies matches unprojected lattice results for QCD string tension.","key_machinery":"The maximal center gauge (MCG) functional $R_{MCG}$ of equation (1), which measures how close gauge-transformed link variables are to center elements, together with center projection of links onto the nearest center element $\\pm 1$ for SU(2). The new element is the statistical treatment of the ensemble of local maxima: the paper works with the distribution of $R_{MCG}$ over random gauge copies rather than with a single best copy, restricts the maximization to the Gaussian core of that distribution, and uses the near-linearity of $\\sigma_{cp}(R_{MCG})$ to select copies that reproduce the unprojected string tension.","core_discovery":"The paper's central claim is that the information needed for the string tension is contained in the ensemble of local maxima of the MCG functional, quantified by $R_{MCG} = \\sum_x \\sum_\\mu |\\mathrm{Tr}[{}^g U_\\mu(x)]|^2$, and that it must be extracted by maximizing only within the approximately Gaussian distribution of $R_{MCG}$ values over random gauge copies. With 100 random copies per configuration on a $24^4$ lattice, the paper finds Gaussian distributions for $2.3 \\le \\beta \\le 2.6$, with small deviations at $\\beta=2.6$ and large ones at $\\beta=2.7$. Because the center-projected string tension $\\sigma_{cp}$ is nearly linear in $R_{MCG}$, taking the two or three highest-$R_{MCG}$ copies within the Gaussian subset reproduces the unprojected string tensions from the literature, whereas unrestricted maximization underestimates them. For $\\beta=2.6$ the ensembles are symmetrized and for $\\beta=2.7$ copies with reversed $R$-dependence of the potential are excluded, after which the remaining copies agree with the other couplings.","pith_inferences":["A testable extension would be to apply the same Gaussian-restricted selection rule to SU(3) gauge theory, where center projection uses the three nontrivial center elements and the distribution of local maxima may take a different shape.","The non-Gaussian tails at $\\beta=2.6$ and $\\beta=2.7$ could be cross-checked with an independent vortex detector; if the excluded copies still contain thick vortices, the exclusion criterion would be a selection effect rather than a physical diagnosis.","The linear $\\sigma_{cp}$--$R_{MCG}$ relation suggests an extrapolation prescription: read the string tension at the Gaussian tip from a linear fit, which would reduce the number of random copies needed per configuration."],"forward_implications":["Unrestricted MCG maximization systematically lowers the center-projected string tension, so the global maximum of the gauge functional is the wrong target for vortex detection.","Taking the two or three highest-$R_{MCG}$ copies inside the Gaussian subset reproduces the unprojected string tension for $2.3 \\le \\beta \\le 2.6$, while taking more than about seven copies overestimates it.","The nearly linear relation between $\\sigma_{cp}$ and $R_{MCG}$ turns the Gribov-copy ensemble into a controlled statistical sample instead of an ambiguity to be minimized away.","At $\\beta=2.6$ and $\\beta=2.7$, the projected potential itself flags bad copies: when the large-$R$ string tension falls below the small-$R$ value, the copy has mis-detected vortices and can be excluded."],"supporting_citations":[{"why":"Analytically shows that trivial maximization of the MCG functional cannot accurately yield the correct physics, motivating the restricted-maximization prescription.","marker":"[12]"},{"why":"Numerically demonstrate that increasing the gauge-functional value reduces the string tension, establishing the problem that the ensemble-based approach addresses.","marker":"[13,14]"},{"why":"Supplies the symmetrization procedure applied to the ensembles of local maxima for beta=2.6.","marker":"[19]"},{"why":"Previously reported the same nearly linear relation for Creutz ratios, which the paper uses as supporting evidence for the linear dependence on the gauge functional.","marker":"[15]"},{"why":"Provide the unprojected string-tension values from the literature that the center-projected results are compared against.","marker":"[8,16–18]"},{"why":"Established the original MCG center-projection method with random gauge copies whose string-tension predictions the paper builds on.","marker":"[6,7]"}],"fun_headline_variants":["String tension from local maxima in Gaussian gauge window","Gaussian gauge copies pin down QCD string tension","Center vortices validated by Gaussian choice of copies","Top gauge copies in Gaussian window match string tension","Maximize only Gaussian-spread copies to extract string tension"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the non-Gaussian deviations seen at $\\beta=2.6$ and $\\beta=2.7$ are artifacts of incorrect vortex detection rather than physical effects; this is what justifies symmetrizing the $\\beta=2.6$ ensembles and dropping the $\\beta=2.7$ copies that show an anomalous potential.","fun_headline_variants_meta":{"raw":{"variants":["String tension from local maxima in Gaussian gauge window","Gaussian gauge copies pin down QCD string tension","Center vortices validated by Gaussian choice of copies","Top gauge copies in Gaussian window match string tension","Maximize only Gaussian-spread copies to extract string tension"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00026,"raw_usage":{"total_tokens":1569,"prompt_tokens":901,"completion_tokens":668,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":594}},"tokens_in":517,"tokens_out":668,"duration_ms":6230,"temperature":1.0,"reasoning_tokens":594,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:37:20.036365+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the analysis on a larger lattice or with a different lattice action and check whether the Gaussian-restricted prescription reproduces the unprojected string tension without any symmetrization or copy exclusion; if agreement at $\\beta=2.7$ requires removing the non-Gaussian tails, the claim that the random-copy ensemble contains the full string-tension information would be selection-dependent rather than measured.","supporting_citations":[{"cited_title":"Center projection vortices in continuum Yang-Mills theory // Nucl.Phys","cited_arxiv_id":null,"evidence_quote":"Analytically shows that trivial maximization of the MCG functional cannot accurately yield the correct physics, motivating the restricted-maximization prescription."},{"cited_title":"Searches for Exotic Hadrons at GlueX","cited_arxiv_id":"1908.09711","evidence_quote":"Supplies the symmetrization procedure applied to the ensembles of local maxima for beta=2.6."},{"cited_title":"A first analysis of the ensemble of local maxima of maximal center gauge","cited_arxiv_id":"2301.05903","evidence_quote":"Previously reported the same nearly linear relation for Creutz ratios, which the paper uses as supporting evidence for the linear dependence on the gauge functional."}],"review_version":1}