Pith. Sign in

REVIEW

Confining Solution of the Dyson-Schwinger Equations in Coulomb Gauge

High Energy Physics - Theory · Submitted 21 Dec 2006 · 13 references

Pith's verdict: A competent numerical companion to the authors' own infrared analysis: it finds the β=1 branch in a full momentum-space solve and gets a linear Coulomb potential, but the exponent extraction has no stability analysis and the potential is algebraically the same statement, not a second result.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-th/0612241 (version 2)
Pith IDQW2NZR3K, Pith's permanent record of this paper and its review
CitationD. Epple, H. Reinhardt, W. Schleifenbaum. "Confining Solution of the Dyson-Schwinger Equations in Coulomb Gauge." arXiv:hep-th/0612241 (2006). https://pith.science/paper/QW2NZR3K
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The Coulomb-gauge equations of Yang-Mills theory admit a second, confining solution whose static quark potential rises strictly linearly.

The reading

This paper asks whether the integral equations that follow from minimizing the vacuum energy of Yang-Mills theory in Coulomb gauge have more than one consistent infrared solution, and whether any of them actually confines. An earlier numerical treatment found only the branch with infrared exponent $beta approx 0.796$; a separate analytic study had argued a second branch exists at $beta = 1$. The authors improve the numerical machinery enough to land on that second branch, and report that it produces a static quark potential rising strictly linearly with separation, together with a running coupling that approaches a finite infrared fixed point. The point of interest is that confinement here is a property of the same coupled equations that determine the gluon and ghost propagators, rather than a separate ingredient bolted on afterwards.

What carries the argument

Three ingredients carry the argument. First, the horizon condition $d^{-1}(mu=0)=0$, imported from Zwanziger, forces the ghost exponent $beta$ positive and is what makes the gluon energy infrared divergent. Second, the sum rule $2beta = alpha + 1$, which follows from the non-renormalization of the ghost-gluon vertex, ties the gluon exponent $alpha$ to the ghost exponent $beta$; combined with the gap equation's infrared limit $lim_{kto 0}(omega(k)-chi(k)) = c$, it is what allows $beta_2=1$ to produce a strictly linear potential rather than a merely rising one. Third, the numerical scheme: Chebyshev expansion on a logarithmic momentum grid, double-logarithmic representation for divergent functions, analytic power-law and logarithmic asymptotics used as fit forms whose exponents and coefficients are still determined by nonlinear least squares, Gauss-Legendre quadrature after a logarithmic change of radial variable, and relaxed iteration with $r approx 0.1$. The design point is that the asymptotics are used as fitting templates, not imposed, so the numerically extracted exponents serve as a check on the analytic predictions.

Extended reading notes

Core claim

The central claim is that the one-loop Dyson-Schwinger equations of the Hamiltonian (variational) approach to Yang-Mills theory in Coulomb gauge possess a second numerical solution, with infrared exponent $beta_2 = 1$, in addition to the previously found $beta_1 approx 0.796$. The authors solve the coupled equations for the ghost form factor $d(k)$, the gluon energy $omega(k)$ and the scalar curvature $chi(k)$ over the whole momentum range, imposing the horizon condition $d^{-1}(0)=0$ and the renormalization choices $bar{chi}(1)=0$, $c=0$, with $bar{omega}(1)$ left as a free parameter (their displayed solution uses $bar{omega}(1)=1.6$). In the infrared, $omega$ and $chi$ are found to share the same divergent power law, consistent with the analytic expectation, and $omega(k)-chi(k)$ tends to a constant. Feeding the solution into the Coulomb term of the Hamiltonian for two static colour charges gives a potential that is linear in $r$ to within an estimated error under one percent, and the ghost-gluon running coupling $alpha(k) = frac{2}{3pi},k,d^2(k),omega^{-1}(k)$ is monotonic with an infrared fixed point matching the analytically predicted $alpha(0) = 16pi/(3N_c)$.

Load-bearing premise

The horizon condition $d^{-1}(0)=0$ is imposed as an input from outside rather than derived from the variational principle, and without it the ghost propagator is not enhanced, the exponent is not forced positive, and the confining solution does not follow.

Editorial extensions

If this is right

  • The variational Hamiltonian approach in Coulomb gauge reproduces a linearly rising static quark potential from its own gap equations, so confinement need not be added as an external assumption about the infrared behaviour of the propagators.
  • The ghost-gluon running coupling has a finite infrared limit $alpha(0)=16pi/(3N_c)$, so the coupling does not diverge in the infrared in this scheme, which constrains how the beta function can behave at low momentum.
  • The Coulomb string tension extracted from the linear fit can be compared with lattice determinations, which the authors note are a factor of roughly 1.5 to 3 larger than the Wilson-loop string tension, in line with the analytic statement that the Coulomb tension is an upper bound.
  • Because the infrared exponents and ultraviolet behaviour are reported to be insensitive to the remaining renormalization constants once the horizon condition is adopted, the qualitative confining picture is claimed to be robust against those choices, with changes confined to intermediate momenta.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the $beta=1$ branch is a genuine feature of the untruncated equations rather than of the one-loop truncation and the fitted asymptotic templates, then the earlier numerical search that found only $beta approx 0.796$ was limited by its discretization and not by the physics; a convergence study in grid density and fit range would settle which it is.
  • The reported insensitivity to $bar{omega}(1)$ suggests a one-parameter family of confining solutions; whether the Coulomb string tension varies across that family is a concrete question the paper's own figures could answer but does not quantify.
  • The linearity claim is a statement about the Coulomb potential alone, so it is testable against lattice Coulomb-gauge measurements of the same quantity at intermediate separations, where the variational truncation is least controlled.
  • Since the horizon condition is what drives the whole chain, a useful extension would be to solve the same truncated system with $d^{-1}(0)$ held at small nonzero values and watch the exponent and the potential's slope move, turning an input assumption into a dial.
Original abstract

The Dyson-Schwinger equations arising from minimizing the vacuum energy density in the Hamiltonian approach to Yang-Mills theory in Coulomb gauge are solved numerically. A new solution is presented which gives rise to a strictly linearly rising static quark potential and whose existence was previously observed in the infrared analysis of the Dyson-Schwinger equations. For the new solution we also present the static quark potential and calculate the running coupling constant from the ghost-gluon vertex.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,964 input tokens, $0.0001 reported.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.07 · human participants 0.01 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.19 · compares with prior results 0.64 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

The Dyson-Schwinger equations arising from minimizing the vacuum energy density in the Hamiltonian approach to Yang-Mills theory in Coulomb gauge are solved numerically. A new solution is presented which gives rise to a strictly linearly rising static quark potential and whose existence was previously observed in the infrared analysis of the Dyson-Schwinger equations. For the new solution we also present the static quark potential and calculate the running coupling constant from the ghost-gluon vertex.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The Dyson-Schwinger equations arising from minimizing the vacuum energy density in the Hamiltonian approach to Yang-Mills theory in Coulomb gauge are solved numerically. method 0.82 0.92 0.08 asserted 1.00
s1 A new solution is presented which gives rise to a strictly linearly rising static quark potential and whose existence was previously observed in the infrared analysis of the Dyson-Schwinger equations. result 0.72 0.74 0.75 asserted 1.00
s2 For the new solution we also present the static quark potential and calculate the running coupling constant from the ghost-gluon vertex. contribution 0.51 0.94 0.15 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.92 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1: Left panel: the gluon energy
Figure 2. shows the dependency of ω(k) and d(k) upon the renormalization constant ω(µ). We find that this constant does not change the infrared exponents and the ultraviolet behaviour but does have impact on the prefactors A, B of the power laws for ω(k), χ(k) and d(k). Also, there are some quantitative changes in the intermediate momentum range
Figure 3. FIG. 3: (left) The running coupling

Figures reproduced from arXiv:hep-th/0612241 by the authors.

References
How later papers used it

49 later papers cite it

Evidence: some. 16 of the citing papers share no author with it, and 2 of those carry a confident label from Pith's citation classifier.

Pith read 18 of the 124 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its data or values 2 of 5 papers with this label · 2 by the same authors

By the same authors

How the literature uses it

47 arXiv papers cite this in 121 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 79 sentences without a confident label

Builds on it 2 papers
Uses its method or tool 3 papers
Uses its data or values 5 papers
  • (Within our approximation, this is the only place where the actual form of the vacuum wave functional enters.) The exact analytical expressions are unknown, but we have reliable data from recent variational calculations [ref] , which are also in decent agreement with high precision lattice measurements [ref] .

    On the temporal Wilson loop in the Hamiltonian approach in Coulomb gauge · Yang-Mills theory in $(3+1)$ dimensions · 0.94

  • Fig. 1 shows the result for the dielectric function ([ref]) obtained from the solutions to the Dyson-Schwinger equations reported in [ref] .

    Dielectric function of the QCD vacuum · The dielectric function · 0.91

  • All numerical plots displayed in this section were calculated by D.Epple [ref] .

    Nonperturbative aspects of Yang-Mills theory · Full numerical solutions · 0.85

Compares with it 2 papers
Cites it as support 4 papers
  • Solutions of the Dyson-Schwinger equations for the three-gluon and four-gluon couplings [ref] and phenomenological studies [ref] of QCD couplings based on physical observables such as $τ$ decay [ref] and the Bjorken sum rule [ref] , show that the QCD $β$ function vanishes and $α_s(Q^2)$ become constant at small virtuality; i.e., effective charges develop an “infrared fixed point." Recent lattice simulations [ref] and nonperturbative analyses [ref] have also indicated an infrared fixed point for QCD.

    AdS/CFT and Light-Front QCD · Introduction · 1.00

  • Our variational approach has given a quite decent description of the infrared sector of Yang-Mills theory as, for example, a linearly rising non-Abelian Coulomb potential [ref] , an infrared diverging gluon energy (expressing confinement) [ref] in accord with lattice data [ref] , an infrared finite running coupling constant [ref] , a perimeter law for the 't Hooft loop [ref] , an area law for the Wilson loop [ref] and a dielectric function of the Yang-Mills vacuum in accord with the bag model picture [ref] .

    Quark Sector of the QCD Groundstate in Coulomb Gauge · Introduction · 0.97

  • It is precisely this picture which is supported by the variational calculation in the Hamiltonian approach to Yang-Mills theory in Coulomb gauge developed in Refs. [ref] .

    A covariant variational approach to Yang-Mills Theory · Gaussian trial measure · 1.00

Cites it as background 21 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.