{"id":"cca80940-522d-464d-a568-b0804f766028","arxiv_id":"2412.12124","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"The paper claims that at leading order in 1/g, the Yang-Mills gluon propagator is a contact term proportional to the delta function, signaling a gluon condensate, but the derivation contains a likely incorrect group-integral step and fits the mass scale by hand.","lead":"By applying their 'effective locality' formalism to pure Yang-Mills theory, the authors compute the strongly coupled gluon propagator and find it proportional to a delta function, implying no propagation and a gluon condensate. The result is meant to support a mass-gap and confinement mechanism, but the derivation rests on an unjustified group-integration identity and on tuning the coupling to match lattice data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (49) is an incorrect O(N) group integral: standard Haar integration gives (vol/N) δ_ij Σ_k 1/ξ_k, and the δ(ξ_k) factor that makes C1 nonzero only emerges after pairing ±ξ with a 2/N prefactor, so Eq. (53) is unsupported.","rationale":"The reader's weakest_assumption correctly identifies Eq. (49) as the load-bearing step, and my independent calculation confirms that the identity is false as stated. The O(N) average of O_{ki}O_{kj} is a textbook result: for any fixed k, ∫ dO O_{ki}O_{kj} = (vol/N)δ_ij. Summing over k with D^{-1}_{kk} = 1/ξ_k gives (vol/N)δ_ij Σ_k 1/ξ_k. This contains no delta function. The delta in Eq. (49) could only arise after the ξ integration, and even then the paired-spectrum calculation yields a prefactor 2/N relative to the paper's claim. Thus the derivation of the contact-term propagator, Eq. (53), is unsupported. I also considered whether the numerical comparison could independently support the claim, but the procedure in Section IIID sets g_R to 10 or 4 and extracts √Δ from the lattice value, so it is a fit, not a prediction. The absence of machine-checked proofs or independent verification further raises correctness risk. Therefore the reader's REJECT verdict is appropriate; I see no reason to change it. The concern is not about disagreement with the physics consensus but about an internally incorrect mathematical step. A corrected group integral might yield a similar qualitative picture, but the paper as presented does not establish it.","tokens_in":17380,"tokens_out":7713,"duration_ms":80912,"concrete_test":"Evaluate Eq. (49) for N=2 using the explicit 2×2 block parameterization of Eq. (34). With O(θ) = [[cos θ, -sin θ], [sin θ, cos θ]] and Haar measure dθ, compute ∫_0^{2π} dθ Σ_{k=1}^2 O_{ki}O_{kj}/ξ_k for i=j=1. The result is π(1/ξ_1 + 1/ξ_2), whereas Eq. (49) gives -iπ vol(O_2) δ(ξ_1) (with vol(O_2)=4π). The discrepancy is immediate. Independently, use ∫ dO O_{ki}O_{kj} = (vol/N) δ_ij and sum over the paired spectrum with the iε prescription to confirm that the delta factor carries an extra 2/N compared to Eq. (49). Then recompute C1/C2 in Eqs. (54)-(55) with this corrected prefactor and check whether the coefficient in Eq. (53) changes by 2/N.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result Eq. (53) depends on Eq. (49), which asserts ∫ dO Σ_k (tO)_{ik} D^{-1}_{kl} O_{lj} = -iπ δ_ij vol(O_N) Σ_{k=1}^{N/2} δ(ξ_k). Standard Haar integration for O(N) gives ∫ dO O_{ki}O_{kj} = (vol(O_N)/N) δ_ij, independent of k. Therefore the left side of Eq. (49) equals (vol(O_N)/N) δ_ij Σ_{k=1}^N 1/(ξ_k+iε), a smooth function of the eigenvalues, not a sum of delta functions. Only after combining the paired spectrum (ξ_p, -ξ_p) and using 1/(x+iε) = v.p.(1/x) - iπδ(x) does a delta appear: each pair contributes [1/(ξ_p+iε)+1/(-ξ_p+iε)] = -2πi δ(ξ_p), so the correct distribution is -(2πi vol(O_N)/N) δ_ij Σ_{p=1}^{N/2} δ(ξ_p), which differs from Eq. (49) by a factor 2/N. For N=32 this is a factor 1/16. This is not a minor prefactor: without the δ(ξ_1) factor, the integral C1 in Eq. (54) would not have the form assumed, and the √Δ behavior of Eq. (53) rests on an unjustified delta. The subsequent numerical comparison in Section IIID, where g_R is chosen and √Δ is extracted from the lattice value, cannot compensate for an invalid mathematical step. The derivation should be redone with the correct O(N) integral; the qualitative contact-term might survive with a modified prefactor, but Eq. (53) as written is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper claims to compute the non-perturbative gluon propagator in strongly coupled Yang-Mills theory starting from the effective-locality generating functional constructed in Ref. [8]. After writing the field-strength tensor as a real symmetric 32x32 matrix M(x), the authors perform an orthogonal-group integration and obtain a leading-order contact term proportional to sqrt(Delta) delta^(4)(x-y) delta_ab g_munu, which is interpreted as a gluon condensate of mass dimension two controlled by a 'meshing parameter' Delta. The paper also presents a random-matrix estimate of the ratio of the two relevant integrals and a numerical comparison in which the meshing parameter is matched to a lattice value of the Yang-Mills scale.","tokens_in":17816,"tokens_out":13144,"duration_ms":139754,"significance":"If the derivation were valid, a pure contact-term gluon propagator at strong coupling would be a striking result with potential implications for the mass gap and confinement. The paper is clearly organized and it is useful that the authors identify the technical bottleneck, namely the Vandermonde integral and the need for random-matrix bounds. However, the central mathematical step, Eq. (49), is incorrect, and the numerical comparison is a fit rather than a test. The manuscript does not provide machine-checked proofs, reproducible code, parameter-free derivations, or independent falsifiable predictions; its main output rests on an unsupported group-integral identity.","major_comments":[{"comment":"Equation (49) is the load-bearing step of the derivation, but the claimed O(N) integration is not correct. For the diagonal D^{-1} used in Eq. (48), the left-hand side is sum_k xi_k^{-1} O_{ki}O_{kj}. The normalized Haar average of O_{ki}O_{kj} is delta_{ij}/N for each k, so the full integral equals (vol(O_N)/N) delta_{ij} sum_{k=1}^N 1/(xi_k+i epsilon), not -i pi delta_{ij} vol(O_N) sum_{k=1}^{N/2} delta(xi_k). With the paired spectrum xi_{N-i+1}=-xi_i this standard result becomes -(2 pi i vol(O_N)/N) delta_{ij} sum_{p=1}^{N/2} delta(xi_p), which differs from Eq. (49) by a factor 2/N and by the way the eigenvalue sum is handled. This is not a cosmetic prefactor: the delta(xi_1) factor in Eq. (54), which produces the nonzero integral C1 and the sqrt(Delta) dependence in Eq. (53), is an artifact of the incorrect identity. As written, Eq. (53) is unsupported.","section":"Eq. (49), Section III.B.2"},{"comment":"The numerical agreement with the Yang-Mills scale is obtained by tuning. In Eq. (71) the authors choose g_R=10 and extract 1/sqrt(Delta) = (230 MeV)^2, and in the following paragraph they state that the same value is obtained at g_R ~ 4 when Eq. (66) is used. Since g_R is a free input and the lattice value of Lambda_YM is the target of the comparison, the proximity of the numbers is a consequence of this tuning, not a prediction. The relation Delta Lambda^4_YM >= 1 invoked after Eq. (52) is also assumed rather than derived. The paper should present these numbers as an illustration, not as independent confirmation of the result.","section":"Section III.D, Eqs. (71) and (66)"},{"comment":"The random-matrix estimate that fixes the numerical coefficient is uncontrolled. The step leading to Eq. (65) replaces (Theta_k^2 - Theta_l^2)^2 by (Theta_k - Theta_l)^2 in C1 and C2, and then uses Wigner's semicircle law to obtain |C1/C2| ~ (1/pi) sqrt(2/N). This replacement is not justified by the gaussian suppression, and the rigorous bound in Eq. (64) is far too wide to serve as a check. Since Eq. (66), and therefore the g_R ~ 4 extraction in Section III.D, depends on this coefficient, the numerical prefactor in the final result is not established.","section":"Section III.B.5, Eq. (65)"}],"minor_comments":[{"comment":"There are several typographical errors: 'wether' should be 'whether', 'conurations' should be 'configurations', and 'untractable' should be 'intractable'.","section":"Conclusion and Section III.B.4"},{"comment":"After the rescaling below Eq. (49), the same symbol xi_i is retained for the dimensionless integration variables; this makes Eqs. (50)-(55) difficult to follow, especially the argument of the exponential and the delta(xi_1) factor.","section":"Section III.B.3"},{"comment":"The replacement of delta^(4)(x-y) by delta^(4)(X-Y) with X = Delta \\hat{X} is introduced without a precise definition of the rescaling, so the relationship between Eq. (53) and Eq. (70) is ambiguous and should be clarified.","section":"Section III.B.7, Eq. (70)"}],"recommendation":"reject","confidential_remarks":"The central mathematical identity in Eq. (49) is incorrect, and the numerical comparison is a fit. The paper would need a complete re-derivation of the O(N) integral and a fresh analysis of the resulting eigenvalue integrals before the claims could be assessed; I do not see a path to acceptance in the current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe short version: this paper claims a contact-term gluon propagator at strong coupling in pure Yang-Mills, but the central identity it rests on, Eq. (49), is off by a factor of N/2. Standard Haar integration for O(N) gives ∫dO O_ki O_lj = (vol/N) δ_kl δ_ij, so the left side of (49) is (vol/N) δ_ij Σ_k 1/(ξ_k+iε), not −iπ vol δ_ij Σ δ(ξ_k). Only after pairing ± eigenvalues and using 1/(x+iε) = v.p. − iπδ do you get delta functions, and the coefficient is then −2πi vol/N. For N=32 that is a factor of 16 in the numerator, which flows directly into C1 and Eq. (53). So the main result as written is unsupported.\n\nThat said, I don't think this is a throwaway. The calculation of ⟨A⟩=0 is careful and shows a real technique for working in the effective-locality formalism. The random-matrix bounds for C1 and C2 are clever, and the crude M → R approximation in Section III.B.6 reproduces the same √Δ δ(x−y) structure, suggesting the qualitative picture is not an artifact of the flawed group integral. The physical message—gluons don't propagate at strong coupling and form a dimension-two condensate—agrees with other work, e.g., Kondo's, so the paper is on the right side of the literature.\n\nThe other soft spots are real but minor by comparison. The numerical section tunes g_R and extracts √Δ from the lattice value, so the 'prediction' of Λ_YM is more of a fit with two free parameters. Also, the whole framework rests on the generating functional from Ref. [8], which hasn't been independently verified. These are fixable with a clearer statement of what is input and what is output.\n\nBottom line: the qualitative conclusion may survive a corrected O(N) integral, but Eq. (53) in this draft is wrong. I would not accept it as is. If the authors fix the group integral and the prefactor, and are more careful about predicted versus fitted quantities, this could become publishable. As it stands, I'd recommend against sending it to review until the math is corrected.","headline":"The paper's contact-term gluon propagator rests on a wrong O(N) group integral (Eq. 49, missing a 2/N factor), so the headline result is unsupported, though the qualitative picture might survive a corrected calculation.","tokens_in":18361,"tokens_out":18240,"would_cite":false,"duration_ms":172868,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Cy"],"model":"deepseek-v4-flash","headline":"At leading order in the strong-coupling expansion, the Yang-Mills gluon propagator becomes a contact term proportional to $\\delta^{(4)}(x-y)$, signalling that gluons do not propagate and suggesting a mass-dimension-two gluon condensate.","keywords":["Yang-Mills","effective locality","gluon propagator","strong coupling","functional methods","random matrix theory","gluon condensate","mass gap"],"falsifier":"Evaluate the orthogonal-group integral in Eq. (49) numerically at $N=32$ or exactly for small $N$ with the paper's measure: if the delta factors $\\delta(\\xi_k)$ do not appear, the leading-order delta-function propagator is not established. A lattice measurement of the strong-coupling gluon propagator near $x=y$ would also test whether the contact amplitude scales as $\\sqrt{\\Delta}$ as claimed.","tokens_in":17092,"feed_emoji":"⚛️","tokens_out":9397,"duration_ms":100171,"temperature":0.7,"pith_summary":"Effective locality is the property that, once all virtual gauge-field exchanges are summed, non-abelian Green's functions reduce to local contact interactions. This paper applies that property to pure Yang-Mills theory in Minkowski spacetime and asks what the gluon propagator looks like when the coupling is strong. Working at leading order in $1/g$, the authors derive a propagator proportional to $\\delta^{(4)}(x-y)\\delta_{ab}g_{\\mu\\nu}$: a contact term with no spacetime propagation. They interpret this as evidence for a mass-dimension-two gluon condensate whose scale is set by a meshing parameter $\\Delta$, and they show the result is consistent with the Yang-Mills scale and with existing estimates of a gluon mass. If the derivation holds, effective locality gives a calculable, non-perturbative route to the mass gap and confinement.","feed_headline":"Gluon propagator is a delta function at strong coupling","feed_subtitle":"Leading 1/g term kills propagation and signals a mass-squared gluon condensate.","key_machinery":"The carrying object is the effective-locality gluon generating functional $Z_{YM}[j]$, Eq. (13), obtained through Halpern's field-strength representation. The functional integral over the antisymmetric tensor fields $\\chi^a_{\\mu\\nu}$ is re-expressed in terms of the $32\\times 32$ real symmetric matrix $M=\\frac{i}{8}\\chi^a\\otimes T^a$, whose spectrum comes in opposite pairs; the measure image theorem converts $d[\\chi]$ into the Vandermonde-weighted eigenvalue measure times the Haar measure on the orthogonal group. The decisive identity is the orthogonal-group average in Eq. (49), which yields a factor $\\delta(\\xi_k)$ for each paired eigenvalue; after analytic continuation and use of the semicircle law to bound $C_1,C_2$, the final contact amplitude is finite and proportional to $\\sqrt{\\Delta}$.","core_discovery":"The paper's central result is that at leading order in the strong-coupling expansion the non-perturbative gluon propagator is $$\\langle T A^a_\\mu(x) A^b_\\nu(y)\\rangle = \\frac{\\pi N}{2g}\\frac{C_1}{C_2}\\sqrt{\\frac{\\$\\Delta$}{4N_c}}\\,\\delta_{ab}g_{\\mu\\nu}\\$delta^{{(4)}}$(x-y) + O($g^{{-2}}$),$$ with $C_1/C_2$ a finite, computable ratio of Vandermonde integrals and $\\Delta$ the spacetime meshing parameter. Because the only spacetime dependence is a delta function, the authors conclude that gluonic degrees of freedom do not propagate in this regime, and because the colour and Lorentz structure is $\\delta_{ab}g_{\\mu\\nu}$, they identify the amplitude as a gluon condensate of mass dimension two. The associated one-point function $\\langle A^a_\\mu(x)\\rangle$ vanishes at the same order. The derivation uses a strong-coupling simplification of the effective-locality generating functional, followed by a random-matrix reduction and an analytic continuation that leaves a finite constant ratio.","pith_inferences":["If the delta-function form is taken at face value, it implies that the strong-coupling vacuum cannot resolve points closer than the meshing scale $\\Delta$; this is a concrete reinterpretation that could be probed by seeing whether the same $\\Delta$ controls higher $n$-point functions.","The random-matrix technology used here is portable: the same spectral and orthogonal-group decomposition could be applied to the effective-locality quark sector, where one would expect a comparable condensate scale to emerge from the same $\\Delta$.","A finite-$N$ numerical evaluation of Eq. (49) would settle whether the delta factor is a genuine property or an artefact of the analytic continuation; this test is within reach because $N=32$ is small enough for high-precision quadrature.","The paper's numerical match, $1/\\sqrt{\\Delta}\\approx (230\\,\\mathrm{MeV})^2$ at $g_R\\simeq 4$, invites comparison in other gauge choices; a gauge-independent condensate would require showing that the contact coefficient is independent of the $\\zeta$-gauge fixing used in Eq. (7)."],"forward_implications":["Gluon propagation is absent at leading order in the strong-coupling regime, and the first correction enters only at $O(g^{-2})$.","The propagator's colour-Lorentz structure $\\delta_{ab}g_{\\mu\\nu}$ identifies the contact term with a condensate of mass dimension two, the same object as the Landau-gauge gluon and ghost mass.","The amplitude vanishes with the meshing parameter $\\Delta$, recovering the perturbative short-distance behaviour expected from asymptotic freedom.","The vanishing one-point function $\\langle A^a_\\mu(x)\\rangle=0$ confirms that the strong-coupling vacuum carries no preferred Lorentz or colour direction.","The method gives a systematic strong-coupling expansion, so subleading orders could reveal residual propagation effects, possibly of a glueball-like kind."],"supporting_citations":[{"why":"Introduces effective locality for QCD and states the gauge-invariant summation property that the paper transfers to the pure Yang-Mills sector.","marker":"[1]"},{"why":"Derives the effective-locality gluon generating functional in Eq. (13) that is the starting point of every calculation here.","marker":"[8]"},{"why":"Supplies the standard functional-manipulation techniques used to pass from the Yang-Mills action to the generating functional.","marker":"[9]"},{"why":"Provides the mass-dimension-two condensate interpretation linking the delta-function propagator to mass gap and confinement.","marker":"[10]"},{"why":"Supplies the field-strength representation used to linearize $F^2$ and introduce the $\\chi$ fields.","marker":"[11]"},{"why":"Gives the analytic continuation used to define the convergent integrals $C_1$ and $C_2$ after the $\\xi\\to\\sqrt{i}\\Theta$ rotation.","marker":"[14]"},{"why":"Supplies the random-matrix normalization constants and the semicircle law used to bound $C_1/C_2$.","marker":"[15]"},{"why":"Justifies the measure image theorem that turns the functional $\\chi$ integration into spectral and orthogonal-group integrations.","marker":"[16]"},{"why":"Gives the volume of the orthogonal group appearing in the pivotal integration identity (49).","marker":"[18]"}],"fun_headline_variants":["No gluon travel: propagator collapses to a delta function","Strong coupling turns gluons into a delta condensate","Gluon propagator goes local: delta at strong coupling","Delta-function gluon: no propagation at strong coupling","Strong coupling: gluon propagator is a delta spike"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result hinges on the integration formula in Eq. (49), which claims that averaging the inverse of the diagonalized matrix over all rotations leaves a delta function at each eigenvalue, because without that delta factor the propagator would not become a contact term.","fun_headline_variants_meta":{"raw":{"variants":["No gluon travel: propagator collapses to a delta function","Strong coupling turns gluons into a delta condensate","Gluon propagator goes local: delta at strong coupling","Delta-function gluon: no propagation at strong coupling","Strong coupling: gluon propagator is a delta spike"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000613,"raw_usage":{"total_tokens":2805,"prompt_tokens":858,"completion_tokens":1947,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":474,"completion_tokens_details":{"reasoning_tokens":1866}},"tokens_in":474,"tokens_out":1947,"duration_ms":15509,"temperature":1.0,"reasoning_tokens":1866,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:24:43.115074+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the orthogonal-group integral in Eq. (49) numerically at $N=32$ or exactly for small $N$ with the paper's measure: if the delta factors $\\delta(\\xi_k)$ do not appear, the leading-order delta-function propagator is not established. A lattice measurement of the strong-coupling gluon propagator near $x=y$ would also test whether the contact amplitude scales as $\\sqrt{\\Delta}$ as claimed.","supporting_citations":[{"cited_title":"Gauge Invariant Summation of All QCD Virtual Gluon Exchanges","cited_arxiv_id":null,"evidence_quote":"Introduces effective locality for QCD and states the gauge-invariant summation property that the paper transfers to the pure Yang-Mills sector."},{"cited_title":"Effective Locality in the pure gluon sector","cited_arxiv_id":null,"evidence_quote":"Derives the effective-locality gluon generating functional in Eq. (13) that is the starting point of every calculation here."},{"cited_title":"Basics of Functional Methods and Eikonal Models ; Editions Fronti` eres: Paris, France, 1990","cited_arxiv_id":null,"evidence_quote":"Supplies the standard functional-manipulation techniques used to pass from the Yang-Mills action to the generating functional."},{"cited_title":"Vacuum condensate of mass dimension 2 as the origin of mass gap and quark confinement","cited_arxiv_id":null,"evidence_quote":"Provides the mass-dimension-two condensate interpretation linking the delta-function propagator to mass gap and confinement."},{"cited_title":"Field–strength formulation of quantum chromodynamics","cited_arxiv_id":null,"evidence_quote":"Supplies the field-strength representation used to linearize $F^2$ and introduce the $\\chi$ fields."},{"cited_title":"Non–Perturbative QCD Amplitudes at Quenched and Eikonal Approximations","cited_arxiv_id":null,"evidence_quote":"Gives the analytic continuation used to define the convergent integrals $C_1$ and $C_2$ after the $\\xi\\to\\sqrt{i}\\Theta$ rotation."},{"cited_title":"Random Matrices; Academic Press: Cambridge, MA, USA, 1967","cited_arxiv_id":null,"evidence_quote":"Supplies the random-matrix normalization constants and the semicircle law used to bound $C_1/C_2$."},{"cited_title":"Johnson and M.L","cited_arxiv_id":null,"evidence_quote":"Justifies the measure image theorem that turns the functional $\\chi$ integration into spectral and orthogonal-group integrations."},{"cited_title":"Volumes of orthogonal groups and unitary groups","cited_arxiv_id":"1509.00537","evidence_quote":"Gives the volume of the orthogonal group appearing in the pivotal integration identity (49)."}],"review_version":1}