Research supported by the Swedish Science Research Council (VR) .
<ph f="cmbx">Spherical transform and Jacobi polynomials on root systems of type BC</ph>

### Genkai Zhang

Department of Mathematics, Chalmers University of Technology and Goteborg University, S-412 96 Goteborg, Sweden E-mail address : genkai@math.chalmers.se
• Abstract. Let $R$  be a root system of type BC in $\mathfrak{a}={\mathbb{R}}^{r}$  of general positive multiplicity. We introduce certain canonical weight function on ${\mathbb{R}}^{r}$  which in the case of symmetric domains corresponds to the integral kernel of the Berezin transform. We compute its spherical transform and prove certain Bernstein-Sato type formula. This generalizes earlier work of the author. Associated to the weight functions there are Heckman-Opdam orthogonal polynomials of Jacobi type on the compact torus, after a change of variables they form an orthogonal system on the non-compact space $\mathfrak{a}$  . We consider their spherical transform and prove that they are the Macdonald-Koornwinder polynomials multiplied by the spherical transform of the canonical weight function.
For rank one case this was proved earlier by Koornwinder.

Introduction

The Gaussian functions and the Hermite polynomials play an important role, in Fourier transform on Euclidean spaces; the Hermite polynomials diagonalize the harmonic oscillator and the corresponding Hermite type function diagonalize the Fourier transform, which make the Plancherel theory more transparent. The generalization of the Plancherel theory to any non-compact symmetric symmetric spaces has been studied intensively and in general there are still no general theory generalizing corresponding the results of Hermite polynomials and their Fourier transform, above all, no concrete orthogonal systems of functions with explicit spherical transforms are constructed. In the present paper we will establish such a theory for root system of type BC. Now associated with any root system in ${\mathbb{R}}^{r}$  there are two kinds of remarkable orthogonal polynomials, namely the Heckman-Opdam orthogonal polynomials giving the spectral decomposition of algebra of the Weyl group invariant polynomials of the Cherednik operators acting on certain ${L}^{2}$  -space on the compact space ${\mathbb{T}}^{r}$  , and the Macdonald-Koornwinder polynomials orthogonal with respect to certain weight functions on ${\mathbb{R}}^{r}$  defined as a product of Gamma functions. Part of the product is in fact the Harish-Chandra Plancherel measure for the spectral decomposition of the algebra of the Weyl group invariant polynomials of the Cherednik operators on the non-compact space ${\mathbb{R}}^{r}$  . In the present paper we will introduce certain canonical weight function ${f}_{-2\nu }\left(x\right)$  for root system of type BC, and we compute its spherical transform $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)$  . The product of $\stackrel{~}{{f}_{-2\nu }}\left(\lambda {\right)}^{2}$  with the the Harish-Chandra measure $|c\left(\lambda \right){|}^{-2}$  gives precisely weight functions in the Macdonald-Koornwinder orthogonality relation; we prove further, roughly speaking, that the Jacobi polynomials on the compact space, after some change of variables, and multiplied by the function ${f}_{-2\nu }$  gives orthogonal system of the ${L}^{2}$  -space on the non-compact space, and that their spherical transform are of the form $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)p\left(\lambda \right)$  where $p\left(\lambda \right)$  are the Macdonald-Koornwinder polynomials.
Thus the function ${f}_{-2\nu }$  plays the role of the Gaussian functions in the Fourier analysis where as the Jacobi and Macdonald-Koornwinder polynomials plays the role of the Hermite polynomials and their Fourier transforms (up to a multiple of the Gaussian). In rank one case this has been proved earlier by Koornwinder [5.
We remark that we are studying here only the spherical transform of Weyl group invariant functions, yet we use the tools of the Cherednik operators acting on general no-invariant functions which makes the computation much easier. Some of our results can easily be generated to that setup.
Our results are motivated by the study of Berezin transform on bounded symmetric domains $G/K$  ([13and [14). Berezin transform appears naturally in the study of branching rule of holomorphic representations and some results here have been proved in the symmetric space case. The Berezin transform is convolution operator with the kernel $h\left(z,z{\right)}^{\nu }$  where $h\left(z,w{\right)}^{-\nu }$  is the reproducing kernel for a weighted Bergman space. In the symmetric space case our function ${f}_{-2\nu }$  is just the kernel $h\left(z,z{\right)}^{\nu }$  written in terms of geodesic coordinates. The Berezin transform is invariant under $G$  and is bounded on the ${L}^{2}$  -space on $G/K$  and is thus a function of a system of generators of invariant differential operators on $G/K$  , also called the symbol of the Berezin transform. The symbol is thus the spherical transform of the function $h\left(z,z{\right)}^{\nu }$  . It is realized in [12the orthogonality relation for the continuous dual Hahn ${}_{3}{F}_{2}$  are Clebsch-Gordan coefficients for the tensor product of holomorphic and anti-holomorphic representations and that the orthogonality measure is actually the product of the Harish-Chandra measure $|c\left(\lambda \right){|}^{-2}$  with the spectral symbol of Berezin transform, which has been computed earlier by Unterberger and Upmeier [8. The Macdonald-Koornwinder polynomials are the multi-variable generalization of the Wilson's ${}_{4}{F}_{3}$  hypergeometric orthogonal polynomials, and the orthogonal relation for the above ${}_{3}{F}_{2}$  orthogonal polynomials can be viewed as the degenerate case of the ${}_{4}{F}_{3}$  -orthogonality relation [9. In this sense the orthogonality relation of the Macdonald-Koornwinder polynomials can be viewed as the most general case, our result in this paper concludes that the relation can be obtained from the more well-understood orthogonality of the Jacobi type polynomials and the Plancherel formula, by introducing the weight function.
The paper is organized as follows. In Section 1 we recall the Plancherel formula for Opdam-Cherednik transform. In Section 2 we find certain Bernstein-Sato type formula for the so-called canonical function ${f}_{\delta }$  . Its spherical transform is computed in Section 3. Finally we compute the spherical transform of certain Jacobi-type polynomials in Sections 4 and 5.
It is our belief that most results in this paper can be generalized to general root systems, that will provide a theory of Hermite-type functions on non-compact space and thus bring together the orthogonal polynomials on compact torus and spherical transform on non-compact space.
I would like to thank Professors Jacques Faraut, Toshiyuki Kobayashi, Erik Opdam and Harald Upmeier for some helpful discussions.

1 Spherical transform and Plancherel formula

Let $\mathfrak{a}={\mathbb{R}}^{r}$  be an Euclidean space with inner product $\left(\cdot ,\cdot \right)$  and let $R\subset {\mathfrak{a}}^{*}$  be a root system of type BC. A prototype of such root systems is the restricted root system of a bounded symmetric domain. We use some familiar notation, in order to be consistent with the notation in the context of bounded symmetric domain and also in the context of Plancherel formula for general root system [7. We fix an orthonormal basis $\left\{{\xi }_{j}{\right\}}_{1}^{r}$  of $\mathfrak{a}$  and a dual basis $\left\{{\varepsilon }_{j}\right\}$  of ${\mathfrak{a}}^{*}$  , i.e., ${\varepsilon }_{j}\left({\xi }_{k}\right)={\delta }_{jk}$  , so that the root system is of the form $R=\left\{2{\varepsilon }_{j};j=1,\cdots ,r\right\}\cup \left\{4{\varepsilon }_{j};j=1,\cdots ,r\right\}\cup \left\{2\left({\varepsilon }_{j}±{\varepsilon }_{k}\right);j\ne k=1,\cdots ,r\right\}$  , with respective multiplicities $\left({k}_{1},{k}_{2},{k}_{3}\right)$  satisfying $2\left({k}_{1},{k}_{2},{k}_{3}\right)=\left(2b,\iota ,a\right).$  We assume that $\iota ,a,b>0$  . We order the roots so that ${\varepsilon }_{1}>{\varepsilon }_{2}>\cdots >{\varepsilon }_{r}>0$  .
The Weyl group is then $W={S}_{r}×{\mathbb{Z}}_{2}^{r}$  consisting of signed permutation o $\left\{{\xi }_{j}\right\}$  . Let $\rho =\frac{1}{2}{\sum }_{\alpha \in {R}^{+}}{k}_{\alpha }\alpha$  be the half sum of positive roots. Then $\rho ={\sum }_{j=1}^{r}{\rho }_{j}{\varepsilon }_{j}={\sum }_{j=1}^{r}\left(\iota +b+\left(r-j\right)a\right){\varepsilon }_{j}.$  We recall briefly in this section the Plancherel formula for the spherical transform associated to the root system $R$  developed by Heckman and Opdam ([4, [3and [7).
Let $d\mu$  be the measure $d\mu \left(t\right)=d{\mu }_{k}\left(t\right)={\prod }_{\alpha \in {R}_{+}}|2sinh\left(\frac{1}{2}\alpha \left(t\right){|}^{2{k}_{\alpha }}dt$  on $\mathfrak{a}$  , and let ${L}^{2}\left(\mathfrak{a}\right)={L}^{2}\left(\mathfrak{a},d\mu \right)$  and ${L}^{2}\left(\mathfrak{a}{\right)}^{W}$  be the corresponding ${L}^{2}$  space and its subspace of $W$  -invariant functions.
Let
 $\begin{array}{cccc}{D}_{j}& ={\partial }_{j}-a{\sum }_{i
be the Cherednik operators acting on functions $f\left(t\right)$  on ${\mathfrak{a}}^{\mathbb{C}}$  , where we identify a function $f\left(t\right)$  on ${\mathfrak{a}}^{\mathbb{C}}$  with $f\left({t}_{1},\cdots ,{t}_{r}\right)$  , for $t={t}_{1}{\xi }_{1}+\cdots +{t}_{r}{\xi }_{r}$  . Here ${s}_{ij},{\sigma }_{ij},{\sigma }_{i}$  are the elements in the Weyl group, ${s}_{ij}=\left(ij\right)$  being the permutation of ${\xi }_{i}$  and ${\xi }_{j}$  , ${\sigma }_{ij}$  the signed permutation, ${\sigma }_{ij}\left({\xi }_{i}\right)=-{\xi }_{j},{\sigma }_{ij}\left({\xi }_{j}\right)=-{\xi }_{i}$  , ${\sigma }_{i}$  the reflection ${\sigma }_{i}\left({\xi }_{i}\right)=-{\xi }_{i}$  , and all mapping ${\xi }_{k}\to {\xi }_{k}$  for $k\ne i,j$  . For later purpose we will also rewrite it as
 $\begin{array}{cccccc}{D}_{j}& ={\partial }_{j}-a{\sum }_{i
The operators $\left\{{D}_{j}\right\}$  are then commuting, and the decomposition of ${L}^{2}\left(\mathfrak{a}\right)$  with respect to the eigenfunctions of $\left\{{D}_{j}\right\}$  is given by Opdam [7in terms of the so-called Cherednik-Opdam, formulated in terms of decomposing representations of a Hecke algebra. We will mostly be concerned with the decomposition of ${L}^{2}\left(\mathfrak{a}{\right)}^{W}$  under $W$  -invariant polynomials of ${D}_{j},j=1,\cdots ,r$  , which is given by spherical transform in terms of the Heckman-Opdam theory of hypergeometric functions ([4and [3); we shall however use the operators ${D}_{j}$  (in place of their symmetric polynomials) to compute the spherical transform of symmetric functions.
For $\lambda \in \left({\mathfrak{a}}^{*}{\right)}^{\mathbb{C}}$  let ${\phi }_{\lambda }$  be as in [7the spherical function. In particular
 $\begin{array}{c}p\left({D}_{1},\cdots ,{D}_{r}\right){\phi }_{\lambda }=p\left(\lambda \left({\xi }_{1}\right),\cdots ,\lambda \left({\xi }_{r}\right)\right){\phi }_{\lambda },\end{array}$ (1.1)
for any $W$  -invariant polynomial $p$  . The spherical transform of a function $f\in {L}^{2}\left(\mathfrak{a}{\right)}^{W}$  is $\stackrel{~}{f}\left(\lambda \right)={\int }_{\mathfrak{a}}f\left(t\right){\phi }_{\lambda }\left(t\right)d\mu \left(t\right).$  The corresponding Plancherel measure is given by $d\stackrel{~}{\mu }\left(\lambda \right)=\frac{\left(2\pi {\right)}^{-r}{c}_{0}^{2}}{c\left(\lambda \right)c\left(-\lambda \right)}d\lambda ,$  where $c\left(\lambda \right){=}^{r}{\prod }_{j=1}{\frac{\Gamma \left({\lambda }_{j}+b\right)\Gamma \left(2{\lambda }_{j}\right)}{\Gamma \left({\lambda }_{j}+b+\frac{\iota }{2}\right)\Gamma \left(2{\lambda }_{j}+2b\right)}}^{r}{\prod }_{1\le j  and ${c}_{0}=c\left(\rho \right){=}^{r}{\prod }_{j=1}{\frac{\Gamma \left({\rho }_{j}+b+1\right)\Gamma \left(2{\rho }_{j}+1\right)}{\Gamma \left({\rho }_{j}+b+\frac{\iota }{2}+1\right)\Gamma \left(2{\rho }_{j}+b+1\right)}}^{r}{\prod }_{1\le j  Namely, we have
 $\begin{array}{c}{\int }_{\mathfrak{a}}|f\left(t\right){|}^{2}d\mu \left(t\right)={\int }_{i{\mathfrak{a}}^{*}}|\stackrel{~}{f}\left(\lambda \right){|}^{2}d\stackrel{~}{\mu }\left(\lambda \right).\end{array}$ (1.2)

2 Bernstein-Sato type formula for the function ${f}_{\delta }$

We define the weight function ${f}_{\delta }$  on $\mathfrak{a}$  by ${f}_{\delta }\left(t\right){=}^{r}{\prod }_{j=1}{cosh}^{\delta }{t}_{j}.$  Motivated by the Berezin transform ([8, [10, [13) we call ${f}_{\delta }$  the canonical function. In the case of bounded symmetric domains ${f}_{\delta }$  is the integral kernel of the Berezin transform considered as a convolution operator defining the so-called canonical representations, and is the analogue of the Gaussian functions in the Euclidean space, see loc. cit..
In this section we prove the following Bernstein-Sato type formula:
Theorem 2.1. There is a $W$  -invariant polynomials of the Cherednik operators ${D}_{j}$  mapping the canonical function ${f}_{\delta }$  to ${f}_{\delta -2}$  ; more precisely, we have
 $\begin{array}{c}{}^{r}{\prod }_{j=1}\left({D}_{j}^{2}-\left(\delta +\rho \left({\xi }_{1}\right){\right)}^{2}\right){f}_{\delta }{=}^{r}{\prod }_{j=1}\left(\delta +a\left(j-1\right)\right)\left(1-\delta -\iota -a\left(r-j\right)\right){f}_{\delta -2}\end{array}$ (2.1)
The proof of it will be divided into two technical lemmas.
Lemma 2.2. The following formula holds ${}^{j}{\prod }_{l=1}\left({D}_{l}+\delta +\rho \left({\xi }_{1}\right)\right){f}_{\delta }\left(t\right){=}^{j}{\prod }_{l=1}\left(\delta +a\left(l-1\right)\right){f}_{\delta }\left(t{\right)}^{j}{\prod }_{l=1}\left(1+tanh{t}_{l}\right)$
• Proof. We prove the lemma by induction on $j$  . First we have, since ${f}_{\delta }$  is Weyl group invariant, ${D}_{1}{f}_{\delta }=\delta tanh{t}_{1}{f}_{\delta }-\rho \left({\xi }_{1}\right){f}_{\delta },$  where we use $\frac{d}{dt}cosht=tanhtcosht$  . Rewriting, $\left({D}_{1}+\delta +\rho \left({\xi }_{1}\right)\right){f}_{\delta }=\delta \left(1+tanh{t}_{1}\right){f}_{\delta },$  which is the claim for $j=1$  . Assume that the equality is true for ${\prod }_{l=1}^{j-1}\left({D}_{l}+\left(\delta +\rho \left({\xi }_{1}\right)\right)\right)$  .
We consider it for $j$  in place of $j-1$  . We need to compute the operator ${D}_{j}+\delta +\rho \left({\xi }_{1}\right)$  on the function ${f}_{\delta }{\prod }_{l=1}^{j-1}\left(1+tanh{t}_{l}\right)$  . This function is invariant under the permutations of the first $j-1$  coordinates and the signed permutations of the last $r-j$  coordinates.
Thus  $\begin{array}{cccccccccc}& {D}_{j}{f}_{\delta }{\prod }_{l\le j-1}\left(1+tanh{t}_{l}\right)& & =\delta tanh{t}_{j}{f}_{\delta }{\prod }_{l\le j-1}\left(1+tanh{t}_{l}\right)& & -a{f}_{\delta }{\sum }_{i
Using the formulas $tanhx±tanhy=\frac{sinh\left(x±y\right)}{coshxcoshy}$  and ${e}^{x}=\left(1+tanhx\right)coshx$  we see that the sum of the $i$  th terms in the two summations is, apart form the factor $a{f}_{\delta }{\prod }_{l\le j-1,l\ne i}\left(1+tanh{t}_{l}\right)$  ,  $\begin{array}{cccc}& \frac{{e}^{{t}_{i}+{t}_{j}}}{{e}^{{t}_{i}+{t}_{j}}-{e}^{-\left({t}_{i}+{t}_{j}\right)}}\left(tanh{t}_{i}+tanh{t}_{j}\right)-\frac{{e}^{{t}_{i}-{t}_{j}}}{{e}^{{t}_{i}-{t}_{j}}-{e}^{-\left({t}_{i}-{t}_{j}\right)}}\left(tanh{t}_{i}-tanh{t}_{j}\right)& & =2\left(1+tanh{t}_{i}\right)tanh{t}_{j},\end{array}$
and thus  $\begin{array}{cccc}& {D}_{j}{f}_{\delta }{\prod }_{l\le j-1}\left(1+tanh{t}_{l}\right)& & =\left(\left(\delta +a\left(j-1\right)\right)tanh{t}_{j}-\rho \left({\xi }_{j}\right)\right){f}_{\delta }{\prod }_{l\le j-1}\left(1+tanh{t}_{l}\right).\end{array}$
It follows then that  $\begin{array}{cccc}& \left({D}_{j}+\delta +\rho \left({\xi }_{1}\right)\right){f}_{\delta }{\prod }_{l\le j-1}\left(1+tanh{t}_{l}\right)& & =\left(\delta +a\left(j-1\right)\right){f}_{\delta }{\prod }_{l\le j}\left(1+tanh{t}_{l}\right)+\left(\rho \left({\xi }_{1}\right)-\rho \left({\xi }_{j}\right)-a\left(j-1\right)\right){f}_{\delta }{\prod }_{l\le j-1}\left(1+tanh{t}_{l}\right).\end{array}$
But $\rho \left({\xi }_{j}\right)=\iota +b+a\left(r-j\right)$  so the second term vanishes and this completes the proof.
Lemma 2.3. The following formula holds
 $\begin{array}{cccc}& {\prod }_{l\ge j}\left({D}_{l}-\left(\delta +\rho \left({\xi }_{1}\right)\right)\right){f}_{\delta -1}{e}^{{t}_{1}+\cdots +{t}_{r}}& & ={\prod }_{l\ge j}\left(1-\delta -\iota -a\left(r-l\right)\right){f}_{\delta -1}{e}^{{t}_{1}+\cdots +{t}_{r}}{\prod }_{l\ge j}\left(1-tanh{t}_{l}\right)\end{array}$
Accepting temporarily the Lemma, we prove Theorem  2.1 . We write the LHS as ${}^{r}{\prod }_{l=1}\left({D}_{l}-\left(\delta +\rho \left({\xi }_{1}\right)\right){\right)}^{r}{\prod }_{l=1}\left({D}_{l}+\left(\delta +\rho \left({\xi }_{1}\right)\right)\right){f}_{\delta }$  Taking $j=r$  in Lemma  2.2 we see that ${\prod }_{l=1}^{r}\left({D}_{l}+\left(\delta +\rho \left({\xi }_{1}\right)\right)\right){f}_{\delta }$  is, disregarding the constant, given by, ${{f}_{\delta }}^{r}{\prod }_{l=1}\left(1+tanh{t}_{j}\right)={f}_{\delta -1}{e}^{{t}_{1}+\cdots +{t}_{r}}.$  Theorem  2.1 then follows from Lemma  2.3 for $j=1$  , by using the identity ${f}_{\delta -1}{{e}^{{t}_{1}+\cdots +{t}_{r}}}^{r}{\prod }_{j=1}\left(1-tanh{t}_{j}\right)={f}_{\delta -2}.$  It remains to proof the Lemma.
• Proof. We first compute ${D}_{r}$  and use induction backward. We have  $\begin{array}{cccccccccc}& {D}_{r}{f}_{\delta -1}{e}^{{t}_{1}+\cdots +{t}_{r}}& & ={f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}+\left(\delta -1\right){f}_{\delta -1}tanh{t}_{r}{\prod }_{l\le r}{e}^{{t}_{l}}& & +a{f}_{\delta -1}{\sum }_{i
The third term can be simplified as $a\left(r-1\right){f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}},$  and the next two terms are $\iota tanh{t}_{r}{f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}+\left(\iota +2b\right){f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}.$  So that the previous formula is then  $\begin{array}{cccc}& \left(\delta -1+\iota \right)tanh{t}_{r}{f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}+\left(1+\iota +2b+a\left(r-1\right)-\rho \left({\xi }_{r}\right)\right){f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}& & =\left(\delta -1+\iota \right)tanh{t}_{r}{f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}+\left(1+b+\rho \left({\xi }_{1}\right)-\rho \left({\xi }_{r}\right)\right){f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}\end{array}$
and consequently  $\begin{array}{cccccc}& \left({D}_{r}-\left({\rho }_{1}+\delta \left({\xi }_{1}\right)\right)\right){f}_{\delta -1}{e}^{{t}_{1}+\cdots +{t}_{r}}& & =\left(\delta -1+\iota \right)tanh{t}_{r}{f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}+\left(1+b-\delta -\rho \left({\xi }_{r}\right)\right){f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}}& & =\left(\delta -1+\iota \right)\left(tanh{t}_{r}-1\right){f}_{\delta -1}{\prod }_{l\le r}{e}^{{t}_{l}},\end{array}$
since $1+b-\delta -\rho \left({\xi }_{r}\right)=-\left(\delta -1+\iota \right)$  . This is the Lemma for $j=r$  .
Assume the lemma is true for the action of ${\prod }_{l\ge j+1}\left({D}_{l}-\left(\delta +\rho \left({\xi }_{1}\right)\right)\right)$  . We compute $\left({D}_{j}-\left(\delta +\rho \left({\xi }_{1}\right)\right)\right){f}_{\delta -1}{\prod }_{l\ge j+1}\left(1-tanh{t}_{l}\right){\prod }_{l\le r}{e}^{{t}_{l}}$  and find it is $-\left(1-tanh{t}_{j}\right)\left(a\left(r-j\right)+\iota +\delta -1\right){f}_{\delta -1}{\prod }_{l\ge j}\left(1-tanh{t}_{l}\right){\prod }_{l\le r}{e}^{{t}_{l}}$  by a straightforward yet tedious computation.
Remark 2.4. Consider the double Hecke algebra generated by the Cherednik operators and the group ring $\mathbb{C}\left[W\right]$  and by $\mathbb{C}\left[P\right]$  of polynomials ${e}^{{\varepsilon }_{j}}$  as multiplication operators on $\mathbb{C}\left[P\right]$  . Some refinements of the above computation then gives certain commutation relations of the Cherednik operators with ${e}^{{\varepsilon }_{1}+\cdots +{\varepsilon }_{r}}$  , which might be of independent interests. This Theorem (or its equivalent version under the spherical transform, see the proof of Theorem 3.2) can possibly be obtained also by using intertwining relations of the Cherednik-Opdam transform with the double Hecke algebras proved by Cherednik [1, it would involves nevertheless many intriguing computations.

3 Spherical transform of the function ${f}_{\delta }$

In this section we shall use Bernstein-Sato type formula for the function ${f}_{-2\nu }$  to derive a recursive formula for its spherical transform $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)$  (for sufficiently large $\nu$  ).
By using a limit formula we then derive a product and, consequently, a Gamma function formula for ${\stackrel{~}{f}}_{-2\nu }\left(\lambda \right)$  .
Denote ${\Gamma }_{a}\left(\sigma \right){=}^{r}{\prod }_{j=1}\Gamma \left(\sigma -\frac{a}{2}\left(j-1\right)\right),$  which in the symmetric domain case is the Gindikin Gamma function; it will be used to simplify certain product formulas.
We compute first some normalization constant, which corresponds to the spherical transform of ${f}_{-2\nu }$  at $\lambda =\rho$  .
Lemma 3.1. Let $\nu >\iota +b+a\left(r-1\right)$  . The integral ${N}_{\nu }={\int }_{\mathfrak{a}}{f}_{-2\nu }\left(t\right)d\mu \left(t\right)$  is given by
 $\begin{array}{cccc}{N}_{\nu }& ={2}^{r\left(2\iota +2b+a\left(r-1\right)}r!{\Gamma }_{a}\left(\frac{\iota +1+2b+a\left(r-1\right)}{2}\right)& & ×\frac{{\Gamma }_{a}\left(\nu -\left(\frac{a}{2}\left(r-1\right)+\iota +b\right)\right)}{{\Gamma }_{a}\left(\nu +\frac{1-\iota }{2}\right)}{\prod }_{1\le i
• Proof. This follows from the Opdam proof of the Macdonald's conjecture. By symmetry we need only to integrate over all $s=\left({s}_{1},\cdots ,{s}_{r}\right)$  with ${s}_{j}\ge 0$  . Making the change of variables ${z}_{j}={tanh}^{2}{t}_{j}$  , we find first that  $\begin{array}{cc}{N}_{\nu }& ={2}^{\left(2\iota +2b+a\left(r-1\right)\right)r}{\int }_{\left[0,1{\right]}^{r}}|{\prod }_{i
which is evaluated by the beta integral (see [6,Ex.7,Sect.10,Chapt.VII) that  $\begin{array}{cccc}& {\int }_{\left[0,1{\right]}^{r}}|{\prod }_{i
In particular it follows that $\frac{{N}_{\nu +1}}{{N}_{\nu }}{=}^{r}{\prod }_{j=1}\frac{\nu -\left(\frac{a}{2}\left(r-1\right)+\iota +b\right)-\frac{a}{2}\left(j-1\right)}{\nu +\frac{1-\iota }{2}-\frac{a}{2}\left(j-1\right)}.$  Changing $j-1$  to $r-j$  and using $\rho \left({\xi }_{1}\right)=\iota +b+a\left(r-1\right)$  , we can rewrite it as the following
 $\begin{array}{c}\frac{{N}_{\nu +1}}{{N}_{\nu }}{=}^{r}{\prod }_{j=1}\frac{\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)-\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu +\frac{1-\iota }{2}-\frac{a}{2}\left(j-1\right)}.\end{array}$ (3.1)
Theorem 3.2. Let $\nu >\iota +b+a\left(r-1\right)$  . The spherical transform of ${f}_{-2\nu }$  is given by $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)={{N}_{\nu }}^{r}{\prod }_{j=1}\frac{\Gamma \left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)±\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)}{\Gamma \left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)±\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)\right)},\lambda \in i{\mathfrak{a}}^{*}.$
The following result is elementary and we omit its proof; it is proved in [11in the case when the root system corresponds to a bounded symmetric domain.
Lemma 3.3. Let $\phi$  be a bounded and continuous function on $\mathfrak{a}$  . Then ${lim}_{\nu \to \infty }\frac{1}{{N}_{\nu }}{\int }_{\mathfrak{a}}{f}_{-2\nu }\left(t\right)\phi \left(t\right)d\mu \left(t\right)=\phi \left(0\right)$
We prove now Theorem 3.2.
• Proof. We write ${\beta }_{\nu }\left(\lambda \right)=\frac{\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)}{{N}_{\nu }}$  and compute it in terms of ${\beta }_{\nu +1}\left(\lambda \right)$  . We note first that it is well-defined for $\lambda \in i{\mathfrak{a}}^{*}$  . Indeed the spherical function ${\phi }_{\lambda }\left(t\right)$  is a bounded function [7, and the function ${f}_{-2\nu }$  is in ${L}^{1}\left({\mathfrak{a}}^{*},d\mu \right)$  by Lemma 3.1. We perform the spherical transform on the identity ( 2.1 ) with $\delta =-2\nu$  . Using ( 1.1 ) we have, ${}^{r}{\prod }_{j=1}\left(\lambda \left({\xi }_{j}{\right)}^{2}-\left(-2\nu +\rho \left({\xi }_{1}\right){\right)}^{2}\right){\beta }_{\nu }\left(\lambda \right){=}^{r}{\prod }_{j=1}\left(-2\nu +a\left(j-1\right)\right)\left(1+2\nu -\iota -a\left(r-j\right)\right){\beta }_{\nu +1}\left(\lambda \right),$  and  $\begin{array}{cccc}{\beta }_{\nu }\left(\lambda \right)& {=}^{r}{\prod }_{j=1}\frac{\left(-2\nu +a\left(j-1\right)\right)\left(1+2\nu -\iota -a\left(r-j\right)\right)}{\left(\lambda \left({\xi }_{j}{\right)}^{2}-\left(-2\nu +\rho \left({\xi }_{1}\right){\right)}^{2}\right)}{\beta }_{\nu +1}\left(\lambda \right)& & {=}^{r}{\prod }_{j=1}\frac{\left(\nu -\frac{a}{2}\left(j-1\right)\right)\left(\nu +\frac{1-\iota }{2}-\frac{a}{2}\left(r-j\right)\right)}{\left(\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)\left(\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)}{\beta }_{\nu +1}\left(\lambda \right)\end{array}$
We write further the denominator as  $\begin{array}{cccc}& {}^{r}{\prod }_{j=1}\left(\nu -\frac{a}{2}\left(j-1\right)\right)\left(\nu +\frac{1-\iota }{2}-\frac{a}{2}\left(r-j\right)\right)& & {=}^{r}{\prod }_{j=1}\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)\left(\nu +\frac{1-\iota }{2}-\frac{a}{2}\left(j-1\right)\right).\end{array}$
Thus ${\beta }_{\nu }\left(\lambda \right)={\frac{{N}_{\nu +1}}{{N}_{\nu }}}^{r}{\prod }_{j=1}\frac{\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)\right)\left(\nu +\frac{1-\iota }{2}-\frac{a}{2}\left(j-1\right)\right)}{\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)-\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)}{\beta }_{\nu +1}\left(\lambda \right).$  Using ( 3.1 ) this becomes  $\begin{array}{cccccc}{\beta }_{\nu }\left(\lambda \right)& ={\beta }_{\nu +1}\left(\lambda {\right)}^{r}{\prod }_{j=1}\frac{\left(\nu +\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)\right)\left(\nu +\frac{1}{2}\rho \left({\xi }_{1}\right)-\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)\right)}{\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)+\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)\left(\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)-\frac{1}{2}\lambda \left({\xi }_{j}\right)\right)}& & ={\beta }_{\nu +1}\left(\lambda {\right)}^{r}{\prod }_{j=1}\left(1+\frac{\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)\left(1-\frac{\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)& & ×{\left(1+\frac{\frac{1}{2}\lambda \left({\xi }_{j}\right)}{\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)}^{-1}{\left(1-\frac{\frac{1}{2}\lambda \left({\xi }_{j}\right)}{\nu -\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)}^{-1}.\end{array}$
Iterating we find  $\begin{array}{cccc}{\beta }_{\nu }\left(\lambda \right)& ={\beta }_{\nu +k}\left(\lambda {\right)}^{k}{{\prod }_{l=1}}^{r}{\prod }_{j=1}\left(1+\frac{\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)\left(1-\frac{\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)& & ×{\left(1+\frac{\frac{1}{2}\lambda \left({\xi }_{j}\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)}^{-1}{\left(1-\frac{\frac{1}{2}\lambda \left({\xi }_{j}\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)}^{-1}\end{array}$
However $\frac{1}{{N}_{\nu +k}}{b}_{\nu +k}\left(\lambda \right)\to 1$  as $k\to \infty$  by Lemma 3.3 since the function ${\phi }_{\lambda }\left(t\right)$  is a continuous bounded function for $\lambda \in i{\mathfrak{a}}^{*}$  . We get thus  $\begin{array}{cccc}{\beta }_{\nu }\left(\lambda \right)& {=}^{\infty }{{\prod }_{l=1}}^{r}{\prod }_{j=1}\left(1+\frac{\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)\left(1-\frac{\frac{1}{2}\left(\iota +b+a\left(j-1\right)\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)& & ×{\left(1+\frac{\frac{1}{2}\lambda \left({\xi }_{j}\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)}^{-1}{\left(1-\frac{\frac{1}{2}\lambda \left({\xi }_{j}\right)}{\nu +l-1-\frac{1}{2}\rho \left({\xi }_{1}\right)}\right)}^{-1}.\end{array}$
Our result follows rewriting the infinite product in terms of the Gamma function, using $\frac{\Gamma \left(A\right)\Gamma \left(B\right)}{\Gamma \left(A+C\right)\Gamma \left(B-C\right)}{=}^{\infty }{\prod }_{s=0}\left(1+\frac{C}{A+s}\right)\left(1-\frac{C}{B+s}\right);$  see [2,p.5.
Remark 3.4. The function $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)$  for $\nu$  as above has exponential decay and is a bounded function for $\lambda \in i{\mathfrak{a}}^{*}$  , and holomorphic in a strip around $i{\mathfrak{a}}^{*}$  . We can then define the Berezin transform on ${L}^{2}\left(\mathfrak{a},d\mu {\right)}^{W}$  by the inverse spherical transform, namely ${B}_{\nu }F\left(t\right)=\int \stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)\stackrel{~}{F}\left(\lambda \right){\phi }_{\lambda }\left(t\right)|c\left(\lambda \right){|}^{-2}d\lambda .$  ${B}_{\nu }$  is then a bounded and positive operator and it has an integral kernel $B\left(t,s\right)$  . In particular $B\left(t,0\right)={f}_{-2\nu }\left(t\right)$  . It follows then from the selfadjoint property of $B$  that $B\left(x,y\right)={f}_{-2\nu }\left(t\right){f}_{-2\nu }\left(s\right)L\left(t,s\right)$  . It would be interesting to find a series expansion for the kernel $L\left(t,s\right)$  in terms of the Jack symmetric polynomials (in the variables $\left({tanh}^{2}{t}_{1},\cdots ,{tanh}^{2}{t}_{r}\right)$  ). In the case of bounded symmetric domain, this is indeed possible, some degenerate cases of the expansion has been studied in [14.

4 Spherical transform of a class of functions

We shall compute the spherical transform of a class of functions of the form ${f}_{-2\nu }\left(t\right)p\left({x}_{1}^{2},\cdots ,{x}_{r}^{2}\right),{x}_{j}=tanh{t}_{j},j=1,\cdots ,r,$  where $p\left({x}_{1}^{2},\cdots ,{x}_{r}^{2}\right)$  is a symmetric polynomial in $\left({x}_{1},\cdots ,{x}_{r}\right)$  . We will prove that they are of the form $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right)q\left(\lambda \right)$  where $q$  is a symmetric polynomial in $\left({\lambda }_{1}^{2},\cdots ,{\lambda }_{r}^{2}\right)$  .
Eventually we will specify $p$  to be Jacobi-type polynomials in the next section.
We denote ${\mathcal{D}}_{j}$  the conjugation of ${D}_{j}$  by the canonical function ${f}_{\delta }$  , ${\mathcal{D}}_{j}={\mathcal{D}}_{j}^{\left(\delta \right)}={f}_{-\delta }{D}_{j}{f}_{\delta },j=1,\cdots ,r.$  The operator ${\mathcal{D}}_{j}^{\left(\delta \right)}$  under the change of variables ${x}_{j}=tanh{t}_{j}$  has the form
 $\begin{array}{cccccc}{\mathcal{D}}_{j}& =\delta {x}_{j}+\left(1-{x}_{j}^{2}\right){\partial }_{j}-\frac{a}{2}{\sum }_{i
We fix now the Bruhat ordering of the symmetric group ${S}_{r}$  of $r$  elements with generators $\left(1,2\right),\left(2,3\right),...,\left(r-1,r\right)$  . The corresponding partial ordering on the set ${\mathbb{N}}^{r}$  of all $r$  -tuples of non-negative integers is then defined by the following: For any $\eta$  let ${\eta }^{*}$  be the unique partition in the ${S}_{r}$  -orbit of $\eta$  , and ${w}^{\eta }$  the unique element in ${S}_{r}$  with shortest length such that $\eta ={w}^{\eta }{\eta }^{*}$  . Then we define $\eta \le \mu$  if ${\eta }^{*}\le {\mu }^{*}$  , and if ${w}^{\eta }\le {w}^{\xi }$  when ${\eta }^{*}={\mu }^{*}$  .
The following lemma can be proved by direct computation, it is similar to the known results that the actions of the operators ${D}_{j}$  on $\mathbb{C}\left[P\right]$  are upper triangular.
Lemma 4.1. The action of the operators ${\mathcal{D}}_{j}$  on the monomials ${x}^{\eta }={x}_{1}^{{\eta }_{1}}\cdots {x}_{r}^{{\eta }_{r}}$  , $\eta \in {\mathbb{N}}^{r}$  , preserves the order in the sense that ${\mathcal{D}}_{j}{x}^{\eta }={a}_{{\eta }^{j},\eta }{x}^{{\eta }^{j}}+{\sum }_{\zeta <{\eta }^{j}}{a}_{\zeta ,\eta }{x}^{\xi },$  where ${\eta }^{j}=\eta +\left(0,\cdots ,0,1,0,\cdots ,0\right)$  and where the leading coefficient is ${a}_{{\eta }^{j},\eta }=\delta +a#\left\{i{\eta }_{j}\right\}+\frac{\iota }{2}\left(1-\left(-1{\right)}^{{\eta }_{j}}\right).$
Using the above lemma we can compute ${\mathcal{D}}_{1}^{2{\eta }_{1}}\cdots {\mathcal{D}}_{r}^{2{\eta }_{r}}1$  successively and get
 $\begin{array}{c}{\mathcal{D}}_{1}^{2{\eta }_{1}}\cdots {\mathcal{D}}_{r}^{2{\eta }_{r}}1={d}_{\eta }{x}^{2\eta }+\text{lower terms}.\end{array}$ (4.1)
where
 $\begin{array}{c}{d}_{\eta }{=}^{r}{{\prod }_{j=0}}^{{\eta }_{j}-1}{\prod }_{k=0}\left(\delta +\left(r-j\right)-2k\right)\left(\delta +\left(r-j\right)a+\iota -1-2k\right)\end{array}$ (4.2)
and where the “lower terms” here stands for a linear combination of ${x}_{1}^{{\mu }_{1}}\cdots {x}_{r}^{{\mu }_{r}}$  with $\mu \prec 2\eta$  . Denoting ${m}_{\eta }$  the symmetric power sum
 $\begin{array}{c}{m}_{\eta }\left({y}_{1},\cdots ,{y}_{r}\right)={\sum }_{\zeta \in {S}_{r}\eta }{y}_{1}^{{\zeta }_{1}}\cdots {y}_{r}^{{\zeta }_{r}}\end{array}$ (4.3)
and taking summation of ( 4.1 ) over all permutations $\zeta \in {S}_{r}\eta$  of $\eta$  , we have,
Lemma 4.2. The operators ${m}_{\eta }\left({\mathcal{D}}_{1}^{2},\cdots ,{\mathcal{D}}_{r}^{2}\right)$  acting on the constant monomial $1$  gives ${m}_{\eta }\left({\mathcal{D}}_{1}^{2},\cdots {\mathcal{D}}_{r}^{2}\right)1={d}_{\eta }{m}_{\eta }\left({x}_{1}^{2},\cdots {x}_{r}^{2}\right)+\text{lower terms},$  where the “lower terms” here stands for a linear combination of symmetric power sums ${m}_{{\eta }^{\prime }}$  of ${x}_{1}^{2},\cdots ,{x}_{r}^{2}$  with ${\eta }^{\prime }\prec \eta$  .
Proposition 4.3. Let $\nu >\iota +b+a\left(r-1\right)$  and let $\eta :{\eta }_{1}\ge {\eta }_{2}\cdots \ge {\eta }_{r}\ge 0$  be a partition. The spherical transforms of a function ${f}_{-2\nu }\left(t\right){m}_{\eta }\left({tanh}^{2}{t}_{1},\cdots ,{tanh}^{2}{t}_{r}\right)$  is given by $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right){l}_{\eta }\left(\lambda \right),$  and ${l}_{\eta }\left(\lambda \right)$  is $W$  -invariant polynomial of $\lambda$  with leading term ${d}_{\eta }^{-1}{m}_{\eta }\left(\lambda \left({\xi }_{1}{\right)}^{2},\cdots ,\lambda \left({\xi }_{r}{\right)}^{2}\right)$  .
• Proof. We apply Lemma 4.2 for $\delta =-2\nu$  . Observe that ${d}_{\eta }\ne 0$  for any partition $\eta$  . It follows that for each partition $\eta$  , ${m}_{\eta }\left({x}_{1}^{2},\cdots ,{x}_{r}^{2}\right)={d}_{2\eta }^{-1}{m}_{\eta }\left({\mathcal{D}}_{1}^{2},\cdots {\mathcal{D}}_{r}^{2}\right)1+{\sum }_{{\eta }^{\prime }\prec \eta }{c}_{{\eta }^{\prime },\eta }{m}_{{\eta }^{\prime }}\left({\mathcal{D}}_{1}^{2},\cdots {\mathcal{D}}_{r}^{2}\right)1.$  In other words ${f}_{-2\nu }{m}_{\eta }\left({x}_{1}^{2},\cdots ,{x}_{r}^{2}\right)={d}_{2\eta }^{-1}{m}_{\eta }\left({D}_{1}^{2},\cdots ,{D}_{r}^{2}\right){f}_{-2\nu }+{\sum }_{{\eta }^{\prime }\prec \eta }{c}_{{\eta }^{\prime },\eta }{m}_{{\eta }^{\prime }}\left({D}_{1}^{2},\cdots ,{D}_{r}^{2}\right){f}_{-2\nu }.$  For each $\zeta$  , the spherical transform of ${m}_{\zeta }\left({D}_{1}^{2},\cdots ,{D}_{r}^{2}\right){f}_{-2\nu }$  is given by, in view of ( 1.1 ), ${m}_{\zeta }\left(\lambda \left({\xi }_{1}{\right)}^{2},\cdots ,\lambda \left({\xi }_{r}{\right)}^{2}\right)\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right).$  Our result then follows from the previous formula.

5 Jacobi type functions and Macdonald-Koornwinder polynomials

We recall briefly the Heckman-Opdam theory on Jacobi polynomials. We will then construct certain functions of Jacobi type and prove that their spherical transforms are Macdonald-Koornwinder type polynomials.
Let $\nu >\iota +b+a\left(r-1\right)$  . Consider now the same root system $R=\left\{2{\varepsilon }_{j};j=1,\cdots ,r\right\}\cup \left\{4{\varepsilon }_{j};j=1,\cdots ,r\right\}\cup \left\{2{\varepsilon }_{j}±{\varepsilon }_{k};j\ne k=1,\cdots ,r\right\}$  with the same ordering and with respective multiplicities ${k}^{\left(\nu \right)}=\left({k}_{1}^{\left(\nu \right)},{k}_{2}^{\left(\nu \right)},{k}_{3}^{\left(\nu \right)}\right)$  , such that $2{k}_{2}^{\left(\nu \right)}=2\left(2\nu -\left(1+\iota +b+a\left(r-1\right)\right)\right)+1,2\left({k}_{1}^{\left(\nu \right)}+{k}_{2}^{\left(\nu \right)}\right)=\iota +2b,2{k}_{3}^{\left(\nu \right)}=a.$  (Note that one of the multiplicities, ${k}_{1}^{\left(\nu \right)}$  , is negative, however they still satisfy the condition ${k}_{\alpha }+{k}_{\frac{\alpha }{2}}\ge 0$  , $\forall \alpha \in R$  , see [7.) The corresponding dual coroot lattice is $\frac{1}{2}\left(\mathbb{Z}{\varepsilon }_{1}+\cdots +\mathbb{Z}{\varepsilon }_{r}\right)=\left(\frac{1}{2}\mathbb{Z}{\right)}^{r}$  and weight lattice is $P=2\left(\mathbb{Z}{\varepsilon }_{1}+\cdots +\mathbb{Z}{\varepsilon }_{r}\right)=\left(2\mathbb{Z}{\right)}^{r}$  . The set of dominant weights is ${P}^{+}=\left\{2\left({\eta }_{1}{\varepsilon }_{1}+\cdots +{\eta }_{r}{\varepsilon }_{r}\right);\eta =\left({\eta }_{1},\cdots ,{\eta }_{r}\right)\in {\mathbb{N}}^{r},{\eta }_{1}\ge \cdots \ge {\eta }_{r}\ge 0\right\},$  and for simplicity we will index the elements in ${P}^{+}$  by the partitions $\eta \in {\mathbb{N}}^{r}$  .
The polynomial algebra $\mathbb{C}\left[P\right]$  of integral weights is then the polynomial algebra generated by ${e}^{±2{\varepsilon }_{j}},j=1,\cdots ,r$  . Consider the inner product in the space $\mathbb{C}\left[P\right]$  ,
 $\begin{array}{c}\left(f,g{\right)}_{{k}^{\left(\nu \right)}}={\int }_{\left[0,\pi {\right]}^{r}}f\left(s\right)\overline{g\left(s\right)}|{\delta }_{{k}^{\left(\nu \right)}}\left(s\right)|ds,\end{array}$ (5.1)
where ${\delta }_{{k}^{\left(\nu \right)}}\left(s\right)={\prod }_{\alpha \in R}|2sin\frac{1}{2}\alpha \left(s\right)\right){|}^{2{k}_{\alpha }^{\left(\nu \right)}}$  .
We recall first the Heckman-Opdam theory of Jacobi polynomials, see [4, [3, [7. The partial ordering on the weight lattice $P$  defined in terms by the ordering of the roots will be denoted (with some abuse of notation) also by $\le$  .
Lemma 5.1. For each $\eta =\left({\eta }_{1},\cdots ,{\eta }_{r}\right)\in {\mathbb{N}}^{r}$  , ${\eta }_{1}\ge \cdots \ge {\eta }_{r}\ge 0$  , there exists a unique polynomial ${P}_{\nu ,\eta }$  on $T$  such that ${P}_{\nu ,\eta }={p}_{\eta }^{W}+{\sum }_{{\eta }^{\prime }<\eta }{c}_{{\eta }^{\prime },\eta }{p}_{\eta }^{W}$  and $\left({P}_{\nu ,\eta },{p}_{\eta }^{W}{\right)}_{{k}^{\left(\nu \right)}}=0,$  where ${p}_{\eta }^{W}={\sum }_{w\in W}w\left({e}^{2\left({\eta }_{1}{\varepsilon }_{1}+\cdots +{\eta }_{r}{\varepsilon }_{r}\right)}\right)$  is the Weyl group orbit sum of the power function. The polynomials $\left\{{P}_{\nu ,\eta }{\right\}}_{\eta }$  forms an orthogonal basis for ${L}^{2}\left(T,|{\delta }_{{k}^{\left(\nu \right)}}\left(s\right)|ds{\right)}^{W}$  .
The inner product $〈{P}_{\nu ,\eta },{P}_{\nu ,\eta }{〉}_{{k}^{\left(\nu \right)}}$  is explicitly computed in [4and [7.
The Jacobi polynomials on $T$  are symmetric with respect to the Weyl group $W$  , thus they are symmetric polynomials of the functions $\frac{{e}^{2i{s}_{j}}+{e}^{-2i{s}_{j}}}{2}=cos2{s}_{j}$  , which can be further identified as symmetric polynomials of $\left({x}_{1}^{2},\cdots ,{x}_{r}^{2}\right)$  , with
 $\begin{array}{c}{x}_{j}:=sin{s}_{j}\in \left[-1,1\right],s={s}_{1}{\xi }_{1}+\cdots +{s}_{r}{\xi }_{r}\in T=i{\mathbb{R}}^{r}/\left(i\pi \mathbb{Z}{\right)}^{r}.\end{array}$ (5.2)
With some abuse of notations we denote it also by ${P}_{\nu ,\eta }\left({x}_{1},...,{x}_{r}\right)$  . We rewrite the characterization of ${P}_{\nu ,\eta }$  in terms of the variables $\left({x}_{1},...,{x}_{r}\right)$  .
The Jacobi polynomial ${P}_{\nu ,\eta }\left({x}_{1},...,{x}_{r}\right)$  is then characterized as the unique polynomials so that ${P}_{\nu ,\eta }\left({x}_{1},...,{x}_{r}\right)={2}^{2{\eta }_{1}+\cdot \cdot \cdot +2{\eta }_{r}}{m}_{\eta }\left({x}_{1}^{2},...,{x}_{r}^{2}\right)+{\sum }_{\zeta <\eta }{c}_{\eta ,\zeta }{m}_{\zeta }\left({x}_{1}^{2},...,{x}_{r}^{2}\right)$  and that $\left({P}_{\nu ,\eta },{m}_{\zeta }{\right)}_{{k}^{\left(\nu \right)}}=0,\zeta <\eta$  where $\left(f,g{\right)}_{{k}^{\left(\nu \right)}}$  is the inner product ( 5.1 ) after the change of variables ( 5.2 ),
 $\begin{array}{c}\begin{array}{cccc}\left(f,g{\right)}_{{k}^{\left(\nu \right)}}& ={2}^{\left(2\iota +2b+a\left(r-1\right)\right)r}{\int }_{\left[-1,1{\right]}^{r}}f\left(x\right)\overline{g\left(x\right)}& & ×{\prod }_{i (5.3)
We now define the Jacobi type function on $\mathfrak{a}$  , ${H}_{\nu ,\eta }\left(t\right):={f}_{-2\nu }\left(t\right){P}_{\nu ,\eta }\left(tanh{t}_{1},\cdots ,tanh{t}_{r}\right)$
Lemma 5.2. The functions $\left\{{H}_{\nu ,\eta }\left(t\right){\right\}}_{\eta }$  form an orthogonal basis for the space ${L}^{2}\left(\mathfrak{a}{\right)}^{W}$  and $〈{H}_{\nu ,\eta },{H}_{\nu ,\eta }{〉}_{{L}^{2}\left(\mathfrak{a}\right)}=〈{P}_{\nu ,\eta },{P}_{\nu ,\eta }{〉}_{{k}^{\left(\nu \right)}}$
• Proof. Changing of variables $t\in \mathfrak{a}\to x\in \left(-1,1{\right)}^{r}$  , ${x}_{j}:=tanh{t}_{j}$  in the definition of inner product in ${L}^{2}\left(\mathfrak{a}\right)$  , we get the above formula. The remainding claim follows from (the symmetric version of ) Weierstrass approximation theorem.
Theorem 5.3. The spherical transforms of the Jacobi type function ${H}_{\nu ,\eta }\left(t\right)$  on $\mathfrak{a}$  are given by $\stackrel{~}{{f}_{-2\nu }}\left(\lambda \right){q}_{\eta }\left(\lambda \right),$  where ${q}_{\eta }\left(\lambda \right)$  are constant multiples of the Macdonald-Koornwinder polynomials, and they satisfy the orthogonality relation ${\int }_{{\mathfrak{a}}^{*}}\stackrel{~}{{f}_{-2\nu }}\left(\lambda {\right)}^{2}{q}_{\eta }\left(\lambda \right){q}_{{\eta }^{\prime }}\left(\lambda \right)|c\left(\lambda \right){|}^{-2}d\lambda =〈{H}_{\nu ,\eta },{H}_{\nu ,\eta }{〉}_{{L}^{2}\left(\mathfrak{a}\right)}=〈{P}_{\nu ,\eta },{P}_{\nu ,{\eta }^{\prime }}{〉}_{{k}^{\left(\nu \right)}}={\delta }_{\eta ,{\eta }^{\prime }}{P}_{\nu ,\eta },{P}_{\nu ,\eta }{〉}_{{k}^{\left(\nu \right)}}.$
• Proof. That the spherical transform is of the given form follows from Proposition 4.3, the orthogonality relation is by the Plancherel formula ( 1.2 ), and that ${q}_{{\eta }^{\prime }}\left(\lambda \right)$  are the Macdonald-Koornwinder polynomials is by the uniqueness of the latter [9.
References

1. I. Cherednik, Inverse Harish-Chandra transform and difference operators, Internat. Math. Res. Notices (1997), no. 15, 733–750. MR MR1470375 (99d:22018)
2. A. Erdelyi et al, Higher transdental functions, vol. 1, McGraw-Hill, New York Toronto London, 1953.
3. G. J. Heckman, Root systems and hypergeometric functions. II, Compositio Math. 64 (1987), no. 3, 353–373.
4. G. J. Heckman and E. M. Opdam, Root systems and hypergeometric functions. I, Compositio Math. 64 (1987), no. 3, 329–352.
5. T. H. Koornwinder, Special orthogonal polynomial systems mapped onto each other by the Fourier-Jacobi transform, Orthogonal polynomials and applications (Bar-le-Duc, 1984), Lecture Notes in Math., vol. 1171, Springer, Berlin, 1985, pp. 174–183. MR MR838982 (87g:33007)
6. I. G. Macdonald, Symmetric functions and Hall polynomials, second ed., Clarendon Press, Oxford, 1995.
7. E. Opdam, Harmonic analysis for certain representations of graded Hecke algebras, Acta Math. 175 (1995), no. 1, 75–121.
8. A. Unterberger and H. Upmeier, The Berezin transform and invariant differential operators, Comm. Math. Phys. 164 (1994), 563–597.
9. J. F. van Diejen, Properties of some families of hypergeometric orthogonal polynomials in several variables, Trans. Amer. Math. Soc. 351 (1999), no. 1, 233–270. MR MR1433128 (99c:33020)
10. G. van Dijk and M. Pevzner, Berezin kernels and tube domains, J. Funct. Anal. 181 (2001), 189–208.
11. G. Zhang, Berezin transform on compact hermitian symmetric spaces, Manuscripta Math. 97 (1998), no. 3, 371–388.
12. , Invariant differential operators on hermitian symmetric spaces and their eigenvalues, Israel J. Math. 119 (2000), 157–185.
13. , Berezin transform on real bounded symmetric domains, Trans. Amer. Math. Soc. 353 (2001), 3769–3787.
14. , Branching coefficients of holomorphic representations and Segal-Bargmann transform, J. Funct. Anal. 195 (2002), 306–349.

Department of Mathematics, Chalmers University of Technology and Goteborg University, S-412 96 Goteborg, Sweden E-mail address : genkai@math.chalmers.se