Fusion Rules for Holomorphic Discrete Series

• Holomorphic discrete series are unitary irreducible representations with lowest weight k > 1/2, realized in weighted Bergman spaces.
• Fusion rules decompose tensor products as D⁺(k1)⊗D⁺(k2) = ⊕ₘ D⁺(k1+k2+m) with each summand appearing once, crucial for quantum and functional analysis.
• Explicit intertwining operators extract fusion components, linking representation theory to orthogonal polynomials, special functions, and quantum channels.

Holomorphic discrete series representations are fundamental objects in the representation theory of non-compact semisimple Lie groups, notably SU(1,1) and its relatives. The set of rules governing the decomposition of tensor products of such representations—commonly referred to as fusion rules—dictates how composite representation spaces break up into irreducible constituents. These rules play a key role in mathematical physics, particularly in the analysis of systems with SU(1,1) symmetry, quantum optics, and in the theory of special functions.

1. Structure of Holomorphic Discrete Series for SU(1,1)

A holomorphic discrete series representation of SU(1,1) is a unitary irreducible representation (UIR) denoted $D^+(k)$, labeled by the lowest weight $k > \tfrac{1}{2}$, with $k$ typically taking values in $\{1, \tfrac{3}{2}, 2, \tfrac{5}{2}, \dots\}$ for SU(1,1) (Gazeau et al., 4 Apr 2025). These series can also be realized as weighted Bergman spaces on the unit disk in appropriate function-theoretic models (Haastrecht, 2024).

The representations exhibit rich orthogonality properties and their characters, when restricted to the maximal compact subgroup U(1), possess a tractable and explicit form: $$\chi_k\bigl(e^{i\theta\sigma_3/2}\bigr) = \frac{e^{-i(2k-1)\theta/2}}{2i\sin(\theta/2)}$$ These character formulas underpin the derivation of fusion rules and ensure the orthogonality of distinct holomorphic discrete series components.

2. Tensor Product Decomposition and Fusion Rules

The tensor product decomposition of two positive holomorphic discrete series representations $D^+(k_1)\otimes D^+(k_2)$ follows a closed, universal fusion rule: $$D^+(k_1) \otimes D^+(k_2) = \bigoplus_{m=0}^\infty D^+(k_1 + k_2 + m)$$ A necessary and sufficient condition for $D^+(k_3)$ to appear in the decomposition is $k_3 - (k_1 + k_2) \in \mathbb{N}_0$, subject to the constraint $k_j > \tfrac{1}{2}$ for all $j$ (Gazeau et al., 4 Apr 2025, Haastrecht, 2024). Each allowed summand appears with multiplicity one.

An equivalent algebraic formulation expresses the decomposition as: $$D^+(k_1)\otimes D^+(k_2) = \sum_{k_3 \ge k_1+k_2,\, k_3-k_1-k_2\in \mathbb{N}_0} D^+(k_3)$$ This closed form, confirmed by explicit character expansions and Fourier analysis, is a cornerstone for analyzing composite quantum systems invariant under SU(1,1).

Input Representations	Allowed Output $D^+(k_3)$	Multiplicity
$D^+(k_1)\otimes D^+(k_2)$	$k_3 = k_1+k_2+m$, $m=0,1,2,\dots$	$1$

3. Generalizations: Higher Tensor Powers and Stratified Model

Fusion rules generalize to higher ($n$-fold) tensor products of holomorphic discrete series. For the scalar discrete series of universal covering $\widetilde{SL_2(\mathbb{R})}$, the diagonal restriction of the outer tensor product $$\pi_{\lambda_1}\boxtimes\dots\boxtimes\pi_{\lambda_n}$$ decomposes as (Labriet, 2022): $$\bigotimes_{i=1}^n \pi_{\lambda_i} \simeq \bigoplus_{k=0}^\infty m_n(k;\Lambda) \cdot \pi_{|\Lambda| + 2k}$$ where $|\Lambda| = \sum_{i=1}^n \lambda_i$ and the multiplicities are given by binomial coefficients: $m_n(k;\Lambda) = \binom{n + k - 2}{n-2}$ This enhances the one-dimensional multiplicity seen in rank-one cases and reflects the structure of invariant orthogonal polynomials on the “stratification space” $D_{n-1}$ associated with the tensor product (Labriet, 2022). For $n=2$ the multiplicity reduces to one, in agreement with SU(1,1) results above.

The stratified model provides a geometric functional model for these decompositions, relating representation-theoretic branching to orthogonal polynomials on a simplex.

4. Explicit Construction of Intertwining Operators

The decomposition is explicitly implemented through a sequence of intertwining operators. For $SU(1,1)$, orthogonal projections $P_n$ (sometimes called “transvectants”) extract the $D^+(k_1 + k_2 + n)$ summand from $D^+(k_1)\otimes D^+(k_2)$ (Haastrecht, 2024): $$P_n F(\zeta) = C_{k_1,k_2,n} \sum_{j=0}^n (-1)^j \frac{(k_1)_j (k_2)_{n-j}}{j!(n-j)!} \left. \partial_z^j \partial_w^{n-j} F(z,w)\right|_{z=w=\zeta}$$ The normalizing constant $C_{k_1, k_2, n}$ is determined by the isometry condition for projections. This operator is both an explicit intertwiner and serves as a Kraus operator for defining quantum channels associated to the components.

For higher tensor powers, symmetry breaking and holographic intertwiners are constructed from orthogonal projections onto the spaces of orthogonal polynomials $Pol_k(D_{n-1})$ on the simplex—showing a deep link between representation theory and multivariate special functions (Labriet, 2022).

5. Quantum Channels and Asymptotic Behavior

Associated to each fusion component is a canonical equivariant quantum channel, implemented by the projection $P_n$. For $A$ a bounded operator on $\mathcal{H}_{\ell}$, the quantum channel to the $n$-th summand is: $\Phi_n(A) = P_n (A \otimes I_m) P_n^*$

$\Phi_n$ is completely positive, weak-$*$ continuous, and trace-preserving up to a known scaling factor: $$\operatorname{Tr}[\Phi_n(A)] = \frac{\ell+m+2n-1}{\ell-1} \operatorname{Tr}(A)$$

As one of the weights (e.g., $m$) tends to infinity, the quantum trace of functional calculus on fusion components converges to classical Berezin–Toeplitz quantization on the unit disk, described via generalized Husimi functions and combinations of Berezin transforms (Haastrecht, 2024): $$\lim_{m\to\infty} \operatorname{Tr}[\psi(\Phi_n(A))] = \int_D \psi(H^{(\ell)}(A)(z))\, d\iota(z)$$ This establishes a connection between the representation-theoretic fusion rules and the semiclassical quantization of phase space.

6. Orthogonality, Completeness, and Character Formulas

Key orthogonality relations and character identities underpin the fusion rules:

• Character orthogonality for SU(1,1) holomorphic discrete series:

$$\frac{1}{4\pi}\int_0^{4\pi} \sin^2\frac{\theta}{2}\; \chi_k(e^{i\theta\sigma_3/2})\;\overline{\chi_{k'}(e^{i\theta\sigma_3/2})}\, d\theta = \delta_{kk'}$$

• Character expansions yield the component structure:

$$\chi_{k_1}(h)\chi_{k_2}(h) = \sum_{m=0}^\infty \chi_{k_1 + k_2 + m}(h)$$

so that only $k_3 = k_1 + k_2 + m$, $m \in \mathbb{N}_0$, emerge in the product decomposition.

For higher tensor powers, completeness and orthogonality of the polynomial bases (e.g., via Jacobi or Gegenbauer systems) guarantee the distinctness and independence of all constructed fusion components (Labriet, 2022).

7. Connections to Branching Laws and Special Functions

The branching and fusion rules for holomorphic discrete series are closely related to the explicit theory of symmetry breaking (restriction to subgroups) and the structure of orthogonal polynomials on symmetric cones and simplices (Labriet, 2022). In higher rank, the decomposition is governed by the combinatorics of orthogonal polynomials, and the explicit functional-analytic realizations of intertwiners link representation theory to advanced special function theory.

This framework not only underpins the explicit spectral analysis of composite systems with SU(1,1) symmetry, but also informs the study of Rankin–Cohen type operators, quantum channels, and geometric quantization.

---

Holomorphic discrete series fusion rules for SU(1,1) and its extensions are thus governed by explicit closed formulas, underpinning a wide range of algebraic, analytic, and geometric structures in representation theory and mathematical physics (Gazeau et al., 4 Apr 2025, Haastrecht, 2024, Labriet, 2022).