Racah problems for the oscillator algebra, the Lie
algebra sl(n), and multivariate Krawtchouk polynomials

Nicolas Crampé, Wouter van de Vijver, Luc Vinet

► To cite this version:

Nicolas Crampé, Wouter van de Vijver, Luc Vinet. Racah problems for the oscillator algebra, the Lie algebra sl(n), and multivariate Krawtchouk polynomials. Annales Henri Poincaré, 2020, 21, pp.3939-3971. 10.1007/s00023-020-00972-8. hal-02885779

HAL Id: hal-02885779

https://hal.science/hal-02885779v1

Submitted on 10 Mar 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RACAH PROBLEMS FOR THE OSCILLATOR ALGEBRA,
THE LIE ALGEBRA sl<sub>n</sub>,
AND MULTIVARIATE KRAWTCHOUK POLYNOMIALS

NICOLAS CRAMPÉ, WOUTER VAN DE VIJVER, AND LUC VINET

ABSTRACT. The oscillator Racah algebra $\mathcal{R}_n(h)$ is realized by the intermediate Casimir operators arising in the multifold tensor product of the oscillator algebra $h$. An embedding of the Lie algebra $sl_{n-1}$ into $\mathcal{R}_n(h)$ is presented. It relates the representation theory of the two algebras. We establish the connection between recoupling coefficients for $h$ and matrix elements of $sl_n$-representations which are both expressed in terms of multivariate Krawtchouk polynomials of Griffiths type.

CONTENTS

<table><tr><td>1. Introduction</td><td>1</td></tr><tr><td>2. The oscillator algebra h and the Racah algebra</td><td>4</td></tr><tr><td>3. Embedding of sl<sub>n-1</sub> into R<sub>n</sub>(h)</td><td>7</td></tr><tr><td>3.1. Embedding of sl<sub>2</sub> into R<sub>3</sub>(h)</td><td>7</td></tr><tr><td>3.2. Embedding of sl<sub>3</sub> into R<sub>4</sub>(h)</td><td>8</td></tr><tr><td>3.3. Embedding of sl<sub>n-1</sub> into R<sub>n</sub>(h)</td><td>9</td></tr><tr><td>3.4. Labelling Abelian algebras</td><td>12</td></tr><tr><td>4. Connection between recoupling coefficients for h and sl<sub>n-1</sub>-representations</td><td>14</td></tr><tr><td>4.1. sl<sub>2</sub> and the Krawtchouk polynomials</td><td>15</td></tr><tr><td>4.2. sl<sub>n-1</sub> and the multivariate Krawtchouk polynomials</td><td>16</td></tr><tr><td>4.3. 6j- and 9j-symbols</td><td>19</td></tr><tr><td>4.4. Automorphisms of R<sub>n</sub>(h) and sl<sub>n-1</sub></td><td>20</td></tr><tr><td>5. Conclusion</td><td>25</td></tr><tr><td>6. Acknowledgements</td><td>25</td></tr><tr><td>Appendix A. Calculation of overlap coefficients</td><td>25</td></tr><tr><td>References</td><td>27</td></tr></table>

1. INTRODUCTION

This paper studies the oscillator Racah algebra $\mathcal{R}_n(h)$ viewed as the centralizer of the diagonal action of the oscillator algebra $h$ [42] in the n-fold tensor product of its universal algebra. We shall find that it admits an embedding of $sl_{n-1}$. Building upon that result, we shall connect the facts that the multivariate Krawtchouk

<footnote>Date: September 30, 2019.</footnote>
polynomials of Griffiths arise as $3(n-1)j$ symbols of $\mathfrak{h}$ as well as matrix elements of the restriction to the group $O(n+1)$ of the symmetric representations of $SU(n+1)$.

There is growing interest in Racah algebras. These are, in particular, identified in the framework of Racah problems where one looks at the recouplings of tensor products of certain Lie algebras. We shall denote by $n$ the number of factors. The cases with $n=3$ for the Lie algebra $su(2)$ (or $su(1,1)$), the quantum algebra $U_q(\mathfrak{sl}_2)$ and the Lie superalgebra $osp(1|2)$ have first been examined. They have led respectively to the (universal versions of the) Racah algebra $R(3)$ [26, 19, 21] the Askey-Wilson algebra $AW(3)$ [27, 30] and the Bannai-Ito algebra $BI(3)$ [20]. In this picture, where there is an implicit map from the abstract Racah algebra onto the centralizer of the diagonal action of say, $su(2)$, $U_q(\mathfrak{sl}_2)$ or $osp(1|2)$ on their triple product, the images of the three generators of the Racah algebra are expressed in terms of the intermediate Casimir elements. The representations of these algebras encompass the bispectral properties of the orthogonal polynomials bearing the same name that are essentially the Racah or $6j$- coefficients of the corresponding algebras whose triple tensor products are considered. In fact this is how the $AW(3)$ was first identified [49] through its realization in terms of the recurrence and $q$-difference operators of the Askey-Wilson polynomials. These Racah algebras have arisen in numerous contexts. They have appeared as symmetry algebras of superintegrable models [22, 11], are featuring centrally in aspects of algebraic combinatorics [45] and are related to the Leonard pairs [44]. They have been related to Double Affine Hecke Algebras (DAHA) [37, 38, 43] and degenerate cases [23]. Algebras over three strands such as the Temperley-Lieb or Brauer ones that arise in Schur-Weyl duality have been shown to be quotients of Racah algebras [6, 5]. Isomorphisms with certain Kauffman-Skein algebras have been established [3, 4]. Howe duality could be used to relate different presentations [17, 18, 16]. Truncated reflection algebras attached to $U_q(\mathfrak{sl}_2)$, to loop algebra or to the Yangian of $sl_2$ have also been found [1, 2, 7] to lead to $AW(3)$ or Racah algebras. Finally, the description of the Bannai-Ito and Askey-Wilson algebras has been cast recently in the framework of the universal $R$-matrix [8, 9]. This offers sufficient cause already to warrant the exploration of the Racah algebra associated to the oscillator algebra.

The study of Racah algebras as centralizers of $n$-fold tensor products with $n$ larger than 3 has been pursued [12, 13, 10, 40]. The recoupling coefficients in these instances are orthogonal polynomials in many variables. In the case of the generalized Racah algebra for example, bases for representations are obtained by diagonalizing the generators of different maximal Abelian subalgebras [12] and the connection coefficients between two such bases are given in terms of multivariate Racah polynomials of the Tratnik type [46, 25]. Given that the Racah polynomials sit at the top of the finite part of the $q=1$ Askey scheme [34, 35], these Tratnik polynomials provide, through specializations and limits, multivariable extensions of all the finite families of orthogonal polynomials in parallel with what occurs in the univariate situation.

We shall here examine the oscillator Racah algebra $R_n(\mathfrak{h})$ for arbitrary $n$ thereby exploring the structure that encodes the properties of the $3nj$-symbols of the oscillator algebra $\mathfrak{h}$. As will be seen, these are given in terms of multivariate Krawtchouk polynomials. Historically, the $3j$- and $6j$- coefficients were obtained in [47, 33] and found to be given both in terms of univariate Krawtchouk polynomials. Looking at the $9j$- symbols of $\mathfrak{h}$, Zhedanov made the observation [50] that these involve
polynomials in two variables orthogonal with respect to the trinomial distribution and depending on one more parameter than the Tratnik ones. We shall indicate below how this generalizes.

Regarding multivariate Krawtchouk polynomials, it is worth recalling and clar-
ifying the following points. Polynomials in $n$ variables that are orthogonal with
respect to the multinomial distribution were introduced by Griffiths in 1971 using
a generating function [28]; these polynomials involve $\frac{1}{2}n(n-1)$ parameters. For a
review see [14]. The specialization to the Krawtchouk family of the Racah poly-
nomials in $n$ variables introduced by Tratnik in 1991 and mentioned before, yields
polynomials also orthogonal with respect to the multinomial distribution but de-
pending in this case on $n$ parameters only (in addition to the maximal degree $N$
[25]. These two sets are hence not the same and their relation remained unclear
for some time largely because of their intricate parametrizations. The bivariate
Krawtchouk polynomials of Griffiths were rediscovered in 2008 in connection with
a probabilistic model and as limits of the $9j$- symbols of $su(2)$ [29]; they were called
Rahman polynomials for a while. Slightly before, Mizukawa and Tanaka [39] had
related the Griffiths polynomials to character algebras and provided an explicit
formula in terms of Gel'fand-Aomoto hypergeometric series.

Of special relevance to the present article is the group theoretical interpretation
of the multivariable Krawtchouk polynomials of Griffiths that was given by Genest,
Vinet and Zhedanov in [24] where they observe that these polynomials arise in the
matrix elements of the representations of the orthogonal group O(n + 1) that act
on the energy eigenspaces of the isotropic (n + 1)–dimensional harmonic oscilla-
tor. In other words, they have shown that the matrix elements of the restriction to
O(n + 1) of the symmetric representations of SU(n + 1) are expressed in terms of the
Krawtchouk polynomials of Griffiths; the parameters of the polynomials are thus
interpreted as the 1/2n(n − 1) parameters, for instance the Euler angles, that specify
rotations in (n + 1) dimensions. This cogent picture has allowed for a complete
characterization on algebraic grounds of the Griffiths polynomials (recurrence rela-
tions, difference equations, generating function etc.) using the covariance properties
of the oscillator creation and annihilation operators under O(n + 1). Furthermore
this approach clarified the connection between the Griffiths and Tratnik classes of
Krawtchouk polynomials by making explicit that the latter is simply a special case
of the former. For example, in the bivariate case (n = 3), while the Griffiths poly-
nomials with 3 parameters correspond to a general rotation in three dimensions,
the Tratnik ones with 2 parameters, arise from rotations that are only products of
two planar rotations about perpendicular axes. Related to this group theoretical
interpretation is the work of Iliev and Terwilliger [31, 32] (see also [41]) where the
Krawtchouk polynomials appear as overlap coefficients between basis elements for
two modules of sl<sub>n+1</sub>(C) with the basis elements for the representation spaces de-
fined as eigenvectors of two Cartan subalgebras related by an anti-automorphism
specified by the parameters.

The embedding of sl<sub>n−1</sub> into R<sub>n</sub>(h) that we shall construct will provide, besides
its intrinsic algebraic interest, a connection between these two manifestations of the
multivariate Krawtchouk polynomials of Griffiths in matrix elements of represen-
tations and in recoupling coefficients.
The paper is organized as follows. In section 2 we introduce the Racah algebra for the oscillator algebra. We also exhibit some properties and find a number of commutation relations that are needed to prove the main theorem of this paper given in the following section. In section 3 we show how to embed the special linear Lie algebra sl$_{n-1}$ into the Racah algebra for the oscillator algebra. We then study Abelian subalgebras of the Racah algebra for the oscillator algebra related to Cartan algebras of the special linear algebra. These are called labelling Abelian algebras and will be the main tool for section 4. In this section we connect the representation theories of sl$_{n-1}$ and of the Racah algebra for the oscillator algebra. We show how multivariate Krawtchouk polynomials both of Tratnik type and of Griffiths type appear as overlap oefficients between bases of irreducible representations diagonalized by the labelling Abelian algebras. We focus briefly on the relation with the 6j- and 9j-symbols. We finish this section by constructing for a number of overlaps an isomorphism of sl$_{n-1}$ and the corresponding rotation matrix. We explain the link of these rotation matrices with the multivariate Krawtchouk polynomials of Griffiths type. A brief conclusion follows. Appendix A records for reference how the overlaps between representation eigenbases associated to equivalent sl$_{2}$ Cartan generators are obtained in terms of univariate Krawtchouk polynomials.

## 2. THE OSCILLATOR ALGEBRA $\mathfrak{h}$ AND THE RACAH ALGEBRA

The oscillator algebra $\mathfrak{h}$ is the Lie algebra generated by four elements $A_{\pm}$, $A_0$ and a central element $a$ with following defining relations:

$$ (1) \qquad [A_{-}, A_{+}] = a, \qquad [A_{0}, A_{\pm}] = \pm A_{\pm}. $$

The Casimir element $Q$ is contained in the universal enveloping algebra $U(\mathfrak{h})$ and is given by:

$$ (2) \qquad Q := aA_0 - A_+A_- $$

We define the elements of $U(\mathfrak{h})^{\otimes n}$ for $1 \le k \le n$

$$ (3) \qquad A_{0,k} := 1^{\otimes(k-1)} \otimes A_0 \otimes 1^{\otimes(n-k)}, $$

$$ (4) \qquad A_{\pm,k} := 1^{\otimes(k-1)} \otimes A_{\pm} \otimes 1^{\otimes(n-k)}, $$

$$ (5) \qquad a_k := 1^{\otimes(k-1)} \otimes a \otimes 1^{\otimes(n-k)}. $$

and for any subset non-empty $K \subset [n] := \{1, \dots, n\}$

$$ (6) \qquad A_{0,K} := \sum_{k \in K} A_{0,k}, \quad A_{\pm,K} := \sum_{k \in K} A_{\pm,k}, \quad a_K := \sum_{k \in K} a_k. $$

We denote the Lie algebra (isomorphic to $\mathfrak{h}$) generated by the operators $A_{0,K}$, $A_{\pm,K}$ and $a_K$ by $\mathfrak{h}_K$. The Casimir element of this algebra is $Q_K$:

$$ (7) \qquad Q_K := a_K A_{0,K} - A_{+,K} A_{-,K}. $$

The operators $Q_K$ will define the algebra of interest of this article.

**Definition 2.1.** We define the oscillator Racah algebra $\mathcal{R}_n(\mathfrak{h})$ to be the subalgebra of $U(\mathfrak{h})^{\otimes n}$ generated by the elements of the set $\{Q_K \mid K \subset [n] \text{ and } K \neq \emptyset\}$.

**Example 2.2.** The easiest non-trivial example is given for $n = 3$. Then a generator is constructed for every non-empty $K \in [3] = \{1, 2, 3\}$. The set of generators are given by

$$ \{Q_1, Q_2, Q_3, Q_{12}, Q_{13}, Q_{23}, Q_{123}\}. $$
**Proposition 2.3.** The following relations hold in $\mathcal{U}(\mathfrak{h})^{\otimes n}$, for any non-empty $K \subset [n]$:

$$[Q_K, A_{0,[n]}) = 0, \quad [Q_K, A_{\pm, [n]}) = 0, \quad [Q_K, a_{[n]}) = 0.$$

Then $\mathcal{R}_n(\mathfrak{h})$ generates a subalgebra of the centralizer of the oscillator algebra $\mathfrak{h}_{[n]}$ in $\mathcal{U}(\mathfrak{h})^{\otimes n}$.

We wish to find the defining commutation relations obeyed by the generators $Q_K$ of $\mathcal{R}_n(\mathfrak{h})$. First we want to point out this lemma:

**Lemma 2.4.** Let $\{K_p\}_{p=1...k}$ be a set of $k$ disjoint subsets of $[n]$. Define $K_B := \cup_{q \in B} K_q$ with $B \subset [k]$. Consider the following map

$$\theta: \mathcal{R}_k(\mathfrak{h}) \to \mathcal{R}_n(\mathfrak{h}): Q_B \mapsto Q_{KB}.$$

This is an injective morphism. We denote its image by $\mathcal{R}_k^{K_1, \dots, K_k}(\mathfrak{h})$. This algebra is isomorphic to $\mathcal{R}_k(\mathfrak{h})$.

**Example 2.5.** Consider a partition of the set $\{1, 2, 3, 4\}$. For example take $K_1 = \{1\}$, $K_2 = \{2, 4\}$ and $K_3 = \{3\}$. Then we have the following injective morphism of $\mathcal{R}_3(\mathfrak{h})$ into $\mathcal{R}_4(\mathfrak{h})$:

$$\begin{align*}
\theta(Q_1) &= Q_{K_1} = Q_1, & \theta(Q_{12}) &= Q_{K_1 K_2} = Q_{124}, \\
\theta(Q_2) &= Q_{K_2} = Q_{24}, & \theta(Q_{13}) &= Q_{K_1 K_3} = Q_{13}, \\
\theta(Q_3) &= Q_{K_3} = Q_3, & \theta(Q_{23}) &= Q_{K_2 K_3} = Q_{234}, \\
\theta(Q_{123}) &= Q_{K_1 K_2 K_3} = Q_{1234}.
\end{align*}$$

Here we introduced the shortened notation $KL := K \cup L$ for sets $K$ and $L$.

*Proof.* We repeat the strategy in [12, section 4.2] and generalize to the $n$-fold tensor product space. In formulas (6) we constructed an algebra $\mathfrak{h}_K$ isomorphic to $\mathfrak{h}$ acting on the components of tensor product whose indices are in $K$. Consider the algebra generated by the union of $\mathfrak{h}_{K_p}$.

$$\mathcal{U}(\mathfrak{h})^{\otimes k} \cong \langle \mathfrak{h}_{K_p} \rangle_{p \in [k]}.$$

The isomorphism is defined on the generators by $1^{\otimes(p-1)} \otimes X \otimes 1^{\otimes(k-p)} \to X_{K_p}$ where $X$ is one of the generators $A_{\pm}$, $A_0$ or $a$. Inside the algebra $\mathcal{U}(\mathfrak{h})^{\otimes k}$ we find $\mathcal{R}_k(\mathfrak{h})$ generated by the operators $Q_B := Q_{i_1 i_2 \dots i_l}$ with $B \subset [k]$. Their images under this isomorphism are $Q_{KB} := Q_{K_{i_1} K_{i_2} \dots K_{i_l}}$. Hence the operators $Q_{KB}$ with $B \subset [k]$ generate an algebra isomorphic to $\mathcal{R}_k(\mathfrak{h})$ inside $\mathcal{R}_n(\mathfrak{h})$. $\square$

Using the strategy of this proof, we can always replace indices by sets in any relation given. For example, consider the following relation:

$$Q_{123} = Q_{12} + Q_{13} + Q_{23} - Q_1 - Q_2 - Q_3.$$

This can be found by straightforward calculation. Using Lemma 2.4 we automatically have for three disjoint sets $K$, $L$ and $M$

$$(8) \qquad Q_{KLM} = Q_{KL} + Q_{KM} + Q_{LM} - Q_K - Q_L - Q_M.$$

This gives us a number of linear dependencies between the generators $Q_A$. From these dependencies one can prove the following:
Let us generate an algebra from the set {$Q_K | \emptyset \neq K \subset [n]$} using the Lie bracket instead of the ordinary multiplication on $U(\mathfrak{h})^{\otimes n}$. Denote it by $\mathcal{R}_n(\mathfrak{h}, [])$). From the relations above we conclude that $\mathcal{R}_n(\mathfrak{h}, [])$ is generated as a vector space by the operators $Q_K$ and their commutators $[Q_K, Q_L]$ over the field $\mathbb{R}(a_1, \dots, a_n)$. We perform some reductions. First the $Q_i$ and $Q_{ij}$ generate all $Q_K$ by Lemma 2.6. The commutators $[Q_{ij}, Q_{jk}]$ can be written as a linear combination of $[Q_{1j}, Q_{jk}]$ by formula (11). By formula (9) we have $[Q_{1j}, Q_{jk}] = [Q_{1k}, Q_{kj}]$ so we require that $j < k$. It follows then that as a vector space $\mathcal{R}_n(\mathfrak{h}, [])$ is generated by the following set:

$$ (15) \qquad \{Q_i | i \in [n]\} \cup \{Q_{ij} | 1 \le i < j \le n\} \cup \{[Q_{1j}, Q_{jk}] | 1 < j < k \le n\}. $$

We will prove later on that this is a basis for $\mathcal{R}_n(\mathfrak{h}, [])$ as a vector space over the field $\mathbb{R}(a_1, \dots, a_n)$. Moreover, we will prove that the set of equalities (12), (13) and (11) together with Lemma 2.7 and Lemma 2.6 exhausts all commutation relations of $\mathcal{R}_n(\mathfrak{h}, [])$.

### 3. EMBEDDING OF $sl_{n-1}$ INTO $\mathcal{R}_n(\mathfrak{h})$

In this section we study the relationship between the special linear Lie algebra $sl_{n-1}$ and $\mathcal{R}_n(\mathfrak{h})$. The Lie algebra $sl_{n-1}$ is generated by the following set of elements: {$e_{kk+1}, e_{k+1k}, h_k | 1 \le k \le n-2$}. By applying Serre's theorem on $sl_{n-1}$ we are guaranteed of a full set of relations which we give here. To this end we introduce the Cartan matrix:

$$ A_{ij} = \begin{cases} 2 & \text{if } i = j \\ -1 & \text{if } j = i \pm 1 \\ 0 & \text{if } |j - i| > 1 \end{cases}. $$

Here are the Chevalley-Serre relations:

$$ (16) \qquad [h_i, h_j] = 0, $$

$$ (17) \qquad [e_{ii+1}, e_{j+1j}] = \delta_{ij}h_i, $$

$$ (18) \qquad [h_i, e_{jj+1}] = A_{ij}e_{jj+1}, \quad [h_i, e_{j+1j}] = -A_{ij}e_{j+1j}, $$

$$ (19) \qquad \mathrm{ad}(e_{ii+1})^{1-A_{ij}}(e_{jj+1}) = 0, \qquad \text{if } i \neq j $$

$$ (20) \qquad \mathrm{ad}(e_{i+1i})^{1-A_{ij}}(e_{j+1j}) = 0. \qquad \text{if } i \neq j $$

The operator $\mathrm{ad}$ is the adjoint action: $\mathrm{ad}(x)(y) := [x, y]$. The set {$h_k | 1 \le k \le n-2$} generates the Cartan algebra of $sl_{n-1}$. When we consider the $sl_2$ case, there are three generators {$e_{12}, e_{21}, h_1$} with following relations:

$$ (21) \qquad [e_{12}, e_{21}] = h_1, \quad [h_1, e_{12}] = 2e_{12}, \quad [h_1, e_{21}] = -2e_{21}. $$

From the Chevalley-Serre relations we see that every triple {$e_{ii+1}, e_{i+1i}, h_i$} generates a copy of $sl_2$.

### 3.1. Embedding of $sl_2$ into $\mathcal{R}_3(\mathfrak{h})$.
Consider the algebra $\mathcal{R}_3(\mathfrak{h})$ and consider the adjoint action of $Q_{12}$ on $\mathcal{R}_3(\mathfrak{h})$. We want to find its eigenspaces. The eigenspace with eigenvalue 0 is a five-dimensional space generated by the central elements {$Q_1, Q_2, Q_3, Q_{123}$} and the operator $Q_{12}$. The eigenvectors with nonzero eigenvalue
are

$$ (22) \qquad \begin{aligned} e_{12} &:= \lambda_1([Q_{12}, [Q_{12}, Q_{23}]] + (a_1 + a_2)[Q_{12}, Q_{23}]), \\ e_{21} &:= \lambda_1([Q_{12}, [Q_{12}, Q_{23}]] - (a_1 + a_2)[Q_{12}, Q_{23}]). \end{aligned} $$

We introduced the number

$$ \lambda_1 = \frac{1}{\sqrt{4a_1a_2a_3(a_1+a_2)^2(a_1+a_2+a_3)}}. $$

One checks easily using relation (14) that:

$$ (23) \qquad \begin{aligned} [Q_{12}, e_{12}] &= (a_1 + a_2)e_{12}, \\ [Q_{12}, e_{21}] &= -(a_1 + a_2)e_{21}. \end{aligned} $$

We also define the following operator

$$ (24) \qquad h_1 := \frac{2Q_{12}}{a_1 + a_2} - \frac{Q_{123}}{a_1 + a_2 + a_3} - \frac{Q_1}{a_1} - \frac{Q_2}{a_2} + \frac{Q_3}{a_3}. $$

We have the following proposition:

**Proposition 3.1.** The operators $e_{12}$, $e_{21}$ and $h_1$ satisfy the commutation relations of $sl_2$.

$$ [e_{12}, e_{21}] = h_1, \quad [h_1, e_{12}] = 2e_{12}, \quad [h_1, e_{21}] = -2e_{21} $$

*Proof.* One checks through straightforward calculation. $\square$

<figure><img src="image_2.png" /></figure>

By Proposition 3.1 we have a map of $sl_2$ into $\mathcal{R}_3(\mathfrak{h})$. Because $sl_2$ is simple, the kernel of this map is either trivial of equal to the whole algebra $sl_2$. Clearly, $e_{12}$, $e_{21}$ and $h_1$ are different from 0 so the kernel must be trivial. This map must therefore be injective and we have indeed an embedding of $sl_2$ into $\mathcal{R}_3(\mathfrak{h})$.

**3.2. Embedding of $sl_3$ into $\mathcal{R}_4(\mathfrak{h})$.** Consider the algebra $\mathcal{R}_4(\mathfrak{h})$. We want to find the common eigenvectors of $Q_{12}$ and $Q_{123}$. By Lemma 2.7 the operators $e_{12}$, $e_{21}$ and $h_1$ commute with $Q_{123}$ and are therefore eigenvectors of both $Q_{12}$ and $Q_{123}$. We can find another set of eigenvectors using Lemma 2.4. Consider the operators $e_{12}$, $e_{21}$ and $h_1$ expressed in the operators $Q_{12}$, $Q_{23}$ and $a_{12}$ and replace the indices as follows by Lemma 2.4: $1 \to \{1, 2\}$, $2 \to 3$ and $3 \to 4$. We find the following operators:

$$ e_{23} := \lambda_2([Q_{123}, [Q_{123}, Q_{34}]] + (a_1 + a_2 + a_3)[Q_{123}, Q_{34}]), $$

$$ e_{32} := \lambda_2([Q_{123}, [Q_{123}, Q_{34}]] - (a_1 + a_2 + a_3)[Q_{123}, Q_{34}]), $$

$$ h_2 := \frac{2Q_{123}}{a_1 + a_2 + a_3} - \frac{Q_{1234}}{a_1 + a_2 + a_3 + a_4} - \frac{Q_{12}}{a_1 + a_2} - \frac{Q_3}{a_3} + \frac{Q_4}{a_4}. $$

The number $\lambda_2$ is given by

$$ \lambda_2 = \frac{1}{\sqrt{4(a_1 + a_2)a_3a_4(a_1 + a_2 + a_3)^2(a_1 + a_2 + a_3 + a_4)}}. $$

By Lemma 2.7 the operators $e_{23}$, $e_{32}$ and $h_2$ commute with $Q_{12}$. They are also eigenvectors of $Q_{123}$ by formula 23:

$$ (25) \qquad \begin{aligned} [Q_{123}, e_{23}] &= (a_1 + a_2 + a_3)e_{23}, \\ [Q_{123}, e_{32}] &= -(a_1 + a_2 + a_3)e_{32}. \end{aligned} $$

The operators $e_{23}$, $e_{32}$ and $h_2$ satisfy the $sl_2$ relations:

$$ [e_{23}, e_{32}] = h_2, \quad [h_2, e_{23}] = 2e_{23}, \quad [h_2, e_{32}] = -2e_{32}. $$
We have the following claim:

**Proposition 3.2.** The operators {$e_{12}, e_{21}, e_{23}, e_{32}, h_1, h_2$} satisfy the commutation relations of $sl_3$.

*Proof.* We already know that both {$e_{12}, e_{21}, h_1$} and {$e_{23}, e_{32}, h_2$} satisfy the $sl_2$ relations. It is also easy to show by Lemma 2.7 that

$$[h_1, h_2] = 0.$$

By straightforward calculation

$$[e_{12}, e_{32}] = 0,$$

$$[e_{21}, e_{23}] = 0.$$

By now we verified relations (16) and (17). By explicit calculation using the definition of $h_1$ and $h_2$ and formula (23) and (25):

$$[h_1, e_{23}] = \left[ -\frac{Q_{123}}{a_1 + a_2 + a_3}, e_{23} \right] = -e_{23},$$

$$[h_1, e_{32}] = \left[ -\frac{Q_{123}}{a_1 + a_2 + a_3}, e_{32} \right] = e_{32},$$

$$[h_2, e_{12}] = \left[ -\frac{Q_{12}}{a_1 + a_2}, e_{12} \right] = -e_{12},$$

$$[h_2, e_{21}] = \left[ -\frac{Q_{12}}{a_1 + a_2}, e_{21} \right] = e_{21}.$$

Relation (18) is also satisfied. We only need check (19) and (20):

$$[e_{12}, [e_{12}, e_{23}]] = 0,$$

$$[e_{23}, [e_{23}, e_{12}]] = 0,$$

$$[e_{21}, [e_{21}, e_{32}]] = 0,$$

$$[e_{32}, [e_{32}, e_{21}]] = 0.$$

This is done by straightforward computation. By Serre's Theorem these are a complete set of defining relations and this concludes the proof. □

By Proposition 3.2 we have a map of $sl_3$ into $R_4(h)$. Because $sl_3$ is simple, the kernel of this map is either trivial or equal to the whole algebra $sl_3$. The generators we used in Proposition 3.2 are different from 0. This map must therefore be injective and we have indeed an embedding of $sl_3$ into $R_4(h)$.

**3.3. Embedding of $sl_{n-1}$ into $R_n(h)$**. As in the previous two sections we construct common eigenvectors for the adjoint action of the Abelian subalgebra $\mathcal{Y} = \{Q_{[k]} | 2 \le k \le n-1\}$. Observe that by formula (6) we have $a_B = \sum_{i \in B} a_i$. To construct the eigenvectors we use Lemma 2.4 in the following way: Take $e_{12}, e_{21}$ and $h_1$ given by formulas (22) and (24) and replace 1 by $[k]$, 2 by $k+1$ and 3 by $k+2$. We obtain the following elements:

$$e_{kk+1} := \lambda_k([Q_{[k+1]}, [Q_{[k+1]}, Q_{k+1k+2}]] + a_{[k+1]}[Q_{[k+1]}, Q_{k+1k+2}]),$$

$$e_{k+1k} := \lambda_k([Q_{[k+1]}, [Q_{[k+1]}, Q_{k+1k+2}]] - a_{[k+1]}[Q_{[k+1]}, Q_{k+1k+2}]),$$

$$h_k := \frac{2Q_{[k+1]}}{a_{[k+1]}} - \frac{Q_{[k+2]}}{a_{[k+2]}} - \frac{Q_{[k]}}{a_{[k]}} - \frac{Q_{k+1}}{a_{k+1}} + \frac{Q_{k+2}}{a_{k+2}}.$$
We introduced the element

$$ \lambda_k = \frac{1}{\sqrt{4a_{[k]}a_{k+1}a_{k+2}a_{[k+1]}^2a_{[k+2]}}} $$

We check that these operators are indeed eigenvectors of the Abelian algebra $\mathcal{Y}$. By Lemma 2.7 the operator $Q_{[l]}$ commutes with $e_{kk+1}$ and $e_{k+1k}$ if $l \neq k+1$ and with $h_k$ for all $k$. If $l = k+1$ we have

$$ \begin{aligned} [Q_{[k+1]}, e_{kk+1}] &= \lambda_k([Q_{[k+1]}, [Q_{[k+1]}, [Q_{[k+1]}, Q_{k+1k+2}]]] + a_{[k+1]}[Q_{[k+1]}, [Q_{[k+1]}, Q_{k+1k+2}])) \\ &= \lambda_k(a_{[k]}^2[Q_{[k+1]}, Q_{k+1k+2}] + a_{[k+1]}[Q_{[k+1]}, [Q_{[k+1]}, Q_{k+1k+2}])) \\ &= a_{[k]}e_{kk+1}. \end{aligned} $$

We used formula (14). Similarly one can show that $[Q_{[k+1]}, e_{k+1k}] = -a_{[k]}e_{k+1k}$. We are now ready to prove the following theorem.

**Theorem 3.3.** The set of operators $\{e_{ii+1}, e_{ii+1}, h_i | i \in [n-1]\} \subset \mathcal{R}_n(\mathfrak{h})$ generate an algebra isomorphic to $\mathfrak{sl}_{n-1}$ for the Lie bracket.

*Proof.* We prove this statement by induction. The cases $n=2$ and $n=3$ have already be obtained in Propositions 3.1 and 3.2. Now assume that $\{e_{ii+1}, e_{ii+1}, h_i | i \in [n-2]\}$ generates $\mathfrak{sl}_{n-2}$. To go to $\mathfrak{sl}_{n-1}$ we add three new operators: $e_{n-2n-1}$, $e_{n-1n-2}$ and $h_{n-2}$. We introduce a morphism using Lemma 2.4:

$$ \sigma : \mathcal{R}_{n-1}(\mathfrak{h}) \to \mathcal{R}_n(\mathfrak{h}) $$

by mapping the indices $1 \to \{1,2\}$ and $i \to i+1$ for every $i > 1$. One should notice that $\sigma(e_{k-1k}) = e_{kk+1}$, $\sigma(e_{kk-1}) = e_{k+1k}$ and $\sigma(h_{k-1}) = h_k$. It maps $\mathfrak{sl}_{n-2}$ as a subalgebra of $\mathcal{R}_{n-1}(\mathfrak{h})$ into $\mathcal{R}_n(\mathfrak{h})$. We will use this map a few times.

We check the relations of $\mathfrak{sl}_n$ using Serre's Theorem. By Lemma 2.4 and Proposition 3.1 we know that $e_{kk+1}$, $e_{k+1k}$ and $h_k$ generate a $\mathfrak{sl}_2$ algebra for every $k$:

$$ \begin{aligned} [h_k, e_{kk+1}] &= 2e_{kk+1}, \\ [h_k, e_{k+1k}] &= -2e_{k+1k}, \\ [e_{kk+1}, e_{k+1k}] &= h_k. \end{aligned} $$

We also have by Lemma 2.7:

$$ [h_k, h_l] = 0. $$

Relation (16) is satisfied so the set $\{h_k | k = 1...n-1\}$ plays the role Cartan algebra of $\mathfrak{sl}_{n-1}$. The following relation we need to check is $[e_{kk+1}, e_{l+1l}] = 0$ if $l \neq k$. If $l \notin \{k-1, k, k+1\}$ this is true by Lemma 2.7. Otherwise set $l = k+1$. Then we can use the map $\sigma$:

$$ [e_{kk+1}, e_{k+2k+1}] = \sigma([e_{k-1k}, e_{k+1k}]) = \sigma(0) = 0. $$

Here we used the relations of the algebra $\mathfrak{sl}_{n-2}$. The proof for $l=k-1$ is analogous. Relation (17) is satisfied. The next commutators we need to calculate are $[h_k, e_{jj+1}]$
and $[h_k, e_{j+1j}]$:

$$
\begin{align*}
[h_k, e_{jj+1}] &= \left[ \frac{2Q_{[k+1]}}{a_{[k+1]}} - \frac{Q_{[k+2]}}{a_{[k+2]}} - \frac{Q_{[k]}}{a_{[k]}}, e_{jj+1} \right] \\
&= 2\delta_{k,j}e_{jj+1} - \delta_{k+2,j+1}e_{jj+1} - \delta_{k,j+1}e_{jj+1}, \\
[h_k, e_{j+1j}] &= \left[ \frac{2Q_{[k+1]}}{a_{[k+1]}} - \frac{Q_{[k+2]}}{a_{[k+2]}} - \frac{Q_{[k]}}{a_{[k]}}, e_{j+1j} \right] \\
&= -2\delta_{k,j}e_{j+1j} + \delta_{k+2,j+1}e_{j+1j} + \delta_{k,j+1}e_{j+1j}.
\end{align*}
$$

If $k \notin \{j-1, j, j+1\}$ then $[h_k, e_{jj+1}] = 0$ and $[h_k, e_{j+1j}] = 0$. Otherwise, we have

$$
\begin{align*}
[h_{j-1}, e_{jj+1}] &= -e_{jj+1}, \\
[h_{j+1}, e_{jj+1}] &= -e_{jj+1}, \\
[h_{j-1}, e_{j+1j}] &= e_{j+1j}, \\
[h_{j+1}, e_{j+1j}] &= e_{j+1j}.
\end{align*}
$$

We have verified relation (18). Finally we need to check the relations (19) and (20).
Specifically, we need to check the following:

$$
(26) \qquad [e_{jj+1}, e_{kk+1}] = 0, \quad \text{if } k \notin \{j-1, j+1\},
$$

$$
(27) \qquad [e_{j+1j}, e_{k+1k}] = 0, \quad \text{if } k \notin \{j-1, j+1\},
$$

$$
(28) \qquad [e_{jj+1}, [e_{jj+1}, e_{kk+1}]] = 0, \quad \text{if } k \in \{j-1, j+1\},
$$

$$
(29) \qquad [e_{j+1j}, [e_{j+1j}, e_{k+1k}]] = 0, \quad \text{if } k \in \{j-1, j+1\}.
$$

The relations (26) and (27) follow by Lemma 2.7. The relations (28) and (29) can
be proven as follows.

$$
[e_{jj+1}, [e_{jj+1}, e_{kk+1}]] = \sigma([e_{j-1j}, [e_{j-1j}, e_{k-1k}]] = \sigma(0) = 0,
$$

$$
[e_{j+1j}, [e_{j+1j}, e_{k+1k}]] = \sigma([e_{jj-1}, [e_{jj-1}, e_{kk-1}]] = \sigma(0) = 0.
$$

Here we used the sl$_{n-2}$ relations. We have shown that the algebra generated by
${e_{ii+1}, e_{ii+1}, h_i | i \in [n-1]}$ is homomorphic to sl$_{n-1}$. To show that it is in fact
isomorphic we need to prove that the map from sl$_{n-1}$ into R$_n$($h$) is injective. This
is seen by the same argument as before based on the simplicity of sl$_{n-1}$ and thus
we have an embedding
$\square$

We state the following corollary.

**Corollary 3.4.** As Lie algebras over the field $\mathbb{K} := \mathbb{R}(a_1, \dots, a_n)$ we have the following isomorphism:

$$
\mathcal{R}_n(\mathfrak{h}, \mathbb{Q}) \cong \mathfrak{sl}_{n-1} \oplus \mathbb{K}^{n+1}.
$$

Proof. The set (15) generates $\mathcal{R}_n(\mathfrak{h}, \mathbb{Q})$ as a vector space over $\mathbb{K} := \mathbb{R}(a_1, \dots, a_n)$
but we do not know yet if this is a basis. We have $n+1$ central elements, $\frac{n(n-1)}{2}$
elements of the form $Q_{ij}$ and $\frac{(n-1)(n-2)}{2}$ elements of the form $[Q_{1j}, Q_{jk}]$. We have
therefore

$$
(30) \quad \dim(\mathcal{R}_n(\mathfrak{h}, \mathbb{Q})) \le n+1 + \frac{n(n-1)}{2} + \frac{(n-1)(n-2)}{2} = n^2 - n + 1.
$$

By Theorem 3.3 we know there exists a subalgebra of $\mathcal{R}_n(\mathfrak{h}, \mathbb{Q})$ isomorphic with $\mathfrak{sl}_{n-1}$. Denote this Lie algebra by $\mathfrak{s}$. The dimension of $\mathfrak{s}$ equals $(n-1)^2 - 1$. We
add to this algebra s the vector space generated by the central elements

$$Z := \langle \{Q_i | i \in [n]\} \cup \{Q_{[n]}\} \rangle.$$

First notice that $s \cap Z = \{0\}$. The algebra s is simple so it does not contain central elements. Second we see that $\{Q_i | i \in [n]\} \cup \{Q_{[n]}\}$ is a basis for Z. If this was not the case, we could find a linear combination of the central elements equal to 0. We show that this is not possible. The $Q_i$ are linearly independent by construction. Assume we have the linear combination:

$$Q_{[n]} = \sum_{i=1}^{n} \lambda_i Q_i.$$

Act with $\mu^* \otimes \bigotimes_{i=1}^{n-1} 1$ and take the commutator with $Q_{23}$. We find

$$0 = [Q_{23}, Q_{[n+1]}] = \left[ Q_{23}, \lambda_1 Q_{12} + \sum_{i=2}^{n} \lambda_i Q_{i+1} \right] \\ = \lambda_1 [Q_{23}, Q_{12}].$$

The element $[Q_{23}, Q_{12}]$ does not equal 0 so $\lambda_1 = 0$. Similarly one can show that all $\lambda_i = 0$ concluding that $\{Q_i | i \in [n]\} \cup \{Q_{[n]}\}$ is indeed a basis for Z. The dimension of Z equals $n+1$. We have

$$s \oplus Z \subset \mathcal{R}_n(h, \emptyset).$$

From this follows that

$$n^2 - n + 1 = \dim(s \oplus Z) \leq \dim(\mathcal{R}_n(h, \emptyset)).$$

We already have an upper bound by inequality (30). Therefore it follows that $\dim(\mathcal{R}_n(h, \emptyset)) = n^2 - n + 1$ and hence

$$\mathcal{R}_n(h, \emptyset) = s \oplus Z \cong s l_{n-1} \oplus \mathbb{K}^{n+1}.$$

This concludes the proof.

□

Theorem 3.3 and Corollary 3.4 both follow from the set of equations calculated in section 2. As a consequence of the isomorphism found in Corollary 3.4 we give an alternative definition of $\mathcal{R}_n(h, \emptyset)$. It is the algebra defined over $\mathbb{R}(a_1, \dots, a_n)$ with generators $Q_i$ and $Q_{ij}$ satisfying the set of equalities (11), (12) and (13), together with Lemma 2.7 and Lemma 2.6.

**3.4. Labelling Abelian algebras.** In the proof of Theorem 3.3 the elements of the Cartan algebra of $sl_{n-1}$ are linear combinations of the central elements and $\{Q_{[k]} | 2 \leq k \leq n-1\}$. It is possible to construct different Cartan algebras if one starts from a different set of generators.

**Definition 3.5.** A set of non-empty subsets of $[n]$ $\mathcal{A}$ is a maximal non-intersecting/ nested set if it satisfies the following properties:

* Every pair of sets in $\mathcal{A}$ is either disjoint or one is included in the other.

* Maximality: It is not possible to add another subset of $[n]$ to $\mathcal{A}$ without contradicting the first property.

By Lemma 2.7 the corresponding generators $Q_A$ with $A \in \mathcal{A}$ will be a set of commuting operators. To each maximally non-intersecting/nested set $\mathcal{A}$ we can associate a graph $T_A$. Let every vertex of the graph represent a set in $\mathcal{A}$. There is an edge between to vertices A and B if either set is included in the other but there is no set $C \in \mathcal{A}$ such that $A \subset C \subset B$ or $B \subset C \subset A$.
**Lemma 3.6.** The graph $T_A$ related to the maximal non-intersecting/nested set $\mathcal{A}$ is a perfect binary tree.

*Proof.* The set $[n]$ is always included in $\mathcal{A}$ by maximality. The set $[n]$ includes every set in $\mathcal{A}$ and will be the root of the tree. The sets with one element are also included in $\mathcal{A}$ by maximality. They are either completely included or disjoint with every set in $\mathcal{A}$. These sets do not have any proper subsets so they cannot have any children and are therefore leaves of the tree. Every other vertex has two children. Assume this is not the case. Let $\mathcal{A}$ have exactly one child $B$. Then we can add $A \setminus B$ to $\mathcal{A}$ contradicting the maximality of $\mathcal{A}$. Assume $\mathcal{A}$ has more than two children. Let $B_1, B_2$ and $B_3$ three sets that are children of $\mathcal{A}$. We can add the set $B_1 \cup B_2$ to $\mathcal{A}$ contradicting the maximality of $\mathcal{A}$. This concludes the proof. $\square$

We have $n$ indices so we have $n$ leaves. A perfect binary tree with $n$ leaves has $n-2$ interior vertices. The generators related to the interior vertices constitute an Abelian algebra with $n-2$ generators.

**Definition 3.7.** Let $\mathcal{A}$ be a maximally non-intersecting/nested set of $[n]$. Then we define the labelling Abelian algebra associated to $\mathcal{A}$ to be:

$$ \mathcal{Y}_{\mathcal{A}} = \{Q_A \mid A \in \mathcal{A} \text{ and } 1 < |A| < n\}. $$

We exclude the elements $Q_i$ and $Q_{[n]}$ as they are central. For each maximally non-intersecting/nested set $\mathcal{A}$ we have constructed a tree $T_A$ and an algebra $\mathcal{Y}_A$. Often we will represent the algebra $\mathcal{Y}_A$ by the tree $T_A$. These trees are similar to the coupling trees introduced in [48].

**Example 3.8.** Consider $\mathcal{R}_3(\mathfrak{h})$. We have three indices 1, 2 and 3. This gives us three possible trees and therefore three different labelling Abelian algebras:

<figure><img src="image_3.png" /></figure>

**Example 3.9.** Consider $\mathcal{R}_4(\mathfrak{h})$. We have four indices 1, 2, 3 and 4 so we have trees with four leaves. We consider two examples:

<figure><img src="image_4.png" /></figure>
We are now in a position to explain the relation between these labelling Abelian algebras and the Cartan algebra of sl$_{n-1}$ inside R$_n(h)$. Assume we have fixed a labelling Abelian algebra $\mathcal{Y}_A$. For every set $A \in \mathcal{A}$ that is not the root or a leaf, we will construct three generators $e_A$, $f_A$ and $h_A$. To do this consider the tree $T_A$. Focus on the subtree consisting of the vertex related to the set A: The children of A which we will call K and L, the parent B of A and the child of B differing from A denoted by M. We have two possibilities:

<figure><img src="image_5.png" /></figure>

We define $e_A$ and $f_A$ as follows:

$$ (31) \qquad \begin{aligned} e_A &:= \lambda_A([Q_A, [Q_A, Q_{LM}]] + a_A[Q_A, Q_{LM}]), \\ f_A &:= \lambda_A([Q_A, [Q_A, Q_{LM}]] - a_A[Q_A, Q_{LM}]). \end{aligned} $$

with

$$ \lambda_A = \frac{1}{\sqrt{4a_K a_L a_M a_A^2 a_B}}. $$

Observe that it is possible to replace $Q_{LM}$ by $Q_{KM}$ in the definitions of $e_A$ and $f_A$ leading to different but equally correct expressions for $e_A$ and $f_A$. To avoid this ambiguity we always choose L to be the right child of A and M the child of B different from A. This, however, fixes an ordering of the leaves. Otherwise it is not possible to speak of a left and right child. Together with the Cartan elements $h_A = [e_A, f_A]$ the set $\{e_A, f_A, h_A\}$ generates the Lie algebra sl$_{n-1}$. This can be proven by repeating the arguments in Theorem 3.3. The maximally non-intersecting/nested set used in Theorem 3.3 is $\{[k] | 1 < k < n\}$. To find the relation between the Cartan algebra $\{h_A | A \in \mathcal{A} \text{ and } 1 < |A| < n\}$ and the labelling Abelian algebra $\mathcal{Y}_A$ we consider explicitly $h_A$:

$$ h_A = \frac{2Q_A}{a_A} - \frac{Q_B}{a_B} - \frac{Q_K}{a_K} - \frac{Q_L}{a_L} + \frac{Q_M}{a_M}. $$

Every Cartan element can be written as a linear combination of elements of the labelling Abelian algebra $\mathcal{Y}_A$ and central elements $Q_i$ and $Q_{[n]}$.

In the next section we will study the representation theory behind $R_n(h)$ and its relation to the Lie algebra sl$_{n-1}$. The relation between Cartan algebras and labelling Abelian algebras will play an important role there.

## 4. CONNECTION BETWEEN RECOUPLING COEFFICIENTS FOR $\mathfrak{h}$ AND sl$_{n-1}$-REPRESENTATIONS

Assume we have a finite dimensional irreducible representation $V$ of $R_n(h)$. Because of corollary 3.4 this is also an irreducible representation for sl$_{n-1}$. In this section we study bases diagonalized by different labelling Abelian algebras or equivalently Cartan algebras of sl$_{n-1}$ and their connection coefficients. We first study the rank one case.
4.1. sl₂ and the Krawtchouk polynomials. Consider the following two Abelian algebras of R₃(h):

<figure><img src="image_6.png" /></figure>

Let {$\psi_k$} be an eigenbasis of $Q_{12}$ and {$\phi_s$} be an eigenbasis of $Q_{23}$. The indices $k$ and $s$ run from 0 to $N$ with $\dim(V) = N + 1$. We are interested in the overlap coefficients $B_{ks}$ between these bases:

$$ (32) \qquad \phi_s = \sum_{k=0}^{N} B_{sk} \psi_k. $$

To study these overlap coefficients we construct the Lie algebra sl₂ inside R₃(h). Corresponding to $Y_{A_1}$, we have the following generators. See also formula (22):

$$ \begin{align*} e &:= \lambda([Q_{12}, [Q_{12}, Q_{23}]] + a_{12}[Q_{12}, Q_{23}]), \\ f &:= \lambda([Q_{12}, [Q_{12}, Q_{23}]] - a_{12}[Q_{12}, Q_{23}]), \\ h &:= \frac{2Q_{12}}{a_{12}} - \frac{Q_{123}}{a_{123}} - \frac{Q_1}{a_1} - \frac{Q_2}{a_2} + \frac{Q_3}{a_3} \end{align*} $$

with $\lambda = \sqrt{4a_1a_2a_3a_{12}^2a_{123}^{-1}}$. The sl₂ Lie algebra related to $Y_{A_2}$ is given by replacing indices $(1 \to 2 \to 3 \to 1)$ in $e, f$ and $h$ given above. We denote these elements by $\tilde{e}, \tilde{f}$ and $\tilde{h}$. The elements $h$ and $\tilde{h}$ diagonalize the bases {$\psi_k$} and {$\phi_s$} respectively.

$$ h\psi_k = \mu_k\psi_k, \qquad \tilde{h}\phi_s = \nu_s\phi_s. $$

The map $\tilde{\gamma}$ is an automorphism of sl₂ so we can write $\tilde{h}$ as a linear combination of $h, e$ and $f$.

$$ \tilde{h} = R_h h + R_e e + R_f f $$

with

$$ \begin{align*} R_h &= -\frac{a_1 a_2 - a_1 a_3 + a_2^2 + a_2 a_3}{a_{12} a_{23}}, \\ R_e &= R_f = 2 \frac{\sqrt{a_1 a_2 a_3 a_{123}}}{a_{12} a_{23}}. \end{align*} $$

Observe that $R_e R_f + R_h^2 = 1$. It is a classical result that the overlap coefficients between two bases related by an inner automorphism of sl₂ are univariate Krawtchouk polynomials [24, 36]. We have therefore relegated the calculations to Appendix A. This is the result:

$$ B_{sk} \sim K_k \left( \frac{\nu_s + N}{2} ; \frac{1 - R_h}{2} ; N \right). $$

The Krawtchouk polynomials depend on 2 values: $R_h$ and $N$. The number $N$ is equal by definition to $\dim(V) - 1$. The dimension of an irreducible representation
of $sl_2$ can be found by considering its Casimir:

$$ (33) \qquad \begin{aligned} C &:= \frac{h^2}{2} + ef + fe \\ &= \frac{1}{2} \left( \frac{Q_1}{a_1} + \frac{Q_2}{a_2} + \frac{Q_3}{a_3} - \frac{Q_{123}}{a_{123}} \right) \left( \frac{Q_1}{a_1} + \frac{Q_2}{a_2} + \frac{Q_3}{a_3} - \frac{Q_{123}}{a_{123}} - 2 \right). \end{aligned} $$

If we act with the Casimir on $\psi_k$ we find:

$$ C\psi_k = \frac{N^2 + 2N}{2}\psi_k. $$

From the action of the Casimir $C$ of $sl_2$ or equivalently the central elements $Q_1, Q_2, Q_3$ and $Q_{123}$ we are able to discern the dimension of the representation and hence the number $N$. For the remainder of the article we will write the dependence on $R_h$ and $C$ explicitly: $B_{sk}(R_h, C)$.

### 4.2. $sl_{n-1}$ and the multivariate Krawtchouk polynomials.
Let $Y_{A_1}$ and $Y_{A_2}$ be two labelling Abelian algebras of $R_n(h)$. Additionally we demand that $A_1$ and $A_2$ differ by only one element:

$$ A_1 \setminus A_2 = \{G_1\}, \quad A_2 \setminus A_1 = \{G_2\}. $$

Let $\{\psi_{\vec{k}}\}$ be diagonalized by the labelling Abelian algebra $Y_{A_1}$ and $\{\psi_{\vec{s}}\}$ be diagonalized by the labelling Abelian algebra $Y_{A_2}$. Let $A_1 = \{A_i | i = 1...n-2\}$, then we have

$$ Q_{A_i} \psi_{\vec{k}} = \mu_{k_i}^{A_i} \psi_{\vec{k}}. $$

For a specific index $j$ it must be that $A_j = G_1$. We want to find the overlap coefficients between the bases $\{\psi_{\vec{k}}\}$ and $\{\phi_{\vec{s}}\}$.

$$ \phi_{\vec{s}} = \sum_{\vec{k}} B_{\vec{s}\vec{k}} \psi_{\vec{k}}. $$

The basis vector $\phi_{\vec{s}}$ is a common eigenvector of the operators $Q_{A_i}$, $i \neq j$ with eigenvalues $\mu_{s_i}^{A_i}$. The basis $\{\psi_{\vec{k}}\}$ also consists of eigenvectors of the operators $Q_{A_i}$, $i \neq j$. The vector $\phi_{\vec{s}}$ must therefore be written as a linear combination of common eigenvectors of $Q_{A_i}$, $i \neq j$ with the same eigenvalues $\mu_{s_i}^{A_i}$, $i \neq j$. It follows that $B_{\vec{s}\vec{k}} = 0$ if $k_i \neq s_i$ for some $i \neq j$. The overlap coefficient can be written as

$$ B_{\vec{s}\vec{k}} = B_{s_j k_j} \prod_{i \neq j} \delta_{s_i k_i} $$

where $\delta_{s_i k_i}$ is the Kronecker delta. To find $B_{s_j k_j}$ consider the common eigenspace $T = \{v \in V | Q_{A_i} v = \mu_{s_i}^{A_i} v \text{ for all } i \neq j\}$. Both $Q_{G_1}$ and $Q_{G_2}$ commute with each $Q_{A_i}$ so they preserve the common eigenspace $T$. In fact $Q_{G_1}$ and $Q_{G_2}$ lie in an algebra isomorphic to $R_3(h)$ that preserves $T$. Let $K = G_1 \setminus G_2$, $L = G_1 \cap G_2$ and $M = G_2 \setminus G_1$. By Lemma 2.4 the algebra $R_3^{K,L,M}(h)$ generated by $\{Q_K, Q_L, Q_M, Q_{KL}, Q_{LM}, Q_{KLM}\}$ is isomorphic to $R_3(h)$. It preserves $T$ as each generator commutes with $Q_{A_i}$ with $i \neq j$. In fact the sets $K, L, M$ and $K \cup L \cup M$ are in $A_1 \cap A_2$. We conclude that $T$ is a representation of $R_3(h)$ with one basis $\{\psi_{\vec{k}}\} \cap T$ diagonalized by $Q_{G_1}$ and the other $\{\phi_{\vec{s}}\} \cap T$ by $Q_{G_2}$. We are basically in the situation discussed in the previous paragraph and represented by the following two trees:
<figure><img src="image_7.png" /></figure>

This means that the overlap coefficients are given by

$$B_{s_j k_j} = B_{s_j k_j} (R_h^{K,L,M}, C_{K,L,M}).$$

The elements $R_h^{K,L,M}$ and $C_{K,L,M}$ are obtained by using the isomorphism between $\mathcal{R}_3(\mathfrak{h})$ and $\mathcal{R}_3^{K,L,M}(\mathfrak{h})$ obtained by replacing the indices 1, 2 and 3 by $K$, $L$ and $M$:

$$R_h^{K,L,M} = - \frac{a_K a_L - a_K a_M + a_L^2 + a_L a_M}{a_K L a_L M},$$

$$C^{K,L,M} = \frac{1}{2} \left( \frac{Q_K}{a_K} + \frac{Q_L}{a_L} + \frac{Q_M}{a_M} - \frac{Q_{KLM}}{a_{KLM}} \right) \left( \frac{Q_K}{a_K} + \frac{Q_L}{a_L} + \frac{Q_M}{a_M} - \frac{Q_{KLM}}{a_{KLM}} - 2 \right).$$

We have shown that the overlap coefficients between two bases diagonalized by labelling Abelian algebras differing by one generators are Krawtchouk polynomials. Let us remove the condition that the labelling Abelian algebras need to differ by one generator. Then the strategy to find the overlap coefficients is to find a series of intermediate bases in such a way that each intermediate basis differs by one generator with the next intermediate basis. For example take the labelling Abelian algebras $\langle Q_{12}, Q_{34} \rangle$ and $\langle Q_{13}, Q_{24} \rangle$ in $\mathcal{R}_4(\mathfrak{h})$. Then we can find a series of intermediate bases:

$$\langle Q_{12}, Q_{34} \rangle \rightarrow \langle Q_{12}, Q_{123} \rangle \rightarrow \langle Q_{13}, Q_{123} \rangle \rightarrow \langle Q_{13}, Q_{24} \rangle.$$

Each step gives us Krawtchouk polynomials as overlaps. More specifically let $\{\psi_{k_1 k_2}^1\}$, $\{\psi_{k_1 l_2}^2\}$, $\{\psi_{s_1 l_2}^3\}$ and $\{\phi_{s_1 s_2}\}$ be the bases diagonalized by each step in this chain. Then the overlap coefficients become:

$$\begin{align*}
\phi_{s_1 s_2} &= \sum_{l_2} B_{s_2 l_2} (R_h^{13,2,4}, C^{13,2,4}) \psi_{s_1 l_2}^3 \\
&= \sum_{l_2 k_1} B_{s_2 l_2} (R_h^{13,2,4}, C^{13,2,4}) B_{s_1 k_1} (R_h^{1,2,3}, C^{1,2,3}) \psi_{k_1 l_2}^2 \\
&= \sum_{l_2 k_1 k_2} B_{s_2 l_2} (R_h^{13,2,4}, C^{13,2,4}) B_{s_1 k_1} (R_h^{1,2,3}, C^{1,2,3}) B_{l_2 k_2} (R_h^{4,3,12}, C^{4,3,12}) \psi_{k_1 k_2}^1.
\end{align*}$$

Here we conclude that the overlap coefficients are

$$(34) \quad B_{\vec{s} \vec{k}} = \sum_{l_2} B_{s_2 l_2} (R_h^{13,2,4}, C^{13,2,4}) B_{s_1 k_1} (R_h^{1,2,3}, C^{1,2,3}) B_{l_2 k_2} (R_h^{4,3,12}, C^{4,3,12}).$$

This gives us a method to calculate connection coefficients between any pair of bases diagonalized by labelling Abelian algebras. We do need to check if there is always a path between two labelling Abelian algebras. To this end we introduce the recoupling graph of $\mathcal{R}_n(\mathfrak{h})$. Let every labelling Abelian algebra be represented by a vertex. Two vertices are connected by an edge if they only differ by one generator.

**Proposition 4.1.** The recoupling graph of $\mathcal{R}_n(\mathfrak{h})$ is connected and its diameter is bounded by $\frac{(n-1)(n-2)}{2}$.
*Proof.* We use an argument applied on the binary trees related to labelling Abelian algebras. On these trees we have two operations:

* A twist: interchanging two children (and their descendant trees) of a vertex. This does not change the labelling Abelian algebra so on the connection graph we stay on the same vertex.

<figure><img src="image_8.png" /></figure>

* A swap: Moving a vertex' child to the other edge. This changes the labelling Abelian algebra by one generator. On the connection graph we move along an edge to a new vertex.

<figure><img src="image_9.png" /></figure>

Because twisting does not change the labelling Abelian algebra, you stay on the same vertex of the connection graph. Every swap on the other hand is related to a step on the connection graph. The proof of connectedness in $\frac{(n-2)(n-1)}{2}$ steps uses induction:

Let $n = 3$ and take any two trees with three leaves. It is easy to see that it takes a single swap combined with a number of twists to get from one tree to the other.

Assume we have proved it for $n-1$. We take two trees with $n$ leaves. We call these trees the initial tree and final tree. In the initial tree there is at least one vertex whose children are leaves. Assume these leaves are labelled $a$ and $b$. Remove those leaves and give the parent vertex the label $a$. We now have a tree with $n-1$ labelled leaves. We removed the label $b$. Remove this label $b$ also from the final tree together with its parent. The final tree now also has $n-1$ labelled leaves. By induction it takes $\frac{(n-2)(n-3)}{2}$ swaps and a number of twists to change the initial tree to the final tree. Add the leaf with label $b$ and its parent again to the final tree where it was removed. In the initial tree we add two leaves to the leaf with label $a$. We remove this label $a$ and add label $a$ and $b$ to the leaves. Only leaf $b$ is in the wrong place. There are at most $n-2$ vertices in between leaf $b$ and were it needs to be in the final tree. It requires n-2 swaps and a number of twists to move leaf $b$ into the right position. The total swaps used to change the initial tree into the final tree is

$$ \frac{(n-2)(n-3)}{2} + n - 2 = \frac{(n-1)(n-2)}{2}. $$

□

**Remark 4.2.** This proof also works for the higher rank Racah algebra for $su(1, 1)$ as in [12] effectively generalizing the connection graph and the proof of connectedness of the connection graph in [12].

We want to conclude this section with a special pair of labelling Abelian algebras of $\mathcal{Y}_1 = \{Q_{[k]} | 2 \le k \le n-1\}$ and $\mathcal{Y}_{n-1} = \{Q_{[2 \dots k]} | 3 \le k \le n\}$. Observe that $\mathcal{Y}_{n-1}$ can be obtained by from $\mathcal{Y}_1$ by cyclicly permuting the indices $i \to i+1$ and $n \to 1$. We can find a path between these two algebras as follows. Define

$$ \mathcal{Y}_l := \{Q_{[k]} | l+1 \le k \le n-1\} \cup \{Q_{[2 \dots k]} | 3 \le k \le l+1\}. $$
These are labelling Abelian algebra with $\mathcal{Y}_{l-1} \setminus \mathcal{Y}_l = Q_{[l]}$ and $\mathcal{Y}_l \setminus \mathcal{Y}_{l-1} = Q_{[2...,l+1]}$. If we determine the overlap associated to each step we find the following connection coefficients:

$$B_{\vec{s}\vec{k}} = \prod_{l=2}^{n-1} B_{s_l k_l} (R_h^{1,[2..l],l+1}, C^{1,[2..l],l+1}).$$

These are the multivariate Krawtchouk polynomials of Tratnik type. They depend on $n-2$ parameters $R_h^{1,[2..l],l+1}$. They are constructed in the same way as the multivariate Racah polynomials were constructed as connections coefficients between labelling Abelian algebras for the Racah algebra in [12]. The number of parameters follows from the number of steps we needed to move through the connection graph. By proposition 4.1 we know that we can find paths up to $\frac{(n-1)(n-2)}{2}$ steps. Further in the paper we will show how to get the multivariate Krawtchouk polynomials of Griffiths type depending on $\frac{(n-1)(n-2)}{2}$ parameters.

### 4.3. 6j- and 9j-symbols.
The 6j-, 9j- and in general the 3nj-symbols for the oscillator algebra $\mathfrak{h}$ can be cast into the framework presented in this article as they are specific overlaps between recoupled bases. Consider the 6j-symbols. Given the algebra $\mathfrak{h} \oplus \mathfrak{h} \oplus \mathfrak{h}$, the 6j-symbols or Racah coefficients are the connection coefficients for coupling the first two oscillator algebras and the last two:

$$(\mathfrak{h} \oplus \mathfrak{h}) \oplus \mathfrak{h} \to \mathfrak{h} \oplus (\mathfrak{h} \oplus \mathfrak{h}).$$

In our framework this is equivalent with finding the overlap between the bases diagonalized by $Q_{12}$ and $Q_{23}$. The 6j-symbols are therefore up to a normalization equal to $B_{sk}(R_h, C)$. One could ask where the 6 $j$'s of the 6j-symbols are. The indices $s$ and $k$ are related to the eigenvalues of $Q_{12}$ and $Q_{23}$ so they are related to $j_{12}$ and $j_{23}$. The number $R_h$ is independent of the representation used. The remaining four $j$'s are hidden in the Casimir $C$. According to formula (2) the Casimir depends on $Q_1, Q_2, Q_3$ and $Q_{123}$ giving the remaining $j_1, j_2, j_3$ and $j_{123}$.

A similar analysis can be given for the 9j-symbols. The 9j-symbols were already identified as bivariate Krawtchouk polynomials in [50]. They are obtained by considering the algebra $\mathfrak{h} \oplus \mathfrak{h} \oplus \mathfrak{h} \oplus \mathfrak{h}$ and the following two bases: the first basis obtained by coupling first factor with the second and the third with the fourth and the second basis defined by coupling the first factor with the third and the second with the fourth. i.e.

$$(1 \oplus 2) \oplus (3 \oplus 4) \to (1 \oplus 3) \oplus (2 \oplus 4).$$

In our framework this is equivalent with finding the overlap coefficients between the bases diagonalized by $\langle Q_{12}, Q_{34} \rangle$ and $\langle Q_{13}, Q_{24} \rangle$. We already calculated the connection coefficients in formula (34). We repeat the solution here:

$$B_{\vec{s}\vec{k}} = \sum_{l_2} B_{s_2 l_2} (R_h^{13,2,4}, C^{13,2,4}) B_{s_1 k_1} (R_h^{1,2,3}, C^{1,2,3}) B_{l_2 k_2} (R_h^{4,3,12}, C^{4,3,12}).$$

These numbers are multivariate Krawtchouk polynomials depending on three parameter $R_h^{13,2,4}, R_h^{1,2,3}$ and $R_h^{4,3,12}$. These numbers $B_{\vec{s}\vec{k}}$ are up to normalization equal to the 9j-symbols. The 9 $j$'s are found in the following way: $\vec{s}$ and $\vec{k}$ are related to the eigenvalues of $\{Q_{12}, Q_{34}\}$ and $\{Q_{13}, Q_{24}\}$ respectively. The others are found by considering the three Casimir elements appearing in the formula. We
can add the following generators:

$$ \{Q_1, Q_2, Q_3, Q_4, Q_{123}, Q_{1234}\}. $$

This gives a total of 10 numbers. This is one too many. The summation runs over $l_1$ which appears not only as an index but also in $C^{1,2,3}$. This number which is related to the eigenvalues of $Q_{123}$ can be considered as being summed away. This leaves us with 9 generators related to the 9 j's in the 9j-symbols. A similar analysis can be done for any $3nj$-symbol.

### 4.4. Automorphisms of $\mathcal{R}_n(h)$ and $sl_{n-1}$.
For each labelling Abelian algebra of $\mathcal{R}_n(h)$ we are able to construct a set of operators that generate $sl_{n-1}$. By Corollary 3.4 these sets of operators must generate the same algebra. This leads to automorphisms of $sl_{n-1}$. In this section we will give a few examples. Additionally we will give the group element of Lie group $SO_{n-1}$ corresponding to each automorphism.

**Example 4.3.** Take the algebra $\mathcal{R}_3(h)$ and construct $sl_2$ from the labelling Abelian algebra $\{Q_{12}\}$: As per formula (31):

$$ e := \lambda([Q_{12}, [Q_{12}, Q_{23}]] + a_{12}[Q_{12}, Q_{23}]), $$

$$ f := \lambda([Q_{12}, [Q_{12}, Q_{23}]] - a_{12}[Q_{12}, Q_{23}]), $$

$$ h := [e, f]. $$

Consider the permutation $(1 \leftrightarrow 2)$. This can be represented by a twist of two branches of our tree as in Proposition 4.1:

<figure><img src="image_10.png" /></figure>

Keeping in mind that $[Q_{12}, Q_{13}] = -[Q_{12}, Q_{23}]$ by formula (9), the images of the generators under this permutation are $\tilde{h} = h, \tilde{e} = -e$ and $\tilde{f} = -f$. The automorphism of $SL_2$ is constructed as follows. Let

$$ (35) \qquad U \begin{pmatrix} \frac{h}{2} & f \\ e & -\frac{h}{2} \end{pmatrix} U^{-1} = \begin{pmatrix} \frac{h}{2} & \tilde{f} \\ \tilde{e} & -\frac{h}{2} \end{pmatrix}. $$

We want to solve for $U$ with $U \in SL_2$. Conjugation by $U$ is the related automorphism of $SL_2$. In this case

$$ U_{(12)} = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}. $$

**Example 4.4.** Next we construct $sl_2$ from the labelling Abelian algebra $\{Q_{23}\}$. We are constructing an isomorphism corresponding to these two trees:

<figure><img src="image_11.png" /></figure>
This is a swap as in Proposition 4.1. We obtain the second set of generators by permuting the permutation $(1 \to 2 \to 3 \to 1)$. The transformed generators $\tilde{h}$, $\tilde{e}$ and $\tilde{f}$ are linear combinations of the original generators $h$, $e$ and $f$:

$$ \tilde{h} = 2 \frac{\sqrt{a_1 a_2 a_3 a_{123}}}{a_{12} a_{23}} (e + f) + \frac{a_1 a_2 + a_2 a_3 + a_2^2 - a_1 a_3}{a_{12} a_{23}} h, $$

$$ \tilde{e} = \frac{a_1 a_3}{a_{12} a_{23}} e - \frac{a_2 a_{123}}{a_{12} a_{23}} f - \frac{\sqrt{a_1 a_2 a_3 a_{123}}}{a_{12} a_{23}} h, $$

$$ \tilde{f} = - \frac{a_1 a_3}{a_{12} a_{23}} e + \frac{a_2 a_{123}}{a_{12} a_{23}} f - \frac{\sqrt{a_1 a_2 a_3 a_{123}}}{a_{12} a_{23}} h. $$

This set of equalities gives an automorphism of sl$_2$. We find the corresponding group element of SO$_2$ by solving the following matrix equation:

$$ (36) \qquad U \begin{pmatrix} \frac{\tilde{h}}{2} & f \\ e & -\frac{\tilde{h}}{2} \end{pmatrix} U^{-1} = \begin{pmatrix} \frac{\tilde{h}}{2} & \tilde{f} \\ \tilde{e} & -\frac{\tilde{h}}{2} \end{pmatrix}. $$

Solving for $U$ gives

$$ U_{(123)} = \begin{pmatrix} \sqrt{\frac{a_1 a_3}{a_{12} a_{23}}} & -\sqrt{\frac{a_2 a_{123}}{a_{12} a_{23}}} \\ \sqrt{\frac{a_2 a_{123}}{a_{12} a_{23}}} & \sqrt{\frac{a_1 a_3}{a_{12} a_{23}}} \end{pmatrix}. $$

A straightforward calculation shows that $U_{(123)}$ is indeed an orthogonal matrix. It represents a rotation with angle $\theta_{1,2,3}$ defined by

$$ \cos(\theta_{1,2,3}) = \sqrt{\frac{a_1 a_3}{a_{12} a_{23}}}. $$

**Example 4.5.** Consider in the algebra $\mathcal{R}_4(\mathfrak{h})$ the Lie algebra sl$_3$ constructed from the labelling Abelian algebra $\{Q_{12}, Q_{123}\}$.

$$ e_{12} := \lambda_1([Q_{12}, [Q_{12}, Q_{23}]] + a_{12}[Q_{12}, Q_{23}]), $$

$$ e_{21} := \lambda_1([Q_{12}, [Q_{12}, Q_{23}]] - a_{12}[Q_{12}, Q_{23}]), $$

$$ h_1 := [e_{12}, e_{21}], $$

$$ e_{23} := \lambda_2([Q_{123}, [Q_{123}, Q_{34}]] + a_{123}[Q_{123}, Q_{34}]), $$

$$ e_{32} := \lambda_2([Q_{123}, [Q_{123}, Q_{34}]] - a_{123}[Q_{123}, Q_{34}]), $$

$$ h_2 := [e_{23}, e_{32}], $$

$$ e_{13} := [e_{12}, e_{23}], $$

$$ e_{31} := [e_{32}, e_{21}]. $$

When performing the permutation $(1 \to 2 \to 3 \to 1)$ we find new generators. These generators are eigenvectors of the labelling Abelian algebra $\{Q_{23}, Q_{123}\}$. We will not explicitly express the new generators as linear combinations but we will give $U \in SL_3$. To find $U$ we solve the following equation

$$ (37) \qquad U \begin{pmatrix} \frac{2\tilde{h}_1 + \tilde{h}_2}{3} & e_{21} & e_{31} \\ e_{12} & \frac{-\tilde{h}_1 + \tilde{h}_2}{3} & e_{32} \\ e_{13} & e_{23} & \frac{-\tilde{h}_1 - 2\tilde{h}_2}{3} \end{pmatrix} U^{-1} = \begin{pmatrix} \frac{2\tilde{h}_1 + \tilde{h}_2}{3} & \tilde{e}_{21} & \tilde{e}_{31} \\ \tilde{e}_{12} & \frac{-\tilde{h}_1 + \tilde{h}_2}{3} & \tilde{e}_{32} \\ \tilde{e}_{13} & \tilde{e}_{23} & \frac{-\tilde{h}_1 - 2\tilde{h}_2}{3} \end{pmatrix}. $$
For the given permutation we find

$$ U_{(123)} = \begin{pmatrix} \sqrt{\frac{a_1 a_3}{a_{12} a_{23}}} & -\sqrt{\frac{a_2 a_{123}}{a_{12} a_{23}}} & 0 \\ \sqrt{\frac{a_2 a_{123}}{a_{12} a_{23}}} & \sqrt{\frac{a_1 a_3}{a_{12} a_{23}}} & 0 \\ 0 & 0 & 1 \end{pmatrix}. $$

The 2 × 2 in the upper right corner is exactly $U_{(123)}$ in the $n = 2$ case. If $U_{(123)}$ acts on a three dimensional space it represents a planar rotation over an angle $\theta_{1,2,3}$. We will denote this matrix alternatively by $R_{x_1 x_2}(\theta_{1,2,3}) := U_{(123)}$. The $x_1 x_2$ index represents the plane that is being rotated. This overlap can be represented by a swap as in Proposition 4.1:

<figure><img src="image_12.png" /></figure>

**Example 4.6.** Consider in the algebra $R_4(\mathfrak{h})$ again the Lie algebra $sl_3$ constructed from the labelling Abelian algebra $Y_1 := \{Q_{12}, Q_{123}\}$. The second labelling Abelian algebra we consider is $Y_2 := \{Q_{12}, Q_{34}\}$. As before we can represent this overlap by a swap as in Proposition 4.1.

<figure><img src="image_13.png" /></figure>

The generators corresponding to the second Abelian algebra are:

$$ \tilde{e}_{12} := \lambda_{12}([Q_{12}, [Q_{12}, Q_{234}]] + a_{12}[Q_{12}, Q_{234}]), $$

$$ \tilde{e}_{21} := \lambda_{12}([Q_{12}, [Q_{12}, Q_{234}]] - a_{12}[Q_{12}, Q_{234}]), $$

$$ \tilde{h}_1 := [e_{12}, e_{21}], $$

$$ \tilde{e}_{23} := \lambda_{34}([Q_{34}, [Q_{34}, Q_{124}]] + a_{34}[Q_{34}, Q_{124}]), $$

$$ \tilde{e}_{32} := \lambda_{34}([Q_{34}, [Q_{34}, Q_{124}]] - a_{34}[Q_{34}, Q_{124}]), $$

$$ \tilde{h}_2 := [e_{23}, e_{32}], $$

$$ \tilde{e}_{13} := [e_{12}, e_{23}], $$

$$ \tilde{e}_{31} := [e_{32}, e_{21}]. $$
Solving equation (37) results in

$$ U_{y_1 y_2} := \begin{pmatrix} 1 & 0 & 0 \\ 0 & \sqrt{\frac{a_3 a_{1234}}{a_{123} a_{34}}} & -\sqrt{\frac{a_4 a_{12}}{a_{123} a_{34}}} \\ 0 & \sqrt{\frac{a_4 a_{12}}{a_{123} a_{34}}} & \sqrt{\frac{a_3 a_{1234}}{a_{123} a_{34}}} \end{pmatrix}. $$

This is again an orthogonal matrix. Similarly to the previous example one can see that the lower-right $2 \times 2$ matrix is a rotation matrix. The matrix $U_{y_1 y_2}$ represents a planar rotation matrx over an angle $\pi/2 - \theta_{12,3,4}$ but in a different plane. We denote this matrix by

$$ R_{x_2 x_3} (\theta_{12,3,4} - \pi/2). $$

Given two labelling Abelian algebras, we are able to construct an isomorphism of $sl_3$ and its corresponding rotation matrix. The previous examples show that if one labelling Abelian algebra is obtained by a swap on the other, the resulting rotation matrices are planar. Twists on the other hand give reflections. This leads to an alternative way to construct the rotation matrix of an isomorphism related to two labelling Abelian algebras. We choose a path in the connection graph between the vertices related to the labelling Abelian algebras. Every step along the path in the connection graph can be represented by a swap as in Proposition 4.1. Each of these swaps leads to a planar rotation. The final rotation matrix will be the product of each planar rotation found along the path. For example, consider these two labelling Abelian algebras: $\langle Q_{12}, Q_{34} \rangle$ and $\langle Q_{13}, Q_{24} \rangle$. The intermediate bases are the following:

$$ \langle Q_{12}, Q_{34} \rangle \to \langle Q_{12}, Q_{123} \rangle \to \langle Q_{13}, Q_{123} \rangle \to \langle Q_{13}, Q_{24} \rangle. $$

For each step we construct the corresponding rotation matrix. This results in

$$ R := R_{x_2 x_3}(\pi/2 - \theta_{12,3,4}) R_{x_1 x_2}(-\theta_{1,2,3}) R_{x_2 x_3}(\theta_{13,2,4} - \pi/2). $$

It is known that every rotation matrix $R \in SO_3$ can be written as the following product of planar rotations:

$$ R_{x_2 x_3}(\theta_1) R_{x_1 x_2}(\theta_2) R_{x_2 x_3}(\theta_3). $$

The angles $\theta_1, \theta_2$ and $\theta_3$ are the so-called Euler angles. By choosing the right path through the connection graph it is possible to give the decomposition of the rotation matrix in planar matrices. The Euler angles also show up in the overlap coefficients between the given bases. The overlap coefficients (34) are multivariate Krawtchouk polynomials of Griffiths type depending on three parameters $R_h^{13,2,4}$, $R_h^{1,2,3}$ and $R_h^{4,3,12}$. The relationship between these parameters and the angles is the following:

$$ R_h^{K,L,M} = \cos(2\theta_{K,L,M}). $$

This relation between multivariate Krawtchouk polynomials and rotations had been discussed earlier in [24].

We generalize the previous analysis to any $n$. To any pair of labelling Abelian algebras we are able to construct overlap coefficients and an isomorphism of $sl_{n-1}$ represented by a rotation matrix in $SO_{n-1}$. Any rotation matrix in $SO_{n-1}$ can be written as the product of $\frac{(n-1)(n-2)}{2}$ planar rotations. In the same way the overlap coefficients are multivariable Krawtchouk polynomials of Griffiths type depending on $\frac{(n-1)(n-2)}{2}$ parameters. By Proposition 4.1 it takes at most $\frac{(n-1)(n-2)}{2}$ steps
through the connection graph. Each step provides us with a parameter and a planar rotation matrix leading to the right number of parameters and planar rotations. We will showcase this with an example that also shows the link with Krawtchouk polynomials of Tratnik type:

**Example 4.7.** Consider the labelling Abelian algebras

$$ \mathcal{Y}_{initial} := \{Q_{[k]} | 1 < k < n\}, \quad \mathcal{Y}_{final} := \{Q_{[k...n]} | 1 < k < n\}. $$

Between these two labelling Abelian algebras, it will take $\frac{(n-1)(n-2)}{2}$ steps to find the overlap.

**Step 1:** Replace $Q_{12}$ by $Q_{23}$ in $\mathcal{Y}_{initial}$. The overlap coefficients are univariate Krawtchouk polynomials:

$$ B_{s_1^1 s_1^2}(R_h^{1,2,3}, C^{1,2,3}). $$

The vector index $s^1$ labels the basis diagonalized by $\mathcal{Y}_{initial}$. The index $s^2$ labels the basis diagonalized by the new labelling Abelian algebra. The related rotation matrix is given by

$$ R_{x_1 x_2}(\theta_{1,2,3}). $$

**Step 2:** We perform two steps: Replace $Q_{123}$ by $Q_{234}$ and then $Q_{23}$ by $Q_{34}$. The overlap coefficients are bivariate Krawtchouk polynomials of Tratnik type:

$$ B_{s_1^2 s_1^3}(R_h^{1,23,4}, C^{1,23,4}) B_{s_2^2 s_2^3}(R_h^{2,3,4}, C^{2,3,4}). $$

The related rotation matrix is given by

$$ R_{x_1 x_2}(\theta_{2,3,4}) R_{x_2 x_3}(\theta_{1,23,4}). $$

**Step $k-1$:** Consider the first $k-1$ elements in the labelling Abelian algebra:

$$ \{Q_{k-1k}, Q_{[k-2..k]}, \dots, Q_{[2..k]}, Q_{[k]}\}. $$

We perform $k-1$ steps. Replace $Q_{[k]}$ by $Q_{[2..k+1]}$. Next replace $Q_{[2..k]}$ by $Q_{[3..k+1]}$ and so on until $Q_{k-1k}$ is replaced by $Q_{kk+1}$. These $k-1$ steps lead to overlap coefficients that are a product of $k-1$ univariate Krawtchouk polynomials. These are more specifically multivariate Krawtchouk polynomials of Tratnik type:

$$ \prod_{l=1}^{k-1} B_{s_l^{k-1} s_l^k}(R_h^{l, [l+1..k], k+1}, C^{l, [l+1..k], k+1}). $$

The related rotation matrix is given by

$$ \prod_{l=1}^{k-1} R_{x_{k-l} x_{k+1-l}}(\theta_{l, [l+1..k], k+1}). $$

Combining all $n-2$ steps we come to the conclusion that the overlap coefficients are given by

$$ B_{s^1 s^{n-1}} = \sum_{m=2}^{n-2} \sum_{s_l^m} \prod_{k=2}^{n-1} \prod_{l=1}^{k-1} B_{s_l^{k-1} s_l^k}(R_h^{l, [l+1..k], k+1}, C^{l, [l+1..k], k+1}). $$

The overlap coefficients are multivariate Krawtchouk polynomials of Griffiths type which are themselves a sum of products of Krawtchouk polynomials of Tratnik type. The overlap coefficients depend on $\frac{(n-1)(n-2)}{2}$ parameters of the form $R_h^{l, [l+1..k], k+1}$.
The corresponding rotation matrix is given by a product of $(\frac{n-1}{2})(\frac{n-2}{2})$ rotation matrices.

$$ \prod_{k=2}^{n-1} \prod_{l=1}^{k-1} R_{x_{k+1-l}x_{k-l}}(\theta_{l,[l+1..k],k+1}) $$

## 5. CONCLUSION

We introduced the oscillator Racah algebra $R_n(h)$. We have shown how to embed $sl_{n-1}$ into the algebra $R_n(h)$. This connects the representation theory for both algebras. In finite irreducible representations we considered bases diagonalized by labelling Abelian algebras of $R_n(h)$. The overlap coefficients between a pair of bases are shown to be multivariate Krawtchouk polynomials of Tratnik or Griffiths type. Isomorphisms of $sl_{n-1}$ related to pairs of labelling Abelian algebras and their corresponding Lie group elements were constructed and their link to the overlap coefficients is explained. This has provided an explanation as to why the recoupling coefficients of the oscillator algebra and the matrix elements of the restrictions to $O(n-1)$ of symmetric representations of $SL(n-1)$ are generically given in terms of the multivariate Krawtchouk polynomials of Griffiths.

## 6. ACKNOWLEDGEMENTS

NC is partially supported by Agence National de la Recherche Projet AHA ANR-18-CE40-0001 and is gratefully holding a CRM-Simons professorship. WVDV thanks the Fonds Professor Frans Wuytack for supporting his research. He is also grateful for the hospitality offered by him at the CRM during his stay. The research of LV is supported in part by a discovery grant of the Natural Science and Engineering Research Council (NSERC) of Canada.

## APPENDIX A. CALCULATION OF OVERLAP COEFFICIENTS

Let $V$ be a finite dimensional representation of $sl_2$ and $\tilde{h}$ an automorphism of $sl_2$. The element $h$ is a Cartan generator of $sl_2$. Let $\{\psi_k\}$ be an eigenbasis of $h$ and $\{\phi_s\}$ be an eigenbasis for $\tilde{h}$. The indices $k$ and $s$ run from 0 to $N$ with $\dim(V) = N+1$. We are interested in the overlap coefficients $B_{ks}$ between these bases:

$$ (38) \qquad \phi_s = \sum_{k=0}^{N} B_{sk} \psi_k. $$

$$ h\psi_k = \mu_k \psi_k \qquad \tilde{h}\phi_s = \nu_s \phi_s. $$

The algebra $sl_2$ has algebra relations $[h, e] = 2e$ and $[h, f] = -2f$ with $e$ is the raising operator and $f$ the lowering operator on $\{\psi_k\}$:

$$ e\psi_k = e_{kk+1}\psi_{k+1} \qquad f\psi_k = f_{kk-1}\psi_{k-1} $$

and $\mu_k = \mu_0 + 2k$. From the algebra relation $[e, f] = h$ it follows that

$$ f_{kk-1}e_{k-1k} - e_{kk+1}f_{k+1k} = \mu_k. $$

Let $A_k := e_{kk-1}f_{k-1k}$. Then we have

$$ A_k - A_{k+1} = 2k + \mu_0. $$

From this we find

$$ A_k = -k(k-1) - \mu_0 k - \Omega $$
with $\Omega \in \mathbb{R}$. We express $\tilde{h}$ as a linear combination of $h, e$ and $f$.

$$ \tilde{h} = R_h h + R_e e + R_f f $$

with $R_e R_f + R_h^2 = 1$. We have set up everything we need to find the overlap coefficients. Let the operator $\tilde{h}$ act on both sides of equality (38).

$$ \tilde{h}\psi_s = \sum_{k=0}^{N} B_{sk}(R_h h + R_e e + R_f f)\psi_k. $$

This gives

$$ \nu_s \phi_s = \sum_{k=0}^{N} B_{sk}(R_h \mu_k \psi_k + R_e e_{kk+1} \psi_{k+1} + R_f f_{kk-1} \psi_{k-1}). $$

We expand the left hand side into the basis $\psi_k$ and we gather the terms on the right hand side:

$$ \sum_{k=0}^{N} \nu_s B_{sk} \psi_k = \sum_{k=0}^{N} (B_{sk} R_h \mu_k + B_{sk-1} R_e e_{k-1k} + B_{sk+1} R_f f_{k+1k}) \psi_k. $$

From this we find the recurrence relation

$$ \nu_s B_{sk} = B_{sk+1} R_f f_{k+1k} + B_{sk} R_h \mu_k + B_{sk-1} R_e e_{k-1k}. $$

We want to recognize this recurrence relation as one of the family of orthogonal polynomials. Let

$$ \tilde{B}_{ks} = \left( \prod_{t=2}^{k} f_{tt-1} R_f \right) B_{sk} $$

to find

$$ \nu_s \tilde{B}_{sk} = \tilde{B}_{sk+1} + R_h \mu_k \tilde{B}_{sk} + R_e R_f e_{k-1k} f_{k-1k} \tilde{B}_{sk-1}. $$

We write the coefficients as polynomials in $x = \nu_s$:

$$ (39) \qquad x \tilde{B}_k(x) = \tilde{B}_{k+1}(x) + R_h \mu_k \tilde{B}_k(x) + R_e R_f A_k \tilde{B}_{k-1}(x). $$

We want to compare this with the recurrence relation of the normalized Krawtchouk polynomials as defined in [35]:

$$ xp_n(x) = p_{n+1} + (n(1-2r) + rN)p_n(x) + r(1-r)n(N+1-n)p_{n-1}(x) $$

with $n = 0, 1, \dots, N$. Let $x = \alpha y + \beta$ and introduce $q_n(y) = p_n(\alpha y + \beta)/\alpha^n$. The polynomial $q_n(x)$ satisfies the following recurrence relation:

$$ yq_n(y) = q_{n+1} + \frac{n(1-2r)+rN-\beta}{\alpha}q_n(y) + \frac{r(1-r)n(N+1-n)}{\alpha^2}q_{n-1}(y). $$

We retrieve equation (39) if we set

$$ \alpha = \frac{1}{2}, \quad r = \frac{1-R_h}{2}, \quad \beta = \frac{N}{2}, \quad k = n, \quad \Omega = 0, \quad \mu_0 = -N. $$
We explicitly write down the polynomials $\tilde{B}_k(x)$.

$$ \begin{aligned} \tilde{B}_k(x) &= 2^k p_k \left( \frac{x+N}{2} \right) \\ &= (-N)_k (1-R_h)^k K_k \left( \frac{x+N}{2}; \frac{1-R_h}{2}, N \right) \\ &= (-N)_k (1-R_h)^k {}_2F_1 \left( -k, -\frac{x+N}{N^2} | \frac{2}{1-R_h} \right). \end{aligned} $$

The overlap coefficients are the Krawtchouk polynomials $K_k(x)$ (defined in the last line of the equation above) up to a normalization factor.

$$ B_{sk} = \left( \prod_{t=2}^{k} f_{tt-1} R_f \right) (-N)_k (1-R_h)^k K_k \left( \frac{\nu_s + N}{2}; \frac{1-R_h}{2}, N \right). $$

REFERENCES

[1] P. Baseilhac. Deformed Dolan-Grady relations in quantum integrable models. *Nucl.Phys. B* 709; 491-521, 2005.

[2] P. Baseilhac and N. Crampe. FRT presentation of classical Askey-Wilson algebras. *Lett. Math. Phys.* 1-21,2019.

[3] D. Bullock and J. H. Przytycki. Multiplicative structure of Kauffman bracket skein module quantization. *Proc. Amer. Math. Soc.* 128; 923-931, 2000.

[4] J. Cooke. Kaufman Skein Algebras and Quantum Teichmüller Spaces via Factorisation Homology. arXiv: 1811.09293, 2018.

[5] N. Crampe, L. Frappat and L. Vinet. Centralizers of the superalgebra osp(1|2): the Brauer algebra as a quotient of the Bannai-Ito algebra. *J. Phys. A* (to appear), and arXiv:1906.03936.

[6] N. Crampe, L. Poulain d'Andecy and L. Vinet. Temperley-Lieb, Brauer and Racah algebras and other centralizers of su(2). arXiv:1905.06346.

[7] N. Crampe, E. Ragoucy, L. Vinet, and A. S. Zhedanov. Truncation of the reflection algebra and the Hahn algebra. *J.Phys. A: Math. Theor.* 52, 35LT01, 2019.

[8] N. Crampe, J. Gaboriaud, L. Vinet and M. Zaimi. Revisiting the Askey-Wilson algebra with the universal R-matrix of $U_q(su(2))$. arXiv: 1908.04806, 2019.

[9] N. Crampe, L. Vinet and M. Zaimi. Bannai-Ito algebras and the universal R-matrix of osp(1|2). to appear.

[10] H. De Bie, H. De Clerq, and W. van de Vijver. The higher rank q-deformed Bannai-Ito and Askey-Wilson algebra. *Commun. Math. Phys.*, 2019.

[11] H. De Bie, V. X. Genest, S. Tsujimoto, L. Vinet, and A. Zhedanov. The Bannai-Ito algebra and some applications. *J. Phys.: Conf. Ser.*, 597:012001, 2015.

[12] H. De Bie, V.X. Genest, L. Vinet, W. van de Vijver. A higher rank Racah algebra and the $(\mathbb{Z}_2)^n$ Laplace-Dunkl operator. *J. Phys. A: Math. Theor.* 51 025203 (20pp), 2018.

[13] H. De Bie, V. X. Genest, L. Vinet, The $\mathbb{Z}_2^n$ Dirac-Dunkl operator and a higher rank Bannai-Ito algebra. *Adv. Math.* 303: 390-414, 2016.

[14] P. Diaconis and R. Griffiths. An introduction to multivariate Krawtchouk polynomials and their applications. *Journal Statistical Planning and Inference*, 154: 39-53, 2014.

[15] L. Dolan and M. Grady. Conserved charges from self-duality. *Phys. Rev. D* **25** 1587?1604, 1982.

[16] L. Frappat, J. Gaboriaud, E. Ragoucy and L. Vinet. The dual pair $(U_q(su(1, 1), o_q^{1/2}(2n)))$, q-oscillators and Askey-Wilson algebras. arXiv: 1908.04277, 2019.

[17] J. Gaboriaud, L. Vinet, S. Vinet and A. Zhedanov. The Racah algebra as a commutant and Howe duality. *J. Phys. A: Math. Theor.* 51 50LT01, 2018.

[18] J. Gaboriaud, L. Vinet, S. Vinet and A. Zhedanov. The dual pair $\mathrm{Pin}(2n) \times \mathrm{osp}(1|2)$, the Dirac equation and the Bannai-Ito algebra. *Nucl. Phys. B.* 937, 226-239, 2018.

[19] V. X. Genest, L. Vinet, and A. Zhedanov. Superintegrability in two dimensions and the Racah-Wilson algebra. *Lett. Math. Phys*, 104: 931-952, 2014.

[20] V. X. Genest, L. Vinet, and A. Zhedanov. The Bannai-Ito polynomials as Racah coefficients of the $sl_{-1}(2)$ algebra. *Proc. Am. Math. Soc*, 142: 1545-1560, 2014.
[21] V. X. Genest, L. Vinet, and A. Zhedanov. The equitable Racah algebra from three $su(1, 1)$ algebras. *J. Phys. A: Math. Theor.*, 47:025203, 2014.

[22] V. X. Genest, L. Vinet, and A. Zhedanov. The Racah algebra and superintegrable models. *J. Phys.:Conf. Ser.*, 512 :012011, 2014.

[23] V. X. Genest, L. Vinet, and A. Zhedanov. The non-symmetric Wilson polynomials are the Bannai-Ito polynomials. *Proc. Amer. Math. Soc.*, 144 :5217-5226, 2015.

[24] V. X. Genest, L. Vinet, and A. Zhedanov. The multivariate Krawtchouk polynomials as matrix elements of the rotation group representations on oscillator states. *J. Phys. A: Math. Theor.*, 46 :505203, 2013.

[25] J. Geronimo and P. Iliev. Bispectrality of multivariable Racah-Wilson polynomials. *Constr. Approx.* 31, 417-457, 2010.

[26] Ya. A. Granovskii and A. S. Zhedanov. Nature of the symmetry group of the 6j-symbol. *Sov. Phys. JETP*, 67:1982-1985, 1988.

[27] Ya. A. Granovskii, and A.S. Zhedanov. Hidden Symmetry of the Racah and Clebsch-Gordan Problems for the Quantum Algebra $sl_q(2)$. arXiv:hep-th/9304138, 1993.

[28] R. Griffiths. Orthogonal polynomials on the multinomial distribution. *Australian journal of Statistics*, 13 (1), 27-25, 1971.

[29] M. R. Hoare and M. Rahman. A probabilistic origin for a new class of bivariate polynomials. *SIGMA*, 4:89-106, 2008.

[30] H.-W. Huang. An embedding of the universal Askey-Wilson algebra into $U_q(sI_2)) \otimes U_q(sI_2)) \otimes U_q(sI_2))$. *Nucl. Phys.B*, 922, 401-434 (2017).

[31] P. Iliev and P. Terwilliger. The Rahman polynomials and the Lie algebra $sl_3(C)$. *Trans. Amer. Math. Soc.* 364, no. 8, 4225-4328, 2012.

[32] P. Iliev. A Lie theoretic interpretation of multivariate hypergeometric. *Compositio Math.*, 148, 991-1002, 2012.

[33] A. U. Klimyk and N. J. Vilenkin. Representation of Lie groups and special functions. *Kluwer Academic*, 1991-1993, Dordrecht.

[34] R. Koekoek and R. F. Swarttouw. "The Askey-scheme of hypergeometric orthogonal polynomials and its q-analogue," Tech. Rep. 98-17, Delft University of Technology, Faculty of Information Technology and Systems, Department of Technical Mathematics and Informationics, 1998.

[35] R. Koekoek, P. A. Lesky, and R. F. Swarttouw. Hypergeometric Orthogonal Polynomials and Their $q$-Analogues. Springer, 2010.

[36] T. Koornwinder. Krawtchouk Polynomials, a Unification of Two Different Group Theoretic Interpretations. *SIAM J. Math. Analysis*, Vol. 13, no.6, 1011-1023, 1982

[37] T. Koornwinder. The Relationship between Zhedanov's Algebra AW(3) and the Double Affine Hecke Algebra in the Rank One Case. SIGMA, 3, 063, 2007.

[38] T. Koornwinder. Zhedanov's Algebra AW(3) and the Double Affine Hecke Algebra in the Rank One Case. II. The Spherical Subalgebra. SIGMA, 4, 052, 2008.

[39] H. Mizukawa and H. Tanaka. $(n+1, m+1)$-hypergeometric functions associated to character algebras. *Proc. Amer. Math. Soc.*, 132: 2613-2618, 2004.

[40] S. Post, and A. Walter, A higher rank extension of the Askey-Wilson Algebra. arXiv:1705.01860.

[41] A. V. Rosenblyum. Spectral analysis of generators of representations of the group $U(3)$. *Theoretical and Mathematical Physics*, 73: 1352-1356, 1987.

[42] R. F. Streater. The representations of the oscillator group. *Commun. Math. Phys.* 4: 217, 1967.

[43] P. Terwilliger. The Universal Askey-Wilson Algebra and DAHA of Type $(C_1^\vee, C_1)$. *SIGMA*, 9:047, 2013.

[44] P. Terwilliger and R. Vidunas. Loonard pairs and the Askey-Wilson relations. J. Algebra Appl. 03; 411-426, 2004.

[45] P. Terwilliger and A. Zitnik. Distance-regular graphs of $q$-Racah type and the universal Askey-Wilson algebra. *J. Combinatorial Theory, Series A*, 125 :98-112, 2014.

[46] M. V. Tratnik. Some multivariable orthogonal polynomials of the Askey tableau-discrete families. *J. Math. Phys.*, 32:2337-2342, 1991.

[47] J. Van der Jeugt. Coupling coefficients for Lie algebra representations and addition formulas for special functions. *J. of Math. Phys.* 38, 2728,1997.
[48] J. Van der Jeugt. $3nj$-coefficients and Orthogonal Polynomials of Hypergeometric type. *Orthogonal polynomials and special functions*, 1817, p.25-92, 2003.

[49] A.S. Zhedanov. Hidden symmetry of the Askey-Wilson polynomials. *Theor. Math. Phys.* 89:1146-1157, 1991.

[50] A. Zhedanov. $9j$-symbols for the oscillator algebra and the Krawtchouk polynomials in two variables. *J. Phys. A: Math. Gen.* **30** 8337, 1997.

INSTITUT DENIS-POISSON CNRS/UMR 7013 - UNIVERSITÉ DE TOURS - UNIVERSITÉ D'ORLÉANS, PARC DE GRAMMONT, 37200 TOURS, FRANCE.

*E-mail address: crampe1977@gmail.com*

DEPARTMENT OF ELECTRONICS AND INFORMATION SYSTEMS, FACULTY OF ENGINEERING AND ARCHITECTURE, GHENT UNIVERSITY, BUILDING S8, KRIJGSLAAN 281, 9000 GHENT, BELGIUM

*E-mail address: Wouter.vandeVijver@UGent.be*

CENTRE DE RECHERCHES MATHÉMATIQUES, UNIVERSITÉ DE MONTRÉAL, P.O. BOX 6128, CENTRE-VILLE STATION, MONTRÉAL, QC H3C 3J7, CANADA

*E-mail address: vinet@crm.umontreal.ca*