Large N Gauge Theory and k-Strings by Shuhang Yang, B.Sc., Tsinghua University , 2007 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE in The Faculty of Graduate Studies (Physics) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) April 2011 c Shuhang Yang, 2011 Abstract We considered the k-antisymmetric representation of U(N) gauge group on two dimensional lattice space and derived the free energy by saddle point approximation in large N limit. k is a large integer comparable with N. Besides Gross-Witten phase transition[1], which happens as the coupling constant changes, we found a new phase transition in the strong coupling system that happens as k changes. The free energy of the weak coupling system is a smooth function of k under continuous limit. We have carefully selected the right saddle point solution among other possible ones. The numerical results match our saddle point calculations. ii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . vi Dedication . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vii 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 2 The U(N) Gauge Theory on The Two Dimensional Lattice Space . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 2.1 Gross-Witten Model . . . . . . . . . . . . . . . . . . . . . . . 3 2.2 k-Antisymmetric Representation . . . . . . . . . . . . . . . . 9 3 The Distribution of The Eigenvalues of The U(N) Gauge Group . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 4 The k-Antisymmetric Representation . . . . . . . . . . . . . 18 4.1 The Generating Function . . . . . . . . . . . . . . . . . . . . 18 4.2 Saddle Point Approximation . . . . . . . . . . . . . . . . . . 20 4.3 The Weak Coupling . . . . . . . . . . . . . . . . . . . . . . . 22 4.4 The Strong Coupling . . . . . . . . . . . . . . . . . . . . . . 25 4.5 Comparison at The Boundry . . . . . . . . . . . . . . . . . . 27 iii Table of Contents 5 Phase Transition . . . . . . . . . . . . . . . . . . . . . . . . . . 28 5.1 The Free Energy of The Weak Coupling System . . . . . . . 28 5.2 The Free Energy of The Strong Coupling System . . . . . . . 32 5.3 The Phase Transition . . . . . . . . . . . . . . . . . . . . . . 33 5.4 Gross-Witten Phase Transition for k-Representation . . . . . 35 6 Remarks and Conclusions . . . . . . . . . . . . . . . . . . . . 43 6.1 Discussion About The Distribution Function . . . . . . . . . 43 6.2 The k-String Tension . . . . . . . . . . . . . . . . . . . . . . 44 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 iv List of Figures 5.1 A=0.8 (p=0.6), Γ+ (θ+ (R)): Green; Γ+ (θ− (R)): Red; Γ− (θ+ (R)): Blue; Γ− (θ− (R)): Black. . . . . . . . . . . . . . . . . . . . . . 5.2 30 A=1 (p=0.3), Γ+ (θ+ (R)): Green; Γ+ (θ− (R)): Red; Γ− (θ+ (R)): Blue; Γ− (θ− (R)): Black. . . . . . . . . . . . . . . . . . . . . . 31 5.3 A=0.8 (p=0.6), Γ+ (R): Red; Γ− (R):Green . . . . . . . . . . 32 5.4 A=1 (p=0.5), Γ+ (R): Red; Γ− (R):Green . . . . . . . . . . . 33 . . . . . . . . . . . . . . . . . . 34 . . . . . . . . . . . . . . . . . . 35 . . . . . . . . . . . . . . . . . . 36 d dR Γ+ (θ+ (R)) d dR Γ+ (θ− (R)) d dR Γ− (θ+ (R)) d dR Γ− (θ− (R)) 5.5 A=0.8 (p=0.6), 5.6 A=0.8 (p=0.6), 5.7 A=0.8 (p=0.6), 5.8 A=0.8 (p=0.6), . . . . . . . . . . . . . . . . . . 37 5.9 A=0.8 (p=0.6), Γ+ (R) . . . . . . . . . . . . . . . . . . . . . . 38 5.10 A=0.8 (p=0.6), Γ− (R) . . . . . . . . . . . . . . . . . . . . . . 39 5.11 p=0.1, Γ1 (R): Red; Γ2 (R):Green . . . . . . . . . . . . . . . . 40 5.12 p=0.2, Γ1 (R): Red; Γ2 (R):Green . . . . . . . . . . . . . . . . 40 1 2.1 , . . . . . . . . . . . . . . . 41 5.14 p=0.5, Γ1 (R): Red; Γ2 (R):Green . . . . . . . . . . . . . . . . 41 5.13 p = Γ1 (R): Red; Γ2 (R):Green ∂ 2 Γw.c ∂λ2 λ=2 − ∂ 2 Γs.c . ∂λ2 λ=2 . . . . . . . . . . . . . . . . . . . . . 42 5.16 Phase Diagram . . . . . . . . . . . . . . . . . . . . . . . . . . 42 5.15 v Acknowledgements I would like to thank Professor Gordon W. Semenoff and Dr. Joanna Karczmarek. Thank my family for their support. vi To My Parents. vii Chapter 1 Introduction Physicists have been trying hard to understand confinement of quarks for about half a century. Although this problem has not been solved with satisfactory, many branches of study in physics have been created and flourished as a result. In 1970s, Kenneth Wilson initiated the method of quantizing the gauge field on space-time lattice in his efforts to demystify quark confinement [2][3]. On the lattice space the minimal interval of space-time is the lattice spacing and there is no need to renormalize the gauge theory. The expectation value of Wilson loop operator is the potential between fermions, and the string tension is the potential energy divided by the spatial length of the loop. Combining his former work on Wilson loop theory, Wilson made the theory of lattice QCD understandable and computable. Another influential method proposed by G. ’tHooft is the large N expansion[6]. Since perturbative expansion fails in the strong interacting systems, ’tHooft expanded the action of SU(N) gauge theory in power of 1/N, where N is large. This new point of view gave encouraging results to explain some phenomena beyond standard model[4]. Large N QCD became an important method in the study of strong interactions of quarks[5]. The U(N) lattice gauge theory is the combination of lattice quantization and large N gauge theory. The theory has been built up and some basic techniques have been developed[7]. In this paper we discussed k-strings, which is the k-representation of the U(N) gauge group on 2 dimensional lattice space. Using the saddle point approximation, we derived the free energy of k-strings. Some numerical results are presented to compare with the saddle point approximation. 1 Chapter 1. Introduction Gross-Witten model is a U(N) gauge theory on two dimensional lattice space[1]. In Gross and Witten’s paper, they showed that the Wilson loop operator of the fundamental representation can be reduced to the loop operator on one plequatte by the U(N) gauge symmetry. In the large N limit, they used the saddle point method to determine the distribution function of the eigenvalues of the U(N) gauge field, then they derived the string tension of the fundamental representation and found a phase transition as the coupling constant changes. We further studied the k-antisymmetric representation of U(N) gauge group, where k is also large and comparable with N. k is a new parameter in our work. We showed that the action and the Wilson loop of the k-antisymmetric representation can be also reduced to one plequatte on lattice space as the case of the fundamental representation. Using the generating function and the saddle point method We calculate the expectation value of the Wilson loop operator and the free energy Γ, which is a function of both the coupling constant λ and the integer k. In the weak coupling system (λ ≤ 2), We find several branches of free energy. None of them is smooth and none could serve as the free energy in the whole range of k N ∈ [0, 1]. The free energy is not simlpy the minimal solution, because we need to choose a proper contour of the integration of the generating function in complex plane. Further work has been done to decide which one of the two is the real saddle point solution. Under strong coupling, there are two branches of free energy from the saddle point equation. We also need to make the choice of which one is the proper saddle solution in different ranges of k N ∈ [0, 1]. A phase transition appears at k N = 12 . This phase transition is easily determined to be of first order, but it changes to be third order at the boundary between weak and strong coupling, where the coupling constant λ = 2. 2 Chapter 2 The U(N) Gauge Theory on The Two Dimensional Lattice Space In this chapter we fisrt review the U(N) gauge theory on the 2-dimensional lattice space and present some results of the fundamental representation done by D. Gross and E. Witten. They showed that by proper gauge choosing, the action and the Wilson loop operator can be reduced into simple forms, which only involve the integarl over one plaquette in the lattice space. Then we further the study on the k-antisymmetric representation. We are going to derive the Wilson loop operator of k-representation and reduce it into a simple form as Gross and Witten did for the fundamental representaion. 2.1 Gross-Witten Model U(N) Gauge thoery on two dimensional lattice is somehow trivial, because we could make the transverse dimension trivial by gauge transformation. This leads to the absence of physics gluons and the result that all physical quantities only depend on the integration over the gauge field on a single plaquette. The dynamical variables are N × N unitary matrices Un,i , which are 3 2.1. Gross-Witten Model functions of n and i. {i, j} is an orthonormal basis of the two dimensional lattice space, and n = n1 i + n2 j is a lattice vector. Un,i transprots a matter field from site n to n + i along i under fundamental representation of U(N). In this paper, we have nothing to do with the matter fields, the fermions. Since the action of the gauge fields can be reduced into one plaquette, the gluons are not physical particles due to the lack of time dependency. But by studying the non-physical gluons - the U(N) gauge fields we can derive the interaction between fermions and further derive the string tension. From Wilson loop theory, we have U † = U −1 = Un+i,−i . n,i (2.1) n,i The Wilson action is defined as 1 Tr g2 S(U ) = U + H.c. , (2.2) plaquette n summing all plaquettes on the lattice space, and U = Un,i Un+i,j Un+i+j,−i Un+j,−j , (2.3) plaquette is a loop orperator over one plaquette. The expectation values of physical observables are < O(U ) >= 1 Z [DU ] exp{S(U )}O(U ), (2.4) where Z is the vacuum-to-vacuum ampliude Z= [DU ] exp{S(U )}, (2.5) and DU is the multiple integration over U(N) group on each lattice with each direction: DU = dUn,k . (2.6) n, k 4 2.1. Gross-Witten Model It also has the property: DU = D(U V ) = D(V U ), (2.7) and dUn,i = dUn,j = 1. (2.8) The expectation value of the Wilson loop operator is WL (g 2 N ) = 1 ZN [DU ] exp{S(U )}Tr U , (2.9) L where L is a closed loop on the lattice space. The calculation of Z and WL can be greatly simplified considering the gauge transformation: Un,k → Vn Un,k V † n+k , k = i or j. (2.10) For arbitrary unitary matrices Vn , the action S(U ) is invariant under this transformation. By properly choosing Vn , we can make the paralell transportation alone transverse dirction trivial Un,i = I, (2.11) for all n. To give such gauge transformation explicitly, we need to find a class of unitary matrices {Vn } to satisfy Un,i = Vn Un,i V † n+i = I, (2.12) for any lattice vector n. To do this we first choose an arbitrary unitary matrix V0 , and from Eqn. (2.12) we have U0,i = V0 U0,i V † = I. i (2.13) Then we get the required unitary matrix on the site of i: Vi = V0 U0,i . (2.14) 5 2.1. Gross-Witten Model Further, at the site i we have Ui,i = Vi Ui,i V † = I, (2.15) V2i = Vi Ui,i . (2.16) 2i and therefore Following this procedure, we can have Vmi = V(m−1)i U(m−1)i,i (2.17) l=m−1 Uli,i , = V0 (2.18) l=0 for arbitrary integer m. Along j direction, we can’t put new restriction at this moment, because it may possibly bring in contradition to Eqn.(2.12). Therefore we leave Vnj , n ∈ Z to be arbitrary U(N) matrices. Based on {Vmi , m ∈ Z} and {Vnj , n ∈ Z}, we can construct the whole class of {Vn }. From Eqn. (2.12) we have Vnj Unj,i V † = I, (2.19) Vnj+i = Vnj Unj,i . (2.20) Vnj+2i = Vnj+i Unj+i,i . (2.21) nj+i so Similarly, Continuing this procedure, we have Vnj+mi =Vnj+(m−1)i Unj+(m−1)i m−1 =Vnj Unj+li,i . (2.22) l=0 Thererfore, We can determine Vn on each lattice site n leaving the arbitrariness of Vmj , which is the remaining gauge freedom. 6 2.1. Gross-Witten Model Under this gauge the action is S(U ) = 1 g2 Tr(Un,j U † n+i,j + H.c.). (2.23) n Originally, the action contains products of gauge fields Un,k on four neighboring lattice sites: n, n + i, n + i + j and n + j and of two different transporting directions: i and j. After the gauge transformation, it contains products of gauge fields only on two neighboring lattices: n and n + j and just one direction: j. However, since the [DU ] is integrating ove all lattice sites, we hope to completely disentangle the products of gauge fields on different lattice sites. This can be done by changing variables: Un+i,j = Wn Un,j . (2.24) Then the action is Z= [dWn ] exp n n 1 Tr(Wn + Wn† ) . g2 (2.25) 2 We notice that the gauge fields Wn in S(W ) are disentangled. So Z = z V /a , where the V is the area of the two dimensional lattice, and a is the spacing. So V /a2 is the number of plaquttes encircled in the loop, and z= dW exp 1 Tr(W + W † ) , g2 (2.26) which is the action on one plaquette and is independent of the lattice vector n. We further make gauge choice Unj,j = I. This can also be done explicitly. We kept the arbitrariness of {Vmj , m ∈ Z} above, but now we need make further restriction to the set of {Vmj , m ∈ Z} to realize the required gauge Unj,j = I. (2.27) In general, we have Un,j = Vn Un,j V † n+j . (2.28) 7 2.1. Gross-Witten Model Still starting from the origin, we have U0,j = V0 U0,j V † = I. (2.29) j So Vj is no longer arbitrary and Vj = V0 U0,j . (2.30) Uj,j = Vj Uj,j V † = I, (2.31) At the site j, we have 2j and so V2j = Vj Uj,j = V0 U0,j Uj,j . (2.32) Iterating this procedure, we have m ∈ Z m−1 Vmj = V0 Ulj,j , (2.33) l=0 and the required gauge is realized. Under this gauge, Wilson loop operator turns to be 1 < Tr WL (g , N ) = N R−1 2 UT i+kj,j k=0 1 >= < Tr N R−1 0 Wli+kj >, (2.34) k=0 l=T −1 where a2 T × R is the area inclosed by the loop L. As we have disentangle the integration of the action into to each lattice point, let’s first look into the integration on a specific lattice site, say dWn . When we look into the gauge field on site n, the gauge fields on other sites Wm do not have any interference. We write R−1 0 Wli+kj = AWn B, (2.35) k=0 l=T −1 where A and B stand for the remaining products of the gauge fields on the sites other than n. For arbitrary unitary matrix V , Tr(Wn +Wn† ) is invariant 8 2.2. k-Antisymmetric Representation under the transformation Wn → V Wn V † . Also dW = d(V W V † ) from Eqn.(2.7). We have 1 N 1 Tr(Wn + Wn† ) Tr(AWn B) g2 1 1 (2.36) = dWn exp 2 Tr(Wn + Wn† ) Tr(AV Wn V † B) N g 1 1 = dV dWn exp 2 Tr(Wn + Wn† ) Tr(AV Wn V † B), N g dWn exp where we insert the integration change the order of integral dV over the U(N) group. We can inter- dV and dV Vij Vkl† = dWn . Further, we have the fact 1 δil δjk . N (2.37) Then dV Ami Vij Wjk Vkl† Blm dV Tr(AV W V † B) = = (2.38) 1 1 Aml Blm Wjj = Tr(AB)TrW. N N By iterating this procedure, we separate out each Wn , and have WL (g 2 , N ) = (w(g 2 , N ))RT , (2.39) where w(g 2 , N ) = 1 z dWn 1 1 TrWn exp 2 Tr(Wn + Wn† ) . N g (2.40) Now, we reduce the Wilson loop operator of fundamental representaion into one plaquette as well as the action. 2.2 k-Antisymmetric Representation We are going to reduce the Wilson loop operator ot k-representation in a similar way. Let’s begin with 2-antisymmetrical representation first. The 9 2.2. k-Antisymmetric Representation fundamental gauge field is still Wn ∈ U (N ), while the representation of the U (N ) gauge group is: (2) Wijkl = 1 [j l] 1 Wi Wk = (Wij Wkl − Wil Wkj ). 2! 2! The 2-antisymmetric representation of U(N) group is N (N −1) 2 (2.41) dimensional. The Wilson loop operator of this representation is: (2) WL (g 2 , N ) = R−1 1 N 0 Tr W (2) li+kj . (2.42) k=0 l=T −1 To make the equations look simple, we label the Ws by integers, as Wα . (2) And now WL = 1 N Tr R×T α=1 (2) Wα . W (2) is a map from Rn ∧ Rn to Rn ∧ Rn , which is a subspace of Rn ⊗ Rn . To be explicit: W (2) (v, w) = (W v, W w), (2.43) where (v, w) ∈ Rn ∧ Rn . To write out the components: (vi , wj ) → (Wil vl , Wjk wk ). (2.44) In the Wislon loop operator, we have the product of a series of W (2) s. We want to separate the trace of the product of the seroes into a product of traces of each one W (2) as we did above for the fundamental representation Let’s first consider a two successive maps: W (2) and Y (2) Y (2) W (2) (v, w) = Y (2) (W v, W w) = (Y W v, Y W w) = (Y W )(2) (v, w), (2.45) for any v, w ∈ Rn . Therefore, we have Y (2) W (2) = (Y W )(2) . (2.46) By induction, we have R×T α=1 (2) R×T Wα(2) = Wα . (2.47) α=1 10 2.2. k-Antisymmetric Representation Since the gauge transformation W1 → V W1 V † , V ∈ U (N ) will not change the Wilson loop operator, we have: (2) WL (g 2 , N ) = 1 Tr(V W1 V † B)(2) N (2.48) where B stands for the product of the other Wα s. To calculate the trace, let insert the integration dV as before, Tr(V W1 V † B)(2) = dV (2) Tr V (2) W1 (V † )(2) B (2) , (2.49) where we have also used Eqn.(2.47). We can put the integarl of U(N) group dV into the expectation < ... >. Eqn. (2.37) is based on Peter-Weyl theorem, which also applies to U(N) group, and we have dV VijR (V R )†kl = 1 δRR δil δjk , dim R (2.50) (R) where R stands for certain representation of the group and Vij is a matix element of R-representation. Therefore, (2) dV Tr V (2) W1 (V † )(2) B (2) = 1 (2) TrW1 TrB (2) n n dim(R ∧ R ) (2.51) Then we have: (2) Tr W1 B (2) = 1 (2) TrW1 TrB (2) . N (N − 1)/2 (2.52) By repeating this procedure, we can separate the trace of the product of W (2) s into product of traces: RT RT Wα(2) Tr TrWα(2) = α=1 . α=1 1 N (N − 1)/2 RT (2.53) Therefore (2) (2) WL (g 2 , N ) = wL (g 2 , N ) RT (2.54) where (2) wL (g 2 , N ) = 1 z dW 1 1 TrW (2) exp 2 Tr(W + W † ) . (2.55) N (N − 1)/2 g 11 2.2. k-Antisymmetric Representation Next we want to derive that for k-antisymmetric representation the expectation value of the Wilson loop operator can be reduced into a single plaquette as well. The k-antisymmetric representaion has its components: j ] [j [j j ...j ] 2 k = Wi1 1 Wij22 ...Wikk , Wi1 i12 ...i k (2.56) where W ∈ U (N ). First, it’s easy to verify the formular: 1 [i1 i2 ...ik ] 1 [j1 j2 ...jk ] 1 1 j2 ...jk ] = (Y W )[j Y W r1 r2 ...rk . k! r1 r2 ...rk k! i1 i2 ...ik k! (2.57) The formular above means Y (k) W (k) = (Y W )(k) . (2.58) By induction, we have (k) Wα(k) = α Wα . (2.59) α So, RT (k) WL (g 2 , N ) = Wα(k) Tr α=1 (k) RT = Tr Wα α=1 = Tr(W1 B)(k) , (2.60) where B is the product of the other Wα s. As before, from the invariant gauge transformation W1 → V W1 V † we have Tr(W1 B)(k) = Tr(V W1 V † B)(k) (k) = Tr V (k) W1 (V † )(k) B (k) We use the similar trick: inserting the integration . (2.61) dV over the U(N) group, 12 2.2. k-Antisymmetric Representation (k) dV Tr V (k) W1 (V † )(k) B (k) = dim(Rn 1 (k) TrW1 TrB (k) , ∧ Rn ∧ ... ∧ Rn ) (2.62) where we also have dim(Rn ∧ Rn ∧ ... ∧ Rn ) = k times N! . k!(N − k)! (2.63) Continue this procedure, we have (k) (k) WL (g 2 , N ) = wL (g 2 , N ) RT , (2.64) where (k) wL (g 2 , N ) = 1 z dW 1 TrW (k) exp N !/k!/(N − k)! 1 Tr(W + W † ) . g2 (2.65) Therefore, for any k ∈ Z, the Wilson loop operator of k-antisymmetric representation is reduced in the same way as in the fundamental representation. 13 Chapter 3 The Distribution of The Eigenvalues of The U(N) Gauge Group In this chapter we present the distribution of eigenvalues of U(N) gauge field derived by Gross and Witten[1], and the results following the distribution function. We will also present the Gross-Witten phase transition in fundamental representaion. Following from the arguement in the last chapter, to determine the free energy of the k-antisymmetric representation ΓAk (g 2 , N ) = − 1 (k) ln WL (g 2 , N ), N (3.1) we only need to consider the action and the Wilson loop operator of the gauge field on one plaquette. The integration dW is over the U(N) group. Since the integrands only depend on the eigenvalues of U(N), we can reduce W into the integral over the eigenvalues of the U(N) group, and diagnolize W by gauge tranformation: V W V † = diag{eiφ1 , eiφ2 , ..., eiφN }, V ∈ U (N ). As U(N) is a compact group, its Haar measure can be reduced to a measure on a subspace of Rn multiplied by a modular function. We have[8], 2π N dφi ∆2 (φi ), dW = 0 (3.2) i=1 14 Chapter 3. The Distribution of The Eigenvalues of The U(N) Gauge Group where (eiφj − eiφi ) ∆(φ) = (3.3) i<j is the Vandermonde deterninant. Then we have the action 2π N z(g 2 , N ) = c 0 i=1 φi − φj 2 exp 2 2 g sin2 dφi i<j N cos(φi ) , (3.4) i=1 where the constant c is to normalize such that z(∞, N ) = 1. We are going to use the steepest-descent method to evaluate the action, so we write the integrand in the form of ef (φ) : 2 2π N z(g 2 , N ) = c dφi exp g2 0 i=1 N cos(φi ) + i=1 i=j φi − φj ln sin . (3.5) 2 We use the steepest-descent method, and the approximated energy is ln z(g 2 , N ) 2 N →∞ N N 1 2 = lim cos(φi ) + N →∞ N 2 g 2 E(λ) = lim i=1 i=j φi − φj ln sin + const. 2 (3.6) The eigrnvalues φi can be approximately determined by the stationarity condition: 2 ∂ φi g2 N cos(φi ) + i=1 i=j φi − φj ln sin =0, 2 (3.7) which is simplified as 1 2 sinφi = λ N cot i=j φi − φj , 2 (3.8) where λ = g 2 N . As N goes to infinity, we treat φi as a value assumed by a continuous function φ at i N: φi = φ(i/N ). (3.9) 15 Chapter 3. The Distribution of The Eigenvalues of The U(N) Gauge Group Then the stationary condition becomes: 2 sinφ(x) = P λ 1 dycot 0 φ(x) − φ(y) , 2 (3.10) where P denotes the principal value of the integration. To find the function φ(x), we need to determine the density of eignevalues ρ(φ) = satisfies φc −φc ρ(φ)dφ = 1 0 dx dx dφ , which = 1 . By solving the equation above, Gross and Witten gave the result in their paper[1]: φ λ 2 2 φ 1/2 −1 λ 1/2 πλ cos 2 ( 2 − sin 2 ) , λ ≤ 2 φc = 2 sin ( 2 ) ρ(φ) = 1 2 λ ≥ 2 φc = π. 2π 1 + λ cos φ , (3.11) In the first case, the density function is gapped, which means it doesn’t cover the whole interval of [0, 2π]. In the second case, the density function is ungapped, and all the values between 0 and 2π can be assumed by φ(x). As the density function has been sovled, the energy can be calculated by replacing the summation N i=1 −E(λ) = by the integral 1 λ2 , 2 λ + 12 ln λ2 − 34 , dφρ(φ). We have λ≥2 (3.12) λ ≤ 2. The expectation value of the Wilson loop operator is λ≥2 λ1 , λ2 ∂ ln z 2 wL (λ) = − = lim wL (g , N ) = N →∞ 2N 2 ∂λ 1 − λ4 , λ ≤ 2, (3.13) and the string tesion σ(λ) = − 1 ln λ, a2 1 ln wL = lim σ(g 2 , N ) = N →∞ a2 1 a2 4 ln 4−λ , λ≥2 (3.14) λ ≤ 2. 16 Chapter 3. The Distribution of The Eigenvalues of The U(N) Gauge Group Gross-Witten phase transition happens at λ = 2. By checking the derivatives of the energy Eqn.(3.12), we can see that this phase transition is of third order. 17 Chapter 4 The k-Antisymmetric Representation As we have reduced the Wilson loop operator of k-antisymmetric representation into one plaquette and have had the distribution function ρ(φ) of the U(N) eigenvalues, we are ready to derive the free energy of k-antisymmetric representation. In this chapter we first present the generating function of Wilson loop operator of k-antisymmetric representation. Then we are going to use the saddle point method to approximate the free energy. From the saddle point equations, we pluge in the distribution function for the weak coupling system. After laborious derivations, we get the paramter k N as a function of t and the approximated free energy as a function of t. Further work has been done to combine them and we have the approximated free energy as a function of k N. 4.1 The Generating Function The calculation for the strong coupling is more straightforward. For the gauge field W ∈ U (N ) we can diagnolize it by gauge transformation. Therefore, without loss gnerality, we can let the gauge field be W = diag{eiφ1 , eiφ2 ...eiφN }. (4.1) 18 4.1. The Generating Function The k-antisymmetric representation of the U(N) group has its components: j ] [j Wi1 1 Wij22 ...Wikk , (4.2) and we have the trace of W (k) as 1 i ] [i TrAk W = Wi1 1 Wii22 ...Wikk k! eiφ1 eiφ2 ...eiφN , = (4.3) i1 <i2 <...<iN Now, we are going to introduce the generating function. First we have the expansion ∞ (−1)n+1 (tW )n , n ln(1 + tW ) = n=1 (4.4) and we have ∞ N Tr ln(1 + tW ) = i=1 n=1 (−1)n+1 iφi n (te ) n N ln(1 + teiφi ). = (4.5) i=1 The generating function is dt eTr ln(1+tW ) = 2πitk+1 N i=1 (1 + teiφi ) 2πitk+1 dt eiφ1 eiφ2 ...eiφk . = (4.6) i1 <i2 <...<ik Then we have the trace evaluated by the generating function TrAk W = dt eTr ln(1+tW ) 2πitk+1 PN = dt e i=1 ln(1+teiφi ) 2πitk+1 , (4.7) 19 4.2. Saddle Point Approximation and the expectation value of the Wilson loop operator is N dt TrAk W = dφi 2πitk+1 i=1 2N N exp cos φi + λ i=1 ln cot2 i<j φi − φj + 2 N ln(1 + teiφi ) . i=1 (4.8) 4.2 Saddle Point Approximation We use the saddle point approximation to evaluate the multiple integration. To do this, we need to put all the integral variables to the exponent, such as tk+1 = e(k+1) ln t . Then the saddle point conditions are: ∂ φi 2N λ N ln cot2 cos φi + i=1 i<j φi − φj 2 N ln(1 + teiφi ) − (k + 1) ln t + = 0, (4.9) i=1 and ∂t 2N λ N ln cot2 cos φi + i=1 i<j φi − φj 2 N ln(1 + teiφi ) − (k + 1) ln t + = 0, (4.10) i=1 From the saddle point conditions we have: 2N − sin(φi ) + λ N i=1 N cot j=i φi − φj iteiφi = 0, + 2 1 + teiφi eiφi k+1 . = iφ t 1 + te i (4.11) (4.12) 20 4.2. Saddle Point Approximation The two equations serve to determine the distribution function of the U(N) eigenvalues for the k-representation and the saddle point t in the complex plane. Since the last term in the first equation is supressed by 1 N, it will not influence the distribution function much. So we can take the distribution function ρ(φ) the same as the one given for the fundamental representation. As N goes to infinity, we can take the continuous limit and replace by N π −π N i=1 dφρ(φ), then the second saddle point equation (4.12) becomes: π RAk (t) = where RAk = k+1 N . −π teiφ , 1 + teiφ (4.13) And we have PN TrAk W = ≈ dφρ(φ) dt e i=1 ln(1+teiφi ) 2πitk+1 π dt N R−π dφρ(φ) ln(1+teiφ )−(k+1) ln t e 2πi π dφρ(φ) ln(1 + teiφ ) − (k + 1) ln tˆ , ≈ exp N (4.14) −π where ρ(φ) and tˆ satisfy the saddle point equations (4.11) and (4.12). Since e−N ΓAk =< TrAk W >, we have the un-normalized free energy: Γ Ak = − dφρ ln(1 + teiφ ) + k+1 ln(t). N (4.15) We are going to combine Eqn.(4.13) and Eqn.(4.15) to determine the free energy ΓAk as a function of k N. To do this, we need to finish the integration in saddle point equations for specific density function ρ(φ). We’ll work on this in the following. 21 4.3. The Weak Coupling 4.3 The Weak Coupling Since the density function is symmetric about the origin: ρ(φ) = ρ(−φ) (4.16) and t = e2θ , from Eqn.(4.13) we have RAk (t) − 1 = 2 π dφρ(φ) −π e2θ eiφ 1 − 2θ iφ 2 1+e e 1 2 + c.c. . (4.17) The integrand is e2θ eiφ 1 − 2θ iφ 2 1+e e 1 2 + c.c. = e4θ − 1 1 + e4θ + 2e2θ cos(φ) 1 2 . (4.18) Let x = sin( φ2 ), and we have: 1 2 e2θ eiφ 1 − 2θ iφ 2 1+e e + c.c. = = 1 2 e4θ − 1 1 + e4θ + 2e2θ − 4e2θ x2 1 sinh 2θ . 2 4 cosh2 θ 1 − x 2 cosh θ (4.19) From the expasion: ∞ 1 1− x2 cosh2 θ = m=0 x2 cosh2 θ m , (4.20) we have RAk 1 − = 2 π sinh 2θ dφρ(φ) 4 cosh2 θ −π ∞ m=0 x2 cosh2 θ m . (4.21) In the weak coupling system, we have the distribution function ρ(φ) = cos φ2 π(2 − 2p) where φc = 2 arcsin λ 2 2 − 2p − sin2 φ , −φc ≤ φ ≤ φc 2 (4.22) and p = wL , which is the Wislon loop operator under the fundamental representation. Weaking coupling means λ ≤ 2, and 22 4.3. The Weak Coupling from Eqn.(3.13) we know wL = 1− λ4 . So we have 2−2p = λ2 . Let A = 2−2p and change variable: sin φ2 = x, we have √ 2 A − x2 dφρ(φ) = √ dx . Aπ −π − A √ π A (4.23) We need to do the integral: √ √ 1 1 A 2 A − x2 2m Am Γ(m − 12 ) m (m − 2 ) × · · · × 2 √ dx x = = A . √ Aπ (m + 1)! πΓ(m + 2) − A (4.24) From Eqn.(4.21) we have sinh 2θ 1 RAk (t) − = 2 4 cosh2 θ ∞ m=0 A cosh2 θ m (m − 12 ) × · · · × 12 . (m + 1)! (4.25) A Further let y = − cosh we have: 2 θ RAk (t) − sinh 2θ 2 1 = ( 1 + y − 1), 2 4 cosh2 θ y (4.26) which is RAk (t) − sinh θ 1 = cosh θ − 2 A cosh2 θ − A . (4.27) Next we are going to derive the free energy for weak coupling system. Similarly, because the distribution function ρ(φ) is even, we have ln(1 + teiφ ) = 1 ln(1 + teiφ ) + c.c., 2 (4.28) where t = e2θ . Then 1 1 ln(1 + teiφ ) + c.c. = ln(1 + e4θ + 2e2θ cos φ). 2 2 (4.29) Let x = sin φ2 , we have 1 1 x2 ln(1 + teiφ ) + c.c. = ln(1 + e2θ ) + ln 1 − 2 2 cosh2 θ . (4.30) Only the second term depends on x, and we have the expansion: ln 1 − x2 cosh2 θ ∞ =− n=1 1 x2n . n cosh2n θ (4.31) 23 4.3. The Weak Coupling Then we need to do a similar integration as above, and from Eqn.(4.24) we have: √ 2 A − x2 x2 dx ln 1 − √ Aπ cosh2 θ − A √ A ∞ − n=1 Further let y = A , cosh2 θ = 1 An (n − 12 ) × · · · × 12 . n cosh2n θ (n + 1)! (4.32) y n (n − 12 ) × · · · × 21 , n (n + 1)! (4.33) and ∞ I(y) = − n=1 then we have 1 k ΓAk = − ln(1 + e2θ ) − I(y) + ln(e2θ ). 2 N (4.34) Now we need to simplify I(y). The first derivation of I(y) is: ∞ y n−1 I (y) = − n=1 (n − 21 ) × · · · × 21 . (n + 1)! (4.35) By adding and multiplying, we get a power series of y which sums up to a simple function, ∞ yn (−y)I (y) + 1 = n=0 (n − 21 ) × · · · × (n + 1)! 1 2 2 = − ( 1 − y − 1). y (4.36) So √ 2 1−y 2 1 I (y) = − 2+ . 2 y y y Integrating the equation above over y, we have √ 2 1−y 2 I(y) = − + 2 ln( 1 − y + 1) + + C, y y (4.37) (4.38) 24 4.4. The Strong Coupling where C is a constant. From (4.33) we know that lim I(y) = 0, (4.39) y→0 and lim y→0 √ 2 1−y 2 − + 2 ln( 1 − y + 1) + + C y y = 1 + 2 ln 2, (4.40) then we have C = −1 − 2 ln 2. From(4.34) we have cosh θ − ΓAk = (2θ sinh θ − cosh θ) 1 cosh2 θ − A + . 2 − ln cosh θ + 4.4 cosh2 θ − A A (4.41) The Strong Coupling In the strong coupling, the density function of the U(N) eigenvalues is: ρ(φ) = 1 (1 + 2p cos(φ)), 2π where now λ ≥ 2 and p = wL = 1 λ. The expectation value of Wilson loop operator is −N ΓAk e dt eN 2πi = (4.42) Rπ −π dφρ(φ) ln(1+teiφ ) tk+1 . (4.43) Consider the integral: π dφρ(φ) ln(1 + teiφ ), (4.44) −π when ||t|| < 1, we can expand the logarithm term ∞ iφ ln(1 + te ) = − n=1 (−1)n (teiφ )n . n (4.45) Only the first term in the summation contributes to the integral and gives π dφ −π 1 (1 + 2p cos(φ))(teiφ ) = pt. 2π (4.46) 25 4.4. The Strong Coupling Therefore, e−N ΓAk = ∞ (N pt)n n=0 n! , Expand eN pt = dt eN pt , 2πi tk+1 ||t|| < 1. (4.47) (N pt)k k! contributes to the and only the term contour integration. We have (N p)k dt (N pt)k = . 2πi k!tk+1 k! Then the un-normalized free energy is 1 ΓAk = − (k ln(N p) − ln(k!)) . N When k is also large, we can further use Stirling’s formula: e−N ΓAk = (4.48) (4.49) ln(k!) = k ln k − k + o(k). (4.50) Then, up to the leading order, we have k k ΓAk = ln( ). N N pe (4.51) When ||t|| > 1, we have 1 ln(1 + teiφ ) = ln(teiφ ) + ln(1 + e−iφ ), t 1 −iφ and expand ln(1 + t e ) in the same way as above. We have (4.52) π p dφρ(φ) ln(1 + teiφ ) = ln(t) + , t −π ||t|| > 1 (4.53) and dt tN N p (N p)N −k t = e , ||t|| > 1. 2πi tk+1 (N − k)! This branch of the un-normalized free energy is 1 ΓAk = − [(N − k) ln(N p) − ln((N − k)!)] . N Use Stirling’s formula again, we have e−N ΓAk = ΓAk = (1 − k ) ln N 1 − Nk pe . When p > 12 , it is weak coupling (λ < 2) . When p < (4.54) (4.55) (4.56) 1 2 the system is strong coupling (λ > 2). We will verify that the two free energies for strong and weak coupling systems match each other at p = 12 , i.e.A = 1(λ = 2). 26 4.5. Comparison at The Boundry 4.5 Comparison at The Boundry Starting from weak coupling system when A = 1, from Eqn.(4.27)we have k 1 = sinh θ cosh θ − cosh2 θ − 1 + . (4.57) N 2 Consider θ as a complex variable, so Nk is a two-valued function of θ: √ sinh2 θ = sinh θ sinh θ(cosh θ − sinh θ) + 21 , k (4.58) = N √ sinh2 θ = − sinh θ sinh θ(cosh θ + sinh θ) + 12 . With t = e2θ , we have 1− k = N t 1 2t , 2. √ sinh2 θ = sinh θ (4.59) √ sinh2 θ = − sinh θ Then let’s check the free energy. When A = 1, from Eqn.(4.41) we have: Γ = (2θ sinh θ − cosh θ) cosh θ − sinh2 θ 1 sinh2 θ + . 2 √ 2 We also have two choices of the sign of sinh θ: √ sinh2 θ = sinh θ −θe−2θ − 21 e−2θ , Γ= √ sinh2 θ = − sinh θ. Γ = θe2θ − 21 e2θ . − ln cosh θ + Pluge in (4.59), we have k 2(1− N ) k 1 − ln N e Γ= Γ = k ln 2k . N eN , (4.60) (4.61) √ sinh2 θ = sinh θ (4.62) √ sinh2 θ = − sinh θ They are the same as the two branches of free energy that we get in the strong coupling system at p = 12 . 27 Chapter 5 Phase Transition 5.1 The Free Energy of The Weak Coupling System In the weak coupling system, we already have the free energy Eqn.(4.41) associated with Eqn.(4.27): R(θ) = sinh θ k = cosh θ − N A 1 cosh2 θ − A + . 2 k N, To solve the free energy as a function of R = (5.1) we must combine the two equations. From the second equation, we could find the inverse function: θ(R) and then pluge it into the first one Γ(θ). Let u = eθ we have: (u + u1 )2 u+ −A= 4 2 1 u − k 2A( 2N ) 1 u− u (5.2) Squaring both sides, it turns 4A − 8A k 1 − N 2 u4 + 16A2 2 k 1 − N 2 − 8A u2 + 8A k = 0. (5.3) N Solve the quadratic equation of u2 , we have u2± (R) = 1 2(1 − R) ± 1 − 2A R − 1 1 − 2A R − 2 2 2 1 2 2 − 4 (1 − R) R. (5.4) 28 5.1. The Free Energy of The Weak Coupling System Then we have θ± (R) = 1 ln u± (R). 2 (5.5) We regard θ as a complex variable, then the free energy Γ(θ) is a two-valued complex function due to cosh2 θ − A. So we have two choices of the free energy: Γ± (θ) = 1 (2θ sinh θ − cosh θ) cosh θ ∓ A − ln cosh θ ± cosh2 θ − A 1 cosh2 θ − A + . 2 Now we have four branches of free energy as a function of R = (5.6) k N: Γ+ (θ+ (R)), Γ+ (θ− (R)), Γ− (θ+ (R)) and Γ− (θ− (R)), see Fig.(5.1) and Fig.(5.2). k N = 1 2, we have u2± ( 12 ) At = 1, and θ± = 0. So we have 1 1 Γ+ (θ+ (R = )) = Γ+ (θ− (R = )), 2 2 (5.7) 1 1 Γ− (θ+ (R = )) = Γ− (θ− (R = )). 2 2 (5.8) and We are going to take the bottom curves as free energy. We have two bottome curves as the moment: one is patched by Γ+ (θ± (R)), and the other one is patched by Γ− (θ± (R)), Γ+ (R) = Γ+ (θ− (R)), R≤ Γ+ (θ+ (R)), R ≥ 12 , Γ (θ (R)), − + R≤ Γ− (θ− (R)), R ≥ 12 , 1 2 (5.9) and Γ− (R) = 1 2 (5.10) see Fig.(5.3) and Fig.(5.4). By checking derivatives, we found that none of these four functions Γ± (θ± (R)) are smooth. see Fig.(5.5), Fig.(5.6), 29 5.1. The Free Energy of The Weak Coupling System Figure 5.1: A=0.8 (p=0.6), Γ+ (θ+ (R)): Green; Γ+ (θ− (R)): Red; Γ− (θ+ (R)): Blue; Γ− (θ− (R)): Black. Fig.(5.7) and Fig.(5.8). However the patched are Γ+ (R) and Γ− (R) are smooth curves, see Fig.(5.9) and Fig.(5.10). Now we need to decide which one of the two is the solution from saddle point approximation. However, the issue cannot be simplily solved by choosing the smaller one. We need to look into the problem and know how each solution actually comes out. First, we note that the free energy ΓAk are multi-valued complex function on t-plane. This is due to the logarithm and the square root terms. The logarithm term contributes a constant to the free energy. Up to a constant, the free energy is 2-valued function, and there are two branch points determined from cosh2 θ = A t± = 2A − 1 ± 2 A(1 − A)i. (5.11) When A < 1, the two branch points are complex numbers and conjugate to each other. To make the free energy single valued, we need draw a branch cut connecting t± . Dr. Karczmarek pointed out that the value of free energy depends on the way we draw the branch cut. The reason is as followed. Let’s first look into some basic facts about the contour choosing and the saddle 30 5.1. The Free Energy of The Weak Coupling System Figure 5.2: A=1 (p=0.3), Γ+ (θ+ (R)): Green; Γ+ (θ− (R)): Red; Γ− (θ+ (R)): Blue; Γ− (θ− (R)): Black. point. First, the integrating contour must encircle the origin based on the generating function method. Secondly, the saddle point should be on the contour of course, and moreover be the summit on the chosen contour. Also, numerical simulation indictes that the saddle point locates on R+ axis and the origin is a highest point of the integrand. Based on the 3 facts, we can continue to discuss the possible ways to draw the branch cut and how to choose contours on the t-plane. To connect the branch points t± , which are conjugate to each other, the branch cut necessarily crosses the R axis. (The branch cuts could also connect the branch points to infinity, but it crosses the real axis at infinity as well.) Let x0 ∈ R+ be the saddle point and t0 be the crossing point of real axis and the branch cut. There are two cases to be considered. One is that t0 ∈ [0, x0 ]. In this case, t0 separates the origin and the saddle point on R asix. The contour must encirle the branch cut, otherwise it is not qualified to be used in saddle point method. The other case is t0 ∈ / [0, x0 ] and t0 is far away from the origin such that we could choose a valid contour without encircling the branch cut. Dr Karczmarek shows that the first kind of branch 31 5.2. The Free Energy of The Strong Coupling System Figure 5.3: A=0.8 (p=0.6), Γ+ (R): Red; Γ− (R):Green cut gives the free energy Γ− (R), and the second gives Γ+ (R). Despite the fact that Γ− (R) is lower than Γ+ (R) as showed in Fig.(5.3) and Fig.(5.4), we take Γ+ (R) as the valid saddle point solution, because the corresponding contour can shrink to the origin without intersecting the branch cut. When 0 < A < 1, Γ+ (R) is a smooth curve, see Fig.(5.9). 5.2 The Free Energy of The Strong Coupling System In the strong coupling system, we are still facing the problem to choose the right free energy among the two branches, see Fig.(5.11) and Fig.(5.12) Γ1 (R) = R ln( R ), pe and Γ2 (R) = (1 − R) ln( where R = 1−R ), pe (5.12) (5.13) k N. 32 5.3. The Phase Transition Figure 5.4: A=1 (p=0.5), Γ+ (R): Red; Γ− (R):Green The numerical result indicates the free energy for strong coupling system to be Γs.c. (R) = R R ln( pe ) 0≤R≤ (1 − R) ln( 1−R pe ). 1 2 (5.14) ≤R≤1 Γs.c. (R) is the bottom curve only when 0 < p < ever, when 1 e < p < 1 2, 1 2 1 e, see Fig.(5.11). How- Γs.c. (R) is not the minimum one in the whole range of R ∈ [0, 1]. The bottom curve is patched by 4 parts, as showed in Fig.(5.13). But the numerical simulation favors Γs.c. (R) in the whole range. The reason is unknown yet. This indicates that the continuous limit N i=1 φi 5.3 →N dφρ(φ) brings in some unknown outgrowth. The Phase Transition Γs.c. (R) has a phase transition at R = 12 , which is obvious from Fig.(5.11), Fig.(5.12) and Fig.(5.13) but not so obvious from Fig.(5.14). The order of the Phase transition can be determined by checking the derivatives of Γ1 33 5.3. The Phase Transition d dR Γ+ (θ+ (R)) Figure 5.5: A=0.8 (p=0.6), and Γ2 at R = 21 , which are the left and right derivatives of Γs.c. at R = 21 . d dR d dR Γ1 (R) = ln R=1/2 Γ2 (R) = − ln R=1/2 Their derivatives at R = 1 2 R p = ln R= 12 1−R p 1 , 2p = − ln R= 21 (5.15) 1 . 2p (5.16) are distinct in general, but both equal zero when p = 21 . So for the case of p = 12 , we need to consider higher order derivatives. d2 dR2 d2 dR2 Γ1 (R) = R= 12 Γ2 (R) = R= 12 1 R = 2, (5.17) R= 12 1 1−R = 2. (5.18) R= 12 The second derivatives match each other at R = 12 . Then we need to further consider the third derivatives: d3 dR3 Γ1 (R) = − R= 12 1 R2 = −4 (5.19) R= 12 34 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.6: A=0.8 (p=0.6), d3 dR3 Γ2 (R) = R= 12 d dR Γ+ (θ− (R)) 1 (1 − R)2 =4 (5.20) R= 21 We see that the third left and right derivatives of Γs.c. at R = 1 2 are not equal, and conclude that it’s a third order phase transition at p = 21 . Therefore, when p < 12 , the phase transition is 1st order, when p = 12 , it is 3rd order. When p > 12 , the system is weak coupling, and the free energy is a smooth function of 5.4 k N. Gross-Witten Phase Transition for k-Representation We have discussed the phase transition that happens when k N changes. The coupling constant is also a parameter of the free energy. In the fundamental representation, Gross-Witten phase transition happens when the coupling constant changes from weak to strong. In the k-antisymmetric representation we have already check that the free energy is continuous at the boundary between weak and strong coupling in Section 4.5. We are going to check 35 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.7: A=0.8 (p=0.6), d dR Γ− (θ+ (R)) the derivatives of the free energy from both sides to see if there is a phase transition. The free energy is also a function of the coupling constant λ. From Eqn.(5.4) and Eqn.(5.5) we know that u± (R, A) and θ± (R, A) are functions of R and A. Since A = 2 − 2p and 1 , λ≥2 λ p = wL (λ) = 1 − λ4 . λ ≤ 2 strong coupling (5.21) weak coupling Therefore, A or p is related to the coupling constant λ. The free energy is function of both R and λ, Γ(R, λ). From Eqn.(5.9) the free energy in the weak coupling system λ ≤ 2) is Γw.c. (R, λ) = Γ (θ (R, λ), λ), + − R≤ Γ+ (θ+ (R, λ), λ), R ≥ 12 , 1 2 (5.22) and from Eqn.(5.14) the free energy for the strong coupling system (λ ≥ 2) 36 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.8: A=0.8 (p=0.6), is Γs.c. (R, λ) = R ln d dR Γ− (θ− (R)) λR e 0≤R≤ 1 2 (5.23) (1 − R) ln λ(1−R) e . 1 2 ≤R≤1 Then we’ll compare the derivatives of Γw.c. (R, λ) and Γs.c. (R, λ). Let’s first consider the case at k N = 12 . From Eqn.(5.4) and Eqn.(5.5) we have θ± (R = 0.5) = 0, (5.24) and from Eqn.(5.6) we know √ √ 1 1 Γw.c. (θ = 0, A) = − (1 − 1 − A) − ln(1 + 1 − A) + . A 2 Since A = λ 2 (5.25) in the weak coupling system, we have Γw.c. (R = 0.5, λ) = − 2 λ 1− 1− λ 2 − ln 1 + 1− λ 2 1 + . (5.26) 2 Then we have ∂Γw.c. 2 (λ) = 2 ∂λ λ 1− 1− λ 2 − 1 , 2λ (5.27) 37 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.9: A=0.8 (p=0.6), Γ+ (R) and ∂Γs.c. 1 (λ) = . ∂λ 2λ (5.28) Then we find ∂Γw.c. ∂λ ∂Γs.c. ∂λ = λ=2 λ=2 1 = . 4 (5.29) So we need to consider the second derivatives: ∂ 2 Γw.c. 4 (λ) = − 3 2 ∂λ λ and We see that 1− 1− λ 2 1 + 2λ 1 − + λ 2 1 , 2λ2 (5.30) ∂ 2 Γs.c. 1 (λ) = − 2 . 2 ∂λ 2λ (5.31) ∂ 2 Γw.c. (λ) → ∞. λ→2 ∂λ2 (5.32) lim So we conclude that the phase transition is of second order at λ = 2 and k N = 21 . In the case k N ∈ (0, 0.5) ∪ (0.5, 1), ∂ 2 Γw.c. ∂λ2 λ=2 does not diverge, but we still have ∂Γw.c. ∂λ = λ=2 ∂Γs.c. ∂λ , (5.33) λ=2 38 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.10: A=0.8 (p=0.6), Γ− (R) and ∂ 2 Γw.c. ∂λ2 = λ=2 ∂ 2 Γs.c. ∂λ2 , (5.34) λ=2 see Fig.(5.15). Therefore, the phase transition is of second order at λ = 2 , k N ∈ (0, 1). Now, we have the whole phase diagram Fig.5.16. 39 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.11: p=0.1, Γ1 (R): Red; Γ2 (R):Green Figure 5.12: p=0.2, Γ1 (R): Red; Γ2 (R):Green 40 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.13: p = 1 2.1 , Γ1 (R): Red; Γ2 (R):Green Figure 5.14: p=0.5, Γ1 (R): Red; Γ2 (R):Green 41 5.4. Gross-Witten Phase Transition for k-Representation Figure 5.15: ∂ 2 Γw.c ∂λ2 λ=2 − ∂ 2 Γs.c . ∂λ2 λ=2 k N ✻ ✲ ✲ ✲ 2nd order 1st order 0.5 ✻ ✻ ✻ ✻3rd order ❄ ✲ ✲ 2nd order ✲ ✲ 2 λ Figure 5.16: Phase Diagram 42 Chapter 6 Remarks and Conclusions 6.1 Discussion About The Distribution Function We stated the the last term of Eqn.(4.11) teiφ 1+teiφ is supressed by 1 N, but this is true only if t is not closed to ei(−φ+π) . However, this is not necessarily true from the saddle point equations Eqn.(4.11) and Eqn.(4.12). Let’s first assume that the distribution function is the same as in fundamental representation, then find the saddle point tˆ. If there is no contradiction, i.e. tˆ is not close to ei(−φ+π) for any φ ∈ [−φc , φc ], we can say that our assumption is reasonable. Let’s start with the weak coupling system. The distribution function of fundamental representation is gapped, i.e. φc < π. Based on the assumed distribution function, the numerical simulation tells us that the saddle point tˆ is positive and real. So the distance between tˆ and {ei(−φ+π) |φ ∈ [−φc , φc ]} is finite. Therefore, in the large N limit, tˆeiφ 1+tˆeiφ is finite and negligible. Next, let’s consider the case in the strong coupling system. The distribution function is ungapped, which means φc = π. All the values on the unit circle in complex plane is assumed by the distribution function. From Eqn.(4.47) and Eqn.(4.54), we get the saddle point k+1 , ||tˆ|| < 1 Np tˆ = k+1p , ||tˆ|| > 1 −1 (6.1) N 43 6.2. The k-String Tension where 0 ≤ p ≤ 1 2. We can see that, when k + 1 = N p or the saddle point tˆ = 1 and tˆeiφ 1+tˆeiφ k+1 N −1 = p diverges at φ = π. This suggests that near φ = π the distribution function ρ(φ) is different from the one of the fundamental representation. The revision of the distribution function is left as an open question. 6.2 The k-String Tension As we have had the expectation value of the Wilson loop operator for kantisymmetric representation of U(N), we can further calculate the k-string tension: σk = − 1 (k) ln wL (g 2 , N ). a2 (6.2) With the normalization factor, the expecation value of Wilson loop operator on one plaquette under k-representation is (k) wL = 1 < TrW (k) > . N !/k!/(N − k)! (6.3) And the k-string tension is σk = − 1 a2 − ln N! k!(N − k)! + ln < TrW (k) > . (6.4) Let’s start with the strong coupling system first, because we have the exact solution under strong coulping and things are more straightforward. From (equation 4.* ) we have under strong coupling (k) wL = k!(N − k)! (N p)k (N − k)!(N p)k = . N! k! N! (6.5) Noting that when k = 1, w(1) = p, which matches (equation 3. *). When k is relatively small, say a finite integer, (N − k)!N k = 1, N →∞ N! lim (6.6) 44 6.2. The k-String Tension so (k) wL ≈ pk for finite k. (6.7) Then we have the ratio of k-string tension σk to the fundamental string tension σ: (k) ln wL σk ≈k = (1) σ ln wL for finite k. (6.8) Both the sine-formular σk sin(kπ/N ) = σ sin(π/N ) (6.9) k(N − k) σk = σ N −1 (6.10) and the Casimir scaling law are satisfied when k is finite and N goes to infinity. We are more concerned about the case when k is comparable with N. From (k) ln wL = ln((N − k)!) − ln(N !) + k ln(N p), (6.11) ans using Stirling’s formular, we have (k) ln wL = (N − k) ln N −k N + k(1 + ln p). (6.12) Now the ration of k-string tension to the fundamental one is (N − k) ln NN−k + k σk =k+ . σ ln p (6.13) The ratio also depends on the coupling constant p. This is not similar as sine-formular nor the Casimir scaling law. Next consider the weak coupling system: σk = 1 a2 −N Γ(k/N, p) + ln k!(N − k)! N! , 0.5 < p ≤ 1. (6.14) Using Stirling’s formular again, we have the k-string tension ration: σk 1 = σ − ln p −N Γ(k/N, p) + N ln N −k N −k − k ln N k . (6.15) 45 Bibliography [1] D. J. Gross and E. Witten, “Possible Third Order Phase Transition In The Large N Lattice Gauge Theory,” Phys. Rev. D 21, 446 (1980). [2] K. G. Wilson, “Quark Confinement,” [3] J. B. Kogut and L. Susskind, “Hamiltonian Formulation Of Wilson’s Lattice Gauge Theories,” Phys. Rev. D 11, 395 (1975). [4] E. Witten, “Current Algebra Theorems For The U(1) Goldstone Boson,” Nucl. Phys. B 156, 269 (1979). [5] A. V. Manohar, “Large N QCD,” arXiv:hep-ph/9802419. [6] Sidney Coleman, “1/N,” Presented at The 1979 International School of Subnuclear Physics: Pointlike Structures Inside and Outside Hadron [7] I. Bars and F. Green, “Complete Integration Of U (N) Lattice Gauge Theory In A Large N Limit,” Phys. Rev. D 20, 3311 (1979). [8] P. Di Francesco, P. H. Ginsparg and J. Zinn-Justin, “2-D Gravity and random matrices,” Phys. Rept. 254, 1 (1995) [arXiv:hep-th/9306153]. 46
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Large N gauge theory and k-strings
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Large N gauge theory and k-strings Yang, Shuhang 2011
pdf
Page Metadata
Item Metadata
Title | Large N gauge theory and k-strings |
Creator |
Yang, Shuhang |
Publisher | University of British Columbia |
Date Issued | 2011 |
Description | We considered the k-antisymmetric representation of U(N) gauge group on two dimensional lattice space and derived the free energy by saddle point approximation in large N limit. k is a large integer comparable with N. Besides Gross-Witten phase transition[1], which happens as the coupling constant changes, we found a new phase transition in the strong coupling system that happens as k changes. The free energy of the weak coupling system is a smooth function of k under continuous limit. We have carefully selected the right saddle point solution among other possible ones. The numerical results match our saddle point calculations. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2011-04-14 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivatives 4.0 International |
IsShownAt | 10.14288/1.0071700 |
URI | http://hdl.handle.net/2429/33648 |
Degree |
Master of Science - MSc |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
GraduationDate | 2011-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/4.0/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 24-ubc_2011_spring_yang_shuhang.pdf [ 314.78kB ]
- Metadata
- JSON: 24-1.0071700.json
- JSON-LD: 24-1.0071700-ld.json
- RDF/XML (Pretty): 24-1.0071700-rdf.xml
- RDF/JSON: 24-1.0071700-rdf.json
- Turtle: 24-1.0071700-turtle.txt
- N-Triples: 24-1.0071700-rdf-ntriples.txt
- Original Record: 24-1.0071700-source.json
- Full Text
- 24-1.0071700-fulltext.txt
- Citation
- 24-1.0071700.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.24.1-0071700/manifest