arxiv: v1 [math-ph] 19 Oct 2007

Transcription

1 Monomal ntegrals on the classcal groups T. Gorn and G. V. López Departamento de Físca, Unversdad de Guadalajara Blvd. Marcelno García Barragan y Calzada Olímpca Guadalajara, Jalsco, Méxco February, 008 arxv: v1 [math-ph] 19 Oct 007 Abstract Ths paper presents a powerfull method to ntegrate general monomals on the classcal groups wth respect to ther nvarant (Haar) measure. The method has frst been appled to the orthogonal group n [J. Math. Phys. 43, 334 (00)], and s here used to obtan smlar ntegraton formulas for the untary and the untary symplectc group. The ntegraton formulas turn out to be of smlar form. They are all recursve, where the recurson parameter s the number of column (row) vectors from whch the elements n the monomal are taen. Ths s an mportant dfference to other ntegraton methods. The ntegraton formulas are easly mplemented n a computer algebra envronment, whch allows to obtan analytcal expressons very effcently. Those expressons contan the matrx dmenson as a free parameter. 1 Introducton Wth the classcal groups, we mean the orthogonal group O(d), the untary group U(d), and the untary symplectc group Sp(d) [1]. They all possess a unque nvarant measure, the Haar measure [], whch s the ntegraton measure commonly used. Integraton formulas for the classcal groups are of nterest n varous felds of mathematcal physcs. A number of dfferent classes of ntegrals have been studed, among those: generatng functons, such as the Harsh-Chandra-Izyson-Zuber ntegral [3, 4], and monomal ntegrals, whch are the concern of the present wor. In some cases, the correspondng ntegral can be solved by character expanson due to Balanten [5] and Balanten and Cassa [6]. Integraton formulas for certan smple monomals have been developed by nuclear physcsts [7, 8]. Ths wor was motvated by the fact that statstcal methods based on the classcal groups are very successful n descrbng certan aspects of nuclear reactons (see [9, 10] and references theren). Later, Mello and Selgman devsed an algebrac method to compute low order monomal ntegrals on U(d) [11], and Samuel solved the problem n full generalty [1]. The result was an explcte formula for arbtrary monomal ntegrals over U(d). The method n [1] s based on the representaton theory of U(d), and we call t the group theoretcal method. Some 5 years ago, t has become clear that group ntegrals on the classcal groups play an mportant role n mesoscopc transport [13, 14, 15], n quantum chaos [16, 17], and more recently n some aspects of quantum nformaton and decoherence, e.g. [18, 19, 0]. These applcatons lead to a renewed nterest n effcent methods for the analytcal calculaton of monomal ntegrals. An nterestng unconventonal approach has been devsed by Prosen et al. [1], and very recently, an mproved nvarant method has been developed for monomals n U(d) [, 3], and also n O(d) [4]. Fnally, Collns and Śnady presented a group theoretcal approach, whch allows to compute monomal ntegrals over all three classcal groups [5]. In [6], a further method has been developped, whch s very dfferent from the prevous ones; the column vector method as t mght be called. It lead to a new type of ntegraton formula for monomal ntegrals on O(d). Ths method s easly mplemented n a computer algebra language. It allows to compute arbtrary monomal ntegrals very effcently, where the result wll always be a ratonal functon n the matrx dmenson of the group. The purpose of the present paper s to apply that method to U(d) and Sp(d). In secton the matrx representatons of the classcal groups are ntroduced, and the notaton for the dfferent group ntegrals used s defned. The general dea of the column vector method, as well as the result for O(d) s revewed n secton 3. In the sectons 4 and 5 the ntegraton formulas for the untary and the untary symplectc group are derved. A summary s gven n secton 6. 0 gorn@ps.mpg.de 1

2 General consderatons In ths secton, we ntroduce the fundamental matrx representatons of the groups under consderaton. Then, we defne the ntegraton measures and specfy some notatonal conventons related to the ntegrals to be calculated..1 The orthogonal group O(d) An element w O(d) s a d dmensonal square matrx wth real entres (w j ), where 1, j d. In addton, w fullflls the followng orthogonalty relatons: w w T = w T w = 1 1, j d : N w w j = δ j. (1) =1 The most general monomal of matrx elements of w s denoted by M 1 (w). The subscrpt 1 wll be replaced by β = (4) n the untary (untary symplectc) case. Here, β s the symmetry parameter of the respectve group. M 1 (w) = p w IrJ r = r=1,j=1 w Mj j M j = q δ Ir δ jjr, () where q s the order of the monomal, and I and J are mult-ndeces of dmenson q. The notaton wth the mult-ndeces I, J s the most common one, but n the present approach the matrx notaton s more convenent. For the ntegraton on O(d), the normalzed Haar measure s used. For the classcal groups, n general, t wll be denoted by σ β, where β s the symmetry parameter ntroduced above. In the present case, the ntegral of an arbtrary monomal s denoted by M = dσ 1 (w) M 1 (w). (3) The normalzaton s such that o = 1, where o s the matrx wth zero elements everywhere. The column vector method has been appled to the orthogonal group, frst [6]. The results obtaned theren are revewed n secton 3. It allows us to ntroduce the general dea of the method, as well as a number of notatonal conventons. The latter wll be useful for the untary and untary symplectc case, also.. The untary group U(d) An element w U(d) s a d dmensonal square matrx wth complex entres (w j ), where 1, j d. In addton, w fullflls the followng orthogonalty condtons: r=1 w w = w w = 1 1, j d : d w wj = δ j. (4) =1 The most general monomal on the untary group depends on matrx elements of w and w. It wll be denoted by M (w): M (w) = p w IrJ r r=1 s=1 q wi s J = s,j=1 w Mj j (w j) Nj M j = p δ Ir δ jjr N j = r=1 q δ I r δ jj r, (5) where I, J and I, J are mult-ndeces of dmenson p and q, respectvely. Agan, whle the notaton wth the mult-ndeces I, J and I, J s the most common one, the matrx notaton wth M and N s more convenent for our purpose. For U(d), the normalzed Haar measure s σ (w), whle the monomal ntegral s denoted by N M = dσ (w) M (w). (6) The ntegraton formula for the untary group s derved n secton 4. s=1

3 .3 The untary symplectc group Sp(d) The untary symplectc group may be defned as the subgroup of U(d) whch s nvarant under the antsymmetrc bllnear form ( ) o Z ½d = : Sp(d) = { w U(d) w ½ d o T Z w = Z. (7) In order to fulfll ths nvarance condton, the matrces w must be of the form ( ) z w w = w, (8) z whth complex d-dmensonal square matrces w and z. For such matrces, the untarty condtons become: w w ν + z z ν = δ,ν z w ν w z = 0, (9) where w and z denote the respectve column vectors of the matrces w and z. Ths parametrzaton wll be used to perform the ntegraton over Sp(d). At ths pont we may note that ( ) ( ) z w Sp(d) : w 1 = w T w = T w z such that w w = ww ½d o =, (10) o ½ d where w s agan n Sp(d). The most general monomal on the untary symplectc group contans matrx elements from four matrces: w, w, z, and z. It s denoted by: M 4 (w) = M 4 (w, z) =,j=1 w Mj j z M j j (w j) Nj (z j) N j, (11) where M, M, N, N are d-dmensonal square matrces wth non-negatve nteger entres. If we collect these matrces n an approprate way n the -matrx, we may equally well wrte: M = ( ) N M N M : M 4 (w) = ( ) N,j=1 w Mj j, N = d,j=1 N j. (1) For Sp(d), we denote the normalzed Haar measure by σ 4 (w, z). The monomal ntegral s denoted by M = dσ 4 (w, z) M 4 (w). (13) 3 Integraton over O(d) For any of the classcal groups, the Haar measure s nvarant under left- and rght-multplcaton wth a fxed group element u, and under tang the nverse. For O(d) these propertes read: dσ 1 (w) f(w) = dσ 1 (w) f(u w) = dσ 1 (w) f(w u) = dσ 1 (w) f(w T ), (14) where f(w) s an arbtrary ntegrable functon of the matrx entres of w. For a monomal ntegral M, as defned n equaton (3), those operatons translate nto correspondng operatons on the nteger matrx M. They can be used to brng the monomal ntegral nto a more convenent form. In other words, we use transposton and/or column permutatons to collect all non-zero elements of M n the frst R d columns, mnmzng R. The general dea s to wrte the ntegral n terms of the full d -dmensonal Eucldean space of all matrx elements of w, and to mplement the restrcton to the group manfold by approprately chosen δ-functons. Snce the monomal M 1 (w) contans matrx elements from only the frst R column vectors of w, the ntegraton may be restrcted to those, as wll be shown wth the help of the fnal result, below. Meanwhle, one may as well assume that R = d. Denotng the column vectors by w 1, w,..., w R, one may wrte: M = dσ 1 (w) 1 (M) = ξ=1 d,r,ξ=1 w M ξ ξ { R dω( w ξ ) = 1 (M) 1 (o), (15) R { δ( w w ν ) dω( w ξ ) dw ξ δ( w ξ 1). (16) w M ξ ξ <ν 3

4 The matrx o n equaton (15) s a d-dmensonal square matrx, whch contans only zeros. The subscrpt n the symbol 1 ndcates the symmetry parameter β, whch s equal to one n the orthogonal case. In equaton (16) the ntegraton regon s the product space of R unt spheres wth constant measure dω( w ξ ). The orthogonalty and the normalzaton s mplemented wth the help of approprately chosen δ-functons. To convnce oneself that the rght hand sde of equaton (15) really yelds the Haar measure, t s suffcent to realze that for R = d the measure used n equaton (16) s ndeed nvarant under transposton as well as left- and rght-multplcaton wth any other fxed group element. Smlar arguments also apply n the case of the untary and the untary symplectc group. To evaluate the ntegral n equaton (16), we ntegrate over the last column vector w R. Ths ntegraton can be done n closed form, and the result s a lnear combnaton of monomals n the matrx elements of the frst R 1 column vectors. Ths means, the resultng expresson wll be of the followng form: 1 (M) = K C(K) N () 1 (M + K) 1 (o) = C(0) N () 1 (o), (17) where K s an nteger matrx wth non-zero elements n the frst R 1 columns only. The sum K runs over a fnte number of such matrces, and for N () 1 (M +K) only the frst R 1 columns (of M and of K) need to be taen nto account. The second relaton n equaton (17) taes care of the proper normalzaton of the measure. The rato between both equatons, yelds the desred recurrence relaton: M (R) = K C(K) M () + K. (18) We wll eventually use the supscrpt of the form M (R) to ndcate that only the frst R column vectors of the respectve matrx are taen nto account. An ntegral, such as the one n equaton (15) wll be called R-vector ntegral. Wth only mnor changes, ths termnology s also appled n the untary and untary symplectc case. The one-vector ntegral The startng pont for the recurrence relaton s the one-vector ntegral. It has frst been computed by Ullah [8]. In our notaton, t reads: m = dσ 1 ( w) w m = ( ) d 1 m/ ( 1 ) m /, (19) where (z) a = Γ(z + a)/γ(z) s the Pochhammer symbol, and m = d m. The ntegral m s non-zero only f all components m are even numbers. The R-vector ntegral The recurrence relaton for the R-vector case has been obtaned n [6]. It relates the desred R-vector ntegral to a lnear combnaton of ()-vector ntegrals. To obtan an explcte result, ths recurrence relaton must be contnued untl R = 1 s reached, where the one-vector result can be nserted. M = χ m:even κ ( m κ ) κ B ( m, κ ; d, ) ( m κ K) cs M () + K, (0) where B(a, b; z 1, z ) = ( ) a b (z 1) b (z 1 ) a b (z 1 z ) a. (1) In ths formula, m s the R-th column vector of M and m s the sum of all ts entres. The prefactor χ m:even s one f m s even, otherwse t s zero. The bnomal of the nteger vectors m and κ s just a short hand for the product of bnomals of ther respectve components: ( ) m ( ) m =. () κ The symbol, κ, denotes a one-vector average, as defned n equaton (19). Due to those two quanttes, the summaton s restrcted to such κ for whch : 0 κ m, κ :even. Analogous to m, κ denotes the sum of all components of the nteger vector κ. The second sum n equaton (0) runs over the d-dmensonal nteger matrx K wth non-negatve entres n the frst R 1 columns, only. The followng vector-multnomal s agan a short hand, whch reads: ( m κ K) = (m κ K 1,..., K, ), (3) 4 κ K

5 where the row sums n K must be equal to m κ. The vector cs denotes the vector of column-sums of K. It has at most non-zero components. Snce our formula contans the one-vector average of cs, the summaton over K should be restrcted to those K, whch have only even column sums. The last term s a ()-vector average, where the matrx entres of K have been added to the matrx entres n the frst R 1 columns of M. As t turns out, the functon B(a, b; z 1, z ) appears agan n the ntegraton formulas for U(d) and Sp(d), although wth slghtly dfferent arguments. Zero column vectors Wth the recurrence formula (0) at hand, t s straght forward to prove that the restrcton of the ntegraton n equaton (16) to the frst R column vectors of w s vald. Ths follows from the fact that the recurrence relaton yelds M = M () f all components of the R th column vector of M are zero. Vanshng ntegrals It s well nown (e.g. reference [10]), that the ntegral M vanshes f any column sum or row sum of M gves an odd number. Ths can be drectly seen from equaton (0), whch requres m to be even n order to yeld a non-zero result. Due to the nvarance propertes of the Haar measure, ths statement apples to any column- or row-vector. 4 Integraton over U(d) For the untary group, we denote the normalzed Haar measure by σ. It has analogous nvarance propertes as σ 1 of the orthogonal group, dscussed n the prevous secton. In the present case, one fnds for a fxed group element u U(d): dσ (w) f(w, w ) = dσ (w) f(u w, w u ) = dσ (w) f(w u, u w ) = dσ (w) f(w, w), (4) where f(w, z) s an analytc ntegrable functon of the matrx entres of w and z. These operatons translate nto correspondng operatons on the nteger matrces M, N, whch leave the monomal ntegral N M, equaton (6), nvarant. In partcular: () smultanous column permutatons: N M = N π M π, () smultaneous row permutatons: N M = πn πm, where π s an arbtrary permutaton matrx, and () conjugate transposton: N M = M T N T. Due to the nvarance under column permutatons, we may assume wthout loss of generalty that the non-zero elements of M and N are all restrcted to the frst R d columns. In other words, ths means that the monomal M (w) contans no matrx elements from column vectors w wth > R. Analogous to the case of the orthogonal group, we agan wrte the ntegral n terms of the flat Eucldean space of all complex matrx elements of w U(d). Then, we mplement the restrcton to the group manfold by approprately chosen δ-functons: N M = dσ (w) (M, N) = ξ=1 d,r,ξ=1 w M ξ ξ { R dω ( w ξ ) (w ξ) N ξ = (M, N) w M ξ ξ (w ξ )N ξ, (5) (o, o) R δ () ( w w ν ). (6) The subscrpt n the symbol stands for the untary case β =. The frst R column vectors of the untary matrx w are denoted by w 1,... w R. The ntegraton regon n equaton (6) s the product space of R unt spheres wth constant measure dω ( w ξ ): dω ( w ξ ) { <ν d (w ξ ) δ( w ξ 1), (7) where we defne the flat measure on the complex plane va z C : d (z) = d(re z)d(imz). The δ-functons n the equatons (6) and (7) mplement the orthogonalty condtons and the normalzaton. Note that the δ-functon n equaton (6) really s the product of two δ-functons: δ () (z) = δ(re z) δ(im z). Snce we assumed that the monomal M (w) contans matrx elements from only the frst R column vectors, we may restrct the ntegraton n (6) to those vectors. As n the orthogonal case, ths does not affect the result of the monomal ntegral, as wll be shown below wth the help of the fnal result. Meanwhle, one may as well assume that R = d. Note that for R = d the ntegraton measure n equaton (6) s ndeed nvarant under the transformatons performed n (4). Ths guarantees that equaton (5) really yelds the Haar measure. 5

6 4.1 The one-vector formula In the one-vector case, the matrces M and N can be replaced by the d-dmensonal vectors m and n. In that case, there are no orthogonalty condtons to obey. We may wrte: n m = dσ (w) (w = N (1) ( m, n) )n N (1) ( o, o) w m N (1) ( m, n) = dω ( w) w m (w )n. (8) As suggested n equaton (7), we ntegrate over the full space R d, whle the normalzaton s mplemented wth the help of a δ-functon. Ths ntroduces an ntegraton constant, denoted by C (d, 1): N (1) ( m, n) = C (d, 1) { The δ-functon s removed as follows: Settng w = u / r we get: N (1) ( m, n) r d+( m+ n)/ 1 = C (d, 1) d (w ) w m (w )n δ ( w 1 ). (9) { d (u ) u m (u ) n δ ( u r ), (30) where m = d m and n = d n. Multplyng both sdes of the equaton wth e r and ntegratng on r from 0 to, the δ-functon dsappears: N (1) ( m, n) Γ ( ) d + m+ n = C (d, 1) d (u ) u m (u )n e u = f(m, n, 0). (31) The Gaussan ntegral f(m, n, α) s defned n equaton (9) n the appendx. The general result for that ntegral s gven n equaton (100). For α = 0, t evaluates to From ths, t follows: N (1) ( m, n) Γ ( d + m+ n f(m, n, 0) = π δ m,n n!. (3) ) = C (d, 1) δ m, n {π m! N (1) ( o, o) Γ(d) = C (d, 1) π d, (33) where δ m, n denotes the product of Kronecer deltas between the components of m and n. Ths fnally leads to n m = δ m, n (d) 1 m m!. (34) 4. The R-vector formula To solve the R-vector ntegral, we start from equaton (6) and separate the ntegraton over the frst R 1 column vectors from the ntegraton over the last clumn vector w R. Let us denote the R-th column vectors of M, N and w wth m, n, and w, respectvely. (M, N) = ( m, n) = ξ=1 { dω ( w ξ ) dω ( w) w M ξ ξ (w ξ) N ξ <ν { δ () ( w w ν ) ( m, n) (35) {w m (w ) n δ () ( w w ). (36) We start by flattenng the ntegraton measure Ω, usng the technque from the one-vector ntegral. Wth d (w ) = d(re w )d(im w ) we may wrte: ( m, n) = C (d, R) { d (w ) w m (w ) δ n ( w w 1) δ () ( w w ). (37) 6

7 The transformaton w = u / r leads to: ( m, n) r d+( m+ n)/ 1 () = C (d, R) { d (u ) u m (u ) δ n ( u u r) δ () ( w u ). (38) Multplyng both sdes wth e r and ntegratng on r from 0 to gves: ( m, n) Γ ( ) { d R m+ n = C (d, R) d (u ) u m (u ) n e u u δ () ( w u ). (39) Usng the Fourer representaton of the one-dmensonal delta functon, we may wrte: dxdy δ () (w) = dxdy e π (x Im w+y Rew) = π e x(w w )+y(w+w ) = Wth the help of ths representaton, we wrte: δ () ( w u ) d (τ ) = Therefore π d (z) π e wz w z = d (z) π e Im(wz). (40) exp ( τ w u τ w ) u. (41) C (d, R) ( m, n) = Γ ( { d ) (τ ) d R m+ n π f(m, n, α ) α = f(m, n, α) = d (u) u m (u ) n e uu e α u α u p = π ( ) m ( ) κ κ! κ=0 ( m κ ) ( n κ ) (α ) m κ α n κ e α α p = mn(m, n). τ w (4) For detals about the computaton of that ntegral, see appendx. In order to put the product of the functons f(m, n, α ) nto a sutable form, we defne addtonal nteger vectors κ and p, where p = mn(m, n ): f(m, n, α ) = π d ( ) m p κ= o ( ) κ {( m κ ) ( ) n κ κ! (α ) m κ α n κ e P α. (43) Wth m and κ, we denote the sum of vector components n m and κ, respectvely. For the expanson of the powers of the coeffcents α and α, we need the nteger matrces K and L. Both are d-dmensonal matrces wth non-zero elements n the frst R 1 columns, only. (α ) m κ = = K ( ) m κ τ w = ( m κ K) K (τ) K w K ( m κ K) (τ ) w K, (44) where 1,... denote the column vectors of K, whle the overbars denote, as usual, the summaton over all vector components. In a completely analogous way, we wrte also: α n κ = L ( n κ L) τ l (w) L, (45) where l 1,... l denote the column vectors of L. Thus, we obtan from equaton (4): ( m, n) = C (d, R) π d ( ) m Γ ( ) d R m+ n p κ= o ( ) κ K,L( m κ K) ( n κ L) {( m κ ) ( ) n κ κ! d (τ ) π (τ ) τ l e P α w K (w ) L. (46) 7

8 Due to the fact that the ntegraton n equaton (6) s restrcted to column vectors w ξ whch are parwse orthogonal, we obtan: α = τ τν w ν w = τ. (47) ν Ths s a crucal pont n our calculaton. It elmnates the remanng matrx elements of w from the exponental, and ( m, n) turns nto a lnear combnaton of smple monomals. The ntegrals over τ are now easly evaluated, usng equaton (31) from the calculaton of the one-vector average: ( m, n) = C (d, R) π d+ ( ) m Γ ( ) d R m+ n p κ= o K,L( m κ K) ( n κ L) ( ) κ {( m κ δ, l! ) ( ) n κ κ! w K (w )L. (48) Now, we nsert ths result nto equaton (35) and wrte for (M, N) and (o, o): (M, N) = C (d, R) π d R+1 ( ) m p {( ) ( ) Γ ( ) ( ) κ m n d R m+ n κ! ( m κ K) ( n κ L) { δ, l! κ= o = C (d, R) π d R+1 ( ) m Γ ( ) d R m+ n { dω ( w ξ ) ξ=1 p κ= o K,L( m κ K) ( n κ L) (o, o) = C (d, R)π d R+1 Γ(d R + 1) Thus we obtan: N M = ( ) m (d R + 1) 1 ( m+ n)/ ( ) κ κ {( m κ { δ, l! κ w M ξ+k ξ ξ ) ( ) n κ κ! K,L (w ξ) N ξ+l ξ { <ν δ( w w ν ) N () (M () + K, N () + L) (49) N () (o, o). (50) p κ= o K,L( m κ K) ( n κ L) κ= o ( ) κ {( m κ ) ( ) n κ κ! { δ, l! N () + L M () + K. (51) Note that due to the multnomal coeffcents ( m κ K) and ( n κ L), for the column sums of K and L, t holds that: = m κ, and also l = n κ. In addton, to obtan a non-zero contrbuton, the column sums of K and L must be the same also: 1 R 1 : = l. Together, both condtons mply that for m n, the whole monomal ntegral gves zero. Therefore: p ( ) ( ) m n N M = δ m, n κ κ B( m, κ; d, R 1) κ κ K,L( m κ K) ( n κ L) l cs cs N () + L M () + K. (5) In ths formula, we use the vector notaton for the product of bnomals and multnomals, as ntroduced n secton 3. We also express the products d κ! and δ, l! n terms of the one-vector averages, gven n equaton (34). In ths way, the coeffcent B( m, κ; d, R 1), defned n equaton (1), comes nto play. The summaton over κ s restrcted to those vectors, for whch : 0 κ mn(m, n ). K and L are d-dmensonal matrces wth non-negatve nteger entres, whch have non-zero elements n ther frst R 1 columns, only. The summaton over K s restrcted to those K, for whch : K = m κ. Smlarly, the summaton over L s restrcted to those L, for whch : L = n κ. In addton, the one-vector average l cs cs mples that the column-sums of K must agree wth the correspondng column sums of L. These column sums are arranged n the frst R 1 components of the d-dmensonal vectors cs and l cs, respectvely. 8

9 Zero column vectors Assume that n the monomal ntegral N M, the R th column of M and N contan both only zeros. Equaton (5) then yelds: N M = N () M (). Ths proves by a smple nducton argument, that the ntegraton space n (6) may ndeed be restrcted to those column vectors, whch contan matrx elements from the monomal M (w). Vanshng ntegrals Denotng monomals n the form: w a 1b 1... w a pb p w α1β 1... w αqβ q (see for nstance, reference [14]), t has been shown that the monomal ntegral w a 1b 1... w a pb p w α1β 1... w αqβ q vanshes, unless p = q and α 1,...,α q s a permutaton of a 1,..., a p and β 1,..., β q s a permutaton of b 1,..., b p. Ths s equvalent to the requrement that all columns sums and row sums of M agree wth the correspondng column sums and row sums of N. In vew of the nvarance of the Haar measure under row- or column permutatons, ths s a trval consequence of the leadng factor δ m, n n equaton (5). 5 Integraton over Sp(d) For the untary symplectc group, we denote the normalzed Haar measure by σ 4. For ( ) z w Sp(d) w = w, let f(w) = f(w, z, w z, z ), (53) where f(w 1, z 1, w, z ) s an arbtrary analytc(?) and ntegrable functon of the matrx elements of w 1, z 1, w, and z. The nvarance propertes of the Haar measure then mply: dσ 4 (w, z) f(w) = dσ 4 (w, z) f(w u) = dσ 4 (w, z) f(u w) = dσ 4 (w, z) f(w ), (54) where u s a arbtary but fxed element of Sp(d). Snce we are nterested n the case, where the functon f s a monomal, we use some partcular permutaton matrces from Sp(d) to obtan nvarance relatons for the monomal ntegral M, under certan transformatons of M. For the followng consderatons, let ( ) ( ) z w v u w = w, u = z u. (55) v (1) For u = π, a d-dmensonal permutaton matrx, and v = o, multplcaton from the left leads to ( ) w π z w w u = π z π w, (56) π such that ( ) N M M = N M = dσ 4 (w) dσ 4 (w) (z π) Mj,j=1 j,j=1 ( w π) M j j w Mj j z M j j (wj) Nj (zj) N j ( ) (zπ) Nj j ( wπ) N j j = ( ) M + N M π 1 N π 1 M π 1 N π 1, (57) where M and N are the sum of all matrx entres n M and N, respectvely. Smlarly, multplcaton wth u from the rght yelds: ( ) M dσ 4 (πz) Mj j ( πw) M j j (πz ) Nj j ( πw ) N j j = ( ) M + N π 1 N π 1 M π 1 N π 1. (58) M,j=1 () For u = o and v = π, multplcaton from the left yelds: M dσ 4 (wπ) Mj,j=1 whereas multplcaton from the rght yelds: M j (zπ) M j j ( ) (w π) Nj j (z π) N j N j = π 1 Mπ 1 Nπ 1 M π 1, (59) ( ) π 1 N π 1 M π 1 N π 1 M. (60) 9

10 (3) Fnally, replacng w wth ts nverse, w, yelds: M dσ 4 (w) ( w j ) Mj (zj )M j ( w j ) Nj (z j ) N j = ( ) M+ N,j=1 ( M T M T N T N T ), (61) where M and N are the sums of all matrx entres of M and N, respectvely. As we wll show below [see equaton (88)], the monomal ntegral M s zero, unless M + M = N + N. Equaton (61) also mples that M s zero, unless M + N = N + M. Both condtons together mply that M 0 M = M = N = N. (6) We can therefore savely gnore the sgn-prefactors n (57), (58), and (61). The results n tem () prove the nvarance of M under synchronous column permutatons [equaton (59)] and synchronous row permutatons [equaton (60)], appled to all submatrces of M. Therefore, the equatons (57) and (58) add only two addtonal operatons: ( ) N M N M The ntegraton ( ) M N M, N ( ) N M N M ( ) N M N. (63) M It s agan possble to wrte the ntegral over Sp(d) n terms of an ntegral over the flat Eucldean space of all complex matrx elements of w and z, mplementng the restrcton to the group manfold by approprately chosen δ-functons: M = dσ 4 (w, z) M 4 (w) = 4 (M) 4 (o), M 4(w) = d,r,ξ=1 w M ξ ξ z M ξ ξ (w ξ )N ξ (z ξ )N ξ. (64) Here, we assume that the monomal M 4 (w) contans no matrx elements from the column vectors w, z wth > R. One may have used the nvarance of M under the above dscussed transformatons n order to mnmze R. The matrx o s of the same shape as M, but all ts entres are equal to zero. 4 (M) = R dω 4 ( w ξ, z ξ ) M 4(w) ξ=1 R δ () ( w w ν + z z ν ) δ () ( z w ν w z ν ). (65) <ν The subscrpt n the symbol 4 s the symmetry parameter β = 4, snce we consder here the untary symplectc case. The frst R column vectors of the matrces w, z are denoted by w 1,... w R, and z 1,... z R, respectvely. The ntegraton regon n equaton (65) s the product space of R unt hyperspheres wth constant measure dω 4 ( w ξ, z ξ ): dω 4 ( w ξ, z ξ ) { d 4 (w ξ, z ξ ) δ( w ξ + z ξ 1) w, z C : d 4 (w, z) = d (w) d (z), (66) where d (z) denotes the flat measure on the complex plane, as already used n equaton (7). The δ-functons n (65) and (66) mplement the orthogonalty condtons and the normalzaton, as defned n equaton (9). Accordng to the assumpton, the monomal M 4 (w) contans matrx elements from the frst R columns of w and z, only. Therefore, we restrct the ntegraton n equaton (65) to those column vectors. The fact that ths does not affect the result of the ntegral, wll be shown only at the end of ths secton, where we have the fnal result, equaton (89) at our dsposal. Note that the ntegraton measure n (65) s ndeed nvarant under the transformatons performed n (54). Ths guarantees that the above constructon really mplements the Haar measure. 5.1 The one-vector formula In the one-vector case, the monomal M 4 (w) contans matrx elements from one column of w only. Ths means that the matrces M, M, N, N have only zero entres n all columns, except for one whch must be the same for all matrces. We denote the column vectors wth non-zero entres by m, m, n, and n and use the followng notaton for the monomal ntegral: n M = n m m = dω 4 ( w, z) w m z m (w ) n (z ) n. (67) 10

11 In order to evaluate that ntegral, we observe that there are no orthogonalty relatons to obey. The ntegral must therefore be equal to the one-vector ntegral q p over the untary group U(d), where p and q are the concatenatons of the vectors m, m and n, n, respectvely. Therefore: n n m m = δ m, n δ m, n (d) 1 m+ m where m and m denote the sum of all components of m and m, respectvely. 5. The R-vector formula m! m!, (68) We am at a recurson formula, whch relates an R-vector ntegral to a lnear combnaton of smpler ()-vector ntegrals. To derve such a relaton, we separate the ntegraton over the R-th column vectors of w and z from the rest. Let us denote the R-th column vectors of M, M, N, N by m, m, n, n. Then we may wrte: 4 (M) = ξ=1 { dω 4 ( w ξ, z ξ ) <ν 4 ( m, m, n, n ) = w M ξ ξ z M ξ ξ (wξ )N ξ (zξ )N ξ { δ () ( w w ν + z z ν ) δ () ( z w ν w z ν ) dω 4 ( w, z) { w m z m (w (z )n )n 4 ( m, m, n, n ) (69) δ () ( w w + z z ) δ () ( z w w z ). We start by flattenng the ntegraton measure Ω 4. Wth equaton (66) and d (z) = d(re z)d(imz), we may wrte: 4 (...) = C 4 (d, R) { d (w ) d (z ) w m z m (z )n (w ) n δ ( w + z 1 ) (70) δ () ( w w + z z ) δ () ( z w w z ). (71) The delta functon for the normalsaton can be removed, usng smlar steps as n the case of the untary group. Ths yelds: C 4 (d, R) 4 (...) = Γ ((d R + 1) + ( m + m + n + n)/) { d (u ) d (v ) u m v m (v ) n (u ) n e u v δ () ( w u + z v ) δ () ( z u w v ). (7) Accordng to equaton (40), the Fourer representatons for the remanng delta functons may be wrtten as δ () [ (w u + z v ) ] = δ () [ (z u w v )] = α = σ w β = d (σ ) π d (τ ) π e Im[σ(w u+z v)] = e Im[τ(zu wv)] = j=1 σ z γ = τ z ε = { d (σ ) π { d (τ ) π Im(αu+βv) e Im(γu εv) e τ w. (73) 11

12 Inserton nto equaton (7), and the exchange of the order of ntegraton gves: 4 (...) = C 4 (d, R) Γ ((d R + 1) + ( m + m + n + n)/) = C 4(d, R) Γ(...) d (u ) u m (u )n e u e Im[u(α+γ)] { d (σ ) d (τ ) π π d (v ) v m (v )n e v Im[v(β ε)] e { d (σ ) d (τ ) f(m, n, α + γ ) f(m, n, β ε ). (74) π π The functon f s the same whch appeared n connecton wth the R-vector ntegral for the untary group, equaton (4), and whch has been computed n the appendx. Insertng the result, equaton (100) and expandng the powers of (α + γ ) and (β ε ), we obtan: f(m, n, α + γ) = π e α+γ p ( m ( ) m κ κ κ=0 ) ( ) n κ! κ m κ p ( ) ( ) f(m, n, β ε) = π e β ε ( ) m κ m n κ! κ =0 κ κ =0 n κ ( ) ( ) m κ n κ (α ) (γ ) m κ α l γ n κ l l l=0 m κ =0 n κ l =0 ( m κ ) ( n κ ) (β ) ( ε ) m κ β l ( ε) n κ l, (76) where we have for the moment suppressed the ndex. In order to completely expand the product of the functons f(m, n, α + γ ) and f(m, n, β ε ), we defne addtonal vector ndeces: p = mn(m, n ), p = mn(m, n ), κ, and κ. The complete expanson then reads: f(m, n, α + γ ) f(m, n, β ε ) = π d e P ( α+γ + β ε ) p κ= o m κ ( ) m κ ( m κ n κ ( m κ = o l= o )( ) p n κ )( n κ l (α ) (γ ) m κ α l κ = o ) m κ ( )( ) ( ) m κ m n κ κ {κ! κ! = o n κ l = o γ n κ l ( m κ )( n κ ) l l (75) (β ) ( ε ) m κ β l ( ε ) n κ l. (77) The argument of the exponental functon can be smplfed, as follows: [ (σ w + τ z )(σνw ν + τν zν) + (σ z τ w )(σνz ν τν wν) ] α + γ + β ε = ν = [ σ σν ( w w ν + z z ν ) + τ σν ( z w ν w z ν ) + σ τν ( w z ν z w ν ) + τ τν ( z z ν + w w ν )] ν = ( σ + τ ), (78) where we have used the orthogonalty relatons n (9). For the expanson of the powers of coeffcents α, γ, etc., we wll need eght matrx ndces, each of dmenson dd. It wll be convenent to combne them nto two bloc-matrces, as follows: ( ) ( ) L K 1 = 1 K 1 L L 1 K 1, K = K L K. (79) 1

13 Then we combne pars of coeffcents, whch lead to monomals of entres from the same submatrx of w: (α ( ε ) ) n κ l = ( ) n κ l ( K 1 )( n κ K ) (σ ) 1 τ K 1,K (β ) γ n κ l = ( K 1)( n κ l K ) (σ) 1 τ K 1,K α l ( ε )m κ = ( ) m κ L 1,L ( l L 1 )( m κ L ) (γ ) m κ β l = ( l L 1)( m κ L ) L 1,L σ l 1 (τ) l w K1 +K (80) z K1 +K (81) σ l 1 (τ ) l (w )L1 +L (8) (z ) L1 +L, (83) where we eventually convert the subscrpts 1, of the eght matrx ndces nto superscrpts, especally when we refer to partcular matrx elements. The parameters 1 1,, etc., denote the column sums of the matrx ndces K 1, K 1, etc.. Pluggng the complete expanson nto (74), we obtan: 4 ( m, m, n, n ) = m κ = o m κ = o C 4 (d, R) π d m+ m ( ) Γ ( (N R + 1) + m+ m + n+ n n κ n κ l= o ( ) + l l = o ) p p κ= o κ = o ( ) ( m κ m κ ( ) ( ) ( ) ( ) m m n n ( ) κ+ κ κ κ κ κ {κ! κ! ) ( n κ ) ( n κ l l ( K 1 )( K 1 )( l L 1 )( l L 1 ) ( n κ l K )( n κ l K )( m κ L )( m κ L ) K 1,K { d (σ ) d 1 1 (τ ) + l 1 + σ l + + l τ e P ( σ + τ ) π π w K1 +K z K1 +K (σ ) 1 (w )L1 +L (z ) L1 (τ ) l +L ). (84) We can now merge each sum and ts correspondng bnomal from the second lne wth two sums and ther correspondng multnomals from the thrd lne: m κ = o ( ) m κ ( K 1 )( m κ L ) = ( m κ L, K 1 ) K 1,L K 1,L m κ ( m κ ) ( K 1 )( m κ L ) = ( m κ L, K 1 ) = o K 1,L K 1,L n κ ( ) n κ ( l L 1 )( n κ l K l ) = ( n κ L 1, K ) l= o n κ l = o L 1,K ( n κ ) l L 1,K L 1,K ( l L 1)( n κ l K ) = L 1,K ( n κ L 1, K ), (85) where an extended vector-multnomal such as ( m K 1, L ) means that the rows of K 1 and L should be concatenated, to yeld a product of multnomals, each wth d elements. Evaluatng the fnal ntegrals on σ and τ, 13

14 we obtan: 4 ( m, m, n, n ) = C 4(d, R) π d R+1 m+ m ( ) Γ ( ) (N R + 1) + m+ m + n+ n K 1,K ( ) ( m κ m κ L K 1 n κ L K 1 n κ L ) 1 K L 1 K w K1 +K z K1 +K (w ) L1 +L (z ) L1 +l p p κ= o κ = o ( ) K 1 + L 1 ( ) ( ) ( ) ( ) m m n n ( ) κ+ κ κ κ κ κ {κ! κ! δ l1 + l 1, δ l + l, + ( )! ( + )!. (86) where, we have replaced and l wth K 1 and L 1, as a more approprate notaton for the sum over all matrx entres n K 1 and L 1, respectvely. We have also further expanded the vector multnomals, concatenatng the matrces and vectors n the vertcal drecton. Now, we may nsert ths result nto equaton (69) n order to obtan 4 (M) as a lnear combnaton of terms whch contan the (R-1)-vector ntegrals N () 4 (M () +K 1 +K ). Note that the matrx ndces n (86) are chosen n order to ft nto that scheme. Dvdng through the correspondng expressons for 4 (o) and N () 4 (o) we obtan: m+ m Γ ((d R + 1)) ( ) p, p ( ) ( ) ( ) ( ) m m n n κ M = ( ) κ+ κ Γ ((d R + 1) + ( m + m + n + n )/) κ κ κ κ (d) κ+ κ κ κ κ κ, κ = o ( ) ( m κ m κ L K 1 n κ L K 1 n κ L ) 1 K l L ( ) K 1 1 K + L 1 (d) cs 1 m κ+ m κ l 1 cs cs cs K 1,K M () + K 1 + K. (87) The vectors 1 cs, cs, l 1 cs, and l cs are d-dmensonal vectors whch contan n ther frst R 1 components the column sums of K 1 and K n the followng order: 1 cs contans the column sums of K 1 and K 1; cs contans those of K and K ; l 1 cs those of L 1 and L 1; and l cs those of L and L. Note that due to the one vector ntegral n these column sums, one may conclude that = l l + l + l m κ+ m κ = n κ+ n κ m+ m = n+ n. (88) Ths allows us to use agan the functon B(a, b; z 1, z ) defned n equaton (1): p, p ( ) ( ) ( ) ( ) m m n n κ M = δ m+ m, n+ n κ κ κ κ B ( m + m, κ + κ ; d, (R 1)) κ κ κ κ, κ = o ( ) ( m κ m κ L K 1 n κ L K 1 n κ L ) 1 K l L ( ) K 1 1 K + L 1 cs 1 l 1 cs cs M () + K 1 + K. (89) cs K 1,K The frst lne of ths expresson agrees precsely wth the correspondng part of the recurrence relaton for the untary group U(d). The dfferences appear n the second lne, where we have to sum over -bloc matrces K 1 and K, wth the restrcton that the row-sums of the bloc-matrces n the followng vector-multnomals must agree wth the correspondng component of the gven d-dmensonal vector. Zero column vectors In the case that the R th column vector of all four matrces M, M, N, and N contan only zeros, the summaton over κ and κ contans only one term, namely κ = κ = o. Smlarly, the summaton over K 1, K also contans only one term, where all matrx entres of K 1 and K are zero. Therefore, we fnd ( ) ( ) N M N () M M = N M = () = M (). (90) N () M () As n the case of O(d) and U(d) n the prevous sectons, ths justfes equaton (65) where we have gnored column vectors of w and z whose entres do not occur n M 4 (w). 14

15 Vanshng ntegrals The monomal ntegral M obvously vanshes f m+ m s dfferent from n+ n. Tang nto account the transformatons from (54), whch do not change the ntegral, ths means that the ntegral vanshes unless the column sums of M fulfll the condton: 1 j d : m j = m d+j, (91) where m j s the column sum of column j of M. An exactly analogous condton holds for the row sums. 6 Conclusons We appled the column vector method, orgnally developed for the orthogonal group n [6], to monomal ntegrals over the untary and the untary symplectc group. As n the orthogonal case, we obtaned recursve ntegraton formulas, where the recurson parameter s the number of non-empty column vectors (these are those matrx columns of the group element, whch contan at least one matrx element from the monomal). In all three cases, the recurson formulas are of a very smlar form and relatvely easy to mplement n a computer algebra envronment. Ths provdes an effcent way to compute such ntegrals analytcally as a functon of the matrx dmenson. The approach presented here s very dfferent from the group theoretcal approach developped n [5], whch also provdes explcte formulas for the analytcal calculaton of arbtrary monomal ntegrals over the classcal groups. We dd not yet compare the effcency of both methods, but we would expect that the result of such a comparson wll depend on the problem at hand. A Gaussan ntegral on the complex plane Here, we compute the Gaussan ntegral from equaton (4). Wrtng for the complex ntegraton varable u = x + y we get: f(m, n, α) = dxdy u m (u ) n e uu e α u α u. (9) Frst note that f(m, n, α) = f(n, m, α), whch mples that t s suffcent to consder the case m n. Then, we use the complex dfferental operators z and z to wrte: f(m, n, α) = ( ) n α m α n dxdy e uu e α u α u z = x y z = x + y. (93) Wth α = α 1 + α, t follows f(m, n, α) = ( ) n m α n α dx e (x α x) dy e (y α 1 y) = ( ) n m α n α π e α α = π ( ) m+n n α (α ) m e α α = π ( ) m+n P(m, n) e α α, (94) where P(m, n) s a polynomal of total order m + n, wth the property that the dfference between the power of the conjugated varable and that of the non-conjugated one s m n, fxed. Hence, for P(m, n) we consder the followng Ansatz: P(m, n) = n =0 C (m,n) (α ) m n+ α = (a ) m n n =0 C (m,n) (α α). (95) For n = 0, the polynomals P(m, 0) are nown, whle for m n > 0, they can be computed by recurson: n = 0 : P(m, 0) = 1 C (m,0) 0 = 1 m > 0 : P(m, n) = e α α n α (α ) m e α α P(m, n + 1) = e α α n α P(m, n) e α α = α P(m, n) + n =0 C (m,n) (m n + ) (α ) m (n+1)+ α. Equatng the coeffcents for correspondng powers, we obtan for m n 0: (m n) C (m,n) C (m,n+1) 0 : = 0 = (m n + ) C (m,n) C (m,n) 1 : 1 n C (m,0) 0 = 1. (97) C n (m,n) : = n + 1 (96) 15

16 Ordered accordng to n and, these coeffcents can be arranged n a two-dmensonal scheme, formng a Pascal trangle. From ths scheme t s easly seen that ( ) C (m,n) n = ( ) (m n + + 1)(m n + + )... m. (98) Therefore, we fnally obtan P(m, n) = n ( ) m! n! (m n + )!! (n )! (α ) m n+ α = =0 n ( ) n m! n! (m )!! (n )! (α ) m α n. (99) As mentoned before, here we have assumed that m n. For m < n by contrast, we use: f(m, n, α) = f (n, m, α) so that n summary, we may wrte: References f(m, n, α) = π ( ) m p ( )! =0 ( m ) ( n =0 ) (α ) m α n e α α [1] H. Weyl, The classcal groups (Prncton Unversty Press, Prncton, 1939). [] A. Haar, Ann. Math. 34, 147 (1933). [3] Harsh-Chandra, Am. J. Math. 79, 87 (1957). [4] C. Itzyson and J.-B. Zuber, J. Math. Phys. 1, 411 (1980). [5] A. B. Balanten, J. Math. Phys. 5, 08 (1984). [6] A. B. Balanten and P. Cassa, J. Math. Phys. 43, 604 (00). [7] N. Ullah and C. E. Porter, Phys. Rev. 13, 948 (1963). [8] N. Ullah, Nucl. Phys. 58, 65 (1964). [9] C. E. Porter (ed.), Statstcal theores of spectra: fluctuatons (Academc, New Yor, 1967). p = mn(m, n). (100) [10] T. A. Brody, J. Flores, J. B. French, P. A. Mello, A. Pandey, and S. S. M. Wong, Rev. Mod. Phys. 53, 385 (1981). [11] P. A. Mello and T. H. Selgman, Nucl. Phys. A 344, 489 (1980). [1] S. Samuel, J. Math. Phys. 1, 695 (1980). [13] P. A. Mello, J. Phys. A: Math. Gen. 3, 4061 (1990). [14] P. W. Brouwer and C. W. J. Beenaer, J. Math. Phys. 37, 4904 (1996). [15] P. A. Mello and N. Kumar, Quantum transport n mesoscopc systems: complexty and statstcal fluctuatons: a maxmum-entropy vewpont (Oxford Unversty Press, Oxford, 004). [16] F. Haae, Quantum sgnatures of chaos (Sprnger, Berln, 1991). [17] M. L. Mehta, Random matrces and the statstcal theory of energy levels (Academc Press, New Yor, 1991). [18] T. Gorn and T. H. Selgman, J. Opt. B: Quantum Semclass. Opt. 4, S386 (00). [19] T. Gorn and T. H. Selgman, Phys. Lett. A 309, 61 (003). [0] R. F. Abreu and R. O. Vallejos, Phys. Rev. A 75, (007). [1] T. Prosen, T. H. Selgman, and H. A. Wedenmüller, J. Math. Phys. 43, 5135 (00). [] S. Aubert and C. S. Lam, J. Math. Phys. 44, 611 (003). 16

17 [3] S. Aubert and C. S. Lam, J. Math. Phys. 45, 3019 (004). [4] D. Braun, J. Phys. A: Math. Gen. 39, (006). [5] B. Collns and P. Śnady, Commun. Math. Phys. 64, 773 (006). [6] T. Gorn, J. Math. Phys. 43, 334 (00). 17