{"id":"faaab93f-48be-4cde-af62-711241d19124","arxiv_id":"2505.16940","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Lattice simulations in the center-symmetric Landau gauge show that the link average and the D33-D88 gluon propagator difference signal center-symmetry breaking across the deconfinement transition.","lead":"This paper tests a new gauge-fixing condition, the center-symmetric Landau gauge, in lattice QCD and measures the gluon propagator below and above the deconfinement transition. The data suggest the average link and the difference between gluon color components are order parameters for the transition.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim rests on gauge-fixed observables, but no Gribov-copy analysis is provided; uncontrolled local maxima of Eq. (5) could bias both the link average and the D33/D88 difference.","rationale":"The central claim is that two gauge-dependent observables act as order parameters. The most fragile step is the gauge-fixing itself: Eq. (5) is maximized with a steepest-descent method, and the paper provides no analysis of how the results depend on the choice of local maximum. This matters because both observables are directly sensitive to the temporal link orientation: the link average is the average of U4, and the D33/D88 split is the difference of color components in the Cartan directions selected by the background. Without a Gribov-copy study, the observed behavior—agreement with the center-symmetric value below Tc and deviation above—could in principle be caused by the algorithm's tendency to converge to particular local maxima. The reader identified the same assumption; the present stress test confirms it and adds that the absence of error bars on the propagator (Fig. 4) makes the D33/D88 part especially hard to evaluate. Since this is a preliminary proceedings report and the concern can be tested with standard methods, the CONDITIONAL verdict remains appropriate; no change is needed.","tokens_in":7017,"tokens_out":19426,"duration_ms":145901,"concrete_test":"Take a subset of 20-50 configurations from each ensemble (64^3×8, T=243 MeV and 64^3×6, T=324 MeV) and repeat the gauge fixing from 10-20 independent random gauge copies per configuration. For each copy, compute the diagonal elements of the link average and the longitudinal gluon propagator components D33 and D88 at the lowest few momenta. If the copy-to-copy spread (standard deviation) of these observables is small compared to the observed deviation from the center-symmetric value above Tc and compared to the D33-D88 difference, the order-parameter claim is robust; if the spread is comparable or larger, the claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 6 concludes that the link average and the D33/D88 difference are order parameters for deconfinement. Both observables are defined after maximizing the gauge-fixing functional (5) with a steepest-descent method (Section 2). The functional is non-convex and admits many local maxima; the algorithm selects one without any Gribov-copy control. Since the functional explicitly contains the center-symmetric background V4 in the temporal term, its global maximum drives U4 toward V4. Below Tc the measured link average (Eqs. 12-14) matches the center-symmetric prediction to within ~0.1%, which could simply reflect the gauge-fixing bias rather than a physical order parameter. Above Tc, the observed deviation and the D33/D88 splitting in Fig. 4 could also be influenced by a sector-dependent distribution of local maxima; no error bars are shown on the propagator, and only one volume and two temperatures are used. Without a Gribov-copy analysis, one cannot rule out that the alleged order-parameter behavior is an artifact of the algorithm.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a lattice implementation of the center-symmetric Landau gauge, defined by maximizing the gauge-fixing functional (5) built from a center-symmetric background field. Using Wilson gauge action ensembles at beta=6.0 on a 64^3 x 8 lattice at T=243 MeV (below the deconfinement temperature) and on a 64^3 x 6 lattice at T=324 MeV (above), the authors compare the link average with analytical predictions from Refs. [13,14] and compute the bare longitudinal gluon propagator for color components 3 and 8. They find that below T_c the link average matches the center-symmetric prediction and that D33 and D88 coincide, while above T_c the link average deviates and D33 and D88 separate. On this basis they conclude in Section 6 that the link average and the difference of the color components of the gluon propagator can serve as order parameters for the deconfinement transition.","tokens_in":7168,"tokens_out":3650,"duration_ms":34587,"significance":"If established, the proposed observables would provide local, gauge-fixed order parameters for the SU(3) deconfinement transition and would strengthen the connection between the center-symmetric Landau gauge formalism and lattice data. The paper has concrete strengths: the analytical predictions from Refs. [13,14] are parameter-free, the lattice data are an independent numerical test of those predictions, the sector-rotation consistency check in Figures 1 and 2 is a useful internal cross-check, and the numerical values in Eqs. (12)-(14) agree with the analytic prediction to about 0.1% below T_c. However, the evidence is preliminary in ways that matter for the central claim: there are only two temperatures and one spatial volume, no continuum extrapolation is attempted, the propagator data in Figure 4 are shown without error bars, and no Gribov-copy analysis is provided for the non-convex gauge-fixing functional. These limitations are partly acknowledged by the authors' characterization of the work as the 'first steps' of a program, but they are not fully reflected in the strength of the wording in Section 6, where the order-parameter interpretation is stated without qualification.","major_comments":[{"comment":"The gauge-fixing functional (5) is non-convex, and the Fourier-accelerated steepest descent method described in the last paragraph of Section 2 selects one local maximum without any Gribov-copy control. Since both the link average (Section 4) and the D33/D88 propagator comparison (Figure 4) are evaluated on gauge-fixed configurations, the close agreement below T_c in Eqs. (12)-(14) may partly reflect the bias of the functional toward the center-symmetric background rather than an independent physical signal, and the D33/D88 splitting above T_c could be influenced by sector-dependent stationary points. The authors should provide a Gribov-copy analysis—for instance, comparing observables obtained from multiple random gauge starts or from several distinct local maxima—or otherwise demonstrate that the reported results are insensitive to the choice of extremum.","section":"Section 2 and Section 5"},{"comment":"The propagator plots in Figure 4 show data points without statistical error bars, yet the central claim that D33=D88 below T_c and that D33 separates from D88 above T_c is made from a visual comparison. Without errors, the reader cannot assess the statistical significance of the splitting, and with only two temperatures and one lattice volume the evidence is too thin to support an order-parameter claim. The authors should report jackknife or bootstrap errors on the propagator points and, ideally, add results at additional volumes or temperatures to show that the effect survives beyond a single ensemble.","section":"Section 5, Figure 4"},{"comment":"The link-average evidence is based on exactly two ensembles: 64^3 x 8 at T=243 MeV and 64^3 x 6 at T=324 MeV. The Monte Carlo history in Figure 3, which is used to argue that the link average tracks the confinement-deconfinement transition, comes from a different simulation taken from Ref. [2], not from the present setup. Finite-volume effects can mix center sectors, as Figure 3 itself shows, and this could bias the averaged link values in an uncontrolled way. The authors should either perform a volume-dependence study or explicitly discuss the expected size of finite-volume corrections before concluding that the link average is an order parameter.","section":"Section 4"}],"minor_comments":[{"comment":"The definition of u0^dagger(mu) in Eq. (6) is written in a compressed notation; please spell out the exponent explicitly in terms of the lattice spacing a, the temperature T = 1/(a N_t), and the generators, so that the lattice implementation is unambiguous.","section":"Equation (6)"},{"comment":"The axis annotations 'aa=3' and 'aa=8' should read 'a=3' and 'a=8' or be replaced by an explicit statement that these are color indices; the vertical axis label D_L(p^2) should also state whether the plotted quantity is the bare or renormalized propagator and in which units.","section":"Figure 4"},{"comment":"Since det<U4> is 1 for an SU(3) matrix, the denominator in Eq. (11) is trivial; please clarify the intended normalization and define all symbols used in the displayed expression.","section":"Equation (11)"},{"comment":"References [19]-[23] contain garbled or incomplete DOI strings; please update them with the correct identifiers or remove the broken URL fragments.","section":"References"},{"comment":"The conclusion that the link average and the D33/D88 difference are order parameters is stated more strongly than the preliminary two-temperature, single-volume, no-error-bar data warrant; consider wording such as 'candidate order parameters' or explicitly list the required checks in the concluding paragraph.","section":"Section 6"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style contribution, and the authors are appropriately careful in describing the results as preliminary. However, the central claim in Section 6 is not backed by the numerical evidence as presented: the missing Gribov-copy analysis and the absence of error bars on the propagator are load-bearing issues, not just presentation concerns. If the venue accepts preliminary studies with clearly qualified conclusions, a shorter revision that softens the order-parameter statement and adds the missing caveats may suffice; otherwise the requested major revisions should be implemented."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is an honest proceedings report that delivers what it promises—a lattice implementation of the center-symmetric Landau gauge and the first gluon-propagator data in that gauge. The D33 = D88 degeneracy below T_c and its splitting above are real, new results. The paper is not yet a full study, and the main caveat is the missing Gribov-copy analysis, but nothing here looks wrong.\n\nThe new content is straightforward. They define the gauge-fixing functional in Eq. (5), maximize it with a Fourier-accelerated steepest descent, and show that the link average in the confined phase sits at the analytically predicted center-symmetric value (the diagonal elements agree to about 0.1%). Above T_c the link average shifts away from that value, and the D33 and D88 longitudinal gluon propagators, which coincide in the confined phase, separate cleanly. The sector-rotation plots in Figs. 1–2 are a nice way to show that the three Z3 sectors map onto each other. For a proceedings this is solid, careful work.\n\nThe soft spots are real but proportionate. First, there is no Gribov-copy analysis. The functional (5) is not convex and the algorithm is steepest descent; different local maxima could give different link averages. Because the functional contains the background field explicitly, the confined-phase agreement might be partly inherited from the gauge condition rather than being an independent check. That possibility is not addressed. Second, the propagator plots in Fig. 4 have no error bars and only two temperatures on a single volume. No continuum extrapolation is attempted. Third, the order-parameter conclusion in Sec. 6 is phrased more strongly than the two-point evidence warrants, though the authors do call it preliminary and point to a companion paper for the link-average study.\n\nThe stress-test worry that the entire deconfined signal could be a gauge-fixing artifact does not hold for the D33/D88 splitting: that quantity is a difference between color components in the same gauge-fixed configuration and is not directly pinned by the background-field term. The magnitude of the split could be affected by Gribov copies, but the existence of the split above T_c is a physical signal. The link average, by contrast, is more exposed to the bias concern. So the right reading is: promising preliminary results, with one systematic check that needs to be done before the order-parameter claims become solid.\n\nWho benefits: lattice QCD groups working on propagators and confinement, and functional-method people comparing with continuum predictions. It deserves a serious referee—a proceedings referee, not a full PRL referee—and should pass with minor revisions (add errors, conduct a Gribov-copy test, soften the conclusion or narrow it to a 'candidate' order parameter). I would take the companion paper more seriously, but this one is worth reading.","headline":"First lattice data for the center-symmetric Landau gauge gluon propagator: promising, honest, but missing Gribov-copy checks and propagator error bars.","tokens_in":7736,"tokens_out":5128,"would_cite":true,"duration_ms":43022,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T25","81V05"],"pacs":["11.15.Ha","12.38.Gc"],"model":"deepseek-v4-flash","headline":"This paper reports lattice evidence that, in the center-symmetric Landau gauge, the normalized link average and the difference between color-3 and color-8 gluon propagators serve as order parameters for the deconfinement transition.","keywords":["center-symmetric Landau gauge","deconfinement transition","gluon propagator","Polyakov loop","lattice QCD","center symmetry","link average","SU(3) gauge theory"],"falsifier":"Perform the same gauge fixing on identical ensembles with many random initial gauge seeds: if the normalized link average or the $D^{33}_L-D^{88}_L$ difference changes between distinct local maxima of $F$ within the same $\\mathbb{Z}_3$ sector by more than the quoted statistical errors, the claimed order-parameter behavior is an artifact of the chosen extremum. A cleaner test is to bracket $T_c$ with several $N_t$ values: if the jump in these observables does not occur at the temperature where the Polyakov loop susceptibility peaks, the identification with the deconfinement transition fails.","tokens_in":6814,"feed_emoji":"⚛️","tokens_out":8987,"duration_ms":70561,"temperature":0.7,"pith_summary":"The paper reports the first lattice computation of the gluon propagator in the center-symmetric Landau gauge and argues that two gauge-fixed observables act as order parameters for the deconfinement transition in pure SU(3) gauge theory. Below the critical temperature, the average temporal link, normalized by the cube root of its determinant, sits at the center-symmetric value and the color-3 and color-8 gluon propagators coincide. Above $T_c$, the link average moves away from that value, configurations cluster into three $\\mathbb{Z}_3$-related sectors, and $D^{33}$ separates from $D^{88}$. The numerical agreement with the analytical predictions in the confined phase is the main evidence, and the same observables are proposed as practical probes of center-symmetry breaking on the lattice.","feed_headline":"Gluon propagator's color split signals deconfinement","feed_subtitle":"Below Tc the link average and D33-D88 match theory; above Tc both deviate — a local probe of deconfinement.","key_machinery":"The object that carries the argument is the gauge-fixing functional $F=\\sum_{x,\\mu}\\operatorname{Re}\\operatorname{Tr}[U_0^\\dagger(\\mu)U_\\mu(x)]$, where $U_0^\\dagger(\\mu)$ encodes the center-symmetric background $e^{i a g \\bar{A}_{0,\\mu}}$ with $\\bar{A}_{0,\\mu}=(T/g)(4\\pi/3)t_7\\delta_{\\mu 0}$. Maximizing $F$ over gauge orbits selects a representative field configuration in each of the three $\\mathbb{Z}_3$ sectors, and $F$ is invariant under the particular center transformation that rotates the Polyakov loop by $2\\pi/3$. The two proposed order parameters are the normalized link average and the color-resolved gluon propagator difference $D^{33}_L-D^{88}_L$; both are predicted analytically to take the center-symmetric value in the confined phase and to deviate in the deconfined phase, and the lattice data are compared with those predictions.","core_discovery":"The central claim is that a lattice version of the center-symmetric Landau gauge, defined by maximizing $\\sum_{x,\\mu}\\operatorname{Re}\\operatorname{Tr}[U_0^\\dagger(\\mu)U_\\mu(x)]$ with the background phase $U_0^\\dagger(\\mu)=\\exp(i a g \\bar{A}_{0,\\mu}\\delta_{\\mu 4})$, can be fixed in a way that makes the deconfinement transition visible directly in gluonic two-point data. The authors compute the link average $\\langle U_4\\rangle/(\\det\\langle U_4\\rangle)^{1/3}$ and the longitudinal gluon propagator $D_L(p^2)$ for color indices 3 and 8 on $64^3\\times N_t$ lattices at $\\beta=6.0$. For $T=243$ MeV (below $T_c$) the link average matches the center-symmetric prediction and $D^{33}=D^{88}$ within errors; for $T=324$ MeV (above $T_c$) the link average deviates and $D^{33}$ and $D^{88}$ split. This is presented as evidence that both quantities can serve as order parameters for the deconfinement phase transition.","pith_inferences":["Beyond the paper: if the same two observables stay well-defined with dynamical quarks, they could serve as diagnostics of the deconfinement crossover, but the explicit breaking of center symmetry by quarks would likely smear the sharp jump seen here.","Beyond the paper: scanning $\\beta$ at fixed $N_t$ to bracket $T_c$ more finely would turn the link-average jump into a quantitative determination of $T_c$; the two temperatures reported here indicate but do not establish this.","Beyond the paper: the lack of a Gribov-copy check means the clearest next test is seed dependence; if different local maxima give different order-parameter values, a selection rule (such as the global maximum) must be specified."],"forward_implications":["The normalized link average and the $D^{33}_L-D^{88}_L$ difference can be computed in ordinary lattice simulations and function as order parameters without relying on the Polyakov loop.","Below $T_c$, matching $D^{33}=D^{88}$ and the link-average value against analytic predictions confirms the center-symmetric background picture; above $T_c$, their deviation is a signal of $\\mathbb{Z}_3$ breaking.","The gluon propagator can be decomposed into longitudinal and transverse parts with the same procedure as standard Landau gauge, making these observables cheap to add to existing analyses.","The three-sector clustering of $(A_4^3,A_4^8)$ above $T_c$ gives a per-configuration view of which $\\mathbb{Z}_3$ sector is selected."],"supporting_citations":[{"why":"Establishes that local order parameters for the transition can be constructed in the center-symmetric Landau gauge, motivating this lattice study.","marker":"[12]"},{"why":"Derives the analytic prediction that the color-3 and color-8 gluon propagators coincide in the symmetric phase.","marker":"[13]"},{"why":"Gives the continuum prediction for the center-symmetric gluon field configuration and link average that the lattice data are compared with.","marker":"[14]"},{"why":"Supplies the Fourier-accelerated steepest-descent method used to maximize the gauge-fixing functional.","marker":"[15]"},{"why":"Provides the standard-Landau-gauge procedure for extracting longitudinal and transverse gluon propagator components that the paper adapts here.","marker":"[1]"},{"why":"Provides the Monte Carlo history near $T_c$ that shows the link-average difference vanishing in the confined phase and nonzero in the deconfined phase.","marker":"[2]"},{"why":"Contains the detailed theoretical framework and updated link-average results underlying this preliminary report.","marker":"[16]"}],"fun_headline_variants":["Center-symmetric Landau gauge exposes deconfinement","Gluon color split reveals deconfinement in lattice QCD","Lattice probe: color 3 and 8 gluon split above Tc","Deconfinement seen in gluon propagator's color asymmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The gauge-fixing algorithm is assumed to land on a representative maximum of the gauge-fixing functional, yet no check is made that different starting points give the same link average and propagators; if the result depends on which maximum is chosen, the order-parameter interpretation above $T_c$ is not robust.","fun_headline_variants_meta":{"raw":{"variants":["Center-symmetric Landau gauge exposes deconfinement","Gluon color split reveals deconfinement in lattice QCD","Lattice probe: color 3 and 8 gluon split above Tc","Deconfinement seen in gluon propagator's color asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00061,"raw_usage":{"total_tokens":2788,"prompt_tokens":842,"completion_tokens":1946,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":458,"completion_tokens_details":{"reasoning_tokens":1870}},"tokens_in":458,"tokens_out":1946,"duration_ms":11390,"temperature":1.0,"reasoning_tokens":1870,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:51:58.107569+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the same gauge fixing on identical ensembles with many random initial gauge seeds: if the normalized link average or the $D^{33}_L-D^{88}_L$ difference changes between distinct local maxima of $F$ within the same $\\mathbb{Z}_3$ sector by more than the quoted statistical errors, the claimed order-parameter behavior is an artifact of the chosen extremum. A cleaner test is to bracket $T_c$ with several $N_t$ values: if the jump in these observables does not occur at the temperature where the Polyakov loop susceptibility peaks, the identification with the deconfinement transition fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that local order parameters for the transition can be constructed in the center-symmetric Landau gauge, motivating this lattice study."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the analytic prediction that the color-3 and color-8 gluon propagators coincide in the symmetric phase."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the continuum prediction for the center-symmetric gluon field configuration and link average that the lattice data are compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Fourier-accelerated steepest-descent method used to maximize the gauge-fixing functional."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard-Landau-gauge procedure for extracting longitudinal and transverse gluon propagator components that the paper adapts here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Monte Carlo history near $T_c$ that shows the link-average difference vanishing in the confined phase and nonzero in the deconfined phase."},{"cited_title":"The center-symmetric Landau gauge meets the lattice","cited_arxiv_id":"2412.07930","evidence_quote":"Contains the detailed theoretical framework and updated link-average results underlying this preliminary report."}],"review_version":1}