arxiv: v1 [math-ph] 1 Nov 2016

Transcription

1 Revisiting SUN) integrals Jean-Bernard Zuber Sorbonne Universités, UPMC Univ Paris 06, UMR 7589, LPTHE, F-75005, Paris, France & CNRS, UMR 7589, LPTHE, F-75005, Paris, France Abstract arxiv: v [math-ph] Nov 206 In this note, I revisit integrals over SUN) of the form DU U ij U ipj p U k l U k nl n. Whe the case p = n is well known, it seems that explicit expressions for p = n N had not appeared in the literature. Simarities and differences, in particular in the large N limit, between the two cases are discussed.

2 Introduction and results In this note, we consider the SUN) integrals DU U i j U ipj p U k l U k nl n, ) with DU the normalized Haar measure, or their generating functions Z p,n J, K) = DU tr KU)) p tr JU )) n, 2) where J and K are arbitrary N N matrices. Such integrals, mainly over the group UN), have been the object of numerous publications in the past, in the context of lattice gauge theories [, 2, 3, 4], or in the large N limit [5, 6], or for their connections with combinatorics [7, 8]. Integrals over SUN) seem to have received less attention, see, however, [2, 9]. Recent work by Rossi and Veneziano [0] has prompted this new investigation. Let us first recall the physical motivations for studying the SUN) integrals ) or 2). Such integrals appear in the context of lattice calculations of baryon spectrum. Indeed, consider a SUN) lattice gauge theory, with link variables denoted U l SUN), Wson action S = β plaquettes tr U P, and lattice averages := l DU l )e S / l DU le S. Following [0], introduce the baryonic Wson loop or book observable B := ɛ i i N ɛ j j N N a= U a) i aj a where the ordered products U a) = l C a) U l, a =, N, stand for N static) quark lines joining two points A and B, a distance r apart, see Fig.. B represents a baryon made of N quarks, created at A and annihated at B. To lowest order in a small β strong coupling) expansion, B = β A ɛ i i N ɛ j j N DV N a= V i aj a, with A the total area of the N sheets of the book, and the last V -integration is carried out on the junction of these sheets and is given by Z N,0. Higher order corrections in β may involve some other integrals Z nn,n. C ) C i B i 3 i 2 j A j 3 j 3) 2 2) V Figure : A book observable for N = 3. C By Z N invariance, it is clear that the above integrals vanish if p n 0 mod N. 3) From a representation theoretic point of view, the number of independent terms in 2), i.e., of independent tensors with the right symmetries in )), is given by the number of invariants in

3 N) p N) n, where N) and N) denote the fundamental N-dimensional representation of SUN) and its complex conjugate. Z n,n J, K) is a well known function of traces of powers of) JK, Weingarten s function []), at least for n < N, and one may collect all Z n,n s into Z W J, K) = n 0 κ 2n n!) 2 Z n,nj, K). 4) For the convenience of the reader, a certain number of known results on these integrals and their generating function are recalled in the Appendix. This case n = p in eqs. ), 2) wl be referred to as the ordinary case. We now turn to the determination of the Z nn,n. First, for n = 0, DU U i j U in j N is given by the only invariant in N N), namely by the totally antisymmetric tensor product, i.e., the determinant of U, equal to in SUN)). Hence DU U i j U in j N = Aɛ i i N ɛ j j N with a constant A determined by contraction with ɛ i i N and use of ɛ i i N U i j U in j N = ɛ j j N det U, hence DU det U = DU = = Aɛ j j N ɛ j j N = AN!, and A = /N!. Thus DU U i j U in j N = N! ɛ i i N ɛ j j N. 5) Equivalently, Z N,0 J, K) = det K. 6) More generally, for n < N, Z Nn,n J, K) is the product of det K by a polynomial of M := JK, invariant under M V MV, for any V SUN), and homogeneous of degree n in M, hence a polynomial that may be expanded on traces of powers of M Z Nn,n = DU tr UK)) Nn tr U J)) n = det K d α t α. 7) α n with a sum over partitions of n denoted α = [ α, 2 α 2,, p αp ], q q α q = n, and with the notations t α := q t αq q, t q := tr JK) q. 8) The coefficients d α are determined through recursion formulae resulting from a contraction of ) with a Kronecker delta: The first coefficients are ready determined δ jk 2 J lk K ji Z Nn,n = N n)n δ Z Nn,n. 9) n = d [] = n = 2 d [2] = N d [ 2 ] = N N n = 3 d [3] = 4 NN ) n = 4 d [4] = 30 NN )N 2) d [,3] = etc. d [,2] = 3N) NN ) 82N) 2 3) N)NN )N 2) d [2 2 ] = d [ 2,2] = 6N)2 4) NN )N 2) d [ 4 ] = N)4 8N) 2 6 N)NN )N 2) 2 d [ 3 ] = N)2 2 NN ) 0) 3N) 2 6) N)NN )N 2)

4 These coefficients are in fact simply related to the analogous coefficients z α in the expansion of the ordinary generating function, see Appendix, eq. A.3). If z α, α n, is written in the form z α = with P α N) a polynomial of N, then d α = P α N) N 2 N 2 ) N 2 n ) 2 ) ) P α N ) N ) N n 2)). 2) This follows from the comparison between the two systems of recursion formulae, see below sect. 2, and may be verified on the first coefficients 0) and A.5). Note that the coefficients d α decrease more slowly than the z α, for fixed n, as N grows. This behavior has consequences on the large N limit of the generating function Z D defined by Z D = κ n n! Z Nn,n = DU e κ tr UK tr U J tr UK) N =: Z D det K. 3) n 0 Note also the different summations in 4) and 3): here κ is a homogeneity parameter for the eigenvalues of M = JK whe in 4), it is κ 2 that plays that role. In both cases, however, we take κ of order N and set κ = N κ. Then whe in the ordinary case, a non trivial limit is obtained by taking the traces t n also of order N, resulting in an exponentation Z W = exp N 2 W W, here one finds that the traces t n have to be kept of order and then Z D = exp NW D, with W D a non trivial function of the t n s. Indeed one finds for the first terms lim N N log Z D det K = κt κ2 2 t2 t 2 ) κ3 3 t3 3t t 2 2t 3 ) κ4 4 t4 6t 2 t 2 2t 2 28t 3 t 5t 4 ) 4) which is just the beginning of a simple formula in the large N limit: W D := lim N N log Z D det K = κ n ) n α q n )! n Catp ) t p ) αp α q )! α n α n p p! in terms of the Catalan numbers Catm) = 2 The recursion relations 2m)! m!m)!. This wl be proved below in sect. 3. Let us return to 9) and look at the action of the differential operator t α det K of 7). Making use of the identities K ji 5) on a typical term det[k] K ji = CofK) ji 6) CofK) ji KX) jl = det K)X 7) one may write for each term in the expansion of Z Nn,n 2 t α det K) = p qα q KJK) q ) jl t αq) det K K ji J lj = p p qα q N )JK) q K ji q s= q 2 α q α q )JK) 2q t α q) det K 2 ) t s JK) q s t αq) det K 8) 3 q<r p qrα q α r JK) qr t αq,r) det K

5 with α q) := α q αq p αp ), α q) := α q αq 2 p αp ), and α q,r) := α q αq r αr p αp ). In this way, 9) yields an overdetermined) system of relations between the coefficients d α at ranks n and n, namely { p d α qα q N )JK) q α n p q s= q 2 α q α q )JK) 2q t α q) detk) 2 = nn n)δ α n d α t α detk). ) t s JK) q s t αq) detk) 9) q<r p } qrα q α r JK) qr t αq,r) detk) Compare these equations with those satisfied by the coefficients z α in the simar expansion of Z n,n J, K) in the ordinary case, see A.4). Their structure is the same, except for changes in the terms of 9) marked in red. In the linear system on the d α, the parameter N in the lhs of A.4) has been changed into N whe the right hand side is multiplied by N n). As a result the solutions of the d α linear system are obtained from those of the z α one by ) for a given n d α = N n)n n ) N ) z α N N. If z α is written as in ), it follows that d α has the form 2), qed. In particular, this relation implies that z α and d α have the same overall sign, namely ) #cyclesα)n. 3 Exponentiation and large N limit The differential equation 9) carries over to the generating function Z D of 3) in the form 2 δ jk Z D = δ N )κ κ 2 ) Z D. 20) J lk K ji κ We write Z D = det K Z D, in which Z D is a function of M := J K invariant under M V MV for any V SUN), and we rewrite the differential equation as 2 ZD M jk N ) Z D = δ N )κ κ 2 ) ZD. 2) M lk M ji M li κ This may be reexpressed as a differential equation wrt the eigenvalues λ i of M generically M is diagonalizable). This is a standard procedure [5] with the result that for any i, i N λ i 2 ZD λ 2 i N ) Z D j i λ j ZD Z D λ j λ i λ j = N )κ κ 2 ) ZD. 22) κ Finally we write Z D = exp NW D, which results in ) λ i N 2 W 2 D λ 2 N 2 WD ) N )N W D N i j i W D λ λ i j W D λ j λ i λ j = N )κ Nκ 2 W D κ. 23) In the large N limit, we rescale κ = N κ, keeping all t n = i λn i of order, and after dropping the subdominant terms, we get for w i := W D the equation λ i w 2 i w i = κ κ 2 W D κ. 24) 4

6 Note that this is in contrast with the ordinary case where the term j i contributes in the large N limit [5].) Now κ is just an homogeneity variable of the λ s, and we may thus substitute κ W D κ for the w i s = j λ j W D λ j Assuming w := i λ iw i known, one finds = j λ jw j. The equation finally reduces to a system of algebraic equations λ i w 2 i w i = κ j λ j w j ). 25) and w is thus the root of λ i w 2 i w i = κ w) 26) w i = 4λ i κ w) 2λ i w = N 2 i 4λi κ w). 27) 2 Equation 26) should be compared with that of the generating function of Catalan numbers Catm) = 2m)! m!m)!, namely Ct) = n 0 Catn)tn, tc 2 t) Ct) = 0. We find that w i = κ w)c κ w)λ i ), hence w i = κ w) n 0 Catn) κ w)λ i ) n 28) and w = i λ i w i = n Catn ) κ w)) n t n. 29) Let f 0 =, f n := ) n Catn )t n for n, y := κw), then 29) reads y = κ n 0 f ny n, whose solution is given by Lagrange formula y = κ w) = n 0 w = n = κ n m n κ n = n κ n d ) n f m z n )! m) n z=0 dz m 0 κ n α n κ n n )! α n n )! n q α q)! q α q! n!f α n! n q α q)! q α q! f α where the multinomial coefficient appears naturally in the expansion of the n -th power, and the n! results from the n-th derivative of z n. Upon integration we get W D = n κ n ) n α q n )! n α n q α q)! q α q! Catp )t p ) αp which establishes 5). Note the simarity of this calculation with the relation between ordinary moments m n and non crossing or free) cumulants f n of a given distribution. The combinatorial or diagrammatical interpretation of 5) remains to be found. Acknowledgements It is a pleasure to thank Gabriele Veneziano for inspiring discussions. p 5

7 References [] D. Weingarten, Asymptotic Behavior of Group Integrals in the Limit of Infinite Rank, J. Math. Phys 9 978) [2] M. Creutz, On invariant integration over SUN), J. Math. Phys ) [3] I. Bars and F. Green, Complete integration of UN) lattice gauge theory in the large N limit, Phys Rev D20 979) ; S. Samuel, UN) integrals, /N and De Wit t Hooft anomalies, J. Math. Phys ) ; I. Bars, UN) integral for the generating functional in lattice gauge theories, J. Math. Phys ) [4] J.-M. Drouffe and J.-B. Zuber, Strong coupling and mean field methods in lattice gauge theories, Physics Reports ) 9 [5] É. Brézin and D. Gross, The external field problem in the large N limit of QCD, Phys. Lett. 97B 980) [6] K.H. O Brien and J.-B. Zuber, A note on UN) integrals in the large N limit, Phys. Lett. 44B 984) [7] B. Collins, Moments and Cumulants of Polynomial random variables on unitary groups, the Itzykson-Zuber integral and free probabity, IMRN ) ; B. Collins and P. Śniady, Integration with respect to the Haar measure of unitary, orthogonal and symplectic group, Comm. Math. Phys ) [8] J. Novak, Complete homogeneous symmetric polynomials in Jucys Murphy elements and the Weingarten function, Banach Center Publications, ) , abs/ ; P. Zinn-Justin, Jucys Murphy elements and Weingarten matrices, Lett. Math. Phys. 9 2) 200), 9 27, [9] J. Carlsson, Integrals over SUN), [0] G.C. Rossi and G. Veneziano, The string-junction picture of multiquark states: an update, Appendix In this Appendix, we recall a certain number of results on the Weingarten s functions Z n,n J, K) and their generating function Z W J, K) of 4). The Z n,n are well known functions of traces of powers of) JK, at least for n < N, see also [2, 3, 4]. For n N, the traces are no longer independent, or equivalently the tensors δ iπl) δ jρk), π, ρ S n, are no longer independent, and we have to deal with a pseudo-inverse of the Gram matrix, see [7, 8]. Consider the integrals on UN) or SUN), this is irrelevant here, assuming n < N) [, 4, 7] DU U i j U inj n U k l U k nl n = τ,σ S n C[σ]) = τ,σ S n n δ ialτa) δ jakτσa) a= Y Young diagr. Y =n 6 χ λ) )) 2 χ λ) [σ]) n! 2 s λ I) A.) p δ iqlτq) δ jqkτσq)

8 where χ λ) [σ]) is the character of the symmetric group S n associated with the Young diagram Y, a function of the class [σ] of σ); thus χ λ) ) is the dimension of that representation; s λ X) is the character of the linear group GLN) associated with Young diagram Y, that is a Schur function when expressed in terms of the eigenvalues of X; s λ I) is thus the dimension of that representation. Finally the coefficient C[σ]) wl be determined below. Alternatively, in terms of generating functions with sources J and K Z n,n J, K) = DUtr KU) n tr JU ) n = n! α C[α]) t α A.2) α n n! where α is the cardinal of class [α] in S n, thus α = p pαp α p!. In A.2), the factor n! comes from the sum over τ and the factor α from that over the elements σ [α]. Thus Z W J, K) := UN) DU exp[κtr KU JU )] = n=0 κ 2n n! z α t α with z α = α C[α]). By the same argument as in sect. 2, the z α satisfy the linear system of recursion relations { p z α qα q NJK) q α n p = nδ q s= q 2 α q α q )JK) 2q t α q) 2 α n z α t α. Explicitly, the first coefficients z α = α C[α]) read ) t s JK) q s t αq) q<r p n = z [] = N n = 2 z [2] = N 2 )N, z [ 2 ] = N 2 ) α n qrα q α r JK) qr t αq,r) } A.3) A.4) n = 3 z [3] = 4 N 2 4)N 2 )N, z [,2] = 3 N 2 4)N 2 ), z [ 3 ] = N 2 2 N 2 4)N 2 )N A.5) n = 4 z [4] = 30 N 2 9)N 2 4)N 2 )N, z [,3] = etc. 3N 2 6) z [2 2 ] =, z N 2 9)N 2 4)N 2 )N 2 [ 2,2] = z [ 4 ] = N 4 8N 2 6 N 2 9)N 2 4)N 2 )N 2 82N 2 3) N 2 9)N 2 4)N 2 )N 2 6N 2 4) N 2 9)N 2 4)N 2 )N, The overall sign of z α is ) #cyclesα) n and its large N behavior z α N 2n#cyclesα). To study the large N limit, we take κ = N κ and t p = Nτ p, with κ and τ p of order. Then Brézin and Gross [5] have shown that W W JK) := lim N N 2 log Z W J, K; N κ) exists and satisfies the coupled equations W W = 2 κ 2 λ i c) 2 N 2N 2 log[ κ 2 λ i c) 2 κ 2 λ j c) 3 2 ] c 4 i i,j { with c = N i κ2 λ i c) 2 for N i κ2 λ i ) 2 2 strong coupling ) 0 for N i κ2 λ i ) 2 2 weak coupling ) 7 A.6) A.7) A.8)

9 The solution has two determinations, in a strong coupling and in a weak coupling phase. Here we are concerned with the strong coupling regime in which we may expand W W = n κ 2n α n w α τ α A.9) with the first terms W W = κ 2 τ κ4 2 τ 2 τ 2 ) 2 κ6 3 2τ 3 3τ τ 2 τ 3 ) κ8 4 24τ 4 48τ 2 τ 2 9τ τ τ 3 5τ 4 ) A.0) and the general term given in [6] w α = ) n 2n 3 q α q)! 2n)! q 2q)! ) αq q!) 2 α q!. A.) 8