Differential Dyson-Schwinger equations for quantum chromodynamics

TL;DR

This paper derives differential Dyson-Schwinger equations for QCD, focusing on two-point functions, and explores their behavior in the large color limit to understand phenomena like chiral symmetry breaking and confinement.

Contribution

It introduces a differential form of Dyson-Schwinger equations for QCD and analyzes their behavior in the 't Hooft limit, providing insights into low-energy phenomena.

Findings

Large N limit describes low-energy QCD behavior

Conditions for confinement are discussed

Chiral symmetry breaking is analyzed

Abstract

Using a technique devised by Bender, Milton and Savage, we derive the Dyson-Schwinger equations for quantum chromodynamics in differential form. We stop our analysis to the two-point functions. The 't~Hooft limit of color number going to infinity is derived showing how these equations can be cast into a treatable even if approximate form. It is seen how this limit gives a sound description of the low-energy behavior of quantum chromodynamics by discussing the dynamical breaking of chiral symmetry and confinement, providing a condition for the latter. This approach exploits a background field technique in quantum field theory.

Equations

$$\mathcal{L}={\mathcal{L}}_{inv}+{\mathcal{L}}_{gf}+{\mathcal{L}}_{FP}$$

$${\mathcal{L}}_{inv}$$

$${\mathcal{L}}_{gf}$$

$${\mathcal{L}}_{FP}$$

$$D_{\mu }\psi$$

$$D_{\mu }c$$

$$Z[j]=\mathcal{N}\int [d\varphi ]e^{iS[\varphi ]-i\int d^{4}xj(x)\varphi (x)}$$

$$Z[j]\to Z[j]⟨e^{i\int d^{4}x\alpha (x)\left(\frac{\delta S}{\delta \varphi (x)}-j(x)\varphi (x)\right)}{⟩}_j,$$

$${⟨\frac{\delta S}{\delta \varphi (x)}⟩}_j=j(x).$$

$${\partial }^2{G}_1$$

$$G_2^{-1}(x-y)$$

$${\partial }^2{G}_1+m^2{G}_1+\lambda G_1^3+3\lambda G_2(0)+\lambda G_3(0,0)=0$$

$$({\partial }^2+m^2)G_2(x-y)+3\lambda G_2(0)G_2(x-y)+3\lambda [G_1(x-y){]}^2{G}_2(x-y)$$

$$+3\lambda G_3(0,x-y)G_1(x)+\lambda G_4(0,0,x-y))={\delta }^4(x-y).$$

$${\partial }^{\mu }{\partial }_{\mu }{A}_{\nu }^a-{\alpha }^{-1}{\partial }_{\nu }(\partial \cdot A^a)+g{f}^{abc}{A}^{b\mu }({\partial }_{\mu }{A}_{\nu }^c-{\partial }_{\nu }{A}_{\mu }^c)+g{f}^{abc}{\partial }^{\mu }(A_{\mu }^b{A}_{\nu }^c)+g^2{f}^{abc}{f}^{cde}{A}^{b\mu }{A}_{\mu }^d{A}_{\nu }^e$$

$$=g{f}^{abc}{\partial }_{\nu }({\bar{c}}^b{c}^c)+g\sum _{q,i}{\bar{q}}^i{\gamma }_{\nu }{T}^a{q}^i+j_{\nu }^a$$

$${\partial }^{\mu }{\partial }_{\mu }{c}^a+g{f}^{abc}{\partial }^{\mu }(A_{\mu }^b{c}^c)={\epsilon }^a$$

$$(i/\partial +g\mathbf{T}\cdot /\mathbf{A}-m_q)q^i={\eta }_q^i.$$

$${\partial }^2{G}_{1\nu }^{(j)a}(x)+g{f}^{abc}(⟨A^{b\mu }{\partial }_{\mu }{A}_{\nu }^c⟩-⟨A^{b\mu }{\partial }_{\nu }{A}_{\mu }^c⟩)Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+g{f}^{abc}{\partial }^{\mu }⟨A_{\mu }^b{A}_{\nu }^c⟩Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]$$

$$+g^2{f}^{abc}{f}^{cde}⟨A^{b\mu }{A}_{\mu }^d{A}_{\nu }^e⟩Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=g{f}^{abc}⟨{\partial }_{\nu }({\bar{c}}^b{c}^c)⟩Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]$$

$$+g\sum _{q,i}⟨{\bar{q}}^i{\gamma }_{\nu }{T}^a{q}^i⟩Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+j_{\nu }^a$$

$${\partial }^2{P}_1^{(\epsilon )a}(x)+g{f}^{abc}{\partial }^{\mu }⟨A_{\mu }^b{c}^c⟩Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]={\epsilon }^a$$

$$(i/\partial -m_q)q_1^{(\eta )i}+g\mathbf{T}\cdot ⟨/\mathbf{A}{q}^i⟩Z_{QCD}^{-1}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]={\eta }_q^i.$$

$$G_{1\nu }^{(j)a}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨A_{\nu }^a(x)⟩$$

$$P_1^{(\epsilon )a}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨c^a(x)⟩$$

$$q_1^{(\eta )i}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨q^i(x)⟩.$$

$$G_{2\nu \kappa }^{(j)ab}(x,x)Z_{YM}[j,\bar{\epsilon },ϵ]+G_{1\nu }^{(j)a}(x)G_{1\kappa }^{(j)b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]=⟨A_{\nu }^a(x)A_{\kappa }^b(x)⟩$$

$$H_2^{(\epsilon )ab}(x,x)Z_{YM}[j,\bar{\epsilon },ϵ]+P_1^{(\epsilon )a}(x)P_1^{(\epsilon )b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]=⟨c^b(x)c^a(x)⟩$$

$${\bar{H}}_2^{(\epsilon )ab}(x,x)Z_{YM}[j,\bar{\epsilon },ϵ]+{\bar{P}}_1^{(\epsilon )a}(x){\bar{P}}_1^{(\epsilon )b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]=⟨{\bar{c}}^b(x){\bar{c}}^a(x)⟩$$

$$P_2^{(\epsilon )ab}(x,x)Z_{YM}[j,\bar{\epsilon },ϵ]+{\bar{P}}_1^{(\epsilon )a}(x)P_1^{(\epsilon )b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]=⟨{\bar{c}}^b(x)c^a(x)⟩$$

$${\partial }_{\mu }{G}_{2\nu \kappa }^{(j)ab}(x,x)Z_{YM}[j,\bar{\epsilon },ϵ]+{\partial }_{\mu }{G}_{1\nu }^{(j)a}(x)G_{1\kappa }^{(j)b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]=⟨{\partial }_{\mu }{A}_{\nu }^a(x)A_{\kappa }^b(x)⟩$$

$$K_{2\nu }^{(\epsilon ,j)ab}(x,x)Z_{YM}[j,\bar{\epsilon },ϵ]+P_1^{(\epsilon )a}(x)G_{1\nu }^{(j)b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]=⟨c^a(x)A_{\nu }^b(x)⟩$$

$$/G_1^{(j)a}(x)q_1^{(\eta )i}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+/W_q^{(\eta )ai}(x,x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨/A^a{q}^i(x)⟩$$

$$S_q^{(\eta )ij}(x,x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+{\bar{q}}_1^{(\eta )i}(x)q_1^{(\eta )j}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨{\bar{q}}^i(x)q^j(x)⟩$$

$$J_{q\mu }^{(\eta )aij}(x,x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+{\bar{q}}_1^{(\eta )i}(x){\gamma }_{\mu }{T}^a{q}_1^{(\eta )j}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨{\bar{q}}^i(x){\gamma }_{\mu }{T}^a{q}^j(x)⟩,$$

$$G_{3\nu \kappa \rho }^{(j)abc}(x,x,x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+G_{2\nu \kappa }^{(j)ab}(x,x)G_{1\rho }^{(j)c}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+$$

$$G_{2\nu \rho }^{(j)ac}(x,x)G_{1\kappa }^{(j)b}(x)Z_{YM}[j,\bar{\epsilon },ϵ]+G_{1\nu }^{(j)a}(x)G_{2\kappa \rho }^{(j)bc}(x,x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]+$$

$$G_{1\nu }^{(j)a}(x)G_{1\kappa }^{(j)b}(x)G_{1\rho }^{(j)c}(x)Z_{QCD}[j,\bar{\epsilon },\epsilon ,{\bar{\eta }}_q,{\eta }_q]=⟨A_{\nu }^a(x)A_{\kappa }^b(x)A_{\rho }^c(x)⟩.$$

$${\partial }^2{G}_{1\nu }^{(j)a}(x)+g{f}^{abc}({\partial }^{\mu }{G}_{2\mu \nu }^{(j)bc}(x,x)+{\partial }^{\mu }{G}_{1\mu }^{(j)b}(x)G_{1\nu }^{(j)c}(x)-{\partial }_{\nu }{G}_{2\mu }^{(j)\mu bc}(x,x)-{\partial }_{\nu }{G}_{1\mu }^{(j)b}(x)G_1^{(j)\mu c}(x))$$

$$+g{f}^{abc}{\partial }^{\mu }{G}_{2\mu \nu }^{(j)bc}(x,x)+g{f}^{abc}{\partial }^{\mu }(G_{1\mu }^{(j)b}(x)G_{1\nu }^{(j)c}(x))$$

$$+g^2{f}^{abc}{f}^{cde}(G_{3\mu \nu }^{(j)\mu bde}(x,x,x)+G_{2\mu \nu }^{(j)bd}(x,x)G_1^{(j)\mu e}(x)$$

$$+G_{2\nu \rho }^{(j)eb}(x,x)G_1^{(j)\rho d}(x)+G_{2\mu \nu }^{(j)de}(x,x)G_1^{(j)\mu b}(x)+$$

$$G_1^{(j)\mu b}(x)G_{1\mu }^{(j)d}(x)G_{1\nu }^{(j)e}(x))=g{f}^{abc}({\partial }_{\nu }{P}_2^{(\epsilon )bc}(x,x)+{\partial }_{\nu }({\bar{P}}_1^{(\epsilon )b}(x)P_1^{(\epsilon )c}(x)))$$

$$+g\sum _{q,i}{J}_{q\nu }^{(\eta )aii}(x,x)+g\sum _{q,i}{\bar{q}}_1^{(\eta )i}(x){\gamma }_{\nu }{T}^a{q}_1^{(\eta )i}(x)+j_{\nu }^a$$

$${\partial }^2{P}_1^{(\epsilon )a}(x)+g{f}^{abc}{\partial }^{\mu }(K_{2\mu }^{(\epsilon ,j)bc}(x,x)+P_1^{(\epsilon )b}(x)G_{1\mu }^{(j)c}(x))={\epsilon }^a$$

$$(i/\partial -m_q)q_1^{(\eta )i}(x)+g\mathbf{T}\cdot /{\mathbf{G}}_1^{(j)}(x)q_1^{(\eta )i}(x)+g\mathbf{T}\cdot /{\mathbf{W}}_q^{(\eta )i}(x,x)q_1^{(\eta )i}(x)={\eta }_q^i.$$

Full text

I Introduction

The main difficulty of quantun chromodynamics (QCD) is that, at low energies, the theory is not amenable to treatment using perturbation techniques. This implies that some non-perturbative methods should be devised to solve them. The most widespread approach is solving the equations of the theory on a large lattice using computer facilities. This permitted to obtain, with a precision of a few percent Fodor:2012gf ; pdg , some relevant observables of the theory. This method improves as the computer resources improve making even more precise the comparison with experiment. Use of numerical techniques is a signal that we miss some sound theoretical approach to compute observables.

A similar situation is seen for the correlation functions of the theory. Studies on the lattice of the gluon and ghost propagators, mostly in the Landau gauge, Bogolubsky:2007ud ; Cucchieri:2007md ; Oliveira:2007px and the spectrum Lucini:2004my ; Chen:2005mg proved that a mass gap appears in a non-Abelian gauge theory without fermions. Theoretical support for these results was presented in Cornwall:1981zr ; Cornwall:2010bk ; Dudal:2008sp ; Frasca:2007uz ; Frasca:2009yp ; Frasca:2015yva providing closed form formulas for the gluon propagator. Quite recently, the set of Dyson-Schwinger equations for this case was solved, for the 1- and 2-point functions, and the spectrum very-well accurately computed both in 3 and 4 dimensions Frasca:2016sky ; Frasca:2015yva ; Frasca:2017slg . Confinement was also proved to be a property of the theory Frasca:2016sky ; Chaichian:2018cyv .

Indeed, the Dyson-Schwinger equations were considered, since the start, the most sensible approach to treat a non-perturbative theory like QCD at low-energies Baker:1976vz ; Eichten:1974et ; Roberts:1994dr and, more recently, Alkofer:2000wg . In any case, the standard technique is to reduce the set of equations, that normally are partial differential equations, to their integral form in momentum space. Some years ago, Bender, Milton and Savage Bender:1999ek proposed to derive the Dyson-Schwinger equations and treat them into differential form. This way to manage these equation was the one used to find the exact solution Frasca:2015yva . This technique appears more general as it permits to work out a solution to a quantum field theory also when a background field is present. This is a rather general situation when a non-trivial solution of the 1-point equation is considered. Such a possibility opens up the opportunity of a complete solution to theories that normally are considered treatable only through perturbation methods. The idea is that, knowing all the correlation functions, a quantum field theory is completely solved.

The aim of this paper is to derive the Dyson-Schwinger equations for QCD in differential form. We obtain them for the 1- and 2-point functions. We show that, in the ’t Hooft limit tHooft:1973alw , the equation can be cast into a treatable form. It appears that a non-local Nambu-Jona-Lasinio model is the proper low-energy limit of the theory Frasca:2016rsi .

The paper is so structured. In Sec. II, we give the main equations and notations. In Sec. III, we presente the technique we will use to obtain the Dyson-Schwinger set of equations. In Sec. IV, we derive the set of equations for the 1- and 2-point correlation functions. In Sec. V, we discuss the ’t Hooft limit. In Sec. VI, we discuss the dynamical breaking of the chiral symmetry and derive a confinement condition. Finally, in Sec. VII, we present the conclusions.

II Basic equations

In QCD, one has the Lagrangian (to fix the notation we assume $(1,-1,-1,-1)$ for the metric signature)

[TABLE]

where ${\cal L}_{{inv}}$ denotes the classical gauge-invariant part, ${\cal L}_{{gf}}$ the gauge-fixing terms and ${\cal L}_{{FP}}$ the Faddeev-Popov (FP) ghost term characteristic of non-Abelian gauge theories Smilga:2001ck

[TABLE]

in the usual notation. We set $\alpha$ for the gauge parameter and $D_{\mu}$ represents the covariant derivative whose explicit forms are given by

[TABLE]

III Dyson-Schwinger equations

It is a rather old result that a tower of equations for the correlation functions of a quantum field theory can be obtained from the equations of motion Dyson:1949ha ; Schwinger:1951ex ; Schwinger:1951hq . The solutions of this set of equations, when performed exactly, solves completely the theory itself. So, it appears like a very powerful non-perturbative approach. Anyway, the set of equations is such that any given equation depends on the correlation functions of higher order making a truncation needed. The set can be obtained using functional techniques as shown in Roberts:1994dr ; Alkofer:2000wg and, more recently, in Bender:1999ek . We are going to describe both the methods to derive them with a simple example while this paper will rely on the latter. It is important to emphasize that, whatever approach one takes, the starting point are always the classical equations of motion satisfying a given variational principle.

III.1 Standard technique

The general starting point is the partition function, for a theory with action $S[\phi]$,

[TABLE]

being $\phi(x)$ a scalar field. The idea is that this functional integral does not change after a reparametrization $\phi(x)\rightarrow\phi(x)+\alpha(x)$, being $\alpha(x)$ an arbitrary function. Then, this gives

[TABLE]

from which, by requiring invariance, we derive the quantum equation of motion

[TABLE]

The procedure can be iterated taking further derivative with respect to $j$ and this will give all the set of Dyson-Schwinger equations for the correlation functions. We will obtain the final result by taking $j=0$ at the end of the computation. We note that from the lhs of eq. (9) one gets the average on the classical equations of motion of the theory that are the starting point of the procedure.

From this point on, the standard technique introduce the vertex functions $\Gamma_{n}(\phi)$ using a Legendre transform. This means that, for a simple $\phi^{4}$ theory, we will get the first two Dyson-Schwinger equations

[TABLE]

Written in this way, the only viable approach is through a Fourier transform to momenta and working with integral equations. This could hide the physical content of the equations and the deep meaning of a possible truncation.

III.2 Bender-Milton-Savage method

Bender-Milton-Savage method Bender:1999ek takes the move from eq. (9), again the average of the classical equations of motion, but does not take any Legendre transform so, vertex functions are not introduced at all. Rather, one always works with higher-order n-point functions $G_{n}(x_{1},x_{2},\ldots,x_{n})$. This permits to preserve the differential structure of the Dyson-Schwinger equations making them particular useful when exact solutions are known. This will give for a $\phi^{4}$ theory Frasca:2015yva

[TABLE]

In this case, we have non-trivial exact solutions for $G_{1}$, being $G_{3}(0,0)=0$ and the the set of Dyson-Schwinger equations becomes treatable without any truncation. This would not have been possible with the Dyson-Schwinger equations written in the integral form as demanded by the standard technique.

IV Correlation functions in QCD

As already stated above, we will use the Bender-Milton-Savage method. In order to get the Dyson-Schwinger equations, we have to start from the classical equations of motion, given by the functional derivatives of the action as in eq. (9), and average on them. These are given by

[TABLE]

Letters $a,b,c,\ldots=1,\ldots,8$ are for gluons and $i,j,k,l,\ldots=1,2,3$ are for colors. The scalar product runs on the former. Flavors are identified with the letter $q=u,d,s,c,t,b$. We fix the gauge to the Landau gauge, $\alpha\rightarrow 0$, and $c,\ \bar{c}$ are the ghost fields. Averaging on the vacuum state and dividing by the partition function $Z_{QCD}[j,\bar{\varepsilon},\varepsilon,\bar{\eta}_{q},\eta_{q}]$ one has

[TABLE]

The one-point functions are given by

[TABLE]

These represent the 1-point function for the gluon, ghost and quark fields respectively. Deriving once with respect to the currents, at the same point because of the averages on the vacuum (see Bender:1999ek ), one has

[TABLE]

being $W^{(\eta)ai}_{q\mu}(x,y)=\delta G^{(j)a}_{1\mu}(x)/\delta\eta_{q}^{i}(y)$, $J_{q\mu}^{(\eta)aij}(x,y)=\gamma_{\mu}T^{a}\delta q_{1}^{(\eta)i}(x)/\delta\eta_{q}^{j}(y)=\gamma_{\mu}T^{a}S^{(\eta)ij}_{q}(x,y)$. Deriving twice one has

[TABLE]

This give us the first set of Schwinger-Dyson equations as

[TABLE]

Setting the currents to zero and noticing that, by translation invariance, is $G_{2}(x,x)=G_{2}(x-x)=G_{2}(0)$, $G_{3}(x,x,x)=G_{3}(0,0)$, $K_{2}(x,x)=K_{2}(0)$, $W^{(\eta)a}(x,x)=W^{(\eta)a}(0)$ and $J_{q\nu}(x,x)=J_{q\nu}(0)$, we get

[TABLE]

The Schwinger-Dyson equation for the two-point functions can be obtained by further deriving eq.(IV). One has

[TABLE]

We have assumed $\langle\bar{q}\bar{q}\rangle=0$, $\langle\bar{q}^{\prime}q\rangle=0$ for $q^{\prime}\neq q$, $Y^{(\varepsilon)b\ldots i}_{nq}(x,y,\ldots,z)=\delta P_{n+1}^{(\varepsilon)b\ldots}(x,y,\ldots)/\delta\eta_{q}^{i}$, $Q^{(\eta)ia}_{\mu}(x-y)=\delta q_{1}^{(eta)i}(x)/\delta j^{a}_{\mu}(y)$, $J_{1q\mu}^{(\eta)aijk}(x,y,z)=\delta J_{q\mu}^{(\eta)aij}(x,y)/\delta\eta^{k}_{q}(z)=\gamma_{\mu}T_{a}\delta S^{(\eta)ij}_{q}(x,y)/\delta\eta^{k}_{q}(z)$, $W^{(\eta)a\ldots i}_{nq\mu}(x,y,\ldots,z)=\delta G^{(j)a\ldots}_{{n+1}\mu}(x,y,\ldots)/\delta{\eta}_{q}^{j}(z)$. Note that, for $n=0$, we omit the index, e.g. $W_{0}=W$. This yields, setting currents to zero and using translation invariance,

[TABLE]

[TABLE]

V ’t Hooft limit

We will give here an approximate solution to the set of equations we obtained in the preceding section. The technique we will use is an iterative one starting from the exact solutions we obtained or the 1- and 2-point functions of the Yang-Mills theory Frasca:2015yva . These correlation functions will be modified by the presence of quarks but, in a first iteration, we assume they are a good non-perturbative approximation working well in the deep infrared limit. In this way, we will be able to derive a confining Nambu-Jona-Lasinio approximation. We start from the set (IV) and take

[TABLE]

For symmetry properties of $P_{2}$ under exchange of indexes the set reduces to

[TABLE]

We have introduced the quark mass matrix

[TABLE]

We note that is degenerate with respect the color index. Now, we write the expected solution in the form in the Landau gauge

[TABLE]

being $\eta_{\nu}^{a}$ some numerical coefficients, $\phi(x)$ a scalar field and $\Delta(x-y)$ the propagator. We will get for the 1-point functions

[TABLE]

In order to understand the behavior of the $\eta-$symbols, we work out the case for SU(2) where the proof of their existence is straightforward. Indeed, in this case, they can be defined as

[TABLE]

that yields

[TABLE]

that implies $\eta_{\mu}^{a}\eta^{a\mu}=3$. This generalizes to SU(N) as

[TABLE]

Similarly, by generalizing the SU(2) case,

[TABLE]

and

[TABLE]

being $g_{\mu\nu}$ the Minkowski metric and $\delta_{\mu\nu}$ the identity tensor. We note that the number of components of $\eta_{\mu}^{a}$ are identical to those of $\gamma_{\mu}T^{a}$, as can be checked already for the aforementioned SU(2) case, as expected and the equation is consistent. This permits to write for the 1-point functions, using eq. (28),

[TABLE]

where, in the first equation, a summation on the repeated index $a$ is implied giving explicitly the term $\eta\cdot\gamma\cdot T$. Therefore, the 1-point function acquires a mass term given by $\mu_{0}^{2}=Ng^{2}\Delta(0)$. Now, we apply eq. (V) to the 2-point function set to obtain

[TABLE]

This set of equations can be solved in the ’t Hooft limit $N\gg 1$, keeping $Ng^{2}$ constant, and $Ng^{2}\gg 1$. We have to show the ordering in $g$ for these equations, as by now we are applying a perturbation technique for a strongly coupled theory. This can be accomplished by a perturbation theory devised in Frasca:2013tma . In order to get the right ordering of terms in a strong coupling limit, we rescale the space-time variables as $x_{\mu}\rightarrow\sqrt{Ng^{2}}x_{\mu}$. E.g. this implies for the first equation in (V)

[TABLE]

Therefore, we realize that the RHS is at least of order $O(1/\sqrt{Ng^{2}})$ and so, negligible with respect to the LHS in the ’t Hooft limit as the scalar field scales accordingly. A similar argument can be applied to eqs. (V). For the quark equation, we will show in a moment that the $\phi$ field, in this approximation, will take a factor $1/\sqrt{Ng^{2}}$. Essentially, the idea, for the gluon field, is that it is of order $Ng^{2}$ while the corrections implied by the quark field are always $\sqrt{Ng^{2}}$ implying the start of an iterative procedure to solve this system of equations. Then, our leading order equations can be written as

[TABLE]

Indeed, for 1-point function we have the non-trivial solution Frasca:2017slg

[TABLE]

being $\mu$ and $\chi$ arbitrary integration constants and $\kappa=\frac{-m^{2}+\sqrt{m^{4}+2Ng^{2}\mu^{4}}}{-m^{2}-\sqrt{m^{4}+2Ng^{2}\mu^{4}}}$. We have set $m^{2}=2Ng^{2}\Delta(0)$ and taken the momenta $p$ so that

[TABLE]

In this way, we are able to see the consistency of our approximation as $\phi\sim O(1/\sqrt{Ng^{2}})$ and then, when applied to the LHS of the given equations, this changes them in higher order contributions. We just observe that the leading order is non-perturbative, due to its dependency on the coupling, and non-trivial and so, it is consistent in the infrared limit. For this reason, we write the second set of equations as

[TABLE]

This set of equations can be solved exactly and yields the self-consistency equations of the theory. In order to present the solutions, we point out that the equation

[TABLE]

has the solution

[TABLE]

Here we have set

[TABLE]

and

[TABLE]

The mass spectrum is given by

[TABLE]

Then,

[TABLE]

This self-consistency equation provides the proper spectrum of a Yang-Mills theory with no fermions Frasca:2017slg , in very close agreement with lattice data. Therefore, we get immediately

[TABLE]

The last equation is an integral equation to be solved iteratively. This will yield the corrections to the gluon propagator.

Similarly, for the quark mass we have to solve the self-consistent system of equations

[TABLE]

given the mass matrix

[TABLE]

that couples both. This set will yield the quark propagator. These equations can be written as

[TABLE]

At this stage, we can act perturbatively by iterating on the leading order solutions. We will see that the solution is consistent as it should. Then, the first approximation for ${\hat{q}}_{1}^{i}(x)=q_{0}^{i}(x)$ is the free particle solution. The second equation will yield the approximation

[TABLE]

that is a non-local Nambu-Jona-Lasinio approximation provided we identify ${\hat{S}}^{ji}_{0q}(y^{\prime}-x)$ with the free Dirac propagator.

VI Dynamical chiral symmetry breaking and confinement in the ’t Hooft limit

We now discuss how dynamical chiral symmetry breaking arises naturally in the ’t  Hooft limit giving also a condition for confinement of quarks. In order to do this we need to evaluate the quark propagator. We consider eq.(V) and we take the Nambu-Jona-Lasinio approximation. One has

[TABLE]

being the quark self-energy

[TABLE]

We realize that the non-local kernel comes out from the gluon propagator $\Delta_{0}(x-y)$. We note that

[TABLE]

Quark self-energy says us that the corrections to the quark propagator can depend on momenta. So, we will have dynamical symmetry breaking until a solution for the quark mass can be found. When the on-shell mass condition fails, varying the coupling, we will have a confined quark. This is a possible definition of quark confinement given by Gribov Gribov:1998kb and Roberts and Williams Roberts:1994dr . Therefore, we compute the gap equation in the limit of very low momenta. This will yield

[TABLE]

where the trace is over flavors, colors and spinor indexes. This yields

[TABLE]

having moved to Euclidean. This equation admits a real solution for the local limit $\Delta_{0}(p)\rightarrow{\rm constant}$ for a critical value of the coupling $g$. This implies a free quark and so, no confinement. In the general case, from the gluon propagator eq. (39), we can write (in the Euclidean limit)

[TABLE]

with the mass spectrum, $m_{n}$, given in eq. (42). We assume no degeneracy between $m_{n}$ and $M_{q}$ and we approximate the gluon propagator in eq.(39) with the correction to the mass $m^{2}$ neglected. This is a fairly good approximation as shown in Frasca:2017slg . So, we have to evaluate

[TABLE]

The integral can be evaluated exactly when a cut-off $\Lambda$ is used, as usual for Nambu-Jona-Lasinio models. This yields

[TABLE]

This equation is amenable to a numerical treatment. It should be solved with the conditions $M_{q}\geq 0$ and $M_{q}\ll\Lambda$ that is, the effective mass of the quark should not exceed the ultraviolet cut-off representing, at least, the boundary of the region where asymptotic freedom starts to set in (generally taken at $\Lambda\approx 1\ {\rm GeV}$). We can normalize this equation to the cut-off $\Lambda$ by introducing the new variables $x=m_{0}/\Lambda$ and $y=M_{q}/\Lambda$ having taken $m_{n}=(2n+1)m_{0}$. The mass $m_{0}$ can be assumed to be that of the $\sigma$ meson or f(500) that we fix to $m_{0}=0.417\ {\rm GeV}$. Then,

[TABLE]

being $\kappa=N_{f}(N^{2}-1)N/8\pi$ and $\alpha_{s}=g^{2}/4\pi$. We note that the cut-off is completely disappeared except for the ratio $m_{q}/\Lambda$ that, for the light quarks, is negligible small. The result for eq.(57), written in implicit form

[TABLE]

and fixing $x=m_{0}/\Lambda$ as said above, is given in fig. 1 where the zeros of the function $\mu(x,y)$ are of interest here.

One can see that, when $\alpha_{s}$ is moving from the asymptotic freedom region, $\alpha_{s}<1$, toward the infrared fixed point generally taken at $\alpha_{s}=\pi$, the number of solutions moves from one simple root, very near 0, where the quark retains its small mass, to 3 roots with the quark gaining a dynamical mass and the root near 0 persisting but becoming negative and so, unphysical. Such a scenario is better appreciated in fig. 2.

This implies that, the chiral symmetry, if quarks would be massless, is broken. It is however interesting to note that there is a critical coupling beyond which the effective mass could overcome the cut-off entering into an unphysical region. This is for $\alpha_{s}\approx 0.5$, given the value of $m_{0}$ we chose.

Indeed, this particular value of the coupling can be evaluate exactly and can give a confinement condition that applies quite generally to light quarks (u,d,s). This holds when the confinement criteria we stated above applies, in agreement with Gribov:1998kb ; Roberts:1994dr , that is the on-shell quark mass fails for the propagator. This is given by assuming that that we evaluate eq, (58) with the condition $y=M_{q}/\Lambda<1$. So, the value $M_{q}=\Lambda$ provides a confinement condition. We will get from eq, (58)

[TABLE]

being $\xi$ given by

[TABLE]

For the values we have chosen, we get $\alpha_{s}\approx 0.5105$ in agreement with our previous qualitative analysis derived from the breaking of the chiral symmetry. So, in ordinary QCD such a minimum exists and the theory is confining.

VII Conclusions

We have derived the set of Dyson-Schwinger equations for QCD for 1- and 2-point correlation functions. We have seen that they can be cast in a treatable form in the ’t Hooft limit. The main aim of this effort is to present a proof of confinement for the theory.

Indeed, in the ’t Hooft limit, we were able to show how the transition from the regime of dynamical breaking of symmetry transits to confinement, deriving a general condition for confinement of light quarks that is satisfied by QCD. We hope in future works to extend this analysis.