Arxiv:0906.4487V1 [Hep-Th] 24 Jun 2009 Wilson Loops

Home , Ordered exponential, Wilson loop

arXiv:0906.4487v1 [hep-th] 24 Jun 2009 ne h ag transformation gauge the Under r omltoso ag hoisadt h ehdof method the mod- to in and 1 lattice role theories the central the gauge a of on play formulations main both which ern The continuum, loops, the Wilson in theories. to and gauge payed inves- other is nonperturbative and attention in QCD used of methods tigations the to duction ouso h isnlosaeeteeyueu nsolving models. in Ana- matrix useful of extremely way. kinds are gauge-invariant various loops the Wilson manifest the through of a reformulated logues in be of formulation loops can lattice Wilson QCD the in theories. role Wil- central gauge effect. a play Aharonov–Bohm loops the son by observable or theory are quantum Abelian loops in in Wilson factors theories. phase gauge essentially non-Abelian are loops Wilson loops Wilson I Lecture classical few a to references with only papers. text myself this In restrict references. I detail contains which [1] textbook ∗ h bla hs atrtasom as transforms factor phase Abelian the lcrncades [email protected] address: Electronic nteelcuentsIgv re eaoia intro- pedagogical brief a give I notes lecture these In o ute tde fti ujc a eomn the recommend can I subject this of studies further For bla hs atri endb h formula the by defined is factor phase Abelian /N expansion. .Dfiiinadbscproperties basic and Definition A. I HS ATR NQED IN FACTORS PHASE II. U speetd h he nrdcoylcue eegvnat given were lectures introductory three The presented. is ASnmes 11.g 11.a 12.38.Aw 11.15.Ha, 11.15.Pg, numbers: PACS 2009. 9–16, February Moscow, [Γ A U µ eaoia nrdcint isnlos atc ag t gauge lattice loops, Wilson to introduction pedagogical A yx ( .INTRODUCTION I. [Γ z ] ) yx −→ g −→ e = ] . g t . . t . re nrdcint isnLosadLarge and Loops Wilson to Introduction Brief A e A ı nttt fTertcladEprmna hsc,Moscow Physics, Experimental and Theoretical of Institute α µ ı ( e ( y z R ) Γ + ) U yx [Γ d 1 e z yx µ ∂ A e ] µ µ α ( − z ( ı z ) α . ) ( x , ) . uiMakeenko Yuri (1) (3) (2) na lcrmgei ed n ocmaepae fa of phases points compare at to and function wave field, electromagnetic an in aalltasotaogsm otu Γ contour some along transport parallel a n analogously and h point the sΓdpnetecp when except Γ-dependent is hrfr h hs atri rnfre steproduct the as transformed is ϕ factor phase the therefore sigfil strength field ishing o o ipycnetdsae teAharonov–Bohm (the spaces connected effect). simply not for rmgei ed oitoueeetoantcfield, electromagnetic introduce To field. tromagnetic h xoeti uttecascl(Euclidean) classical the just is exponent The st erpae ytecvratderivative covariant the by replaced be to is nerlrpeetto 9 o h rpgtro scalar a of path- propagator The the for field. (9) electromagnetic representation external integral an in particle ( h hs atrpasterl faprle transporter parallel a of role the plays factor phase The aefnto ttepoint the at function wave A e scnie unu atcei lsia elec- classical a in particle quantum a consider us Let aefnto ttepoint the at function wave A o h rpgtrw get we propagator the For y G ) ϕ ( † ,y x, × z z h 7hIE itrSho fPhysics, of School Winter ITEP 37th the ( µ µ x ( (0)= er n h 1 the and heory τ ): Z ; )= y .Poaaosi xenlfield external in Propagators B. A x fe utpiainb h hs factor: phase the by multiplication after y = ) µ µ D ∂ µ ϕ z U U µ † 2 1 ϕ →∇ −→ [Γ ( ( [Γ ( t y Z x e ) yx ) 0 yx ) ∞ U ] Russia , F − ] [Γ −→ d µ ϕ g µν /N 1 2 g ∼ . τ . ( t yx R t x . x ( = . 0 e τ epnino QCD of -expansion z ) ] − d ) eti uteisoccur subtleties Certain )). and “ e t A ϕ ∂ g g 1 2 ∗ ∼ ∼ z ı ˙ . . µ α τm N ( µ t 2 µ t . y . ( ( ( − x x t y x ) z 2 )+ı “ ) “ esol rtmake first should we , ϕ sapr ag (van- gauge pure a is ) ı ϕ ϕ stasomdas transformed is ϕ rnfrslk n at one like transforms eA † e † ( ( ( x ( y R x µ x 0 )” ) τ )” ( )” , d x , t . ITEP–TH–21/09 . ) z ˙ yx . µ ( t h result The . ) A action µ ( z ( t )) . fa of (7) (6) (5) (4) (8) (9) ∂ µ 2

The phase diﬀerence depends on (the real part of)

µ µ ıe R + dz Aµ(z) −ıe R − dz Aµ(z) e Γyx e Γyx ıe dzµA (z) ıe dσµν F ıeHS electron beam HΓ µ R µν ¨XXX = e = e = e , (12) ¨¨ XXz ¨ XXX ¨* solenoid XXX interference where the closed contour Γ is composed from Γ+ and ¨ ¨¨ yx ¨ ¨ plane Γ− . It does not depend on the shape of Γ+ and Γ− but H¨source H=06 ¨ xy yx yx s H ¨ depends only on HS — the magnetic flux through the H ¨ Hj ¨¨* solenoid. HH ¨ H ¨ electronH¨ beam III. YANG-MILLS THEORIES screen Modern theories of fundamental interactions are gauge theories. The principle of local gauge invariance was introduced by H. Weyl for the electromagnetic interaction FIG. 1: Scheme of the Aharonov–Bohm experiment. in analogy with general covariance in Einstein’s theory of gravitation. An extension to non-Abelian gauge groups was given by Yang and Mills in 1954 [3]. particle in an external electromagnetic field is due to A crucial role in gauge theories is played by the phase Feynman. factor which is associated with parallel transport in an We can alternatively rewrite Eq. (9) as external gauge field. The phase factors are observable in quantum theory, in contrast to classical theory. This is µ −S [Γyx]+ıe R dz Aµ(z) G(x, y; A) = e free Γyx , (10) analogous to the Aharonov–Bohm effect for the electromagnetic field. XΓyx where we represented the (parametric invariant) integral A. Gauge invariance over dt as the contour integral along the trajectory Γyx over The principle of local gauge invariance deals with the dzµ = dt z˙µ(t) . (11) gauge transformation (g.t.) of a matter field ψ, which is given by The transition amplitude of a quantum particle in a g.t. classical electromagnetic field is the sum over paths of ψ(x) ψ′(x) = Ω(x) ψ(x) . (13) −→ the Abelian phase factor (1). Here Ω(x) G with G being a semisimple Lie group which is called∈ the gauge group (G = SU(3) for QCD). C. Aharonov–Bohm effect Equation (13) demonstrates that ψ belongs to the fundamental representation of G. The unitary gauge group is when Transverse components of the electromagnetic field describe photons. Longitudinal components are related to Ω−1(x) = Ω†(x) , (14) gauging the phase of a wave function, i.e. permit one to compare its values at different space-time points when an while an extension to other Lie groups is straightforward. electron is placed in an external electromagnetic field. Then we have In quantum mechanics, the wave-function phase itself † g.t. ′ † † † is unobservable. Only the phase differences are observ- ψ (x) ψ (x) = ψ (x)Ω (x) . (15) −→ able, e.g. via interference phenomena. The phase differ- ence depends on the value of the phase factor for a given In analogy with QCD, the gauge group G = SU(N) is usually associated with color and the proper index of ψ path Γyx along which the parallel transport is performed. The phase factors are observable in quantum theory, in is called the color index. contrast to classical theory. This is seen in the Aharonov– The gauge transformation (13) of the matter field ψ Bohm effect [2] whose scheme is depicted in Fig. 1. can be compensated by a transformation of the non- Abelian gauge field which belongs to the adjoint rep- Electrons do not pass inside the solenoid where the mag- Aµ netic field is concentrated. Nevertheless, a phase differ- resentation of G: ence arises between the electron beams passing through g.t. (x) ′ (x) the two slits. The interference picture changes with the Aµ −→ Aµ value of the electric current. = Ω(x) (x)Ω†(x)+ıΩ(x) ∂ Ω†(x) . (16) Aµ µ 3

It is convenient to introduce the Hermitian matrix B. Non-Abelian phase factors (Wilson loops) [ (x)]ij = g Aa (x) [ta]ij , (17) Aµ µ Xa To compare phases of wave functions at distinct points, a non-Abelian extension of the parallel transporter is where g is the gauge coupling constant. needed. The proper extension of the Abelian formula (1): The matrices [ta]ij are the generators of G (a = 2 µ 1,...,N 1 for SU(N)) which are normalized such that ı RΓ dz Aµ(z) − U[Γyx] = P e yx , (29) tr tatb = δab, (18) includes the symbol P of path-ordering. where tr is the trace over the matrix indices i and j. Although the matrices µ(z) do not commute, the Quite often another normalization of the generators path-ordered exponential onA the RHS of Eq. (29) is de- ã˜b 1 ab with an extra factor of 1/2, tr t t = 2 δ , is used for fined unambiguously. This is obvious after rewriting the historical reasons, in particular tã = σa/2 for the SU(2) phase factor in an equivalent form group, where σa are the Pauli matrices. This results in µ τ µ ı R dz Aµ(z) ı R dt z˙ (t) Aµ(z(t)) the redefinition of the coupling constant,g ˜2 =2g2. P e Γyx = P e 0 . (30) Equation (17) can be inverted to give The path-ordered exponential in Eq. (29) can be un- a 1 a derstood as Aµ(x) = tr µ(x) t . (19) g A τ Substituting U[Γ ] = [1 + ı dt z˙µ(t) (z(t))] . (31) yx Aµ tY=0 Ω(x)= e ıα(x), (20) we obtain for an infinitesimal α: Using Eq. (11), Eq. (31) can also be written as g.t. adj µ δ µ(x) = α(x) . (21) U[Γ ] = [1 + ı dz (z)] . (32) A ∇µ yx Aµ Here z∈YΓyx adj If the contour Γ is discretized, then the non-Abelian µ α ∂µα ı [ µ, α] (22) yx ∇ ≡ − A phase factor is approximated by is the covariant derivative in the adjoint representation of G, while U[Γyx] fun M µ ψ ∂µψ ı µψ (23) µ zi + zi−1 ∇ ≡ − A = lim 1+ı(zi zi−1) µ , M→∞ − A 2 is that in the fundamental representation. It is evident Yi=1 that (33) adj fun µ B(x) = [ µ ,B(x)] , (24) ∇ ∇ which obviously reproduces (32) in the limit zi−1 zi. where B(x) is a matrix-valued function of x. The non-Abelian phase factor (29) is an element→ of the The QCD action is given in the matrix notation as gauge group G itself, while belongs to the Lie algebra Aµ S , ψ, ψ¯ of G. A Matrices are rearranged in inverse order under Hermi- 4 1 2 tian conjugation: = d x ψγ¯ µ (∂µ ı µ) ψ + mψψ¯ + tr , Z − A 4g2 Fµν U †[Γ ] = U[Γ ] . (34) (25) yx xy where The notation Γyx means the orientation of the contour from x to y, while Γxy denotes the opposite orientation = ∂ ∂ ı [ , ] (26) Fµν µAν − ν Aµ − Aµ Aν from y to x. These two result in opposite orders of mul- is the (Hermitian) matrix of the non-Abelian field tiplication for the matrices in the path-ordered product. strength. The phase factors obey the backtracking (zig-zag) con- This action is invariant under the local gauge transfor- dition mation since U[Γyx] U[Γxy] = 1 . (35) g.t. † µν (x) Ω(x) µν (x)Ω (x) (27) F −→ F The gauge field µ in the discretized phase factor (33) or is chosen at the centerA of the ith interval in order to g.t. δ µν (x) = ı [ µν (x), α(x)] (28) satisfy Eq. (35) at finite discretization. F − F Under the gauge transformation (16) U[Γyx] transfor the infinitesimal gauge transformation. forms as For the Abelian group G = U(1), the formulas recover g.t. those for QED. U[Γ ] Ω(y) U[Γ ]Ω†(x) . (36) yx −→ yx 4

This formula stems from the fact that gauge theory is a nonperturbative regularization of gauge theory and its nontrivial definition beyond perturbation µ g.t. µ ′ [1 + ı dz µ(z)] 1+ıdz µ(z) theory. It provides a nonperturbative quantization of A −→ A µ † gauge fields by a lattice. Lattice gauge theories uses = Ω(z + dz) [1 + ı dz µ(z)] Ω (z) (37) A an analogy between quantum field theory and statistical which can be proven by substituting Eq. (16), so that mechanics and offers a possibility of applying nonper- Ω†(z) and Ω(z) cancel in the definition (32) at the inter- turbative methods: the strong-coupling expansion or the mediate point z. numerical Monte Carlo method to QCD, which provide A consequence of Eq. (36) is that ψ(x), transported by evidence for quark confinement. the matrix U[Γyx] to the point y, transforms under the gauge transformation as ψ(y): IV. LATTICE FORMULATION g.t. U[Γyx] ψ(x) “ψ(y)”. (38) ∼ A. The lattice

Therefore, U[Γyx] is, indeed, a parallel transporter. ¯ It follows from these formulas that ψ(y) U[Γyx] ψ(x) is A lattice approximates continuous space by a discrete gauge invariant: set of points. In Euclidean formulation the lattice is along all four coordinates, while the time is left continuous in g.t. ψ¯(y) U[Γyx] ψ(x) ψ¯(y) U[Γyx] ψ(x) . (39) Hamiltonian approach. −→ The lattice is defined as a set of points of d-dimensional Another consequence of Eq. (36) is that the trace of the Euclidean space with coordinates phase factor for a closed contour Γ is gauge invariant: xµ = nµa , (41) µ µ ı H dz Aµ(z) g.t. ı H dz Aµ(z) tr P e Γ tr P e Γ . (40) −→ where the components of the vector This is quite similar to the Abelian phase factor. The sufficient and necessary condition for the phase nµ = (n1,n2,...,nd) (42) factor to be independent on a local variation of the path is are integer numbers. The points (41) are called the lattice the vanishing of . The formulas of this type are well- Fµν sites. known in differential geometry where parallel transport The dimensional constant a, which equals the distance around a small closed contour determines the curvature. between the neighboring sites, is called the lattice spac- µν in Yang–Mills theory is the proper curvature in an ing. Dimensional quantities are measured in units of a, internalF color space while is the connection. Aµ thereby setting a = 1. A 2d lattice with periodic boundary conditions is depicted in Fig. 2. Sites with the same numbers are iden- A historical remark tified. The spatial size of the lattice has L1 = 6 and L2 = 4. An analogous 4d lattice is called hypercubic. An analog of the phase factor was first introduced by The next concepts are the link and plaquette of a lat- H. Weyl in 1919 [4] in his attempt to describe gravi- tice which are shown in Fig. 3. A link l = x; µ con- tational and electromagnetic interactions of electron on nects two neighboring sites x and x + aµˆ, where{ }µ ˆ is a equal footing. What he did is associated in modern lan- unit vector along the µ-direction (µ = 1,...,d). A pla- guage with the scale rather than the gauge transforma- quette p = x; µ, ν is the elementary square enclosed by tion, i.e. the vector-potential was not multiplied by ı as { } in Eq. (1). This explains the term “gauge invariance” – gauging literally means fixing a scale. 1 ppppppp23456 1 The factor of ı was inserted by London in 1927 [5] after creation of quantum mechanics and the recognition ppppppp that the electromagnetic interaction corresponds to the 7 7 6 freedom of choice of the phase of a wave function and not a ppppppp?- to a scale transformation. 8 x a 8 ppppppp 9 9 Lecture II ppppppp Lattice gauge theories 1 2 3 4 5 6 1

Lattice gauge theories were proposed by Wilson in FIG. 2: Two-dimensional lattice with periodic boundary con- 1974 [6] to explain quark conﬁnement in QCD. Lattice ditions. 5

x + aνˆ x + aµˆ + aνˆ where the integral is along the link x; µ . As a 0, this yields { } → d d U (x) e ıaAµ(x) (46) p µ →

l so that Uµ(x) is expressed via the exponential of the µth xd xd+ aµˆ xd xd+ aµˆ component of the vector potential at the center of the link. Since the path-ordered integral in Eq. (45) depends on FIG. 3: Link (left) and plaquette (right). the orientation, links are oriented. The same link, which connects the points x and x+aµˆ, can be written either as pppppp ppppppppppp x; µ or as x + aµˆ; µ . The orientation is positive for ppppppppppp {µ > 0} in the{ former− case} (i.e. the same as the direction pppppp ppppppppppp ppppppppppp of the coordinate axis) and is negative in the latter case. pppppp ppppppppppp The link variable Uµ(x) is assigned to links with posi- ppppppppppp tive orientations. The U-matrices which are assigned to pppppp ppppppppppp links with negative orientations are given by

† (a) (b) U−µ(x + aµˆ) = Uµ(x) . (47)

This is a one-link analog of Eq. (34). FIG. 4: Description of continuum configurations by lattices. It is clear from the relation (45) between the lattice and continuum gauge variables how to construct lattice phase factors (Wilson loops) – to construct the contours four links in the directions µ and ν. The set of four links from the links of the lattice. which bound the plaquette p is denoted as ∂p. For an infinite lattice, the number of degrees of freedom = (but enumerable). To limit the this number, ∞ C. The Wilson action the lattice has a finite size L1 L2 Ld in all directions. Periodic boundary conditions× ×···× are usually imposed to reduce finite-size effects. An important role in the lattice formulation is played by the phase factor for the simplest closed contour on the lattice: the (oriented) boundary of a plaquette. The pla- B. Matter and gauge fields on the lattice quette variable is composed from the link variables (45) as

Matter field, say a quark field, is attributed to the † † U(∂p) = U (x) U (x + aνˆ) U (x + aµˆ) U (x) . (48) lattice sites, therefore a continuous field ϕ(x) is approxi- ν µ ν µ mated by its values at the lattice sites The link variable transforms under the gauge transformation as ϕ(x) = ϕ . (43) ⇒ x g.t. † Uµ(x) Ω(x + aµˆ) Uµ(x)Ω (x) , (49) The lattice field ϕx is a good approximation of a contin- −→ uous field ϕ (x) when a than the characteristic size of where the matrix Ω(x) is attributed to the lattice sites. a given configuration. A≪ description of continuum field This defines the lattice gauge transformation. configuration by lattices is illustrated by Fig. 4. Lattice The plaquette variable transforms under the lattice (a) is “coarse” and can represent the given continuum gauge transformation as field configuration very roughly, while lattice (b) is “fine” and its spacing is small enough. g.t. U(∂p) Ω(x) U(∂p)Ω†(x) . (50) The gauge field is attributed to the links of the lattice: −→ (x) = U (x) . (44) Its trace over the color indices is gauge invariant: Aµ ⇒ µ g.t. It looks natural since a link is characterized by a coor- tr U(∂p) tr U(∂p) . (51) −→ dinate and a direction – the same as (x). Sometimes Aµ the notation Ux,µ is used as an alternative for Uµ(x) to This invariance is used in constructing an action of a emphasize that it is attributed to links. lattice gauge theory. The simplest Wilson action is The link variable Uµ(x) can be viewed as 1 x+aµˆ µ Slat[U] = 1 Re tr U(∂p) . (52) ı Rx dz Aµ(z) − N Uµ(x) = P e , (45) Xp 6

The summation is over all the elementary plaquettes of Haar measure under multiplication by an arbitrary group the lattice (i.e. over all x, µ, and ν), regardless of their element from the left or from the right: orientations. ′ Since a reversal of the orientation results in complex dU = d(ΩU) = d(UΩ ) , (60) conjugation: guarantees the gauge invariance of the partition function (58). tr U(∂p) reor. tr U †(∂p) = [tr U(∂p)]∗, (53) −→ This invariance of the Haar measure is crucial for the Wilson formulation of lattice gauge theories. the action can be rewritten as the sum over oriented pla- An explicit expression for the Haar measure can be quettes written for the SU(2) gauge group, whose element is parametrized by a unit four-vector a : 1 1 µ S [U] = 1 tr U(∂p) . (54) lat 2 − N U = a I + ı~a~σ a2 =1 , (61) orientX p 4 µ where ~σ are the Pauli matrices. The Haar measure for As a 0, the lattice action (52) becomes the action SU(2) is of continuum→ gauge theory. To show this, we note that 4 2 3 1 (1) 2 U(∂p) exp ıa (x)+ a , (55) dU = daµ δ a 1 , (62) → Fµν O π2 µ − µY=1 where (x) is non-Abelian field strength (26). Fµν 2 In the Abelian theory, the expansion (55) is easily since det U = aµ. found from the Stokes theorem. The commutator of The partition function (58) characterizes vacuum effects in quantum theory. Physical quantities are given µ(x) and ν (x), which arises in the non-Abelian case, AcomplementsA the field strength to the non-Abelian one, by the averages as is ensured by the gauge invariance. −1 −βS[U] F [U] = Z (β) dUµ(x)e F [U] , (63) A transition to the continuum limit is performed as h i Z Yx,µ 4 a→0 1 4 a d x . (56) where F [U] is a gauge-invariant functional of Uµ(x). · · · −→ 2 Z · · · Xp Xµ,ν These become the expectation values in the continuum theory as a 0 and β = N/g2. Expanding the RHS of Eq. (55) in a, we reproduce the → continuum action V. WILSON LOOPS ON A LATTICE a→0 1 4 2 Slat d x tr (x) . (57) −→ 4N Z Fµν Xµ,ν Lattice phase factors are associated with paths drawn on the lattice. To write down the phase factor on the lattice via the D. The Haar measure link variables, let us specify the (lattice) contour C by its initial point x and by the directions (some may be The partition function of a pure lattice gauge theory negative) of the links forming the contour is defined by C = x; µ ,...,µ . (64) { 1 n} −βS[U] Z(β) = dUµ(x) e , (58) The lattice phase factor U(C) is given by Z Yx,µ U(C) = U (x + aµˆ + + aµˆ − ) µn 1 · · · n 1 · · · U (x + aµˆ ) U (x) . (65) where the action is given by Eq. (52). × µ2 1 µ1 This is the analog of a partition function in statistical For the links with a negative direction, it is again conve- mechanics at an inverse temperature β given by nient to use Eq. (47). N A closed contour hasµ ˆ1 + +µ ˆn = 0. The gauge β = . (59) invariant trace of the phase factor· · · for a closed contour is g2 called the Wilson loop. This formula results from comparing Eq. (57) with the The average of the Wilson loop is determined to be gauge-field part of the continuum action (25). 1 W (C) tr U(C) A subtle question is what is the measure dUµ(x) in ≡ N Eq. (58). To preserve the gauge invariance at finite lattice −1 −βS[U] 1 spacing, the integration is over the Haar measure which = Z (β) dUµ(x)e tr U(C). (66) Z N is an invariant group measure. The invariance of the Yx,µ 7

(R, 0) (R,T ) sum over states and ﬁnally we ﬁnd large T W (R T ) e −E0(R)T . (72) ? × −→ 6 - Since nothing in the derivation relies on the lattice, it (0, 0) (0,T ) holds for a rectangular loop in the continuum theory as well. VI. STRONG-COUPLING EXPANSION FIG. 5: Rectangular Wislon loop.

Lattice path integrals can be calculated either by per- This average is often called the Wilson loop average. turbation theory in g2 or by an expansion in β 1/g2. ∝ A very important role is played by the Wilson loop This is called the strong-coupling expansion in analogy of averages for rectangular contours. Such a contour lying the high-temperature expansion in statistical mechanics in the (x, t)-plane is depicted in Fig. 5. This rectangular since β is the analog of inverse temperature. Wilson loop is of the size R T . To perform the strong-coupling expansion, we expand × The Wilson loop average is related for T R to the exponential of the lattice action in β. Then the prob- ≫ the interaction energy of the static (i.e. inﬁnitely heavy) lem is to calculate the integrals over unitary group: quarks, separated by a distance R, by the formula

i1···im,k1···kn i1 im †k1 †kn T ≫R −E (R)T I = dUU U U U , (73) W (R T ) e 0 . (67) j1···jm,l1···ln j1 jm l1 ln × ∝ Z · · · · · · It can be proved in the axial gauge = 0, where A4 where the Haar measure is normalized U4(x) = 1 so that only vertical segments contribute to U(R T ). Denoting × dU = 1 . (74) R Z ı R dz1 A1(z1,...,t) Ψij (t) P e 0 , (68) ≡ h iij The integral (73) is nonvanishing only if n = we have m (mod N), i.e. if n = m + kN with integer k. 1 For the simplest case of m = n = 1, the answer can W (R T ) = trΨ(0)Ψ†(T ) . (69) × N easily be found by using the unitarity of U and the or- thogonality relation Inserting a sum over a complete set of intermediate states k 1 dUU i U † = δiδk . (75) Z j l N l j n n = 1 , (70) | ih | Xn The simplest Wilson loop average is that for the loop we obtain which coincides with the boundary of a plaquette. It is called the plaquette average 1 W (R T ) = Ψ (0) n n Ψ† (T ) × N ij ji Xn 1 W (∂p) = tr U(∂p) . (76) 1 2 N = Ψ (0) n e −EnT , (71) N ij Xn To calculate the plaquette average to order β, it is where En is the energy of the state n . As T , only suﬃcient to retain only the terms (β) in the expansion the ground state with the lowest energy| i survives→ ∞ in the of the exponentials: O

1 ′ 1 dUµ(x) 1+ β Re tr U(∂p ) tr U(∂p) Z N N Yx,µ h Xp′ i W (∂p)= + β2 . (77) 1 ′ O dUµ(x) 1+ β Re tr U(∂p ) Z N Yx,µ h Xp′ i 8

- ′ ∂p ? ? 6 6∂p qqq q -

FIG. 6: Boundaries of plaquettes with opposite orientations.

(a) (b) - - - ? ? ?6 FIG. 8: Lines of force between static quarks. 6 6 6 ------? ? ? ? ? ? ? 6 6 6 6 6 6 the minimal surface is a piece of the plane: ------6 ? ? ? ? ? ? W (R T ) = [W (∂p)]R T (81) 6 6 6 6 6 6 × - to the leading order in β. Contribution of more complicated surfaces, do not ly- FIG. 7: Filling of a loop with plaquettes. ing in the plane of the rectangle, is W (C) βarea and is, therefore, suppressed since their areas are∼ larger than The group integration can then be performed by remem- Amin. bering that VII. AREA LAW AND CONFINEMENT i † k 1 i k dUµ(x) [Uµ(x)] [U (y)] = δxy δµν δ δ (78) Z j ν l N l j The exponential dependence of the Wilson loop aver- at different links. age on the area of the minimal surface (as in Eq. (80)) is Using this property of the group integral in Eq. (77), called the area law. It is customarily assumed that if an we immediately see that the denominator is equal to 1 area law holds for loops of large area in the pure SU(3) (each link is encountered no more than once), while the gauge theory then quarks are confined. In other words, only nonvanishing contribution in the numerator is from there are no physical in or out quark states. This is the plaquette p′, which coincides with p but has the oppo- the essence of Wilson’s| confinementi h | criterion. The argu- site orientation as is depicted in Fig. 6. The boundaries ment is that physical amplitudes (e.g., the polarization of the plaquettes p and p′ have opposite orientations ∂p operator) do not have quark singularities when the Wil- and ∂p′, respectively. son criterion is satisfied. The final answer for the plaquette average is

β W (∂p) = for SU(N) with N 3 , A. Linear potential 2N 2 ≥  (79) β  W (∂p) = for SU(2) .  A justification for the Wilson criterion is based on the 4 relationship (67) between the Wilson loop average and   the potential energy of interaction between static quarks. The result for SU(2) differs by a factor of 1/2 because When the area law tr U(∂p) is real for SU(2) so that the orientation of the large C plaquettes can be ignored. W (C) e −KAmin(C) (82) According to Eq. (78), a nonvanishing result emerges −→ only when plaquettes, arising from the expansion of the holds for large loops, the potential energy is a linear func- exponentials of Eq. (66) in β, completely cover a surface tion of the distance between the quarks: enclosed by the given loop C as is shown in Fig. 7. In this case each link is encountered twice (or never), once in the E(R) = KR . (83) positive direction and once in the negative direction, so that all the group integrals are nonvanishing. Lines of force between static quarks are depicted in The leading order in β corresponds to filling a minimal Fig. 8 for (a) linear and (b) Coulomb interaction poten- surface: tials. For the linear potential the lines of force are contracted into a tube, while they are distributed over the W (C) = [W (∂p)]Amin(C) , (80) whole space for the Coulomb one. The coefficient K is called the string tension because where W (∂p) is given by Eq. (79) and Amin(C) is the area the gluon field between quarks contracts to a tube or (in units of a2) of the minimal surface. For a rectangle string, whose energy is proportional to its length. The 9

2 value of K is the energy of the string per unit length. ln a K 6 This string is stretched with the distance between quarks and prevents them from moving apart to macroscopic strong coupling distances. Equation (81) gives 1 2N 2 1 K = ln = ln 2Ng2 (84) J a2 β a2 J asymptotic scaling for the string tension to the leading order of the strong- J J - coupling expansion. The next orders of the strong- 0 J 1/g2 coupling expansion result in corrections in β to this for- J mula. Conﬁnement holds in the lattice gauge theory to J any order of the strong-coupling expansion. J

FIG. 9: String tension versus 1/g2 B. Asymptotic scaling .

Equation (84) establishes the relationship between lat- For such values of g2, where asymptotic scaling holds, tice spacing a and the coupling g2. Let K equals its the lattice gauge theory has a continuum limit. experimental value The knowledge of the two asymptotic behaviors says K = (400 MeV)2 1 GeV/fm (85) nothing about the behavior of a2K in the intermediate ≈ region of g2 1. There can be either a smooth transition ′ which results from the slope of the Regge trajectory α = between these∼ two regimes or a phase transition. Numer- ′ −2 1/2πK. The slope α = 1 GeV from the ρ – A2 – g ical methods were introduced to study this problem. trajectory. The renormalizability prescribes that variations of a, which plays the role of a lattice cutoff, and of the bare C. Relation to second-order phase transition charge g2 should be made simultaneously in order that K does not change. Given Eq. (84), this procedure calls for a as Continuum limits of a lattice system are reached at g2 . In other words, the lattice spacing is→ large ∞ in the points of second-order phase transitions when the the→ strong-coupling ∞ limit, compared with 1 fm. Such a correlation length becomes infinite in lattice units. coarse lattice cannot describe the continuum limit and, The correlation length is inversely proportional to in particular, the rotational symmetry. ΛQCD: In order to pass to the continuum, the lattice spacing dg2 a should be decreased. Equation (84) shows that a de- ξ Λ−1 = a exp , (87) creases with decreasing g2. However, this formula ceases ∼ QCD Z (g2) to be applicable in the intermediate region of g2 1 and, B therefore, a 1 fm. ∼ where (g2) is the Gell-Mann–Low function. The only To further∼ decrease a, we further decrease g2. While chanceB for the RHS to diverge is to have a zero of (g2) 2 2 B no analytic formulas are available at intermediate values at some fixed point g = g∗. Therefore, the bare coupling 2 2 2 of g , the expected relation between a and g for small should approach the fixed-point value g∗ to describe the g2 is predicted by the known two-loop Gell-Mann–Low continuum. 2 function of QCD. (0) = 0 for a non-Abelian gauge theory so that g∗ =0 For pure SU(3) Yang-Mills, Eq. (84) is replaced at is aB fixed-point value of the coupling constant. Therefore, small g2 by the continuum limit is associated with g2 0. → 102 2 121 1 8π 2 2 K = const e −8π /11g , (86) · a2 11g2 Lecture III where the two-loop Gell-Mann–Low function is used. 1 The exponential dependence of K on 1/g2 is called /N Expansion asymptotic scaling. Asymptotic scaling sets in for some value of 1/g2. The dependence of the string tension on An effective coupling constant of QCD becomes large at 1/g2 is shown in Fig. 9. The strong-coupling formula large distances, so fluctuations of scales of different orders (84) holds for small 1/g2. The asymptotic-scaling for- of magnitude are essential and there is no small param- mula (86) sets in for large 1/g2. Both formulas are not eter. ’t Hooft proposed in 1974 [7] to use the number of applicable in the intermediate region 1/g2 1. colors N of the gauge group SU(N) as such a parameter ∼ 10 and to perform an expansion in 1/N. The motivation Each line, often termed the index line, represents the was the 1/N-expansion in statistical mechanics. Kronecker delta-symbol and has an orientation which is The 1/N-expansion of QCD rearranges perturbation indicated by arrows. This notation is obviously consis- theory in a way consistent with a string picture. The tent with the space-time structure of the propagator that accuracy of large-N QCD is of the order of the ratios of describes a propagation from x to y. ij meson widths to their masses (10–15%). While QCD is Arrows are a result of the fact that the matrix Aµ is simplified in the large-N limit, it is not yet solved. Hermitian and its off-diagonal components are complex conjugate. Double lines appear generically in all models describing VIII. INDEX OR RIBBON GRAPHS matrix fields in contrast to vector (in internal symmetry space) fields, whose propagators are depicted by single In order to describe the 1/N-expansion of QCD, it is lines. convenient to use the matrix-field representation The three-gluon vertex is depicted in the double-line notation as ij a a ij [Aµ(x)] = Aµ(x) [t ] . (88) j1 i1 i1 j1 ...... Xa ...... i2 ...... j3 j2 ...... i3 ...... − ...... The matrix (88) is Hermitian and differs from (17) by g...... ij . . . . The propagator of the matrix field A (x) reads j2 i3 i2 j3 i j i j i j i j i j i j 1 g δ 1 3 δ 2 1 δ 3 2 δ 1 2 δ 2 3 δ 3 1 , (94) Aij (x) Akl(y) = δilδkj δij δkl D (x y) , ∝ − µ ν Gauss − N µν − where the subscripts 1, 2 or 3 refer to each of the three (89) gluons. The relative minus sign arises from the commuta- where we have assumed, as usual, a gauge-fixing to de- tor in the cubic-in-A term in the action (25). The color fine the gluon propagator in perturbation theory. For 1 part is antisymmetric under an interchange of gluons. instance, one has The (momentum space) space-time part 1 δµν γ (p ,p ,p )= δ (p p ) Dµν (x y) = (90) µ1µ2µ3 1 2 3 µ1µ2 1 2 µ3 − 4π2 (x y)2 − − +δµ2µ3 (p2 p3) + δµ1µ3 (p3 p1) (95) − µ1 − µ2 in the Feynman gauge. is also antisymmetric. We consider all three gluons as Equation (89) can be derived immediately from the incoming so their momenta obey p +p +p = 0. The full standard formula 1 2 3 vertex is symmetric as is prescribed by Bose statistics. Aa (x) Ab (y) = δabD (x y) (91) The four-gluon vertex involves six terms – each of them µ ν Gauss µν − is depicted by a cross – which differ by interchanging of multiplying by the generators of the SU(N) gauge group the color indices. We depict the color structure of the according to the definition (88) and using the complete- four-gluon vertex for simplicity in the case when i1 = ness condition j2 = i, i2 = j3 = j, i3 = j4 = k, i4 = j1 = l, but i,j,k,l take on different values. Then only the following term is 2− N 1 1 left: (ta)ij (ta)kl = δilδkj δij δkl for SU(N) . − N l i ... Xa=1 ...... (92) l ...... i ...... g2 , Alternatively, Eq. (89) can be derived directly from a ...... k... j ∝ ...... path integral over the matrix fields...... k j (96) We concentrate only on the structure of diagrams in the index space, i.e. the space of the indices associated where there are no delta-symbols since the color structure with the SU(N) group. We shall not consider, in most is fixed. We pick up only one color structure by equating cases, space-time structures of diagrams which are pre- indices pairwise. scribed by Feynman’s rules. Diagrams of perturbation theory can now be com- Omitting at large N the second term in parentheses pletely rewritten in the double-line notation. The sim- on the RHS of Eq. (89), we depict the propagator by the plest one describing the one-loop correction to the gluon double line propagator is depicted in Fig. 10. The sum over the N i - l indices is associated with the closed index line. The con- Aij (x) Akl(y) δilδkj = . (93) µ ν Gauss ∝ j k tribution of this diagram is g2N 1. In order for the large-N limit∼ to be∼ nontrivial, the bare coupling constant g2 should satisfy 1 Here and in Eq. (95) δµν is due to Euclidean metric and is to be 1 g2 . (97) substituted by −gµν in Minkowski space. ∼ N 11

......

FIG. 10: Double-line representation of a one-loop diagram for FIG. 11: Double-line representation of a four-loop diagram. gluon propagator.

...... This dependence on N is also prescribed by the ...... asymptotic-freedom formula ...... 2 ...... 12π ...... g2 = (98) 11N ln (Λ/Λ ) QCD FIG. 12: Double-line representation of a nonplanar diagram. of the pure SU(N) gauge theory. Thus, the contribution of the diagram in Fig. 10 is of A. ’t Hooft topological expansion order g2N 1 in the large-N limit. The∼ double lines∼ can be viewed as bounding a piece of a plane. These lines represent a two-dimensional object. Let us consider a typical diagram for the gluon propa- In mathematics these double-line graphs are often called gator depicted in Fig. 11. The sum over the N indices is ribbon graphs or fatgraphs. They are connected with associated with each of the four closed index lines, whose Riemann surfaces. number is equal to the number of loops. The contribution of this diagram is g8N 4 1. Diagrams of this∼ type, which∼ can be drawn on a sheet Remark on the U(N) gauge group of paper without crossing any lines, are called planar diagrams. For such diagrams, the addition of a loop in- evitably results in the addition of two three-gluon (or one The double-line representation of the diagrams holds, four-gluon) vertices. A planar diagram with n loops has strictly speaking, only for the U(N) gauge group, whose 2 n closed index lines. It is of order generators 2 2 n2 n2-loop planar diagram g N 1, (102) A a √ A B AB 2 ∼ ∼ T = t ,I/ N , tr T T = δ A =1,...,N . (99) so that all planar diagrams survive in the large-N limit. obey the completeness condition Let us now consider a nonplanar diagram of the type depicted in Fig. 12. The diagram has six three-gluon N 2 vertices but only one closed index line (although it has ij kl T A T A = δilδkj for U(N) . (100) three loops!). The order of this diagram is g6N 2 ∼ ∼ AX=1 1/N . This nonplanar diagram can be drawn without line- Elements of both SU(N) and U(N) can be represented crossing on a surface with one handle (or hole) which in in the form mathematics is called a torus or a surface of genus one. A plane is then equivalent to a sphere and has genus zero. ıB U = e , (101) A general Riemann surface with h holes has genus h. The above evaluations of the order of the diagrams can where B is a general Hermitian matrix for U(N) and a be described by the single formula traceless Hermitian matrix for SU(N). The large-N limit of both the U(N) and SU(N) groups 1 genus genus-h diagram . (103) is the same. ∼ N 2

The expansion in 1/N rearranges perturbation-theory di- IX. PLANAR AND NONPLANAR GRAPHS agrams according to their topology as demonstrated in 1974 by ’t Hooft [7]. It is referred to as the topological The double-line representation of perturbation-theory expansion or the genus expansion. diagrams is very convenient to estimate their orders in Only planar diagrams associated with genus zero sur- 1/N. Each three- or four-gluon vertex contributes a fac- vive in the large-N limit. The problem of summing tor of g or g2, respectively. Each closed index line con- the planar graphs is complicated but simpler than that tributes a factor of N, while g2 1/N. of summing all the graphs, since the number of planar ∼ 12 graphs with n0 vertices grows geometrically at large n0: planar.

# (n ) no of planar graphs constn0 , (104) p 0 ≡ ∼ as shown by Tuttle [8] and Koplik, Neveu, Nussinov [9], B. Topological expansion (continued) while the total number of graphs grows factorially with n0. There is no dependence in Eq. (104) on the number of external lines of a planar graph which is assumed to Equation (103) holds, strictly speaking, only for the be much less than n0. gluon propagator, while the contribution of all planar di- There is a big diﬀerence between the planar diagrams agrams to a connected n-point Green function is gn−2, and the ladder diagrams which describe e+e− elastic scat- which is its natural order in 1/N. The three-gluon∼ Green tering in QED. For the ladder with n rungs, there are n! function is g, the four-gluon one is g2 and so on. The ladder diagrams, but only one of them is planar. This contributions∼ of all planar diagrams are∼ of the same order shows why the number of planar graphs is much smaller 1 in the large-N limit, independently of the number of than the total number of graphs, most of which are non- ∼external lines, for the Wilson loop average

∞ x1 xn−1 µ 1 ıg H dx Aµ(x) n µ1 µ2 µn (n) tr P e Γ = ı dx dx . . . dx Gµ ···µ (x1,...,xn) (105) N I 1 Z 2 Z n 1 n nX=0 Γ x1 x1

with shrunk to a point. Then a bounded piece of a plane will n be topologically equivalent to a sphere with a puncture. (n) g Gµ1···µn (x1,...,xn) tr [Aµ1 (x1) Aµn (xn)] . We draw planar diagrams in a plane with an extended ≡ N h · · · i (106) boundary (boundaries) rather than in a sphere with a The factor of 1/N, which normalizes the trace, provides puncture (punctures). The closed boundary is associated the natural normalization. The ordering along a closed with the trace over the color indices of the multi-point path implies cyclic-ordering in the index space as de- Green function picted in Fig. 13, where we omit the arrows for simplicity. The boundary represents the Wilson loop = a trajec- This diagram has n0 = 10 vertices, n1 = 12 gluon prop- tory of a heavy quark in the fundamental representation. agators, n2 = 4 closed index lines, and B = 1 boundary. The color indices of the external lines are contracted by the Kronecker delta-symbols (represented by the single C. Topological expansion and quark loops lines) in a cyclic order. The extra factor of 1/N arises from the normalization. The order in 1/N of the diagram It is easy to incorporate quarks in the topological ex- in Fig. 13 is 1/N 2 in accord with Eq. (103). ∼ pansion. A quark ﬁeld belongs to the fundamental repre- Analogously, the color indices in Eq. (106) are con- sentation of the gauge group SU(N) and its propagator tracted in the cyclic order. The delta-symbols, which is represented by a single line contract the color indices, are depicted by the single lines. They can be viewed as a boundary of the diagram. The ψ ψ¯ δ = i - j . (107) actual size of the boundary is not essential – it can be i j ∝ ij

The arrow indicates, as usual, the direction of propa-

...... gation of a (complex) ﬁeld ψ. These arrows are often ...... omitted for simplicity...... The diagram for the gluon propagator which involves ...... one quark loop is depicted in Fig. 14a. It involves one ...... quark loop and has no closed index lines so that its order ...... 2 6 2 ...... is g 1/N. The diagram in Fig. 14b is g N ...... ∼ ∼ ∼ ∼ ...... 1/N analogously...... It is evident from this consideration that quark loops ...... are not accompanied by closed index lines. One should ...... add a closed index line for each quark loop in order for ...... a given diagram with L quark loops to have the same double-line representation as for pure gluon diagrams. FIG. 13: Generic double-line index diagram. Therefore, given Eq. (103), diagrams with L quark loops 13

...... associated with a given graph, is open and has B + L ...... boundaries. This surface can be closed by attaching a ...... cap to each boundary. The single lines then become dou- ...... ble lines together with the lines of the boundary of each ...... cap. We have already considered this procedure when deducing Eq. (108) from Eq. (103). (a) (b) The number of faces for a closed Riemann surface constructed in such a manner is n2 + L + B, while the num- FIG. 14: Diagrams for gluon propagator which involve quark ber of edges and vertices are n1 and n0, respectively. loop. Euler’s theorem states that

χ 2 2h = n2 + L + B n1 + n0 . (113) are suppressed at large N by ≡ − − Therefore the RHS of Eq. (112) can be rewritten as 1 L+2·genus L quark loops . (108) N n2−n1+n0−B = N 2−2h−L−2B. (114) ∼ N The single-line representation of the quark loops is sim- We have thus proven that the order in 1/N of a generic ilar to that of the Wilson loop. Such a diagram emerges graph does not depend on its order in the coupling con- in gluon corrections to the vacuum expectation value of stant and is completely expressed via the genus h and the the quark operator number of virtual quark loops L and external boundaries B by 1 ¯ O = ψψ , (109) 1 2h+L+2(B−1) N generic graph . (115) ∼ N where the factor of 1/N is introduced to make it (1) in O the large-N limit. Therefore, the external boundary can For B = 1, we recover Eqs. (103) and (108). be viewed as a single line associated with valence quarks.

E. ’t Hooft versus Veneziano limits D. The proof of topological expansion In QCD there are several species or flavors of quarks To prove Eqs. (103) and its quark counterpart (108), (u-, d-, s- and so on). We denote the number of flavors let us consider a generic diagram in the index space which by Nf and associate a Greek letter α or β with a flavor (3) has n0 three-point vertices (either three-gluon or quark– index of the quark field. (4) The quark propagator then has the Kronecker delta- gluon ones), n four-gluon vertices, n1 propagators (ei- 0 symbol with respect to the flavor indices in addition to ther gluon or quark ones), n2 closed index lines, L virtual quark loops and B external boundaries. Its order in 1/N Eq. (107): is α ¯β αβ ψi ψj δ δij . (116) 1 n(3)+2n(4) n n −n(3)/2−n(4)−B ∝ g 0 0 N 2 N 2 0 0 (110) D E N B ∼ Their contraction results in B as has already been explained. The extra factor of 1/N Nf arises from the extra normalization factor of 1/N in op- δαα = Nf . (117) erators associated with external boundaries. αX=1 The number of propagators and vertices are related by Therefore, an extra factor of Nf corresponds to each (3) (4) closed quark loop for the Nf flavors. 2n1 = 3n0 +4n0 , (111) The limit when Nf is fixed as N is called the since three- and four-point vertices emit three or four ’t Hooft limit. Only valence quarks→ are ∞ then left (the propagators, respectively, and each propagator connects quenched approximation). In order for a meson to decay two vertices. We then rewrite the RHS of (110) as into other mesons built out of quarks, a quark–antiquark pair must be produced out of the vacuum. Consequently, n −n(3)/2−n(4) −B n −n +n −B N 2 0 0 = N 2 1 0 , (112) the ratios of meson widths to their masses are

Γtotal Nf where n = n(3) + n(4) is the total number of vertices. (118) 0 0 0 M ∼ N The exponent on the RHS of Eq. (112) can be expressed via the Euler characteristic χ of a given graph in the ’t Hooft limit. The ratio on the LHS of Eq. (118) of genus h. An appropriate Riemann surface, which is is 10–15% experimentally for the ρ-meson. The hope of 14 solving QCD in the ’t Hooft limit is the hope to describe Here Γ denotes one of the combination of the γ-matrices: QCD with this accuracy. An alternative large-N limit of QCD, when Nf N 1 ∼ Γ = I, γ , γ , ıγ γ , Σ = [γ ,γ ] ,... . (124) as N , was proposed by Veneziano in 1976 [10]. A 5 µ µ 5 µν 2ı µ ν general→ diagram ∞ with L quark loops will contribute

L genus The factorization of the gauge-invariant quark opera- Nf 1 L quark loops , (119) tors holds both in the ’t Hooft and Veneziano limits: ∼ N N 2 since each quark loop results in N . f OΓ1 OΓn = OΓ1 OΓn + (1/(Nf N)) . The quark loops are not suppressed at large N in the h · · · i h i···h i O (125) Veneziano limit The nonfactorized part, which is associated with connected diagrams, is 1/N in the ’t Hooft limit. This N N (120) ∼ f ∼ → ∞ leads, in particular, to the coupling constant of meson– meson interaction of order 1/N. The Veneziano limit is if the diagram is planar. analogous to pure Yang-Mills. It is the Veneziano limit that is related to the hadronic topological expansion in the dual-resonance models. In The factorization can be seen (at all orders of per- the Veneziano limit hadrons can have finite widths ac- turbation theory) from Eq. (115) for the contribution of cording to Eq. (118). a generic connected graph of genus h with B external boundaries which are precisely associated with the Wil- son loops (or the quark operators OΓ). The diagrams X. LARGE-N FACTORIZATION with gluon lines emitted and absorbed by the same operator are products of diagrams having only one boundary. The vacuum expectation values of several colorless or Their contribution is of order one. The diagrams with white operators, which are singlets with respect to the gluon lines emitted and absorbed by two different oper- gauge group, factorize in the large-N limit of QCD (or ators have two boundaries. This proves the factorization property (125) at all orders of perturbation theory. other matrix models) as was first noticed by A. A. Migdal and independently E. Witten in late 1970’s. The large-N factorization can also be verified beyond The simplest gauge-invariant operators in a pure perturbation theory at all orders of the strong-coupling SU(N) gauge theory are the closed Wilson loops expansion in the SU(N) lattice gauge theory. A nonperturbative proof of the factorization was given using µ 1 ıg H dz Aµ(z) quantum equations of motion (the loop equations) [11]. Φ(C)= tr P e C . (121) N They obey the factorization property

Φ(C ) Φ(C ) = Φ(C ) Φ(C ) + N −2 . XI. CONCLUSION h 1 · · · n i h 1 i···h n i O (122) The factorization implies a semiclassical nature of the We have considered in these lecture notes the basic large-N limit of QCD (a saddle point in the path integral features of the methods for nonperturbative studies of for certain variables). gauge theories, which were developed in the second half The factorization property also holds for gauge- of 1970’s – early 1980’s. Their contemporary applica- invariant operators constructed from quarks as in tions in high-energy physics are extremely broad: from Eq. (109). For the case of several ﬂavors Nf , we nor- the scattering of particles at very high energies to the at- malize these quark operators by tempts of constructing a uniﬁed theory of all interactions, 1 including gravity. Some of these issues are considered in OΓ = ψ¯Γψ . (123) other lectures at this School. Nf N

[1] Y. Makeenko, Methods of contemporary gauge theory [8] W.T. Tuttle, Can. J. Math. 14, 21 (1962). (Cambridge Univ. Press, 2002). [9] J. Koplik J., A. Neveu, and S. Nussinov, Nucl. Phys. B [2] Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 123, 109 (1977). [3] C.N. Yang and R.L. Mills, Phys. Rev. 96, 191 (1954). [10] G. Veneziano, Nucl. Phys. B 117, 519 (1976). [4] H. Weyl, Ann. der Phys. 59 v. 10, 101 (1919). [11] Y. M. Makeenko and A. A. Migdal, Phys. Lett. B 88, [5] F. London, Z. Phys. 42, 375 (1927). 135 (1979). [6] K. G. Wilson, Phys. Rev. D 10, 2445 (1974). [7] G. ’t Hooft, Nucl. Phys. B 72, 461 (1974).