# On the role of dynamical quark mass generation in chiral symmetry breaking in QCD

## Abstract

The phenomenon of dynamical quark mass generation is studied in QCD within the framework of a gauge invariant formalism. An exact relationship is established between the equation satisfied by the scalar part of the two-point gauge invariant quark Green’s function and the quark-antiquark bound state equation in the chiral limit. A possible nontrivial solution of the former yields a massless pseudoscalar solution of the bound state equation with vanishing total momentum. The result is also corroborated by the corresponding Ward-Takahashi identity. The problem is explicitly solved in two-dimensional QCD in the large- limit.

###### keywords:

QCD, quark, dynamical mass generation, chiral symmetry breaking, Wilson loop, gauge invariant Green’s function## 1 Introduction

It has been known for a long time, since the advent of the Nambu–Jona-Lasinio model Nambu:1961tp (), that a mechanism of chiral symmetry breaking in theories with fermions as the elementary matter fields is provided by the dynamical mass generation phenomenon. This mechanism was also analyzed in the framework of QED by Baker, Johnson and Lee Baker:1964zza (), but the presence of the axial-vector current anomaly in Abelian sectors of the chiral group did not allow drawing a similar conclusion.

The problem was later studied in the framework of QCD, with the aid of models applied to the Dyson–Schwinger equations Dyson:1949ha (); Schwinger:1951ex (); Alkofer:2000wg (); Fischer:2006ub () in the Coulomb gauge, where the gluon propagator in the rainbow-ladder approximation is considered with its instantaneous expression, with an appropriate confining behavior Finger:1980dw (); LeYaouanc:1983iy (); Adler:1984ri (); Alkofer:1988tc (). The dynamical mass generation phenomenon was confirmed and phenomenological applications to the meson spectrum were considered. An extension of the approximation beyond the rainbow-ladder one preserved the stability of the phenomenon Bender:1996bb (). On the other hand, the use of the Ward-Takahashi identity for the axial-vector current Preparata:1969hg () establishes a general relationship, independent of approximation schemes of the Bethe–Salpeter kernel, between the dynamically generated quark mass and the pseudoscalar Goldstone boson wave function Maris:1997hd ().

In spite of these positive results, the notion of dynamical quark mass remains ambiguous. Contrary to the case of nonconfining theories, the quark mass is not an observable quantity. Furthermore, the quark propagator not being a gauge invariant object, it is difficult to give a physical interpretation to its singularities. In many examples, an infrared cutoff is needed to make it finite.

The aim of this talk is to present results obtained within a gauge invariant formalism, developed in recent years Sazdjian:2007ng (); Sazdjian:2013rva (), where an exact relationship is established between the dynamical quark mass term and the pseudoscalar Goldstone boson wave function in the total zero momentum limit. The explicit resolution of the problem in analytic form in two-dimensional QCD in the large- limit illustrates the advantage of the gauge invariant approach in providing a physical insight into the interpretation of the dynamical mass generation phenomenon.

## 2 Gauge invariant quark Green’s functions

To define gauge invariant quark Green’s functions, one needs to introduce gluon field path-ordered phase factors (Wilson lines) Mandelstam:1968hz (); Nambu:1978bd (), which are the operators of parallel transport of fields from one point to the other. The paths followed by the phase factors might be arbitrarily chosen, however, for the quark Green’s functions polygonal lines seem to correspond to the optimal choice. First, they can be decomposed along a succession of straight line segments. The latter are Lorentz invariant in form and this in turn allows the classification of the polygonal lines according to the number of segments they contain. Second, polygonal lines form a complete set of lines for the description of the present problem, in the sense that no other types of line are needed to complete the description.

The classification adopted for the polygonal lines is applied to the two-point gauge invariant quark Green’s functions (2PGIQGF). We designate by the path-ordered phase factor along the oriented straight line segment going from to . The quark fields, with mass parameter , are assumed to belong to the defining fundamental (-dimensional) representation of the color gauge group . The 2PGIQGF with a polygonal line composed of the succession of segments with junction points is designated by and defined as

(1) | |||||

(The vacuum expectation value is defined in the path integral formalism; spinor indices are omitted and the color indices are implicitly summed.) The simplest such function is , having a phase factor along one straight line segment:

(2) |

Use of the equations of motion of Green’s functions and integrations yield functional relations between the Green’s functions of different classes of polygonal line. It turns out that is the only dynamically independent 2PGIQGF. All other 2PGIQGFs , with , are calculable from through the functional relations. These involve essentially Wilson loops Wilson:1974sk () along polygonal contours and functional derivatives acting on their sides.

The equation satisfied by can be brought to the following form:

(3) | |||||

where the kernels () contain Wilson loop averages along polygonal contours, which are at most -sided, and () 2PGIQGFs and their derivative. The total number of derivatives contained in is , each derivative acting on a different segment. Since the high-index s can be expressed in terms of , Eq. (3) is an integrodifferential equation in , which is the primary unknown quantity to be solved. (Integration symbols on the intermediate variables have been omitted in the right-hand side.) Equation (3) is the analog of the self-energy Dyson–Schwinger equation in the present formalism.

## 3 Bound state wave functions

The previous definitions and approach can also be applied to four-point Green’s functions. We introduce two quark fields (in flavor space) and , with masses and , respectively, and consider gauge invariant four-point Green’s functions made of two quark and two antiquark fields and appropriate phase factors along polygonal lines. Considering a mesonic bound state, with total momentum , one defines corresponding gauge invariant wave functions with a similar classification as for the 2PGIQGFs ():

(4) | |||||

In the simplest case , one has

(5) |

The above wave functions, for all s, describe the same bound state, but differ in their expressions due to their differences in their contents with respect to the phase factor lines.

As in the case of the 2PGIQGFs, functional relations exist between the above wave functions: s, with , are calculable in terms of .

The wave function satisfies two Dirac-type integrodifferential equations, which can be written in the following compact forms:

(6) | |||||

(7) | |||||

The kernels are similar to those met in Eq. (3); they contain Wilson loop averages along polygonal contours with sides; the additional indices 1 or 2 reflect the differences that arise from the presence there of the 2PGIQGFs of quarks 1 and/or 2. The kernels , with indices 1 or 2, are obtained from the s by functional derivation with respect to one of the 2PGIQGFs. The star operations are compact notaions for the various integration operations on the intermediate variables. Since the wave functions are calculable in terms of , the above equations are ultimately wave equations for .

## 4 Chiral symmetry breaking

To analyze the chiral symmetry breaking possibility, we consider the 2PGIQGF . Since the single straight line segment that is present in the path-ordered phase factor depends only on the end points and , we can decompose along two Poincaré invariant parts, made of a vector and of a scalar:

(8) |

(.) In perturbation theory, the scalar component is proportional to the quark mass . When the latter vanishes, also vanishes. Therefore, if the solution of the equation of yields a scalar part that does not vanish in the limit of a massless quark, then one is in the presence of a nonperturbative solution, which may be described as corresponding to the phenomenon of dynamical mass generation. This mass term might have a complicated structure, neither being a constant nor producing a simple pole in the 2PGIQGF.

The equation satisfied by can be extracted from that of [Eq. (3)] by considering the anticommutator of with . One obtains in the massless quark limit

(9) | |||||

This equation is to be compared with that of the wave function in the massless limit:

(10) | |||||

(The quark indices 1 and 2 have been removed, since the quarks are massless.) Equations similar to (9) and (10) with operators acting from the right on the variable could also be written down.

We observe that the equation of is the same as the equation of , provided the general correspondences are also made. Actually, these correspondences are independently justified by comparing the equations each of them satisfy. They represent the general consistency conditions for the similarity of the set of all equations of the above types.

Therefore, if the equation of has, in the massless quark limit, a nonvanishing solution for , then this is also a solution of the bound state equation of for a pseudoscalar state with zero total momentum , since does not depend on . This is possible only if ; this means that the pseudoscalar state is massless and represents a Goldstone boson for chiral symmetry breaking.

A few comments are in order here. First, the above results are exact statements of QCD, although for unrenormalized quantities, since no approximations were made throughout the calculations; all expressions of the kernels are explicitly known in terms of Wilson loop averages and functional derivatives. Second, one may adopt, at the practical level of resolution of the equations, an approximation scheme based mainly on the truncation of the series of kernels according to the number of sides of the polygonal contours of the Wilson loops. There are indications that the series expansions in the right-hand sides of the integrodifferential equations are perturbative with respect to the inverse of the number of sides of polygonal contours. Actually, the first kernels and vanish for symmetry reasons and therefore the leading kernels are and , corresponding to triangular Wilson loops. In any event, to maintain the consistency between the various equations, one should adopt the same type of approximation in all equations. Third, the conclusions that are drawn are gauge invariant, since all quantities under consideration are individually gauge invariant and this property is not altered by the approximation scheme of the truncation type. This is one of the advantages of the gauge invariant formalism.

## 5 Ward-Takahashi identities

Ward–Takahashi identities (WTI) express in compact form relations among Green’s functions that emerge from symmetry properties of a theory. Here we sketch the WTI adapted to the gauge invariant formalism that we are using and corresponding to chiral symmetry.

Considering the case of two different quarks, 1 and 2, we define the axial-vector current and the pseudoscalar density as

and introduce the full axial-vector and pseudoscalar vertex functions

( are the quark spinor indices.) The WTI takes the form

The analysis that follows is then standard. One passes to momentum space, where only two momenta are independent; the total momentum enters in the singularities of the vertex functions coming from the mesonic sector. One first considers the chiral limit , which removes the first term of the right-hand side of Eq. (5) and transforms, in the limit , the contributions of the 2PGIQGFs into a term proportional to the anticommutator of with . If the latter contribution does not vanish in these limits, then it should be compensated by the left-hand side. This is possible only if the axial-vector vertex function possesses a pseudoscalar pole at . The residue of this pole is proportional to the wave function of the state and to the weak decay constant of the Goldstone boson. One thus finds the proportionality relation between the Goldstone boson wave function at and the scalar part of the 2PGIQGF.

In the next step, one considers the case where the quark masses are different from zero and takes the limit of to the physical mass of the would-be Goldstone boson. In this limit, only the two vertex functions of Eq. (5) survive. One then considers an expansion in the pseudoscalar vertex residue in terms of the quark masses, keeping only the leading contribution. Remembering that the Goldstone boson wave function is proportional to the scalar part of the 2PGIQGF, one ends up with the Gell-Mann–Oakes–Renner relation GellMann:1968rz ()

(15) |

where is the quark condensate in the chiral limit per flavor unit.

## 6 Two-dimensional QCD

Two-dimensional QCD in the large- limit is a simplified laboratory for real QCD ’tHooft:1973jz (); ’tHooft:1974hx (). It possesses the main characteristics of confinement and allows the study of the infrared behavior of the theory in a more explicit way. On the other hand, Wilson loop averages can be explicitly calculated: they satisfy the area law, due to the fact that the coupling constant of the theory is dimensionful from the start Kazakov:1980zi (); Kazakov:1980zj (); Bralic:1980ra ().

Equation (3) has also its counterpart in two dimensions. It turns out that the kernels for vanish due to the particular structure of the Wilson loop averages. ( vanishes in general for symmetry reasons.) One remains solely with the kernel . Equation (3) then reduces to a nonlinear nonlocal equation in . The equation can be solved explicitly and analytically by analyzing the singularity structure of the 2PGIQGF Sazdjian:2010ku ().

The solution is infrared finite with singularities in momentum space lying on the positive real axis of (timelike region). They are represented by an infinite number of branch points, characterized by positive masses ( with threshold singularities equal to in . The expression of is, in momentum space,

(16) |

where is the string tension. The masses are greater than the free quark mass and ordered according to increasing values. For massless quarks they remain positive. The masses and the coefficients , the latter being also positive, satisfy, for general , an infinite set of algebraic equations that are solved numerically. Their asymptotic values, for large values of such that , are

(17) |

The masses can be considered as dynamically generated quark masses, since they do not exist in the QCD Lagrangian. The gauge invariant treatment of the problem has displayed two main features: First, the dynamically generated masses are infinite in number. Second, they do not produce simple poles in the quark Green’s functions, but stronger singularities; this fact might prevent quarks from being observed as asymptotic free states.

The appearance of dynamically generated masses is also accompanied by a nonvanishing value of the quark condensate in the chiral limit. Since the quark condensate is a local gauge invariant quantity, it can also be calculated in various gauges. Calculations done in the light cone gauge Zhitnitsky:1985um (); Burkardt:1995eb (), in the axial gauge Li:1986gf () and in the present approach Sazdjian:2010ku (), give numerically the same result.

## 7 Conclusion

Use of a gauge invariant formalism has confirmed the general relationship that exists in QCD between the wave function of the pseudoscalar Goldstone boson of chiral symmetry breaking and the scalar part of the two-point gauge invariant quark Green’s function, in case the latter is nonvanishing in the chiral limit. The relationship also persists under truncation schemes of the exact interaction kernels.

Two-dimensional QCD in the large- limit, provides an explicit illustration of the above phenomenon. The dynamical quark mass generation produces an infinite number of branch point singularities, stronger than simple poles, which are infrared finite and describe the mechanism of chiral symmetry breaking in that theory.

### References

- Y. Nambu, G. Jona-Lasinio, Dynamical model of elementary particles based on an analogy with superconductivity. I, Phys. Rev. 122 (1961) 345.
- M. Baker, K. Johnson, B. W. Lee, Broken symmetries and zero-mass bosons, Phys. Rev. 133 (1964) B209.
- F. J. Dyson, The S matrix in quantum electrodynamics, Phys. Rev. 75 (1949) 1736.
- J. S. Schwinger, On the Green’s functions of quantized fields. 1., Proc. Nat. Acad. Sci. 37 (1951) 452.
- R. Alkofer, L. von Smekal, The infrared behavior of QCD Green’s functions, Phys. Rept. 353 (2001) 281.
- C. S. Fischer, Infrared properties of QCD from Dyson-Schwinger equations, J. Phys. G 32 (2006) R253.
- J. Finger, J. E. Mandula, J. Weyers, The pion in QCD, Phys. Lett. B96 (1980) 367.
- A. Le Yaouanc, L. Oliver, O. Pène, J.-C. Raynal, Spontaneous breaking of chiral symmetry for confining potentials, Phys. Rev. D 29 (1984) 1233.
- S. L. Adler, A. C. Davis, Chiral symmetry breaking in Coulomb gauge QCD, Nucl. Phys. B 244 (1984) 469.
- R. Alkofer, P. A. Amundsen, Chiral symmetry breaking in an instantaneous approximation to Coulomb gauge QCD, Nucl. Phys. B 306 (1988) 305.
- A. Bender, C. D. Roberts, L. von Smekal, Goldstone theorem and diquark confinement beyond rainbow ladder approximation, Phys. Lett. B 380 (1996) 7.
- G. Preparata, W. I. Weisberger, Ultraviolet divergences in radiative corrections to weak decays, Phys. Rev. 175 (1968) 1965.
- P. Maris, C. D. Roberts, P. C. Tandy, Pion mass and decay constant, Phys. Lett. B 420 (1998) 267.
- H. Sazdjian, Integral equation for gauge invariant quark two-point Green’s function in QCD, Phys. Rev. D 77 (2008) 045028.
- H. Sazdjian, Gauge invariant bound state equations for quark-antiquark systems in QCD, Phys. Rev. D 88 (2013) 025034.
- S. Mandelstam, Feynman rules for electromagnetic and Yang-Mills fields from the gauge independent field theoretic formalism, Phys. Rev. 175 (1968) 1580.
- Y. Nambu, QCD and the string model, Phys. Lett. B80 (1979) 372.
- K. G. Wilson, Confinement of quarks, Phys. Rev. D10 (1974) 2445.
- M. Gell-Mann, R. J. Oakes, B. Renner, Behavior of current divergences under , Phys. Rev. 175 (1968) 2195.
- G. ’t Hooft, A planar diagram theory for strong interactions, Nucl. Phys. B 72 (1974) 461.
- G. ’t Hooft, A two-dimensional model for mesons, Nucl. Phys. B 75 (1974) 461.
- V. A. Kazakov, I. K. Kostov, Nonlinear strings in two-dimensional gauge theory, Nucl. Phys. B 176 (1980) 199.
- V. A. Kazakov, Wilson loop average for an arbitrary contour in two-dimensional gauge theory, Nucl. Phys. B 179 (1981) 283.
- N. E. Bralic, Exact computation of loop averages in two-dimensional Yang-Mills theory, Phys. Rev. D 22 (1980) 3090.
- H. Sazdjian, Spectral properties of the gauge invariant quark Green’s function in two-dimensional QCD, Phys. Rev. D 81 (2010) 114008.
- A. R. Zhitnitsky, On chiral symmetry breaking in QCD in two-dimensions (), Phys. Lett. B165 (1985) 405.
- M. Burkardt, Trivial vacua, high orders in perturbation theory and nontrivial condensates, Phys. Rev. D 53 (1996) 933.
- M. Li, Large two-dimensional QCD and chiral symmetry, Phys. Rev. D 34 (1986) 3888.