Arxiv:0705.3572V1 [Math.CA]

Home , Symmetrization

arXiv:0705.3572v1 [math.CA] 24 May 2007 ulda space Euclidean uain(ymti)group (symmetric) mutation rcs hs nre r h sa xoeta ucin fone of functions exponential usual the are entries whose trices, det( emnnso h same the of terminants om.Atsmercmliait xoeta ucin w deno (we functions exponential multivariate E correspondin Antisymmetric the and forms. functions exponential antisymmetr multivariate of metrized ratio a is the character to such according each Moreover, characters, functions. these representa symmetric dimensional are finite algebras Lie of in Characters appear solutions systems. antisymmetric integrable and Symmetric particles. identical eemn ore eisepnin nmliait ymti n an and symmetric multivariate functions; in exponential expansions series Fourier determine ye fatsmercadsmercmliait xoeta fun exponential multivariate symmetric and antisymmetric of types eemn nerlmliait ore transforms; Fourier multivariate integral determine aus hydtriemliait nt ore transforms. Fourier finite multivariate determine they values; emtto group permutation UCIN N ORSODN ORE TRANSFORMS FOURIER CORRESPONDING AND FUNCTIONS λ − nmteaia n hoeia hsc,vr fe eda with deal we often very physics, theoretical and mathematical In h i fti ae st ecieadt td ymtie n ant and symmetrized study to and describe to is paper this of aim The a functions c may (a) we variable, one of functions exponential the of case the in As b functions (b) c functions (c) ucin b r niymti smerc ihrsett eleme to respect with (symmetric) antisymmetric are (b) Functions ( e x 2 ), π AT)YMTI UTVRAEEXPONENTIAL MULTIVARIATE (ANTI)SYMMETRIC i λ λ Abstract. ieFuirtasom.Egnucin fteitga Fo integral the of found. Eigenfunctions transforms. transform Fourier Fourier transform nite Fourier integral such series, of Fourier types corresponding three are There transforms. the conditio functions exponential multivariate boundary operator antisymmetric certain and Laplace satisfying the domains of fundamental eigenfunctions funct are determin exponential functions are as These entries defined whose matrices are of They terminants functions. exponential variate i x ( = j ) i,j n λ =1 1 E λ , E E ymti utvraeepnnilfunctions exponential multivariate Symmetric . E λ − λ n − 2 m edfieadsuysmerzdadatsmerzdmulti- antisymmetrized and symmetrized study and define We − ( λ , . . . , S ( hc r ymti ratsmercwt epc oteper- the to respect with antisymmetric or symmetric are which x ( x n x and ) and ) At)ymtiso ucin a r ecie yawider a by described are (a) functions of (Anti)symmetries . and ) n × n E ), n .KIY N .PATERA J. AND KLIMYK A. E S λ + E arcs( ento fatdtriat e below). see antideterminants of definition (a matrices x λ n + ( m + o xml,sc ucin eciecletosof collections describe functions such example, For . x ( ( = 1. x with ) ( x ,where ), Introduction with ) x 1 x , λ 2 1 x , . . . , ( = m x λ ( = 1 ( = λ , n )aedtriat of determinants are )) x 2 m 1 λ , . . . , x , 1 m , 2 x , . . . , ,admliait fi- multivariate and s, 2 ecrepn Fourier correspond re n m , . . . , re rnfrsare transforms urier oso n variable. one of ions ), s osymmetric To ns. :epnin into expansions s: ncorresponding on λ nsadantide- and ants n i aeafiiestof set finite a take ) cfunctions. ic ∈ in fsemisimple of tions n variable, E R ), elfruafor formula Weyl λ hs functions these ; ore trans- Fourier g + ucin nthe on functions ctions: m ( x h hoyof theory the nie three onsider i r antide- are ) tisymmetric ete by them te ∈ n t fthe of nts E × Z λ − they ; n ( isym- x ma- = ) 2 A. KLIMYK AND J. PATERA group, since exponential functions e2πimx, m Z, of one variable are invariant with respect to shifts x x + k, k Z. These∈ (anti)symmetries are described by → ∈ aff elements of the affine symmetric group Sn which is a product of the group Sn and the group Tn, consisting of shifts in the space En by vectors r = (r1, r2,...,rn), Z aff aff rj . A fundamental domain F (Sn ) of the group Sn is a certain bounded subset∈ of Rn (see subsection 2.2). + The functions Eλ (x) give solutions of the Neumann boundary value problem on − a closure of the fundamental domain F (Sn). The functions Eλ (x) are solutions of the Laplace equation ∆f = µf on the domain F (Sn) vanishing on the boundary ∂F (Sn) of F (Sn). aff Functions on the fundamental domain F (Sn ) can be expanded into series in the functions (a). These expansions are an analogue of the usual Fourier series for functions of one variable. Functions (b) determine an (anti)symmetrized Fourier integral transforms on the fundamental domain F (Sn) of the symmetric group Sn. This domain consists of points x En such that x1 > x2 > > xn. Functions (c) are used to determine∈ (anti)symmetric finite··· (that is, on a finite set) Fourier transforms. These finite Fourier transforms are given on grids consisting of aff points in the fundamental domain F (Sn ). Symmetric and antisymmetric exponential functions are closely related to symmetric and antisymmetric orbit functions defined in [1], [2] and studied in detail in [3] and [4]. In fact, symmetric and antisymmetric exponential functions are connected with orbit functions corresponding to the Coxeter–Dynkin diagram An. Discrete orbit function transforms, corresponding to Coxeter–Dynkin diagrams of low order, were detailly studied and it was shown that they are very useful for applications [5]–[13]. The exposition of the theory of orbit functions in [3] and [4] strongly depends on the theory of Weyl groups, properties of root systems, etc. In this paper we avoid this dependence. We use only the permutation (symmetric) group and properties of determinants and antideterminants. It is well-known that a determinant det(a )n of the n n matrix (a )n is defined as ij i,j=1 × ij i,j=1 det(a )n = (det w)a a a ij i,j=1 1,w(1) 2,w(2) ··· n,w(n) wX∈Sn

= (det w)a a a , (1) w(1),1 w(2),2 ··· w(n),n wX∈Sn where Sn is the symmetric group of n symbols 1, 2,...,n, the set (w(1), w(2),... , w(n)) means the set w(1, 2,...,n), and det w denotes a determinant of the transform w, that is, det w = 1 if w is an even permutation and det w = 1 otherwise. Along with a determinant, we shall use an antideterminant det+ of− the matrix n (aij )i,j=1 which is deﬁned as a sum of all terms, entering to the expression for the corresponding determinant, taken with the sign +,

det+(a )n = a a a = a a a . ij i,j=1 1,w(1) 2,w(2) ··· n,w(n) w(1),1 w(2),2 ··· w(n),n wX∈Sn wX∈Sn (2) Symmetrized and antisymmetrized multivariate polynomials were studied by several authors (see, for example, [14] and [15]). In this paper we investigate symmetric and antisymmetric multivariate exponential functions. (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 3

2. Symmetric and antisymmetric multivariate exponential functions

2.1. Deﬁnition. A symmetric multivariate exponential function of x = (x1, x2, ...,xn) is deﬁned as the function n E+(x) E+ (x) = det+ e2πiλixj λ (λ1,λ2,...,λn) i,j=1 ≡ e2πiλ1x1 e2πiλ1x2 e2πiλ1xn ···  e2πiλ2x1 e2πiλ2x2 e2πiλ2xn  =det+ ···  ············   e2πiλnx1 e2πiλnx2 e2πiλnxn   ···  e2πiλ1xw(1) e2πiλ2xw(2) e2πiλnxw(n) = e2πihλ,wxi, (3) ≡ ··· wX∈Sn wX∈Sn where λ = (λ1, λ2,...,λn) is a set of real numbers, which determines the function + Eλ (x), and λ, x denotes the scalar product in the n-dimensional Euclidean space E , λ, x =h ni λ x . When λ , λ ,...,λ are integers, we denote this set of n h i i=1 i i 1 2 n numbers as mP (m ,m ,...,m ), ≡ 1 2 n n E+ (x) = det+ e2πimixj . (4) m1,m2,...,mn i,j=1 It is seen from the expression (2) for an antideterminant det+ that the symmetric + exponential functions Eλ (x) satisfy the relation

+ 2πi(λ1+λ2+···+λn)a + Eλ (x1 + a, x2 + a,...,xn + a)= e Eλ (x). (5) + This means that it is enough to consider the function Eλ (x) on the hyperplane x + x + + x = b, 1 2 ··· n where b is a ﬁxed number (we denote this hyperplane by ). A transition from one Hb hyperplane b to another c is fulﬁlled by multiplication by a usual exponential function e2πHi|λ|(c−b), where Hλ = λ + λ + + λ , | | 1 2 ··· n E+(x) = e2πi(λ1+λ2+···+λn)(c−b)E+(x) . λ |x∈Hb λ |x∈Hc + It is useful to consider the functions Eλ (x) on the hyperplane 0. For x 0 we have the relation H ∈ H E+ (x)= E+(x). (6) λ1+ν,λ2+ν,...,λn+ν λ It is seen from the expression (2) for an antideterminant det+ that its expression does not change under permutations of rows or under permutations of columns. This means that for any permutation w S we have ∈ n + + + + Eλ (wx)= Eλ (x), Ewλ(x)= Eλ (x). (7) + Therefore, it is enough to consider only symmetric exponential functions Eλ (x) with λ = (λ1, λ2,...,λn) such that λ λ λ . 1 ≥ 2 ≥···≥ n Such λ are called dominant. The set of all dominant λ is denoted by D+. Below, considering symmetric exponential functions E+(x), we assume that λ D . λ ∈ + 4 A. KLIMYK AND J. PATERA

Antisymmetric multivariate exponential functions of x = (x1, x2,...,xn) are deﬁned as the functions n E−(x) E− (x) := det e2πiλixj λ (λ1,λ2,...,λn) i,j=1 ≡ e2πiλ1x1 e2πiλ1x2 e2πiλ1xn ···  e2πiλ2x1 e2πiλ2x2 e2πiλ2xn  det ··· ≡  ············   e2πiλnx1 e2πiλnx2 e2πiλnxn   ···  (det w)e2πiλ1 xw(1) e2πiλ2xw(2) e2πiλnxw(n) = (det w)e2πihλ,wxi, ≡ ··· wX∈Sn wX∈Sn (8) where λ = (λ1, λ2,...,λn) is a set of realnumbers. When λ1, λ2,...,λn are integers, we denote them as m = (m1,m2,...,mn),

− 2πimixj n Em(x) = det e i,j=1 . (9) − It is seen from properties of determinants that the functions Eλ (x) satisfy the relation

− 2πi(λ1+λ2+···+λn)a − Eλ (x1 + a, x2 + a,...,xn + a)= e Eλ (x), (10) that is, it is enough to consider functions E−(x) on some hyperplane . As in the λ Hb case of symmetric exponential functions, a transition from one hyperplane b to − H another c for the function Eλ is fulﬁlled by means of multiplication by a usual H 2πi|λ|(c−b) exponential function e , where λ = λ1 + λ2 + + λn. For points x of the hyperplane we have the relation | | ··· H0 E− (x)= E−(x). (11) λ1+ν,λ2+ν,...,λn+ν λ − It follows from properties of determinants that Eλ (x) = 0 if λ has at least two coinciding numbers or if x has at least two coinciding coordinates. For any permutation w S we receive ∈ n − − − − Ewλ(x) = (det w)Eλ (x), Eλ (wx) = (det w)Eλ (x). (12) − This means that it is enough to consider antisymmetric exponential functions Eλ (x) for λ = (λ1, λ2,...,λn) such that λ > λ > > λ . 1 2 ··· n + Such λ are called strictly dominant. The set of these λ is denoted by D+.

2.2. Aﬃne symmetric group and fundamental domains. We have seen that + the functions Eλ (x) are symmetric with respect to the permutation group Sn, that is, E+(wx) = E+(x), w S . The symmetric exponential functions E+ λ λ ∈ n m with integral m = (m1,m2,...,mn) admit additional symmetries related to the periodicity of the exponential functions e2πiry, r Z, y R. These symmetries are ∈ ∈ described by the discrete group of shifts in the space En by vectors r e + r e + + r e , r Z, 1 1 2 2 ··· n n i ∈ where e, e2,..., en are the unit vectors in directions of the corresponding axes. We denote this group by Tn. Permutations of Sn and shifts of Tn generate a group (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 5

aff aff which is denoted as Sn and is called the affine symmetric group. The group Sn is a semidirect product of its subgroups Sn and Tn, Saff = S T , n n × n −1 where Tn is an invariant subgroup, that is, wtw Tn for w Sn and t Tn. An open connected simply connected set F Rn∈is called a fundamental∈ ∈ domain aff ⊂ for the group Sn (for the group Sn) if it does not contain equivalent points (that is, ′ ′ aff points x and x such that x = wx, where w belongs to Sn or Sn, respectively) and aff if its closure contains at least one point from each Sn -orbit (from each Sn-orbit). aff Rn aff Recall that a Sn -orbit of a point x is the set of points wx, w Sn . + ∈ ∈ It is evident that the set D+ of all points x = (x1, x2,...,xn) such that x > x > > x 1 2 ··· n constitute a fundamental domain for the group Sn (we denote it as F (Sn)). The set of points x = (x , x ,...,x ) D+ such that 1 2 n ∈ + 1 > x > x > > x > 0 1 2 ··· n aff aff constitute a fundamental domain for the affine group Sn (we denote it as F (Sn )). + As we have seen, the functions Eλ (x) are symmetric with respect to the permu- + tation group Sn. This means that it is enough to consider the functions Eλ (x) only + on the closure of the fundamendal domain F (Sn). Values of Eλ on other points are received by using symmetricity. + Symmetricity of functions Em with integral m = (m1,m2,...,mn) with respect aff to the affine symmetric group Sn , E+ (wx + r)= E+ (x), w S , r T , (13) m m ∈ n ∈ n + means that we may consider Em(x) only on the closure of the fundamental domain aff F (Sn ), that is, on the set of points x such that 1 x1 x2 xn 0. + ≥ ≥ ≥ ··· ≥ ≥ Values of Em(x) on other points are obtained by using the relation (13). − The exponential functions Em(x) with integral m = (m1,m2,...,mn) also admit additional symmetries related to the periodicity of the usual exponential functions e2πiry, r Z, y R. These symmetries are described by the affine symmetric group aff ∈ ∈ Sn . We have E−(wx + r) = (det w)E−(x), w S , r T , (14) m m ∈ n ∈ n − that is, it is enough to consider the functions Em(x) only on the closure of the aff − fundamental domain F (Sn ). Values of Em(x) on other points are obtained by using the relation (14). The functions E−(x), λ = (λ , λ ,...,λ ), λ R, are antisymmetric with λ 1 2 n i ∈ respect to the symmetric group Sn, E−(wx) = (det w)E−(x), w S . λ m ∈ n − For this reason, we may consider Eλ only on the fundamental domain F (Sn). 2.3. Properties. Symmetricity and antisymmetricity of symmetric and antisymmetric multivariate exponential functions are main properties of these functions. However, they possess many other interesting properties. + Behavior on boundary. The symmetric and antisymmetric functions Eλ (x) − and Eλ (x) are finite sums of exponential functions. Therefore, they are continuous n functions of x1, x2,...,xn and have continuous derivatives of all orders in R . 6 A. KLIMYK AND J. PATERA

The closure F (Sn) of the fundamental domain F (Sn) without points of F (Sn) is called a boundary of the fundamental domain F (Sn) and is denoted by ∂F (Sn). A point x = (x , x ,...,x ) F (S ) belongs to ∂F (S ) if and only if at least two 1 2 n ∈ n n coordinates xi, xj in x coincide. It is clear that the boundary ∂F (Sn) is composed of points of F (Sn) belonging to the hyperplanes given by the equations x = x , i,j =1, 2,...,n, i = j. i j 6 aff aff Similarly, the boundary ∂F (Sn ) of the fundamental domain F (Sn ) consists of aff aff points of F (Sn ) which do not belong to F (Sn ). A point x = (x1, x2,...,xn) aff aff ∈ F (Sn ) belongs to ∂F (Sn ) if and only if at least two coordinates xi, xj in x coincide or if one of the conditions x1 = 1, xn = 0 is fulfilled. − It follows from properties of determinants that the function Eλ (x) vanishes on the boundary ∂F (Sn), E−(x)=0, for x ∂F (S ). λ ∈ n − Z − This relation is true for Em(x), mi . In this case we also have Em(x) = 0 for aff ∈ points x ∂F (Sn ) such that x1 xn = 1. ∈ − + For symmetric multivariate functions Eλ (x) we have ∂E+(x) λ =0 for x ∂F (S ) ∂n ∈ n where n is the normal to the boundary ∂F (Sn).

Complex conjugation. Let λ = (λ1, λ2,...,λn) be a strictly dominant element, that is, λ > λ > > λ . We have 1 2 ··· n − 2πi((wλ)1x1+···+(wλ)nxn) Eλ (x)= (det w)e , (15) wX∈Sn where (wλ)1, (wλ)2,..., (wλ)n are the coordinates of the point wλ. The element (λ , λ ,...,λ ) is strictly dominant if the element (λ , λ ,... , − n n−1 1 1 2 λn) is strictly dominant. In the group Sn there exists an element w0 such that

w0(λ1, λ2,...,λn) = (λn, λn−1,...,λ1)

It is easy to calculate that the set (λn, λn−1,...,λ1) is obtained from (λ1, λ2,...,λn) by n(n 1)/2 permutations of two neighboring numbers. Clearly, det w0 = 1 if this number− is even and det w = 1 otherwise. Thus, 0 − det w0 =1 for n =4k and n =4k +1, det w = 1 for n =4k 2 and n =4k 1, 0 − − − where k is an integer. It follows from here that in the expressions for the exponential functions E− (x) and E− (x) there are summands (λ1,λ2,...,λn) −(λn,λn−1,...,λ1)

e2πihw0λ,xi = e2πi(λnx1+···+λ1xn) and e−2πi(λnx1+···+λ1xn), (16) respectively, which are complex conjugate to each other. Moreover, the ﬁrst expression is contained with the sign (det w ) in E− (x), that is, the expressions 0 (λ1,λ2,...,λn) (16) are contained in E− (x) and E− (x) with the same sign (λ1,λ2,...,λn) −(λn,λn−1,...,λ1) for n =4k, 4k + 1 and with opposite signs for n =4k 2, 4k 1, k Z. − − ∈ (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 7

Similarly, in the expressions (15) for the function E− (x) and for the (λ1,λ2,...,λn) function E− (x) all other summands are (up to a sign, which depends −(λn,λn−1,...,λ1) on a value of n) pairwise complex conjugate. Therefore,

E− (x)= E− (x) (17) (λ1,λ2,...,λn) −(λn,λn−1,...,λ1) for n =4k, 4k + 1 and

E− (x)= E− (x) (18) (λ1,λ2,...,λn) − −(λn,λn−1,...,λ1) for n =4k 2, 4k 1. According− to (17)− and (18), if (λ , λ ,...,λ )= (λ , λ ,...,λ ), (19) 1 2 n − n n−1 1 − then the function Eλ is real for n = 4k, 4k +1 and pure imaginary for n = 4k 2, 4k 1. Moreover, the right hand side of (15) for this case consists of pairs of− terms− which give sine or cosine functions. It is representable as a sum of cosines of angles if n = 4k, 4k + 1 and as a sum of sines of angles multiplied by i = √ 1 if n =4k 2, 4k 1. − − − + It is proved similarly that for the symmetric exponential functions Eλ (x) and E+ (x) we have the following relation w0λ

E+ (x)= E+ (x). (λ1,λ2,...,λn) −(λn,λn−1,...,λ1) + If (λ1, λ2,...,λn) = (λn, λn−1,...,λ1), then the function Eλ (x) is real. In this + − case Eλ (x) can be represented as a sum of cosines of the corresponding angles. Scaling symmetry. For c R, let cλ = (cλ , cλ ,...,cλ ). Then ∈ 1 2 n − 2πihcwλ,xi 2πihwλ,cxi − Ecλ(x)= (det w)e = (det w)e = Eλ (cx). wX∈W wX∈W

− − The equality Ecλ(x)= Eλ (cx) expresses the scaling symmetry of exponential func- − tions Ecλ(x). The scaling symmetry is true also for symmetric exponential func- + + tions, Ecλ(x)= Eλ (cx). Duality. Due to invariance of the scalar product , with respect to the sym- h· ·i metric group Sn, wµ,wy = µ,y , for x = (x1, x2,...,xn), xi = xj , i = j, we have h i h i 6 6 − − 2πihλ,w 1xi 2πihx,wλi − Eλ (x)= (det w)e = (det w)e = Ex (λ). wX∈W wX∈W This relation expresses the duality of antisymmetric orbit functions. The duality is + + true also for symmetric exponential functions, Eλ (x)= Ex (λ). aﬀ Orthogonality on the fundamental domain F (Sn ). Antisymmetric expo- − + Z nential functions Em with m = (m1,m2,...,mn) D+, mj , are orthogonal aﬀ ∈ ∈ on F (Sn ) with respect to the Euclidean measure,

aff −1 − − ′ F (Sn ) Em(x)Em′ (x)dx = Sn δmm , (20) | | Z aff | | F (Sn ) 8 A. KLIMYK AND J. PATERA where the overbar means complex conjugation, Sn means a number of elements aff | | aff in the set Sn, and F (Sn ) is an area of the fundamental domain F (Sn ). This relation follows from| the equality|

− − ′ Em(x)Em′ (x)dx = Sn δmm ZT | | (where T is the torus in E consisting of points x = (x , x ,...,x ), 0 x < n 1 2 n ≤ i 1), which is a consequence of orthogonality of the exponential functions e2πihµ,xi − (entering into the definition of Em(x)) for different sets µ. T T aff −1 If to assume that an area of is equal to 1, = 1, then F (Sn ) = Sn and formula (20) takes the form | | | | | |

− − ′ Em(x)Em′ (x)dx = δmm . (21) Z aff F (Sn ) In the expression (3) for symmetric exponential functions there can be coinciding + summands. For this reason, for symmetric exponential functions Em with m = (m ,m ,...,m ) D , m Z, the relation (21) is replaced by 1 2 n ∈ + j ∈

+ + ′ Em(x)Em′ (x)dx = Sm δmm . (22) Z aff | | F (Sn ) where Sm is a number of elements in the subgroup Sm of Sn consisting of elements w S | such| that wm = m. ∈ n Orthogonality of symmetric and antisymmetric exponential functions. Let wi (i =1, 2,...,n 1) be the permutation of coordinates xi and xi+1. We create ext aff− aff aff aff the domain F (Sn ) = F (Sn ) wiF (Sn ), where F (Sn ) is the fundamental ∪ + domain for the affine group Saff . Since for m = (m ,m ,...,m ) D , m Z, n 1 2 n ∈ + j ∈ we have E+ (w x)= E+ (x) and E−(w x)= E−(x), then m i m m i − m + − Em(x)Em′ (x)dx =0. (23) Z ext aff F (Sn ) Indeed, due to symmetry and antisymmetry of symmetric and antisymmetric exponential functions, respectively, we have

+ − + − + − Em(x)Em′ (x)dx = Em(x)Em′ (x)dx + Em(x)Em′ (x)dx Z ext aff Z aff Z aff F (Sn ) F (Sn ) wiF (Sn )

+ − + − = Em(x)Em′ (x)dx + Em(x)( Em′ (x))dx =0. Z aff Z aff − F (Sn ) F (Sn ) The relation (23) is a generalization of the orthogonality of the functions sine and cosine on the interval (0, 2π). 2.4. Special cases. The special case of symmetric and antisymmetric exponential 1 functions at (λ1, λ2,...,λn) = 2 (n 1,n 3,..., n +3, n + 1) ρ is of great interest since it is met in the representatin− − theory.− The antisymmet− ric≡ exponential − function Eρ (x) is given by the formula E−(x) = (2i)n(n−1)/2 sin π(x x ). (24) ρ i − j 1≤i

It follows if to represent sin π(xi xj ) in terms of exponential functions, then to fulﬁl multiplication of these functions− and to compare with the expression (8) for − Eρ (x). (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 9

′ Let us set (m1,m2,...,mn) = (n 1,n 2,..., 1, 0) ρ . The antisymmetric − − − ≡ exponential function Eρ′ (x) can be written down in the form of the Vandermonde determinant, n − 2πi(n−i)xj 2πixk −2πixl Eρ′ (x) = det e = (e e ). (25) i,j=1 − k

3. Solutions of the Laplace equation

The Laplace operator on the n-dimensional Euclidean space En in the orthogonal coordinates x = (x1, x2,...,xn) has the form ∂2 ∂2 ∂2 ∆= 2 + 2 + + 2 . ∂x1 ∂x2 ··· ∂xn We take any summand in the expression for symmetric or antisymmetric multivariate exponential function and act upon it by the operator ∆. We get ∆e2πi((w(λ))1x1+···+(w(λ))nxn) = 4π2 λ, λ e2πi((w(λ))1x1+···+(w(λ))nxn), − h i − where λ = (λ1, λ2,...,λn) determines Eλ (x). Since this action of ∆ does not depend on a summand from the expression for symmetric or antisymmetric exponential function, we have ∆E−(x)= 4π2 λ, λ E−(x), ∆E+(x)= 4π2 λ, λ E+(x). (27) λ − h i λ λ − h i λ The formula (27) can be generalized in the following way. Let σk(y1,y2,...,yn) be the k-th elementary symmetric polynomial of degree k, that is, σ (y ,y ,...,y )= y y y . k 1 2 n k1 k2 ··· kn 1≤k1

∂2 ∂2 ∂2 ± 2 k 2 2 2 ± σk ∂x2 , ∂x2 ,..., ∂x2 Eλ (x) = ( 4π ) σk(λ1, λ2,...,λn)Eλ (x). (28) 1 2 n − Note that n diﬀerential equations (28) are algebraically independent. − Thus, antisymmetric exponential functions Eλ (x) are eigenfunctions of the op- ∂2 ∂2 ∂2 erators σk ∂x2 , ∂x2 ,..., ∂x2 , k = 1, 2,...,n, on the fundamental domain F (Sn) 1 2 n of the symmetric group Sn satisfying the boundary condition E−(x)=0 for x ∂F (S ). (29) m ∈ n 10 A. KLIMYK AND J. PATERA

+ Similarly, symmetric exponential functions Eλ (x) are eigenfunctions of the op- ∂2 ∂2 ∂2 erators σk ∂x2 , ∂x2 ,..., ∂x2 , k = 1, 2,...,n, on the fundamental domain F (Sn) 1 2 n satisfying the boundary condition ∂ E+(x) m =0 for x ∂F (S ), ∂ n ∈ n where n is the normal to the boundary ∂F (Sn). That is, they are solutions of the Neumann boundary value problem for the domain F (Sn).

4. Symmetric and antisymmetric Fourier transforms Symmetric and antisymmetric exponential functions determine symmetric and antisymmetric multivariate Fourier transforms which generalize the usual Fouirier transform. As in the case of exponential functions of one variable, (anti)symmetric exponential functions determine three types of Fourier transforms:

± (a) Fourier transforms related to the exponential functions Em(x) with m = (m ,m ,...,m ), m Z (Fourier series); 1 2 n j ∈ (b) Fourier transforms related to E±(x) with λ D ; λ ∈ + (c) symmetrized and antisymmetrized multivariate ﬁnite Fourier transforms.

aff 4.1. Expansions in (anti)symmetric exponential functions on F (Sn ). Let aff f(x) be symmetric (with respect to the affine symmetric group Sn ) continuous function on the n-dimensional Euclidean space En which has continuous derivatives. We may consider this function on the torus T which is a closure of the union of the sets wF (Saff ), w S . The function f(x), as a function on T, can be expanded in n ∈ n exponential functions e2πim1x1 e2πim2x2 e2πimnxn , m Z. We have ··· i ∈ f(x)= c e2πim1x1 e2πim2x2 e2πimnxn , (30) m ··· mXi∈Z where m = (m1,m2,...,mn). Due to symmetry f(wx) = f(x), w Sn, for any w S we have ∈ ∈ n 2πim − x 2πim − x f(wx)= c e2πim1xw(1) e2πimnxw(n) = c e w 1(1) 1 e w 1(n) n m ··· m ··· mXi∈Z mXi∈Z

= c e2πim1x1 e2πimnxn = f(x)= c e2πim1x1 e2πimnxn . wm ··· m ··· mXi∈Z mXi∈Z

Therefore, the coeﬃcients cm satisfy the conditions cwm = cm, w Sn. Collecting in (30) exponential functions at the same c , w S , we obtain∈ the expansion wm ∈ n + f(x)= cmEm(x), (31) mX∈P+ where P = D Z. Thus, any symmetric (with respect to S ) continuous function + + ∩ n f on T which has continuous derivatives (that is, any continuous function on D+ with continuous derivatives) can be expanded in symmetric exponential functions E+(x), m P . m ∈ + (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 11

By the orthogonality relation (22), the coeﬃcients cm in the expansion (31) are determined by the formula

−1 + cm = Sm f(x)Em(x)dx, (32) | | Z aff F (Sn ) where, as before, Sm is a number of elements in the subgroup Sm of Sn consisting of w S such that| wm| = m. Moreover, the Plancherel formula ∈ n 2 −1 2 cm = Sm f(x) dx (33) | | | | Z aff | | mX∈P + F (Sn ) holds, which means that the Hilbert spaces with the appropriate scalar products are isometric. Formula (32) is the symmetrized Fourier transform of the function f(x). Formula (31) gives an inverse transform. Formulas (31) and (32) give the symmetric multivariate Fourier transforms corresponding to the symmetric exponential functions + + Em(x), m P . ∈ − Analogous transforms hold for antisymmetric exponential functions Em(x), m + + Z ∈ P+ D+ . Let f(x) be antisymmetric (with respect to the symmetric group Sn) continuous≡ ∩ function on the n-dimensional torus T, which has continuous derivatives. aﬀ We may consider this function as a function on F (Sn ). Then we have the expansion

− − f(x)= cmEm(x), where cm = f(x)Em(x)dx. (34) Z aff X+ F (Sn ) m∈P+ Moreover, the Plancherel formula holds:

2 2 cm = f(x) dx. (35) | | Z aff | | X+ F (Sn ) m∈P+ 2 aff Let 0(F (Sn )) denote the Hilbert space of functions on the fundamental domain aff L aff aff F (Sn ), which behave on the boundary ∂F (Sn ) of the fundamental domain F (Sn ) − in the same way as the functions Em(x) do. Let

f1,f2 = f1(x)f2(x)dx h i Z aff F (Sn ) be a scalar product in this space. The formulas (34)-(35) show that the set of − + 2 aff exponential functions Em(x), m P+ , form an orthogonal basis of 0(F (Sn )). ext aff aff ∈ aff L Let F (Sn ) = F (Sn ) F (wiSn ) denote the set from section 2. Then we can extend the symmetric andS antisymmetric Fourier transforms to the functions from the Hilbert space 2(F ext(Saff )) with the scalar product L n

f1,f2 = f1(x)f2(x)dx. h i Z ext aff F (Sn ) This transform is of the form + ′ − f(x)= cmEm(x)+ cmEm(x), (36) mX∈P+ X+ m∈P+ where

−1 + ′ − cm = Sm f(x)Em(x)dx, cm = f(x)Em(x)dx. (37) | | Z aff Z aff F (Sn ) F (Sn ) 12 A. KLIMYK AND J. PATERA

+ The corresponding Plancherel formula holds. The functions Em(x), m P+, + ∈ and E−(x), m P , form a complete orthogonal basis of the Hilbert space m ∈ + 2(F ext(Saﬀ )). L n

4.2. Multivariate Fourier transforms on the fundamental domain F (Sn). aff The expansions (31) and (34) of functions on the fundamental domain F (Sn ) are respectively expansions in the symmetric and antisymmetric exponential functions + − Em(x) and Em(x) with integral m = (m1,m2,...,mn). The exponential functions + − Eλ (x) and Eλ (x) with λ lying in the fundamental domain F (Sn) (and not oblig- atory integral) are not invariant (anti-invariant) with respect to the corresponding aff affine symmetric group Sn . They are invariant (anti-invariant) only with respect to the permutation group Sn. A fundamental domain of Sn coincides with the set D+ consisting of the points x such that m >m > >m . For this reason, the + 1 2 ··· n functions E−(x), λ D+, and E+(x), λ D , determine Fourier transforms on λ ∈ + λ ∈ + D+. We began with the usual Fourier transforms on Rn:

f˜(λ)= f(x)e2πihλ,xidx, (38) ZRn

f(x)= f˜(λ)e−2πihλ,xidλ. (39) ZRn

Let the function f(x) be anti-invariant with respect to the symmetric group Sn, that is, f(wx) = (det w)f(x), w Sn. It is easy to check that f˜(λ) is also anti- invariant with respect to the group∈ S . Replace in (38) λ by wλ, w S , multiply n ∈ n both sides by det w, and sum these both side over w Sn. Then instead of (38) we obtain ∈

f˜(λ)= f(x)E−(x)dx, λ D+, (40) Z λ + D+ ∈ where we have taken into account that f(x) is anti-invariant with respect to Sn. Similarly, starting from (39), we obtain the inverse formula:

f(x)= f˜(λ)E−(x)dλ. (41) Z λ D+ For the transforms (40) and (41) the Plancherel formula

f(x) 2dx = f˜(λ) 2dλ Z Z D+ | | D+ | | holds. Formulas (40) and (41) determine the antisymmetric multivariate Fourier transform on the domain F (Sn). Similarly, starting from formulas (38) and (39) we receive the symmetric multivariate Fourier transform on the domain F (Sn):

f˜(λ)= f(x)E+(x)dx, f(x)= f˜(λ)E+(x)dλ. (42) Z λ Z λ D+ D+ The corresponding Plancherel formula holds. (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 13

5. Multivariate antisymmetric and symmetric finite Fourier transforms Along with the integral Fourier transform in one variable there exists a discrete Fourier transform in one variable running over a finite set. Similarly, it is possible to introduce finite multivariate antisymmetric and symmetric Fourier transforms, based on antisymmetric and symmetric exponential functions. We first consider the finite Fourier transform in one variable, which will be used below. Then we expose a general antisymmetric and symmetric Fourier transforms. Under exposition we use the methods developed in [16]. 5.1. Finite Fourier transform. Let us fix a positive integer N and consider the numbers −1/2 emn := N exp(2πimn/N), m,n =1, 2, ,N. (43) N ··· The matrix (emn)m,n=1 is unitary, that is,

emkenk = δmn, ekmekn = δmn. (44) Xk Xk Indeed, according to the formula for a sum of a geometric progression we have ta + ta+1 + + ta+r = (1 t)−1ta(1 tr+1), t =1, ··· − − 6 ta + ta+1 + + ta+r = r +1, t =1. ··· Setting t = exp(2πi(m n)/N), a = 1 and r = N 1, we prove (44). Let f(n) be a function− of n 1, 2 ,N . We− may consider the transform ∈{ ··· } N N f(n)e N −1/2 f(n) exp(2πimn/N)= f˜(m). (45) mn ≡ nX=1 nX=1 N Then due to unitarity of the matrix (emn)m,n=1, we express f(n) as a linear com- bination of conjugates of the functions (43): N f(n)= N −1/2 f˜(m) exp( 2πimn/N). (46) − mX=1 The function f˜(m) is a ﬁnite Fourier transform of f(n). This transform is a linear map. The formula (46) gives an inverse transform. The Plancherel formula N N f˜(m) 2 = f(n) 2 | | | | mX=1 nX=1 holds for transforms (45) and (46). This means that the ﬁnite Fourier transform conserves the norm introduced in the space of functions on 1, 2,...,N . { } 5.2. Antisymmetric multivariate discrete Fourier transforms. We use the discrete exponential function (43), e (s) := N −1/2 exp(2πims), s F 1 , 2 ,..., N−1 , 1 , m Z≥0, (47) m ∈ N ≡{ N N N } ∈ and make a multivariate discrete exponential function by taking a product of n copies of these functions,

em(s) := e (s )e (s ) e (s ) m1 1 m2 2 ··· mn n = N −n/2 exp(2πim s ) exp(2πim s ) exp(2πim s ) (48) 1 1 2 2 ··· n n 14 A. KLIMYK AND J. PATERA where s = (s ,s ,...,s ) F n and m = (m ,m ,...,m ) (Z≥0)n. Now we take 1 2 n ∈ N 1 2 n ∈ these multivariate functions for integers mi such that m1 > m2 > > mn 0 and make an antisymmetrization. As a result, we obtain a finite··· version of≥ the antisymmetric exponential functions (8): − −1/2 n −1/2 −n/2 − E˜m(s) := S det(e (s )) = S N Em(s), (49) | n| mi j i,j=1 | n| where, as before, Sn is the order of the symmetric group Sn. The n-tuples s|in| (49) runs over F n F F (n times). We denote by N ≡ N ×···× N Fˆn the subset of F n consisting of s F n such that N N ∈ N s >s > >s . 1 2 ··· n ˆn aff aff The set FN is a finite subset of the fundamental domain F (Sn ) of the group Sn . ˆn Note that acting by permutations w Sn upon FN we obtain the whole set n ∈ FN without those points which are invariant under some nontrivial permutation w S . Clearly, the function (49) vanishes on the last points. ∈ n Since the discrete exponential functions em(s) satisfy the equality em(s) = e (s), we do not need to consider them for all values m Z≥0. It is enough to m+N ∈ consider them for m 1, 2,...,N . By Dˆ n we denote the set of integer n-tuples ∈{ } N m = (m1,m2,...,mn) such that N m >m > >m > 0. ≥ 1 2 ··· n We need a scalar product in the space of linear combinations of the functions (48). It is natural to give it by the formula n n

′ ′ ′ ′ em(s),em (s) emi (si),em (si) := emi (si)em (si)= δmm (50) h i≡ h i i i iY=1 iY=1 siX∈FN where m ,m′ 1, 2,...,N . Here we used the relation (44). i i ∈{ } ′ ˆ n Proposition 1. For m, m DN the discrete functions (49) satisfy the orthogonality relation ∈

− − − − E˜m(s), E˜m′ (s) = S E˜m(s)E˜m′ (s)= δmm′ , (51) h i | n| sXˆn ∈FN where the scalar product is determined by formula (50). ′ ′ ′ Proof. Since m1 > m2 > > mn > 0 and m1 > m2 > > mn > 0, then due to the deﬁnition of the scalar··· product we have ···

− − − − E˜m(s), E˜m′ (s) = E˜m(s)E˜m′ (s) h i sXn ∈FN n −1 = Sn em (si)em′ (si)= δmm′ , (52) | | w(i) w(i) wX∈Sn iY=1 siX∈FN where (mw(1),mw(2),...,mw(n)) is obtained from (m1,m2,...,mn) by action by − the permutation w S . Since functions E˜m(s) are antisymmetric with respect to ∈ n Sn, then − − − − E˜m(s)E˜m′ (s)= S E˜m(s)E˜m′ (s). | n| sX∈F n sXˆn N ∈FN This proves the proposition. (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 15

ˆn n Let f be a function on FN (or an antisymmetric function on FN ). Then it can be expanded in the functions (49) as

− f(s)= amE˜m(s). (53) mXˆ n ∈DN

The coeﬃcients am are determined by the formula

− am = S f(s)E˜m(s). (54) | n| mXˆn ∈FN ˆ n We have taken into account the facts that numbers of elements in DN and in ˆn FN are the same and that the discrete functions (49) are orthogonal with respect to the scalar product (51). We call expansions (53) and (54) the antisymmetric multivariate discrete Fourier transforms. These expansions can be written in terms − n of the exponential function Em(s) = det(exp(2πimisj ))i,j=1,

−n/2 − −n/2 − f(s)= N amEm(s), am = N S f(s)Em(s). (55) | n| mXˆ n mXˆn ∈DN ∈FN

5.3. Symmetric multivariate discrete Fourier transforms. Let us give a symmetric multivariate discrete Fourier transforms. For this we take the multivariate exponential functions (48) for integers mi such that N m m m 1 ≥ 1 ≥ 2 ≥···≥ n ≥ and make a symmetrization. We obtain a ﬁnite version of the symmetric exponential functions (3),

+ −1/2 + n −1/2 −1/2 + E˜m(s) := S det (e (s )) = S N Em(s), (56) | n| mi j i,j=1 | n| where the discrete functions em(s) are given by (47). The n-tuples s in (56) run over F n F F (n times). We denote by N ≡ N ×···× N F˘n the subset of F n consisting of s = (s ,s ,...,s ) F n such that N N 1 2 n ∈ N s s s . 1 ≥ 2 ≥···≥ n ˘n aﬀ The set FN is a ﬁnite subset of the closure of the fundamental domain F (Sn ). ˘n Note that acting by permutations w Sn upon FN we obtain the whole set n ∈ FN , where each point, having some coordinates mi coinciding, are repeated several times. Namely, a point s is contained Ss times in wF˘n ; w S , where Ss is the | | { N ∈ n} subgroup of Sn consisting of elements w Sn such that ws = s. ˘ n ∈ By DN we denote the set of integer n-tuples m = (m1,m2,...,mn) such that N m m m 1. ≥ 1 ≥ 2 ≥···≥ n ≥ ′ ˘ n Proposition 2. For m, m DN the discrete functions (56) satisfy the orthogonality relation ∈

+ + −1 + + E˜m(s), E˜m′ (s) = S Ss E˜m(s)E˜m′ (s)= Sm δmm′ . (57) h i | n| | | | | sX˘n ∈FM Proof. This proposition is proved in the same way as Proposition 1, but we ˘n ˆn have to take into account the diﬀerence between FM and FM . Due to the deﬁnition 16 A. KLIMYK AND J. PATERA of the scalar product we have

+ + + + E˜m(s), E˜m′ (s) = E˜m(s)E˜m′ (s) h i sXn ∈FN n −1 = Sn Sm em (si)em′ (si)= Sm δmm′ , | | | | w(i) w(i) | | wX∈Sn iY=1 siX∈FN Here we have taken into account that there appear additional summands (with respect to (52)) because some summands on the right hand side of (56) can coincide. + Since functions E˜m(s) are symmetric with respect to Sn, then

+ + −1 + + E˜m(s)E˜m′ (s)= S Ss E˜m(s)E˜m′ (s), | n| | | sX∈F n sX˘n N ∈FN ˘n where we have taken into account that under an action by Sn upon FN a point s appears Ss times in F n . This proves the proposition. | | N ˘n n Let f be a function on FN (or a symmetric function on FN ). Then it can be expanded in functions (56) as

+ f(s)= amE˜m(s). (58) mX˘ n ∈DN

The coeﬃcients am are determined by the formula

−1 −1 + am = S Sm Ss f(s)E˜m(s). (59) | n|| | | | sX˘n ∈FN The expansions (58) and (59) follows from the facts that numbers of elements ˘ n ˘n in DN and in FN are the same and from the orthogonality relation (57). We call expansions (58) and (59) the symmetric multivariate discrete Fourier transforms.

6. Eigenfunctions of (anti)symmetric Fourier transforms

Let Hn(x), n = 0, 1, 2,... , be the well-known Hermite polynomials of one variable. They satisfy the relation ∞ 2πipx −πp2 −m −πx2 e e Hm(√2πp)dp = i e Hm(√2πx) (60) Z−∞ (see, for example, subsection 12.2.4 in [17]). We create polynomials of many variables

Hm(x) H (x , x ,...,x ) := H (x )H (x ) H (x ). (61) ≡ m1,m2,...,mn 1 2 n m1 1 m2 2 ··· mn n The functions −|x|2/2 e Hm(x), mi =0, 1, 2,..., i =1, 2, . . . , n, (62) where x is the length of the vector x, form an orthogonal basis of the Hilbert | | space L2(Rn) with the scalar product f ,f := f (x)f (x)dx, where dx = h 1 2i Rn 1 2 dx1 dx2 dxn. R We make··· symmetrization and antisymmetrization of the functions

−π|x|2 m(x) := e Hm(√2πx) H (ANTI)SYMMETRIC MULTIVARIATE EXPONENTIAL FUNCTIONS 17

(obtained from (62) by replacing x by √2πx) by means of symmetric and antisymmetric multivariate exponential functions:

+ −π|x|2 √ −|m| −π|λ|2 sym √ Eλ (x)e Hm( 2πx) = i e Hm ( 2πλ), (63) ZRn

− −π|x|2 √ −|m| −π|λ|2 anti √ Eλ (x)e Hm( 2πx) = i e Hm ( 2πλ). (64) ZRn sym anti It is easy to see that the polynomials Hm and Hm indeed are symmetric and antisymmetric, respectively, with respect to the group Sn, sym sym anti anti Hm (wλ)= Hm (λ), Hm (wλ) = (det w)Hm (λ), w S . ∈ n sym For this reason, we may consider Hm (λ) for values of λ = (λ1, λ2,...,λn) such anti that λ1 λ2 λn and Hm (λ) for values of λ such that λ1 > λ2 > > λn. ≥ ≥···≥sym ··· The polynomials Hm are of the form sym + n Hm (λ) = det (Hmi (λj ))i,j=1 (65) anti and the polynomials Hm of the form anti n Hm (λ) = det (Hmi (λj ))i,j=1 . (66) anti Moreover, Hm (λ) = 0 if mi = mi+1 for some i =1, 2,...,n 1. For this reason, sym − we may consider the polynomials Hm (λ) for integer n-tuples m such that m1 anti ≥ m m and the polynomials Hm (λ) for integer n-tuples m such that 2 ≥ ··· ≥ n m1 >m2 > >mn. Let us apply··· the symmetric Fourier transform (42) (we denote it as F) to the symmetric functions (65). Taking into account formula (63) we obtain

2 1 2 −π|x| sym √ + −π|x| sym √ F e Hm ( 2πx) = Eλ (x)e Hm ( 2πx)dx Sn ZRn | | −|m| −π|λ|2 sym = i e Hm (√2πλ),

−π|x|2 sym that is, functions e Hm (√2πx) are eigenfunctions of the symmetric Fourier transform F. Since these functions for mi = 0, 1, 2,... , i = 1, 2,...,n, m1 2 Rn ≥ m2 mn, form an orthogonal basis of the Hilbert space L0( ) of functions ≥···≥2 n 2 from L (R ) symmetric with respect to Sn (that is, of the Hilbert space L (D+)), then they constitute a complete set of eigenfunctions of this transform. Thus, this 2 Rn transform has only four eigenvalues i, i, 1, 1 in L0( ). This means that we have F4 = 1. − − Now we apply the antisymmetric Fourier transform (40) (we denote it as F˜) to −π|x|2 anti the antisymmetric function e Hm (√2πx). Taking into account formula (64) we obtain

2 1 2 ˜ −π|x| anti √ − −π|x| anti √ F e Hm ( 2πx) = Eλ (x)e Hm ( 2πx)dx Sn ZRn | | −|m| −π|λ|2 anti = i e Hm (√2πλ),

−π|x|2 anti that is, functions e Hm (√2πx) are eigenfunctions of the transform F˜. Since these functions for mi =0, 1, 2,... ; i =1, 2,...,n, m1 >m2 > >mn 0, form 2 Rn ··· 2 ≥Rn an orthogonal basis of the Hilbert space L−( ) of functions from L ( ) antisymmetric with respect to W , then they constitute a complete set of eigenfunctions 18 A. KLIMYK AND J. PATERA of this transform. Thus, this transform has only four eigenvalues i, i, 1, 1. This − − means that, as in the previous case, we have F˜4 = 1. Acknowledgements. We are grateful for the hospitality extended to A.K. at the Center de Recherches Mathématiques, Universitéde Montréal, during the prepa- ration of this paper. We acknowledge also partial support for this work from the National Science and Engineering Research Council of Canada, MITACS, the MIND Institute of Costa Mesa, California, and Lockheed Martin, Canada.

References [1] Patera J., C-functions of compact semisimple Lie groups as special functions, Proc. Inst. Math. Nat. Acad. Sci. of Ukraine, 30 (2004), 1152-1160. [2] Patera J., Compact simple Lie groups and their C-, S-, and E-transforms, Symmetry, integrability and Geometry: Methods and Applications, 1 (2005), paper 025, 6 pages. [3] Klimyk A. U. and Patera J., Orbit functions, Symmetry, integrability and Geometry: Methods and Applications, 2 (2006), paper 06, 60 pages. [4] Klimyk A. U. and Patera J., Antisymmetric orbit functions, Symmetry, integrability and Geometry: Methods and Applications, 3 (2007), paper 023, 83 pages. [5] Atoyan A. and Patera J., Properties of continuous Fourier extension of the discrete cosine transform and its multidimensional generalization, J. Math. Phys. 45 (2004), 2468–2491. [6] Patera J. and Zaratsyan A., Discrete and continuous cosine transform generalized to Lie groups SU(2) × SU(2) and O(5), J. Math. Phys. 46 (2005), 053514, 17 pages. [7] Patera J. and Zaratsyan A., Discrete and continuous cosine transform generalized to Lie groups SU(2) and G2, J. Math. Phys. 46 (2005), 113506, 25 pages. [8] Kashuba I. and Patera J., Discrete and continuous exponential transforms of simple Lie groups of rank two, J. Phys. A: Math. Theor. 40 (2007), 1751–1774. [9] Atoyan A. and Patera J., Continuous extension of the discrete cosine transform, and its applications to data processing, in Group Theory and Numerical Analysis, CRM Proc. Lecture Notes, Vol. 39, Amer. Math. Soc., Providence, RI, 2005, 1–15. [10] Atoyan A., Patera J., Sahakian V., and Akhperjanian A., Fourier transform method for imaging atmospheric Cherenkov telescopes, Astroparticle Phys. 23 (2005), 79–95. [11] Patera J., Zaratsyan A., and Zhu H.-M., New class of interpolation methods based on dis- cretized Lie group transform, SPIE Electronic Imaging, 6064A-06 (2006), S1. [12] Germain M., Patera J., and Zaratsyan A., Multiresolution analysis of digital images using the continuous extension of discrete group transform, SPIE Electronic Imaging, 6064A-03 (2006), S2. [13] Germain M., Patera J. and Allard Y., Cosine transform generalized to Lie groups SU(2) × SU(2), O(5), and SU(2) × SU(2) × SU(2): application to digital image processing, Proc. SPIE 6065 (2006), 387–395. [14] Karlin S. and McGregor J., Determinants of orthogonal polynomials, Bull. Amer. Math. Soc. 68 (1962), 204-209. [15] Berens H., Schmid H., and Xu Y., Multivariate Gaussian cubature formulas, Arch. Math., 64 (1995), 26-32. [16] Moody R. V. and Patera J., Orthogonality within the families of C-, S-, and E-functions of any compact semisimple Lie group, Symmetry, integrability and Geometry: Methods and Applications, 2 (2006), paper 76, 14 pages. [17] Vilenkin N. Ja. and Klimyk A. U., Representations of Lie Groups and Special Functions, Vol. 2, Kluwer, Dordrecht, 1993.

Bogolyubov Institute for Theoretical Physics, Kiev 03143, Ukraine E-mail address: [email protected]

Centre de Recherches Mathematiques,´ Universite´ de Montreal,´ C.P.6128-Centre ville, Montreal,´ H3C3J7, Quebec,´ Canada E-mail address: [email protected]