\section{The Kauffman skein module}

The \emph{Kauffman skein module} $\mathcal S^A(M)$  of an oriented 3--dimensional manifold $M$ depends on a parameter $A=\mathrm{e}^{\pi\mathrm{i} \hbar}\in \mathbb C-\{0\}$, and is defined as follows: one first considers the vector space freely generated by  all isotopy classes of framed links in the thickened surface $S \times [0,1]$, and then one takes the quotient of this space by two relations:
\begin{itemize}
\item the first relation is the  \emph{skein relation} that
$
[L_1] = A^{-1} [L_0] + A [L_\infty]
$
whenever the three links $L_1$, $L_0$ and $L_\infty\subset S\times [0,1]$ differ only in a little ball where they are as represented on Figure~\ref{fig:SkeinRelation};

\item the second relation states that $[L\cup O] = -(A^2 +A^{-2})[L]$ whenever the knot $O$ is the boundary of a disk $D$ endowed with framing transverse to $D$, and the framed link $L$ is disjoint from the disk $D$.

\end{itemize}

\begin{figure}[htbp]

\SetLabels
( .5 * -.4 ) $L_0$ \\
( .1 * -.4 )  $L_1$\\
(  .9*  -.4) $L_\infty$ \\
\endSetLabels
\centerline{\AffixLabels{\includegraphics{SkeinRelation.pdf}}}
\vskip 15pt
\caption{A Kauffman triple}
\label{fig:SkeinRelation}
\end{figure}

In the special case where $M=S\times [0,1]$ for an oriented surface $S$, we write $\mathcal S^A(S) = \mathcal S^A \bigl (S\times [0,1]\bigr)$ and we note that this module now comes with a natural multiplication. Indeed, if $[L_1]$, $[L_2]\in \mathcal S^A(S)$ are respectively represented by the framed links $L_1$, $L_2$, we can consider their superposition
$$
[L_1] \cdot [L_2] = [L_1' \cup L_2'] \in \mathcal S^A(S)
$$
 represented by the union of the framed link $L_1'\subset S\times[0, \frac12]$ obtained by rescaling $L_1 \subset S\times [0,1]$ and of the framed link $L_2'\subset S \times [\frac12, 1]$ obtained by rescaling $L_2 \subset S\times [0,1]$. This endows the skein module $\mathcal S^A(S)$ with the structure of an algebra, called the  \emph{Kauffman skein algebra} of the oriented surface $S$.

\section{The Witten-Reshetikhin-Turaev topological quantum field theory}
\label{sec:WRTtqft}

We briefly review a few fundamental properties of the Witten-Reshetikhin-Turaev   topological quantum field theory, and refer to
\cite{BHMV3man, BHMVTQFT, LickorishSkein3manInvs, LickorishBook, LickorishHandbook,TuraevBook} for  details and proofs.

The Witten-Reshetikhin-Turaev   topological quantum field theory $Z^A_{\mathrm{WRT}}$ depends on the choice of a primitive $2N$--root of unity $A$, and is defined over the category $\mathcal C$ defined as follows: the objects of $\mathcal C$ are closed oriented surfaces $S$; the morphisms from $S_1$ to $S_2$ are pairs $(M,L)$ where $M$ is a compact oriented  3--manifold with $\partial M = (-S_1) \sqcup S_2$, where $L$ is a framed link in the interior of $M$, and where $M$ is endowed with a $p_1$--structure. The precise definition of a $p_1$--structure can be found in \cite[App.~B]{BHMVTQFT} but, for the purpose of the current article, it suffices to know that it captures  certain homotopic information on the tangent bundle of the manifolds considered.

In particular, $Z^A_{\mathrm{WRT}}$ associates a finite-dimensional vector space $V_S= Z^A_{\mathrm{WRT}}(S)$  to each closed oriented surface $S$, and a linear map $Z^A_{\mathrm{WRT}}{(M,L)} \colon V_{S_1} \to V_{S_2}$ to each morphism $(M,L)$ as above. In addition, the vector space $V_\varnothing$ associated to the empty surface $\varnothing$ comes equipped with a canonical identification with $\mathbb C$; in other words, this vector space is 1--dimensional and contains a preferred basis element that we will denote by 1.

The topological quantum field theory $Z^A_{\mathrm{WRT}}$ satisfies many  properties, in particular those that characterize topological quantum field theories. Most of these features will play no direct  r\^ole in the current article. However the following fact, which is grounded in properties of the quantum group $\mathrm U_q(\mathfrak{sl}_2)$  underlying the construction of $Z^A_{\mathrm{WRT}}$, is crucial for our purposes.

\begin{lem}
\label{lem:WRTmorphSkeinRel}
The linear maps $Z^A_{\mathrm{WRT}}{(M,L)}$ associated to the morphisms of the category $\mathcal C$ satisfy the skein relation that, as linear maps $ V_{S_1} \to V_{S_2}$,
$$
Z^A_{\mathrm{WRT}}{(M,L_1)} =   A^{-1} Z^A_{\mathrm{WRT}}{(M,L_0)} + A \, Z^A_{\mathrm{WRT}}{(M,L_\infty)}
$$
whenever the framed links $L_1$, $L_0$ and $L_\infty$ form a Kauffman triple in the manifold $M$ with $\partial M = (-S_1) \sqcup S_2$, for a fixed $p_1$--structure on $M$.

Also,
$$
Z^A_{\mathrm{WRT}}(M, L \cup O) = -(A^2 + A^{-2}) Z^A_{\mathrm{WRT}}(M,L)
$$
whenever the knot $O$ is the boundary of a disk $D$ endowed with framing transverse to $D$, and the framed link $L$ is disjoint from the disk $D$.    \qed
\end{lem}

A first consequence of Lemma~\ref{lem:WRTmorphSkeinRel} is that the skein algebra $\mathcal S^A(S)$ acts on the space $V_S$, by considering the special case $M= S\times [0,1]$. Indeed, for any framed link $L \subset S \times[0,1]$, the pair $(S \times [0,1] ,L)$ can be seen as a morphism from $S$ to $S$, and therefore induces a linear map $Z^A_{\mathrm{WRT}}(S\times [0,1], L) \colon V_S \to V_S$.

\begin{lem} \label{Lemma:rho}
\label{lem:WRTrepExists}
There exists a unique algebra homomorphism
$$
\rho\colon \mathcal S^A(S) \to \mathrm{End}(V_S)
$$
such that, for every framed link $L$ in $S\times [0,1]$,
$$
\rho \bigl( [L] \bigr) = Z^A_{\mathrm{WRT}}(S\times [0,1], L)
$$
when $S\times [0,1]$ is endowed with the product $p_1$--structure.
\end{lem}
\begin{proof}
Lemma~\ref{lem:WRTmorphSkeinRel} shows that the rule $L \mapsto Z^A_{\mathrm{WRT}}(S\times [0,1], L)$ is compatible with the skein relation, and therefore induces a linear map $\rho\colon \mathcal S^A(S) \to \mathrm{End}(V_S)$.  The multiplication law of $\mathcal S^A(S)$ is defined by the superposition operation, which itself corresponds to the composition of morphisms $(S\times [0,1], L)$ in the category $\mathcal C$. It follows  that $\rho$ is an algebra homomorphism.
\end{proof}

This homomorphism
$
\rho\colon \mathcal S^A(S) \to \mathrm{End}(V_S)
$
is the \emph{Witten-Reshetikhin-Turaev representation} of the skein algebra $\mathcal S^A(S)$.

Another application of Lemma~\ref{lem:WRTmorphSkeinRel} will enable us to perform computations in the space $V_S$.
Consider a 3-manifold $M$ with boundary $\partial M =S$, endowed with a $p_1$--structure. A framed link $L \subset M$ provides a morphism  $(M,L)$ from the empty surface $\varnothing$ to $S$. This provides a linear map $Z^A_{\mathrm{WRT}}(M,L)$ from $V_\varnothing=\mathbb C$ to $V_S$, and in particular specifies an element $Z^A_{\mathrm{WRT}}(M,L)(1) \in V_S$. As above, Lemma~\ref{lem:WRTmorphSkeinRel} shows that the map $L \mapsto Z^A_{\mathrm{WRT}}(M,L)(1)$ defines a linear map $\Phi_M \colon \mathcal S^A(M) \to V_S$, from the skein module of $M$ to the space $V_S$.

\begin{lem} \label{Lemma:Phi}
\label{lem:WRTspaceSpannedBySkeins}
If $M$ is an oriented  $3$--manifold with boundary $\partial M = S$, the linear map $\Phi_M \colon \mathcal S^A(M) \to V_S$ defined by $\Phi_M\bigl([L]\bigr)=Z^A_{\mathrm{WRT}}(M,L)(1)$ is surjective. \qed
\end{lem}

We will use Lemma~\ref{lem:WRTspaceSpannedBySkeins} in the special case where $S$ is connected and where $M$ is a handlebody $H$ with boundary $\partial H = S$. Choose an identification $H\cong\Sigma \times [0,1]$ of this handlebody with the product of the interval with  a compact oriented surface  $\Sigma$ with  boundary.
Select also a trivalent spine for $\Sigma$, namely a trivalent graph $\Gamma$ embedded in the interior of  $\Sigma$ such that $\Sigma$ deformation retracts to $\Gamma$.  We will use this data to describe a basis for~$V_S$.

An \emph{$N$--admissible weight system} assigns to each edge $e$ of $\Gamma$ a non-negative  integer weight $w(e)$ such that the following conditions hold:
\begin{enumerate}

\item at every vertex of $\Gamma$, the weights  of the edges $e_1$, $e_2$, $e_3$ adjacent to this vertex satisfy the triangle inequalities $w(e_1) \leq w(e_2) +w(e_3) $, $w(e_2) \leq w(e_1) +w(e_3) $ and $w(e_3) \leq w(e_1) +w(e_2) $;

\item if $N$ is odd, the weight $w(e)$ of each edge $e$ is even and bounded by $N-2$; in addition, at each vertex, the sum of the weights of the adjacent edges is bounded by $2N-4$;

\item if $N$ is even, the weight $w(e)$ of each edge $e$ is bounded by $\frac N2 -2$; in addition, at each vertex, the sum of the weights of the adjacent edges is even and bounded by $N-4$;

\end{enumerate}
Let $\mathcal W_\Gamma$ denote the (finite) set of all $N$--admissible weight systems for $\Gamma$.

An $N$--admissible edge weight system $w\in \mathcal W_\Gamma$ specifies  an element $\beta_w$ of the skein module $\mathcal S^A(H)$, by replacing  each edge $e$ of $\Gamma$ weighted by $w(e) \leq N-2$ by a copy of the $w(e)$--th Jones-Wenzl  idempotent, represented by a box \SetLabels
\mathrm{e}( .5* .5) $w(e) $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-5pt}{\AffixLabels{\includegraphics{JW.pdf}}}.
 The precise definition of the element $\beta_w\in \mathcal S^A(H)$ associated to the weight system $w\in \mathcal W_\Gamma$ can be found for instance in \cite{LickorishBook} or \cite{TuraevBook}, but we can give a flavor of the construction. For $a\geq 0$ such that $A^{4k}\neq 1$ for every $k$ with $0<k<a$ (namely, for $0\leq a <N$ is $N$ is odd, and for $0 \leq a<  \frac N2$ if $N$ is even), the  Jones-Wenzl  idempotent  \SetLabels
\mathrm{e}( .5* .5) $a $ \\
\endSetLabels
\raisebox{-5pt}{\AffixLabels{\includegraphics{JW.pdf}}}
is a certain formal linear combination of families of disjoint arcs, each with $a$ strands emanating from each end of the box lying on $\Sigma$. For instance,

\begin{align*}\textstyle
\SetLabels
\mathrm{e}( .5* .5) 4 \\
\endSetLabels
\raisebox{-5pt}{\AffixLabels{\includegraphics{JW.pdf}}}
&=\textstyle
\raisebox{-7pt}{\includegraphics{JW40.pdf}}
+
\frac{A^4+A^{-2}}{A^4+A^{-4}}\,
\raisebox{-7pt}{\includegraphics{JW42.pdf}}
+
\frac{A^4+1+A^{-4}}{A^6+A^2+A^{-2}+A^{-6}}\,
\raisebox{14pt}{\rotatebox{180}{\includegraphics{JW41.pdf}}}\\
&\quad\textstyle
+
\frac{A^4+1+A^{-4}}{A^6+A^2+A^{-2}+A^{-6}}\,
\raisebox{-7pt}{\includegraphics{JW41.pdf}}
+
\frac1{A^4+A^{-4}}\,
\raisebox{-7pt}{\includegraphics{JW43.pdf}}\\
&\quad\textstyle
+
\frac1{A^4+A^{-4}}\,
\raisebox{-7pt}{\reflectbox{\includegraphics{JW43.pdf}}}
+
\frac1{A^4+A^{-4}}\,
\raisebox{14pt}{\rotatebox{180}{\includegraphics{JW43.pdf}}}
+
\frac1{A^4+A^{-4}}\,
\raisebox{14pt}{\rotatebox{180}{\reflectbox{\includegraphics{JW43.pdf}}}}\\
&\quad\textstyle
+
\frac1{A^6+A^2+A^{-2}+A^{-6}}\,
\raisebox{-7pt}{\includegraphics{JW46.pdf}}
+
\frac1{A^6+A^2+A^{-2}+A^{-6}}\,
\raisebox{-7pt}{\reflectbox{\includegraphics{JW46.pdf}}}\\
&\quad\textstyle
+
\frac{A^2+A^{-2}}{(A^4+A^{-4})(A^4+1+A^{-4})}\,
\raisebox{-7pt}{\includegraphics{JW44.pdf}}
+
\frac{A^2+A^{-2}}{(A^4+A^{-4})(A^4+1+A^{-4})}\,
\raisebox{-7pt}{\reflectbox{\includegraphics{JW44.pdf}}}\\
&\quad\textstyle
+
\frac1{(A^4+A^{-4})(A^4+1+A^{-4})}\,
\raisebox{-7pt}{\includegraphics{JW45.pdf}}
+
\frac{(A^2+A^{-2})^2}{(A^4+A^{-4})(A^4+1+A^{-4})}\,
\raisebox{-7pt}{\includegraphics{JW47.pdf}}
\end{align*}

\noindent These strands are then connected by disjoint arcs near the vertices of $\Gamma$, using the fact that at each vertex the $w(e)$  add up to an even number and satisfy the triangle inequalities. See Figure~\ref{fig:vertex}, where $w(e_1) =6$, $w(e_2)=4$ and $w(e_3) =4$.

\begin{figure}[htbp]

\SetLabels
\mathrm{e}( .26* .5) $w(e_1) $ \\
( .84* .67) \rotatebox{60}{$w(e_2) $} \\
(.84 * .31)\rotatebox{-60}{$w(e_3) $} \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\centerline{\AffixLabels{\includegraphics{Vertex.pdf}}}

\caption{}
\label{fig:vertex}
\end{figure}

For an $N$--admissible weight system $w\in \mathcal W_\Gamma$, let $\beta_w \in \mathcal S^A(H)$ be associated to $w$ as above, and let $b_w \in V_S$ be the image of $\beta_w$ under the map $\Phi_H \colon \mathcal S^A(H) \to V_S$ of Lemma~\ref{lem:WRTspaceSpannedBySkeins}.

\begin{lem}
\label{lem:WRTbasis}
The subset $\mathcal B_\Gamma = \{ b_w \}_{w\in \mathcal W_\Gamma}$ is a basis for the vector space $V_S$.  \qed
\end{lem}

The weight system space $\mathcal W_\Gamma$ contains a special element $0$, assigning weight $0$ to each edge of $\Gamma$. This provides a prefered element $b_0\in V_S$, represented by the empty skein $[\varnothing ]\in \mathcal S^A(H)$. By definition, $b_0$ is the \emph{vacuum element} of $V_S$.

\section{The Witten-Reshetikhin-Turaev representation is irreducible}

I the rest of the article, $S$ will always denote a connected closed oriented surface.
Consider the Witten-Reshetikhin-Turaev representation  $\rho \colon \mathcal S^A(S) \to \mathrm{End}(V_S)$ of Lemma~\ref{lem:WRTrepExists}.

\begin{thm}
\label{thm:WRTirred}
The Witten-Reshetikhin-Turaev representation  $\rho \colon \mathcal S^A(S) \to \mathrm{End}(V_S)$ is irreducible.
\end{thm}

\begin{proof}
We will split the proof into several steps.
Let $W \subset V_S$ be a non-trivial linear subspace that is invariant under the image $\rho \bigl( \mathcal S^A(S) \bigl)$. We want to show that $W$ is equal to the whole space $V_S$.
For this, we will use a handlebody $H\cong \Sigma\times[0,1]$ bounding the surface $S$, a trivalent spine $\Gamma$ for the surface $\Sigma$, and the basis $\mathcal B_\Gamma = \{ b_w \}_{w\in \mathcal W_\Gamma}$ of Lemma~\ref{lem:WRTbasis}.

By hypothesis, $W$ is non-trivial, and therefore contains a non-trivial element $\sum_{w\in \mathcal W_\Gamma} \alpha_w b_w$ with $\alpha_w \in \mathbb C$.
Our first step is   borrowed from  \cite{RobertsIrred} and \cite{RobertsSkeinMCG}.

\begin{lem}
\label{lem:InvSpaceContainsBasis}
Let $\Gamma$ be a trivalent spine for the surface $\Sigma$.
If $\sum_{w\in \mathcal W_\Gamma} \alpha_w b_w$ is in the invariant subspace $W$, then every basis element $b_w \in \mathcal B_\Gamma$ with non-zero coefficient $\alpha_w\neq 0$ also belongs to $W$.
\end{lem}
\begin{proof}
[Proof of Lemma~\ref{lem:InvSpaceContainsBasis}]
For every edge $e$ of $\Gamma$, there exists a disk $D_e \subset H$ such that $D_e \cap \partial H$ is equal to the boundary $\partial D_e$, and such that $D\cap \Gamma = D\cap e$ consists of a single point; this is an immediate consequence of the fact that  the handlebody $H$ deformation retracts to the graph $\Gamma$. Consider $[\partial D_e] \in \mathcal S^A(S)$ and its image $\rho \bigl( [ \partial D_e] \bigr) \in \mathrm{End}(V_S)$. For every basis element $b_w\in \mathcal B_\Gamma$ associated to the weight system $w\in \mathcal W_\Gamma$, a computation as in \cite[Lemma~14.2]{LickorishBook} shows that
$$
\rho\bigl( [\partial D_e] \bigr) (b_w) = - (A^{2(w(e)+1)} + A^{-2(w(e)+1)}) b_w.
$$

As $A$ is a primitive $2N$--root of unity, the numbers $A^{2(i+1)}+A^{-2(i+1)}$ are non-zero and distinct as $i$ ranges over all admissible $i$, that is, over all $i \in \{0,2,4,\dots, N-3\}$ when $N$ is odd and over all $i \in \{ 0, 1, 2, \ldots, \frac N2 -2\}$ when $N$ is even.

After these observations, the lemma is just a matter of elementary linear algebra.
Consider an element $v=\sum_{w\in \mathcal W_\Gamma} \alpha_w b_w$   of $W$ such that $\alpha_{w_1}\neq 0$. We want to show that $b_{w_1}$ also belongs to $W$

If all the other coefficients $\alpha_w$ where $w\neq w_1$ are equal to 0, then $b_{w_1} = \frac 1 {\alpha_{w_1}} v\in W$ and we are done.

Otherwise,  there exists another weight system $w_2\neq w_1$ with $\alpha_{w_2}\neq 0$. The fact that $w_2\neq w_1$ in $\mathcal W_\Gamma$ means that there exists an edge $e$ of $\Gamma$ such that $w_2(e)\neq w_1(e)$. Then the invariant subspace $W$ also contains the element
\begin{align*}
v' &= \bigl(A^{2(w_2(e)+1)} + A^{-2(w_2(e)+1)}\bigr) v  + \rho\bigl( [\partial D_e] \bigr) (v) \\
&= \sum_{w\in \mathcal W_\Gamma} \bigl (A^{2(w_2(e)+1)} + A^{-2(w_2(e)+1)} -A^{2(w(e)+1)} - A^{-2(w(e)+1)}\bigr)  \alpha_w b_w
\end{align*}
Note that, in the basis $\mathcal B_\Gamma = \{ b_w\}_{w\in \mathcal W_\Gamma}$, the coordinate of $v'$ corresponding to $b_{w_1}$ is still non-zero since $w_2(e)\neq w_1(e)$, but that $v'$  now has one fewer non-zero coordinate than $v$ (because the coordinate of $v'$corresponding to $b_{w_2}$ is equal to  0). We can therefore replace $v$ by $v'\in W$, which is simpler.

Iterating this construction, we eventually reach an element of $W$ that has exactly one non-zero coordinate, corresponding to $b_{w_1}$. This proves that $b_{w_1}$ belongs to $W$, as required.
\end{proof}

We will need to extend  Lemma~\ref{lem:InvSpaceContainsBasis} to a slightly more general framework. Let a \emph{partial spine} for the surface $\Sigma$ be the union $\Gamma$ of a finite family of disjoint trivalent graphs and simple closed curves in $\Sigma$, such that each component of $\Sigma - \Gamma$ contains at least one component of the boundary $\partial \Sigma$. This condition guarantees that $\Gamma$ can be enlarged to a trivalent spine $\hat \Gamma$ for $\Sigma$, by adding a few vertices and edges.

The notion of an $N$--admissible weight system straightforwardly extends to partial spines: such a weight system consists of an $N$--admissible edge weight system on each trivalent graph component of $\Gamma$; and it assigns to each closed curve component $C$ of $\Gamma$ an even weight  $w(C)\in\{0,2,4, \dots, N-3\}$ if $N$ is odd, or a weight $w(C)\in\{0,1,2, \dots, \frac N2- 2 \}$ if $N$ is even. Again, plugging Jones-Wenzl idempotents into the edges and closed curve components of $\Gamma$ associates an element $\beta_w \in \mathcal S^A(H)$ to each $N$--admissible weight system $w\in \mathcal W_\Gamma$. As before,  we denote by $b_w=\Phi_H(\beta_w) \in V_S$  the image of $\beta_w$ under the map $\Phi_H \colon \mathcal S^A(H) \to V_S$ of Lemma~\ref{lem:WRTspaceSpannedBySkeins}, and we define $\mathcal B_\Gamma = \{b_w\}_{w\in \mathcal W_\Gamma}$. The only major difference is that $\mathcal B_\Gamma$ does not necessarily generate the Witten-Reshetikhin-Turaev space $V_S$.

\begin{lem}
\label{lem:InvSpaceContainsBasisPartialSpine}
The statement of Lemma~{\upshape\ref{lem:InvSpaceContainsBasis}} also holds when $\Gamma$ is only a partial spine for $\Sigma$. Namely, if $\Gamma$ is a partial spine for $\Sigma$ and if  $\sum_{w\in \mathcal W_\Gamma} \alpha_w b_w$ belongs to the invariant subspace $W$, then every  element $b_w \in \mathcal B_\Gamma$ with non-zero coefficient $\alpha_w\neq 0$ also belongs to $W$.
\end{lem}

\begin{proof}
Enlarge the partial spine $\Gamma$ to a trivalent spine $\hat \Gamma$ for $\Sigma$, by adding vertices inside of the edges and closed curve components of $\Gamma$ and then adding edges connecting these new vertices as necessary. A weight system $w\in \mathcal W_\Gamma$  determines an $N$--admissible edge weight system for $\hat \Gamma$ as follows: it assigns to each edge of $\hat\Gamma$ that is contained in $\Gamma$ the $w$--weight of the edge or closed curve component of $\Gamma$ that contains it; and it assigns  weight 0 to each of the new edges of $\hat \Gamma - \Gamma$. This defines an inclusion $\mathcal W_\Gamma \subset \mathcal W_{\hat \Gamma}$, and we will use the same letter $w$ to denote the original  $w\in \mathcal W_\Gamma$ and the edge weight system $w\in \mathcal W_{\hat \Gamma}$ for $\hat \Gamma$  that it defines.

We saw that a weight system $w\in \mathcal W_\Gamma$ for $\Gamma$ determines an element $\beta_w \in \mathcal S^A(H)$. Similarly, weighing the edges of $\hat \Gamma$ with $w \in \mathcal W_{\hat \Gamma}$ defines another element $\hat \beta_w \in \mathcal S^A(H)$. It turns out that $\hat \beta_w=\beta_w$. Indeed, this immediately follows from the idempotent property
\SetLabels
\mathrm{e}( .3* .5) $a $ \\
\mathrm{e}( .72* .5) $a $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-5pt}{\AffixLabels{\includegraphics{JW2.pdf}}}
$=$
\SetLabels
\mathrm{e}( .5* .5) $a $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-5pt}{\AffixLabels{\includegraphics{JW.pdf}}}
 of Jones-Wenzl idempotents, which shows that adding vertices inside of the edges and closed curve components of $\Gamma$ does not change the associated element of $\mathcal S^A(H)$; it is immediate that adding weight 0 edges also has no impact. As a consequence, the inclusion $\mathcal W_\Gamma \subset \mathcal W_{\hat \Gamma}$ induces an inclusion $\mathcal B_\Gamma \subset \mathcal B_{\hat \Gamma} \subset V_S$.

With this observation, every element $\sum_{w\in \mathcal W_\Gamma} \alpha_w b_w \in V_S $ can also be written as $\sum_{w\in \mathcal W_{\hat \Gamma}} \hat\alpha_w b_w$, by setting $\hat\alpha_w = \alpha_w$ when $w\in \mathcal W_{ \Gamma}\subset \mathcal W_{\hat\Gamma}$ and $\hat\alpha_w = 0$ when $w\in \mathcal W_{\hat \Gamma}-\mathcal W_\Gamma$. Lemma~\ref{lem:InvSpaceContainsBasisPartialSpine} then immediately follows by applying Lemma~\ref{lem:InvSpaceContainsBasis} to the trivalent spine $\hat \Gamma$.
\end{proof}

 Lemma~\ref{lem:InvSpaceContainsBasis} shows that the invariant subspace $W$ contains at least one  element $b_w\in \mathcal B_\Gamma$ associated to a weight system $w\in \mathcal W_\Gamma$ for a trivalent spine $\Gamma$. Our next goal is to show that $W$ contains the vacuum element $b_0$ corresponding to the zero weight system $0\in \mathcal W_\Gamma$ for all partial spines $\Gamma$. The following definition is designed to  measure progress in this direction.

 The \emph{complexity} of a weight system $w\in \mathcal W_\Gamma$ for a partial spine $\Gamma$  is defined as the triple
 $$
  |w| = \bigl( e(\Gamma), \max(w), n_{\max}(w) \bigr) \in \mathbb N^3
 $$
 where $e(\Gamma)$ is the number of edges of $\Gamma$, $\max(w)$ is the largest weight assigned by $w$ to the edges and closed curve components of $\Gamma$, and where $n_{\max}(w)$ is the number of edges and closed curve components where this maximum is attained (and where $\mathbb N$ denotes the set of non-negative integers). We endow $\mathbb N^3$ with the lexicographic order.

\begin{lem}
\label{lem:DecreaseComplexity}
If the invariant subspace $W$ contains an element $b_w\in \mathcal B_\Gamma$ associated to a non-zero weight system $w\in \mathcal W_\Gamma$ for a partial spine $\Gamma$, then $W$ contains another element $b_{w'} \in \mathcal B_{\Gamma'}$, represented by a partial spine $\Gamma'$ and a weight system $w' \in \mathcal W_{\Gamma'}$, such that  $|w'|<|w|$.
\end{lem}

\begin{proof}
We distinguish cases.

\medskip
\noindent\textsc{Case 1:} \emph{ The weight system $w$ assigns weight $0$ to an edge $e$ of $\Gamma$.}
The admissibility properties of $w$ imply that, if the endpoints of $e$ are not distinct and  correspond to the same vertex of $\Gamma$,  the third edge $e'$ emanating from this vertex has weight $w(e')$ equal to 0. Replacing $e$ by $e'$ if necessary, we can therefore assume that the endpoints of $e$ are distinct. The admissibility condition then shows that, at each of these endpoints, the two other adjacent edges have the same $w$--weight. Let $\Gamma'$ be the partial spine obtained from $\Gamma$ by removing $e$, and combining the edges that meet at each of its end vertices. By the above observation, $w$ induces a weight system $w'\in \mathcal W_{\Gamma'}$, and as in the proof of Lemma~\ref{lem:InvSpaceContainsBasisPartialSpine}, they are represented by the same basis element in $V_S$. By construction, $\Gamma'$ has one fewer edge than $\Gamma$, so that $|w'|<|w|$.

\begin{figure}[htbp]

\centerline{
$\displaystyle
\SetLabels
\mathrm{e}( .5* .5) \rotatebox{0}{\scalebox{.60}{$ w(e_0) $}} \\
( .135 *.85 )  \rotatebox{-60}{\scalebox{.60}{$ w(e_2) $}} \\
( .865* .625) \rotatebox{60}{\scalebox{.60}{$ w(e_3) $}} \\
( .13* .13)  \rotatebox{60}{\scalebox{.60}{$ w(e_1) $}} \\
( .87*.35 )  \rotatebox{-60}{\scalebox{.60}{$ w(e_4) $}}\\
\endSetLabels
\raisebox{-30pt}{\AffixLabels{{\includegraphics{Flip.pdf}}}}
=  \quad  \sum_{w'} \  \begin{Bmatrix} w(e_1) &w(e_2) & w'(e_0') \\ w(e_3) &w(e_4) & w(e_0)  \end{Bmatrix} \
\SetLabels
( .5* .4) \rotatebox{90}{\scalebox{.55}{$ w'(e_0') $}} \\
( .25* .9)  \rotatebox{-30}{\scalebox{.60}{$ w(e_2) $}} \\
( .75*.8 ) \rotatebox{30}{\scalebox{.60}{$ w(e_3) $}} \\
( .26* .07)  \rotatebox{30}{\scalebox{.60}{$ w(e_1) $}} \\
( .75*.16 )  \rotatebox{-30}{\scalebox{.60}{$ w(e_4) $}} \\
\endSetLabels
\raisebox{-36pt}{\AffixLabels{\rotatebox{90}{\includegraphics{Flip.pdf}}}}
$
}
\caption{The Flip Relation}
\label{fig:Flip}
\end{figure}

\medskip
\noindent\textsc{Case 2:} \emph{No edge weight is $0$, and the maximum weight $\max(w)$ is attained on an edge $e_0$ of $\Gamma$ whose endpoints are distinct.}

We apply to $\Gamma$ the classical Flip Relation represented in Figure~\ref{fig:Flip}.  The Flip Move replaces $\Gamma$ by a new partial spine $\Gamma'$ that differs from $\Gamma$  only in the edge $e_0$, and replaces $e_0$ by an edge $e_0'$ that connects differently the four edges meeting $e_0$.  The Flip Relation
$$
\beta_w = \sum_{w'} \begin{Bmatrix} w(e_1) &w(e_2) & w'(e_0') \\ w(e_3) &w(e_4) & w(e_0)  \end{Bmatrix} \beta_{w'}
$$
expresses the element $\beta_w \in \mathcal S^A(H)$ represented by  $w\in \mathcal W_{\Gamma}$ as a linear combination of  elements $\beta_{w'}\in \mathcal S^A(H)$ where $w'$ ranges over all $N$--admissible weight systems for $\Gamma'$ that coincide with $w$ over all edges that are common to $\Gamma$ and $\Gamma'$.

A key feature of this relation are the coefficients $ \begin{Bmatrix} w(e_1) &w(e_2) & w'(e_0') \\ w(e_3) &w(e_4) & w(e_0)  \end{Bmatrix}\in \mathbb C$, known as $6j$--symbols. A precise computation of these $6j$--symbols can be found in \cite{KauffmanLins} or \cite{MasbaumVogel}. The corresponding formula is usually complicated, and expresses a $6j$--symbol as a sum of several terms, each of which is a product of quantum integers and their inverses. However, it is somewhat simpler for the ``smallest'' of the $N$--admissible weight systems $w'$ that are compatible with $w$.

Consider the weight system $w_1'$ that  coincides with $w$ on $\Gamma-e_0 = \Gamma'-e_0'$ and  assigns weight
$$w_1'(e_0') = \max \bigl\{ | w(e_1) - w(e_4)|, |w(e_2) - w(e_3)|\bigr \}$$
to the edge $e_0'$. The formula is specially designed so that $w_1'$ is $N$--admissible. In fact, $w_1'$ is the edge weight that minimizes the weight $w'(e_0')$ among all $N$--admissible weight systems $w'\in \mathcal W_{\Gamma'}$ that coincide with $w$ outside of $e_0'$. We will not need this minimizing property, but  the following other feature of $w'$ is critical for our purposes:  for this specific weight system  $w_1'\in \mathcal W_{\Gamma'}$,  the sum occurring in the formula of \cite{MasbaumVogel} consists of a single term, and expresses the $6j$--symbol  $ \begin{Bmatrix} w(e_1) &w(e_2) & w_1'(e_0') \\ w(e_3) &w(e_4) & w(e_0)  \end{Bmatrix}$ as a product of non-zero quantum integers and their inverses. In particular, this $6j$--symbol is different from 0.

Remembering that $b_w= \Phi_H(\beta_w)\in \mathcal B_{\Gamma}$ and $b_{w'} = \Phi_H(\beta_{w'})\in \mathcal B_{\Gamma'}$ for  the map $\Phi\colon \mathcal S^A(H) \to V_S$ of Lemma~\ref{lem:WRTspaceSpannedBySkeins},
$$
b_w = \sum_{w'} \begin{Bmatrix} w(e_1) &w(e_2) & w'(e_0') \\ w(e_3) &w(e_4) & w(e_0)  \end{Bmatrix} b_{w'}
$$
in the Witten-Reshetikhin-Turaev space $V_S$. By hypothesis, $b_w$ belongs to the invariant subspace $W\subset V_S$. Lemma~\ref{lem:InvSpaceContainsBasisPartialSpine}  then shows that $W$ also contains the element $b_{w_1'} \in \mathcal B_{\Gamma'}$ corresponding to   $w_1'\in \mathcal W_{\Gamma'}$, since its coefficient in the above sum is different from 0.

By our hypothesis that the weights assigned by $w$  to the edges of $\Gamma$ are non-zero and bounded by $w(e_0)=\max(w)$, the weight $w_1'(e_0') $ defined above is strictly less than $\max(w)$. It follows that $w_1'$ has lower complexity $|w_1'|<|w|$ than $w$.

\medskip
\noindent\textsc{Case 3:} \emph{No edge weight is $0$, and the maximum weight $\max(w)$ is attained on an edge $e_0$ of $\Gamma$ whose endpoints are identified.}

Because the endpoints of $e_0$ are identified, we cannot apply a Flip Move at $e_0$. Instead, we will apply such a move at the remaining edge $e_1$ that is adjacent to the vertex corresponding to the two ends of  $e_0$. This gives a new partial spine $\Gamma'$, obtained from $\Gamma$ by replacing the edge $e_1$ by an edge $e_1'$ as in Figure~\ref{fig:Flip2}.

\begin{figure}[htbp]

\centerline{
$\displaystyle
\SetLabels
\mathrm{e}( .43* .49) \rotatebox{0}{\scalebox{.60}{$ w(e_1) $}} \\
( .115 *.83)  \rotatebox{-60}{\scalebox{.60}{$ w(e_3) $}} \\
( .11* .13) \rotatebox{60}{\scalebox{.60}{$ w(e_2) $}} \\
\mathrm{e}( .955* .5)  \rotatebox{90}{\scalebox{.60}{$ w(e_0) $}} \\
\endSetLabels
\raisebox{-30pt}{\AffixLabels{{\includegraphics{Flip2.pdf}}}}
=  \quad  \sum_{w'} \  \begin{Bmatrix} w(e_2) &w(e_3) & w'(e_1') \\ w(e_0) &w(e_0) & w(e_1)  \end{Bmatrix} \
\SetLabels
( .475* .4) \rotatebox{90}{\scalebox{.55}{$ w'(e_1') $}} \\
( .935* .4)  \rotatebox{90}{\scalebox{.60}{$ w(e_0) $}} \\
( .24*.9 ) \rotatebox{-30}{\scalebox{.60}{$ w(e_3) $}} \\
( .25* .06)  \rotatebox{30}{\scalebox{.60}{$ w(e_2) $}} \\
\endSetLabels
\raisebox{-36pt}{\AffixLabels{\rotatebox{0}{\includegraphics{Flip3.pdf}}}}
$
}
\caption{A special case of the Flip Relation}
\label{fig:Flip2}
\end{figure}

We then consider for $\Gamma'$ the $N$--admissible weight system $w_1'$ that assigns weight
$$w_1'(e_1') = \max \bigl \{ w(e_0)-w(e_2), w(e_0)-w(e_3)\bigr  \} $$
to the edge $e_1'$ and coincides with $w$ outside of $e_1'$. As in Case~2, the formula of \cite{MasbaumVogel} shows that the $6j$--symbol occurring as coefficient of $\beta_{w_1'}$ in the Flip Relation of Figure~\ref{fig:Flip2} is non-zero. An application of Lemma~\ref{lem:InvSpaceContainsBasisPartialSpine} again proves that the invariant subspace $W\subset V_S$ contains the element $b_{w_1'}\in \mathcal B_{\Gamma'}$ associated to $w_1' \in \mathcal W_{\Gamma'}$.

Note that $w_1'(e_1') < w(e_0) = \max (w)$, so that $|w_1'| \leq |w|$. Because the inequality is not necessarily strict, we are not quite done yet. However, we can now apply a Flip Move to $\Gamma'$ at the edge $e_0$, and use Case~2  to conclude.

\medskip
\noindent\textsc{Case 4:} \emph{$N$ is odd, and the maximum weight $\max(w)$ is attained on a closed curve component $C$ of $\Gamma$.}

Push this simple closed curve $C\subset \Sigma$ to  $\Sigma \times\{1\} \subset \partial H =S$ to consider it as a partial spine in a thickening of $S$, and let $[C^{S_2}] \in \mathcal S^A(S)$ be defined by assigning weight 2 to this partial spine, namely by plugging a Jones-Wenzl idempotent    \SetLabels
\mathrm{e}( .5* .5) $2 $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-5pt}{\AffixLabels{\includegraphics{JW.pdf}}}
in $C$; see Remark~\ref{rem:ChebyWRT}
 to explain the notation.

\begin{figure}[htbp]

$\SetLabels
( .5* .3) \tiny$ a$ \\
(.5 * .08) \tiny $ 2$ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-17pt}{\AffixLabels{\includegraphics{TorusJW2.pdf}}}
=
\SetLabels
( .5* .19) \tiny$ a\kern -3pt -\kern -3pt 2$ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-17pt}{\AffixLabels{\includegraphics{TorusJW.pdf}}}
+
\SetLabels
( .5* .19) \tiny$ a$ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-17pt}{\AffixLabels{\includegraphics{TorusJW.pdf}}}
+
\SetLabels
( .5* .19) \tiny$ a\kern -3pt +\kern -3pt 2$ \\
( * ) $ $ \\
( * ) $ $ \\
( * ) $ $ \\
\endSetLabels
\raisebox{-17pt}{\AffixLabels{\includegraphics{TorusJW.pdf}}}
$
\caption{A multiplication property for Jones-Wenzl idempotents}
\label{fig:TorusJW}
\end{figure}

If the subspace $W\subset V_S$ contains $b_w \in \mathcal B_\Gamma$, it also contains $\rho\bigl( [C^{S_2}] \bigr)(b_w)$ by invariance of $W$ under the action of $\rho\bigl( \mathcal S^A(S) \bigr)\subset \mathrm{End}(V_S)$. The relation of Figure~\ref{fig:TorusJW}, valid in a solid torus neighborhood of $C$, enables us to compute this element and gives
$$
\rho\bigl( [C^{S_2}] \bigr)(b_w) = b_{w'} + b_w +  b_{w''}
$$
where the weight systems $w'$ and $w''$ for $\Gamma$ coincide with $w$ outside of the closed curve component $C$, and respectively assign weight $w(C)-2$ and $w(C) +2$ to $C$. See \cite[Lemma~14.11]{LickorishBook} for a proof of this relation.

There is a little caveat needed here when $w(C)$ is equal to its maximum possible value $N-3$. Then, $w''$ assigns weight $N-1$ to $C$, and consequently is no longer $N$--admissible. The associated element $\beta_{w''} \in \mathcal S^A(H)$ still makes sense, but its image $\Phi_H (\beta_w)\in V_S$ under the map of Lemma~\ref{lem:WRTspaceSpannedBySkeins} is equal to 0. We consequently set $b_{w''}=0$ in this case.

In all cases, we can apply Lemma~\ref{lem:InvSpaceContainsBasisPartialSpine} to $\rho\bigl( [C^{S_2}] \bigr)(b_w) \in W$, and we conclude that $b_{w'}$ belongs to $W$. By construction, $w'(C) < w(C) = \max(w)$, so that $w'\in \mathcal W_\gamma$ has lower complexity $|w'|<|w|$. This concludes the proof in this case.

\medskip
\noindent\textsc{Case 5:} \emph{$N$ is even, and the maximum weight $\max(w)$ is attained on a closed curve component $C$ of $\Gamma$.}

The proof is almost identical to that of Case~4, except that we do not have to worry about keeping all weights even. Because of this, it suffices to consider the action of the element $[C]\in \mathcal S^A(S)$ represented by the closed curve $C$. Then, a computation similar to that of  Figure~\ref{fig:TorusJW} (see again \cite[Lemma~14.11]{LickorishBook}) gives that
$$
\rho\bigl( [C] \bigr)(b_w) = b_{w'} + b_w +  b_{w''}
$$
where the weight systems $w'$ and $w''$ for $\Gamma$ coincide with $w$ outside of the closed curve component $C$, and respectively assign weight $w(C)-1$ and $w(C)+1$ to $C$ (with $b_{w''}=0$ when $w(C) = \frac N2-2$). An application of Lemma~\ref{lem:InvSpaceContainsBasisPartialSpine} then shows that $W$ contains the basis element $b_{w' } \in \mathcal B_\Gamma$. Again, $w'\in \mathcal W_\gamma$ has lower complexity $|w'|<|w|$, and this concludes the proof in this case.

Since the five cases considered exhaust all possibilities, the proof of Lemma~\ref{lem:DecreaseComplexity} is now complete.
\end{proof}

We are now almost done with the proof of Theorem~\ref{thm:WRTirred}.

Recursively applying Lemma~\ref{lem:DecreaseComplexity}, we eventually reach a partial spine $\Gamma'$ such that the invariant subspace $W$ contains the element $b_0 \in \mathcal B_{\Gamma'}$ associated to the trivial weight system $0 \in \mathcal W_{\Gamma'}$. By definition $b_0$ is also the image of the empty skein $[\varnothing] \in \mathcal S^A(H)$ under the map $\Phi_H$ of Lemma~\ref{lem:WRTspaceSpannedBySkeins}.

Let $L$ be a framed link in the handlebody $H$. Push $L$ into a tubular neighborhood of the boundary $\partial H = S$, so that $L$ defines a skein $[L] \in \mathcal S^A(S)$.

By invariance of $W$ under the under the action of $\rho\bigl( \mathcal S^A(S) \bigr)\subset \mathrm{End}(V_S)$, it contains the element $\rho\bigl( [L] \bigr)(b_0)$. However, we also have, from the definition of the Witten-Reshetikhin-Turaev homomorphism $\rho$, that
$$
\rho\bigl( [L] \bigr)(b_0) = \rho\bigl( [L] \bigr) \bigl( \Phi_H([\varnothing]) \bigr) = \Phi_H\bigl([L\cup \varnothing] \bigr) =\Phi_H\bigl([L] \bigr)
$$
where, in the first two terms, $[L]$ denotes the element of $\mathcal S^A(S)$ represented by $L$, whereas $[L]$ is the element of $\mathcal S^A(H)$ represented by $L$ in the last two terms.

This proves that the invariant subspace $W\subset V_S$ contains the image $\Phi_H\bigl([L] \bigr) $ of every skein $[L] \in \mathcal S^A(H)$. Since these skeins generate $\mathcal S^A(H)$ and since $\Phi_H \colon \mathcal S^A(H) \to V_S$ is surjective by Lemma~\ref{lem:WRTspaceSpannedBySkeins}, this proves that $W$ is equal to the whole space $V_S$.

This completes the proof of Theorem~\ref{thm:WRTirred}.
\end{proof}

\section{The classical shadow of the Witten-Reshetikhin-Turaev representation}

\subsection{Threading polynomials along a framed link}

Given a framed (connected) knot $K$ in a 3--dimensional manifold $M$ and  a polynomial
$
P(x) = \sum_{i=0}^n a_i x^i
$, we can consider  the linear combination
$$
[K^P] = \sum_{i=0}^n a_i\, [K^{(i)}] \in \mathcal S^A(M)
$$
where, for each $i$, $K^{(i)}$ is the framed link obtained by taking $i$ parallel copies of $K$ in the direction indicated by the framing.
More generally, if $L\subset M$ is a framed link with components $K_1$, $K_2$, \dots, $K_l$, define
$$
[L^P] = \sum_{0\leq i_1, \,i_2,\dots, \,i_l \leq n} a_{i_1} a_{i_2} \dots a_{i_l}\, [K_1^{(i_1)} \cup K_2^{(i_2)} \cup \dots \cup K_l^{(i_l)} ] \in \mathcal S^A(M).
$$
By definition,  $[L^P]\in \mathcal S^A(M)$ is obtained by \emph{threading the polynomial $P$ along the framed link $L$}.

We will apply this construction to the (normalized) Chebyshev polynomials of the first and second type.

The \emph{$n$--th Chebyshev polynomial of the first type} $T_n(x)$ is defined by the properties that $T_n(x) = xT_{n-1}(x) - T_{n-2}(x)$, $T_0(x) =2$ and $T_1(x) =x$. The \emph{$n$--th Chebyshev polynomial of the second type} $S_n(x)$ is defined by the same recurrence relation $S_n(x) = xS_{n-1}(x) - S_{n-2}(x)$, the same  initial condition $S_1(x) =x$, but differs in the other initial condition $S_0(x) =1$. The two types of Chebyshev polynomials are related by the property that $T_n(x) = S_n(x) - S_{n-2}(x)$ for every $n$.

\begin{rem}
\label{rem:ChebyWRT}
The Chebyshev polynomials of the second type $S_n(x)$  are ubiquitous in the Witten-Reshetikhin-Turaev theory and, more generally, in the representation theory of the quantum group $\mathrm U_q(\mathfrak{sl}_2)$. In particular,  for each framed link $L$ in a 3--manifold $M$,  the element of $\mathcal S^A(M)$ obtained by plugging the $n$--th Jones-Wenzl idempotent in each component of $L$ is equal to the element $[L^{S_n}] \in \mathcal S^A(M)$. See  \cite[\S 13]{LickorishBook} or \cite[p.~715]{LickorishHandbook}.
\end{rem}

The following facts, which can for instance be found in Lemma 6.3 of \cite{BHMV3man}, are  crucial for our computations.

\begin{lem}
\label{lem:ChebyRelations}
Suppose that $A$ is a primitive $2N$--root of unity with $N$ odd, and let $V_S$ be the Witten-Reshetikhin-Turaev space of the surface $S$.
Let $K$ and $L$ be two disjoint framed links in a $3$--manifold $M$ bounded by $S$. Then, for the homomorphism $\Phi_M \colon \mathcal S^A(M ) \to V_S$ of Lemma~{\upshape\ref{lem:WRTspaceSpannedBySkeins}},
\begin{enumerate}
\item $\Phi_M \bigl( [K^{S_{N-1}} \cup L] \bigr) = 0$;
\item $\Phi_M \bigl( [K^{S_{N-2-n}} \cup L] \bigr)  = \Phi_M \bigl( [K^{S_{n}} \cup L] \bigr) $ for every integer $n$. \qed
\end{enumerate}
\end{lem}

\subsection{The classical shadow of the Witten-Reshetikhin-Turaev representation}
As usual, $S$ is a connected closed oriented surface.
Consider the character variety
$$
\mathcal R_{\mathrm{SL}_2(\mathbb C)}(S) = \{ \text{homomorphisms } r \colon \pi_1(S) \to \mathrm{SL}_2(\mathbb C)\}/\kern -4pt/ \mathrm{SL}_2(\mathbb C)
$$
where $\mathrm{SL}_2(\mathbb C)$ acts on homomorphisms $\pi_1(S) \to \mathrm{SL}_2(\mathbb C)$ by conjugation, and where the double bar indicates that one takes the quotient in the sense of geometric invariant theory. In practice, this means that two homomorphisms $r$, $r'\colon \pi_1(S) \to \mathrm{SL}_2(\mathbb C)$ represent the same point of $\mathcal R_{\mathrm{SL}_2(\mathbb C)}(S)$ if and only if they induce the same trace functions, namely if and only if $\mathrm{Tr}\,r(\gamma) = \mathrm{Tr}\,r'(\gamma) $ for every $\gamma \in \pi_1(S)$.

\begin{thm}[\cite{BonWonSkeinReps1}]
\label{thm:ClassicShadow}
Suppose that $A$ is a primitive $2N$--root of unity with $N$ odd. If $\rho\colon \mathcal S^A(S) \to \mathrm{End}(V)$ is an irreducible representation of the skein algebra $\mathcal S^A(S)$, then there exists a unique character $r _\rho\in \mathcal R_{\mathrm{SL}_2(\mathbb C)}(S)$ such that
$$
\rho \bigl( [K^{T_N}] \bigr) = - \mathrm{Tr}\,r_\rho(K)\,\mathrm{Id}_V
$$
for every framed knot $K \subset S \times [0,1]$. \qed
\end{thm}

By definition, $r _\rho\in \mathcal R_{\mathrm{SL}_2(\mathbb C)}(S)$ is the \emph{classical shadow} of the representation $\rho\colon \mathcal S^A(S) \to \mathrm{End}(V)$. See also \cite{Le} for an alternative approach to the key properties underlying this statement.

\begin{thm}
\label{thm:ClassiShadowWRT}
When $A$ is a primitive $2N$--root of unity with $N$ odd, the classical shadow of the Witten-Reshetikhin-Turaev representation $\rho\colon \mathcal S^A(S) \to \mathrm{End}(V_S)$ is the trivial character $\iota \in  \mathcal R_{\mathrm{SL}_2(\mathbb C)}(S)$, represented by the trivial homomorphism $\pi_1(S) \to \mathrm{SL}_2(\mathbb C)$.
\end{thm}

\begin{proof}
This is a relatively simple consequence of Lemma~\ref{lem:ChebyRelations}. Identify $S\times [0,1]$ to a tubular neighborhood of the boundary $S=\partial M$ in the 3--manifold $M$. To compute $\rho \bigl( [K^{T_N}] \bigr) \in \mathrm{End}(V_S)$ for a framed knot $K \subset S \times [0,1]$, Lemma~\ref{lem:WRTspaceSpannedBySkeins} shows that it suffices to consider its action on those elements of $V_S$ of the form $v= \Phi_M \bigl([L] \bigr) $  for a framed link $L\subset M$. Pushing $L$ away from the neighborhood  $S\times [0,1]$ of $S=\partial M$ in $M$,
\begin{align*}
\rho \bigl( [K^{T_N}] \bigr) (v)& =  \Phi_M\bigl( [K^{T_N} \cup L] \bigr) =    \Phi_M\bigl( [K^{xS_{N-1} - 2S_{N-2}} \cup L] \bigr)  \\
&=  \Phi_M\bigl( [K^{xS_{N-1}} \cup L] \bigr)  -  2\Phi_M\bigl( [K^{S_{N-2}} \cup L] \bigr),
\end{align*}
using the property that $T_n(x) = S_n(x) - S_{n-2}(x) = xS_{n-1}(x) -2S_{n-2}(x)$ for every $n$.

The term $[K^{xS_{N-1}} \cup L] \in \mathcal S^A(M)$ is also equal to $[K^{S_{N-1}}\cup K' \cup L] $ where $K'$ is a push-off of $K$ in the direction given by the framing. Its image $\Phi_M \bigl(  [K^{S_{N-1}}\cup K' \cup L]  \bigr)$ in $V_S$ is therefore equal to 0 by Part~(1) of Lemma~\ref{lem:ChebyRelations}. Similarly, Part~(2) of Lemma~\ref{lem:ChebyRelations} shows that
$$
\Phi_M\bigl( [K^{S_{N-2}} \cup L] \bigr) = \Phi_M\bigl( [K^{S_{0}} \cup L] \bigr) = \Phi_M\bigl( [K^{1} \cup L] \bigr)= \Phi_M\bigl( [ L] \bigr) =v.
$$
(Note that, since $1=x^0$, the skein $[K^1]=[K^{(0)}]$ is represented by 0 copies of the knot $K$, and is therefore trivial.)

Therefore, $\rho \bigl( [K^{T_N}] \bigr) (v) = -2v$ for every $v =  \Phi_M\bigl( [ L] \bigr)  \in V_S$ represented by a framed link $L \subset M$. Since these elements generate $V_S$ by Lemma~\ref{lem:WRTspaceSpannedBySkeins}, it follows that $\rho \bigl( [K^{T_N}] \bigr)=-2\,\mathrm{Id}_{V_S}$.

If $r_\rho$ is the classical shadow of the Witten-Reshetikhin-Turaev representation $\rho\colon \mathcal S^A(S) \to \mathrm{End}(V_S)$, this proves that $\mathrm{Tr}\, r_\rho(K)=2$ for every knot $K \subset S\times[0,1]$. This means that $r _\rho\in \mathcal R_{\mathrm{SL}_2(\mathbb C)}(S)$ is the character represented by the trivial homomorphism.
\end{proof}

\bibliographystyle{amsplain}


\section{Introduction}

Let $G$ be a connected real semi-simple Lie group and $P = MAN$
a minimal parabolic subgroup. Let $H <G$ be a closed
and connected subgroup.  We call $Z=G/H$ {\it real spherical} provided
that there is an open $P$-orbit on $Z$.
In \cite{KS} we have shown that
this condition implies that there are only finitely many $P$-orbits
on $Z$. The purpose of this paper is to explore the representation theoretic
significance of real sphericity.

\par This paper relies in part on the forthcoming article \cite{KKS} on the local structure
of real spherical spaces, which forces us to assume that the Lie algebra
of $H$ is an algebraic subalgebra of the Lie algebra of $G$.

\par For a Harish-Chandra
module $V$ with smooth completion $V^\infty$, we show that
if $G/H$ is real spherical with $PH$ open then
\begin{equation}\label{u-bound} \dim \mathrm{Hom}_H
(V^\infty, \mathbb{C}) \leq \dim (V/\bar\mathfrak{n} V)^{M\cap H}\, ,\end{equation}
with
$\bar\mathfrak{n} ={\rm Lie}(\bar N)$ corresponding to an
opposite parabolic subgroup.

In this context we recall that $V/\bar\mathfrak{n} V$ is  finite dimensional,  a
consequence of the Casselman-Osborne lemma (see \cite{W}, Sect.~3.7).
For symmetric spaces  (which
are real spherical) finite dimensionality of
$\mathrm{Hom}_H (V^\infty, \mathbb{C})$ was originally established by van den Ban
in \cite{vdB2}, Cor.~2.2.
Finally we remark that certain bounds on $\dim \mathrm{Hom}_H
(V^\infty, \mathbb{C})$ were obtained with a different technique by Kobayashi and Oshima in \cite{KO}, Thm.~2.4.

The bound in (\ref{u-bound}) is essentially sharp as equality
is obtained for $H=\bar N$ and generic irreducible representations $V$.
However, a statement is presented in Thm.~\ref{thm1} which
in general can be stronger than (\ref{u-bound}).

The main part of the proof of (\ref{u-bound}) is elementary in the sense that it
only invokes simple methods of ordinary differential equations,
applied to generalized matrix coefficients
on $Z$ (cf.~the proof of the subrepresentation theorem in \cite{W},
Sect.~3.8).  However these methods
typically result in asymptotic expansions only. To prove that a matrix coefficient which
is asymptotically zero (i.e.~of super-exponential decay), is in fact vanishing, we need more
elaborate analytic methods.
This is done in the end of the paper where we adapt some results from
\cite{CM} and \cite{vdB} to show that the relevant system of differential equations
has a regular singularity at infinity.

\par Suppose that $PH$ is open. According to \cite{KKS} there
exists a parabolic subgroup  $Q\supset P$ with a Levi decomposition $Q=LU$
such that $Q\cap H\subset L$ and $L/ (Q\cap H) Z(L)$ is compact.

In case that $V$ is irreducible and $H$-spherical, i.e.
$\mathrm{Hom}_H(V^\infty, \mathbb{C})\neq \{0\}$, we prove
the spherical subrepresentation theorem which asserts that
$V$ is a submodule of an induced module $\mathrm{Ind}_{\bar Q}^G \tau$,
where $\tau$ is an irreducible finite dimensional
representation of $L$ which is
$L\cap H$-spherical. This was established for symmetric spaces by Delorme
in \cite{Patrick}. For the group, it is the subrepresentation theorem of
Casselman~\cite{cass}.

\section{Some structure theory for real spherical spaces}\label{str}

Let $G$ be a real reductive group and $H<G$ a closed subgroup
with finitely many connected components. In what follows Lie algebras are
always denoted by corresponding lower case German letters, i.e. $\mathfrak{g}$ is
the Lie algebra of $G$, $\mathfrak{h}$ the Lie algebra of $H$ etc.

\par Let $P<G$ be a minimal parabolic subgroup of $G$. The homogeneous space
$Z=G/H$ is called {\it real spherical}, provided that $P$ admits open orbits
on $Z$. This means that by replacing $P$ with a conjugate we can obtain
that $PH$ is open, that is
$$\mathfrak{g} = \mathfrak{p} +\mathfrak{h}\, .$$

If $\mathfrak{l}$ is a reductive
Lie algebra, then we denote by $\mathfrak{l}_n$ the sum of all simple non-compact ideals
of $\mathfrak{l}$. We recall the following result from \cite{KKS}.

\begin{prop}\label{lst}
Let $H\subset G$ be algebraic groups over $\mathbb{R}$, and assume
$Z=G/H$ is real spherical. Let $P$ be a minimal parabolic subgroup such that
$PH$ is open.
Then there exists a parabolic subgroup $Q$ such that
\begin{enumerate}
\item\label{eins} $Q\supset P$.
\item\label{zwei} There is a Levi-decomposition $Q=LU$ such that
$$\mathfrak{l}_n\subset  \mathfrak{q} \cap \mathfrak{h} \subset \mathfrak{l}\, ,$$
and in particular, $L/(L\cap H)Z(L)$ is compact
(here $Z(L)$ denotes the center of $L$).
\item\label{vier} $L\cap H$ has finitely many components.
\item\label{drei} $QH=PH$
\end{enumerate}
\end{prop}

Through most of the paper we only need the properties (\ref{eins})-(\ref{zwei}), and for that
it clearly suffices to assume that $Z$ is real spherical and
{\it locally} algebraic,
that is, there are real algebraic groups $H_1\subset G_1$ with
Lie algebras $\mathfrak{h}$ and $\mathfrak{g}$.
We assume throughout that $Z$ is locally algebraic, but emphasize that
this assumption is only used to obtain (\ref{eins})-(\ref{zwei}) above.
 In Remarks \ref{Lcap H-rep} and \ref{second remark}
we need also (\ref{vier}), which is a consequence if
$G$ and $H$ are algebraic since then $L$ is algebraic.
The last property (\ref{drei})
is important for \cite{KKS} but will not be
needed here.

Note that in case $H$ is symmetric, then the minimal $\sigma\theta$-stable
parabolic subgroups (\cite{vdB3} Sect.~2)
satisfy   (\ref{eins})-(\ref{drei}).
Here $\theta$ is a Cartan involution which commutes
with the involution $\sigma$ which defines $\mathfrak{h}$.

\par If $G/H$ is real spherical and $P\subset Q=LU$
is as above, we let $\theta$ be a Cartan involution of $G$
which leaves $L$ stable. The existence of $\theta$ follows
since $L$ is a reductive subgroup of $G$.
Let $\mathfrak{g}=\mathfrak{k}+\mathfrak{s}$ be the Cartan decomposition and
$K\subset G$ the corresponding maximal compact subgroup,
then $K_L=L\cap K$ is maximal compact in $L$.
Let $\mathfrak{a}$ be a maximal abelian subspace in $\mathfrak{l}\cap\mathfrak{s}$.
We may assume $A$ is contained in $L\cap P$,
since this intersection is a minimal parabolic subgroup in $L$.
Let $L=K_L A N_L$  be an Iwasawa decomposition of $L$ and
put $N=N_LU$. Note that $\mathfrak{a}$ is maximal abelian in $\mathfrak{s}$ as well. Let $M$ be the centralizer
of $A$ in $K$. Then $P=MAN$, and it follows from (\ref{zwei}) above that
$\mathfrak{a}=\mathfrak{a}\cap\mathfrak{z}(\mathfrak{l})+\mathfrak{a}\cap\mathfrak{h}.$
Let $\mathfrak{a}_Z \subset \mathfrak{a}\cap \mathfrak{z}(\mathfrak{l})$ be a vector space complement to
$\mathfrak{a}\cap\mathfrak{h}$,
$$\mathfrak{a}=\mathfrak{a}_Z\oplus(\mathfrak{a}\cap\mathfrak{h}),$$
and $\mathfrak{m}_Z\subset \mathfrak{m}$ a subspace such that
$\mathfrak{a}_Z+\mathfrak{m}_Z$ complements $(\mathfrak{a}+\mathfrak{m})\cap\mathfrak{h}$ in $\mathfrak{a}+\mathfrak{m}$,
then we arrive at the direct sum
decomposition
\begin{equation}\label{ldeco} \mathfrak{g} = \mathfrak{h} \oplus \mathfrak{a}_Z \oplus \mathfrak{m}_Z \oplus\mathfrak{u}\, .\end{equation}

Let $L_n$ be the analytic subgroup of $L$ with Lie algebra $\mathfrak{l}_n$.
Since $\mathfrak{l}=\mathfrak{m}+\mathfrak{a}+\mathfrak{l}_n$ and $\mathfrak{l}_n\subset \mathfrak{h}$, and since $M$ meets every component of $L$
(see \cite{KnappBeyond}, Prop. 7.33),
we conclude
\begin{equation*
L=MAL_n
\end{equation*}
and $L_n\subset H$.
For any Lie group $J$ we denote by $J_0$
its identity component.

\begin{lemma}\label{LcapH} There exists
a vector subgroup $A_h\subset MA$ such that
\begin{equation}\label{identity component LcapH}
(L\cap H)_0=(M\cap H)_0A_hL_n.
\end{equation}
Moreover,
if $L\cap H$ has finitely many connected components, then
\begin{equation}\label{all of LcapH}
L\cap H=(M\cap H)A_hL_n,
\end{equation}
In this case $L\cap H$ is a real reductive group and
$L\cap H\cap K$ is a maximal compact subgroup.
\end{lemma}

\begin{proof} We have $L\cap H=(MA\cap H)L_n$.
Since the Lie algebra $\mathfrak{m}+\mathfrak{a}$ is compact, then so is its intersection with $\mathfrak{h}$.
It follows that $(MA\cap H)_0$ is the direct product of a compact group $S$
and a vector group $A_h$. The compact group projects trivially to $A$ along $M$,
hence $S\subset M$
(whereas it need not be the case that $A_h\subset A$). This implies (\ref{identity component LcapH}).
The argument for (\ref{all of LcapH}) is the same. The last statement follows
from (\ref{all of LcapH}).
\end{proof}

\section{Finite multiplicity}

We assume throughout this section that $G/H$ is spherical
with $PH$ open, and that a Cartan involution $\theta$
of $G$ has been chosen as described in the preceding section.
Accordingly we write $P=MAN$. Let $\Sigma^+(\mathfrak{g}, \mathfrak{a})$ be the set
of positive roots of $\mathfrak{a}$ in $\mathfrak{g}$, and let $\bar P=\theta(P)=MA\bar N$ be
the opposite parabolic subgroup.

\par For a Harish-Chandra module $V$ we denote by $V^\infty$ its unique smooth
moderate growth Fr\'echet completion.
Note that  $V^\infty$ is a $G$-Fr\'echet module and we set
$V^{-\infty}:= (V^\infty)'$ for its strong dual. Our goal is to provide a bound for the dimension of
the space of $H$-invariants
$$(V^{-\infty})^H:=  \mathrm{Hom}_H  (V^\infty, \mathbb{C})$$
where $\mathrm{Hom}$ stands for continuous linear homomorphisms.

\par For $v\in V^\infty$ and $\eta$ a continuous $H$-invariant functional on $V^\infty$ we form the
matrix coefficient
$$m_{v,\eta}(g):= \eta (\pi(g^{-1})v) \qquad (g\in G)\, .$$
Note that $m_{v,\eta}$ is a smooth function on $G$, even analytic  for $v\in V$.
We start with a general lemma which we prove later. We shall say that
$f:\mathbb{R}\to\mathbb{C}$ is of {\it super exponential decay} for $t\to\infty$ if
$f(t)=O(e^{\lambda t})$ for all $\lambda\in\mathbb{R}$.

\begin{lemma}\label{lem1} Let $V$ be a Harish-Chandra module.
Let $X\in \mathfrak{a}$ be any element which is strictly anti-dominant with respect to $P$ and
for which $\alpha(X)\in\mathbb{Q}$ for each $\alpha\in\Sigma(\mathfrak{g},\mathfrak{a})$.
Fix $\eta\in (V^{-\infty})^H$. Suppose that
for all $v\in V$ the function
$$\mathbb{R} \ni t\mapsto m_{v,\eta} (\exp(tX)) \in \mathbb{C}$$
is of super exponential decay for $t\to \infty$. Then $\eta=0$.
\end{lemma}

\begin{proof} Let $v\in V$ and set $F_v(t)= m_{v,\eta} (\exp(tX))$.
We may assume that $\alpha(X)$ is an integer for each root $\alpha$.
With the new variable $z=e^{-t}$ we will show in Section \ref{appendix}
that then $F_v$ admits an expansion
\begin{equation}\label{power}
F_v(t)=\sum_{j=1}^N \sum_{k=0}^M  z^{\lambda_j} (\log z)^k f_{j,k}(z),\quad(t\gg 0),
\end{equation}
where $\lambda_j\in \mathbb{C}$ and the $f_{j,k}$ are holomorphic functions in
a neighborhood of $z=0$ in $\mathbb{C}$. The fact that $F_v$ is of superexponential
decay then forces $F_v=0$ (this follows for example
from \cite{W}, Lemma 4.A.1.2).  Then
$F_v(0)=\eta(v)=0$, and hence $\eta=0$ since $v\in V$ was arbitrary.
\end{proof}

\par For a Harish-Chandra module $V$ we recall that $V/\bar\mathfrak{n} V$ is a finite dimensional module
for $\bar P=MA\bar N$ with $\bar N$ acting trivially.
Recall from Proposition \ref{lst} the
parabolic subgroup $Q=LU$ and its
subalgebra $\mathfrak{q} =\mathfrak{l}+\mathfrak{u} \supset \mathfrak{p}$. Let $\bar\mathfrak{q}=\theta(\mathfrak{q})=\mathfrak{l}+\bar\mathfrak{u}$ and
define a subalgebra $\bar\mathfrak{q}_1$ of $\bar\mathfrak{q}$ by
\begin{equation} \label{eq1}
\bar\mathfrak{q}_1:=(\mathfrak{l}\cap\mathfrak{h}) +{\bar\mathfrak{u}} \subset \mathfrak{h} + {\bar\mathfrak{n}}. \end{equation}
Note that $\bar\mathfrak{n}\subset \mathfrak{l}_n+\bar\mathfrak{u} \subset\mathfrak{q}_1$. Hence the
quotient $V/\bar\mathfrak{q}_1 V$ is finite dimensional. Moreover,
it carries a natural action of $L\cap K\cap H$ since $\bar\mathfrak{q}_1$ is $L\cap H$-invariant.

\begin{theorem}\label{thm1}  Let $V$ be a Harish-Chandra module and
$V^\infty$ its unique smooth Fr\'echet completion. Then
$$ \dim \mathrm{Hom}_H  (V^\infty, \mathbb{C}) \leq \dim (V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}\, .$$
In particular {\rm (\ref{u-bound})} from the introduction is valid.
\end{theorem}

\begin{rmk}\label{Lcap H-rep}
If $L\cap H$ has finitely many components then
it follows from Lemma \ref{LcapH} that the trivial action of
$\mathfrak{l}\cap\mathfrak{h}$ on $V/\bar\mathfrak{q}_1 V$ lifts to an action of $L\cap H$ which
agrees with the natural action of $L\cap H\cap K$.
In this case then $(V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}=(V/\bar\mathfrak{q}_1 V)^{L\cap H}$ for this action.
\end{rmk}

\begin{proof}
\par
Set $\Upsilon=(V^{-\infty})^H$. It is almost immediate from the definitions
(see \cite{KSS} eq.~(3.7)) that there exists $\delta\in \mathfrak{a}^*$
such that for all $v\in V^\infty, \eta\in\Upsilon$ there is a constant $C_{v,\eta}>0$
such that
\begin{equation} \label{ibound}
|m_{v,\eta} (a)| \leq C_{v, \eta}  a^{\delta} \qquad (a\in {A^-}) \, .\end{equation}

\par Fix an element $X\in \mathfrak{a}$ as in Lemma \ref{lem1}. After rescaling, we may assume that
$\min_{\alpha \in \Sigma^-(\mathfrak{g}, \mathfrak{a})}  \alpha(X)=1$.
For all $v\in V^\infty$ and $\eta\in\Upsilon$ we define
\begin{equation*}
F_{v,\eta}(t):= m_{v,\eta} (\exp(tX)) \qquad (t\geq 0)\,  .
\end{equation*}

We say that $\Lambda\in \mathbb{R} \cup\{-\infty\}$ bounds $\eta\in\Upsilon$,
provided for all $\Lambda'>\Lambda$ and  all $v\in V$,
\begin{equation} \label{ibound2}
\sup_{t\geq 0} \,e^{-t \Lambda'} |F_{v,\eta} (t)| <\infty \, .
\end{equation}
Let $\Upsilon_\Lambda\subset\Upsilon$ denote the space of elements bounded by $\Lambda$,
then $\Upsilon_{\Lambda_1}  \subset \Upsilon_{\Lambda_2}$ for $\Lambda_1\le\Lambda_2$.
It follows from (\ref{ibound}) that
\begin{equation}\label{Ups}
\Upsilon=\cup_\Lambda\Upsilon_\Lambda,
\end{equation}
and from Lemma \ref{lem1} that
\begin{equation}\label{zero}
\Upsilon_{-\infty}=\cap_\Lambda\Upsilon_\Lambda=\{0\}
\end{equation}

The element $X$ acts on the space
$(V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}$. First note that
$\mathfrak{l}\cap\mathfrak{h}=\mathfrak{l}_n+\left((\mathfrak{m}+\mathfrak{a})\cap\mathfrak{h}\right)$
is normalized by $\mathfrak{a}$. Hence so is $\bar\mathfrak{q}_1$ and thus
$X$ acts on $V/\bar\mathfrak{q}_1 V$.
Likewise, since $\mathfrak{a} =\mathfrak{a}\cap\mathfrak{z}(\mathfrak{l}) + \mathfrak{a}\cap \mathfrak{h}$
we find that $\operatorname{Ad}(l)X=X \mod \mathfrak{l}\cap\mathfrak{h}$ for $l\in L\cap K\cap H$, and hence
$X$ preserves the space of $L\cap K\cap H$-invariant vectors in the quotient.

We denote by $\Xi$ the set of values $-\operatorname{Re}\lambda$
where $\lambda\in\mathbb{C}$ is an eigenvalue for $X$ on
$(V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}$, and write
$\Xi=\{\mu_1,\dots,\mu_l\}$ where
$$\mu_{l+1}:=-\infty<\mu_l<\dots<\mu_1<\mu_0:=+\infty.$$
Let $m_1,\dots, m_l$ denote
the sums of the algebraic multiplicities of the corresponding eigenvalues
$\lambda$. Then
\begin{equation*}
m_1+\dots+m_l=n_0:=\dim (V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}.
\end{equation*}
We shall prove:
\begin{enumerate}
\item\label{one}
$\Upsilon_\Lambda=\Upsilon_{\mu_{k+1}}$ for $\mu_{k+1}\le\Lambda<\mu_{k}$ and
$k=0,\dots,l$
\item\label{two} The codimension
of $\Upsilon_{\mu_{k+1}}$ in $\Upsilon_{\mu_{k}}$ is  at most $m_{k}$ ($k=1,\dots,l$).
\end{enumerate}
It is easily seen that these statements together with (\ref{Ups}) and
(\ref{zero})
imply the theorem. Before proving (\ref{one}) and (\ref{two})
we need some preparations.

Let $\bar w_1, \ldots, \bar w_{n_0}$ be a basis for $(V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}$
and let $B$ denote the corresponding $n_0\times n_0$-matrix defined by
$$X\bar w_j=\sum_{k} b_{jk} \bar w_k$$
for $j=1,\dots,n_0$.
We choose a representative $w_j \in V$ for each $\bar w_j$
and define
$$u_j=Xw_j-\sum_k b_{jk}w_k\in \bar\mathfrak{q}_1 V.$$

We can arrange that  $B$
consists of block matrices $B_1,\dots,B_l$ along the diagonal, such that each
$B_k$ is an $m_k\times m_k$ matrix all of whose eigenvalues
have real part $-\mu_k$.
In the following we shall make the identification
\begin{equation*
\mathbb{C}^{n_0}=\mathbb{C}^{m_1}\times \dots\times \mathbb{C}^{m_l}
\end{equation*}
and write elements $x\in\mathbb{C}^{n_0}$ accordingly
as $x=(x_1,\dots,x_l)$.

For a given $\eta\in\Upsilon$ let
$$F(t)=(F_1(t),\dots,F_l(t))\in \mathbb{C}^{n_0}$$
where $F_k(t)\in\mathbb{C}^{m_k}$ is the $m_k$-tuple with entries $F_{w_j,\eta}(t)$,
corresponding to the $k$-th block of $B$.
Likewise we put
$$
R(t)=(R_1(t),\dots,R_l(t))\in \mathbb{C}^{n_0}$$
where $R_k(t)$ has entries
$F_{u_j,\eta}(t)$.
Then for each $k$,
$$ F_k'(t) = -B_k  F_k(t) - R_k(t). $$
Hence
\begin{equation} \label{eqk}
F_k(t) = e^{-tB_k}c_0- e^{-tB_k}\int_0^t e^{sB_k} R_k(s) \,ds\,
\end{equation}
where  $c_0=F_k(0)$.

\par  Next we recall from \cite{KSS}, proof of Thm.~3.2 (with
$\mathfrak{a}^+,\bar\mathfrak{n}$ interchanged by $\mathfrak{a}^-,\mathfrak{n}$), that the elements in
$\bar\mathfrak{n} V$ satisfy improved bounds as follows.
Assume $\eta\in\Upsilon_\Lambda$
and $\Lambda'>\Lambda$. Then
(\ref{ibound2}) is valid for all $v\in V$, and it follows
for all $ u \in \bar\mathfrak{n} V$ that
\begin{equation}\label{ibound3}
\sup_{t\geq 0} \, e^{-t (\Lambda'-1)} |F_{u,\eta} (t)| <\infty  .
\end{equation}

The bound (\ref{ibound3}) is valid for $u\in \bar\mathfrak{q}_1V$ as well.
To see this we note first that since $\mathfrak{a} =\mathfrak{a}\cap\mathfrak{z}(\mathfrak{l}) + \mathfrak{a}\cap \mathfrak{h}$ we have
\begin{equation}\label{XY}
m_{v,\eta}(\exp(tX))=m_{v,\eta}(\exp(tY)), \qquad(v\in V),
\end{equation}
for some $Y\in\mathfrak{a}$ which centralizes $\mathfrak{l}$.
Then $u= \sum_j X_jv_j +u'$ with $X_j\in \mathfrak{l}\cap\mathfrak{h}$, $v_j\in V$ and
$u' \in \bar\mathfrak{n} V$, and it follows that
$YX_jv_j=X_jYv_j\in \mathfrak{h} V$. Hence $m_{u,\eta}(\exp(tY))=m_{u',\eta}(\exp(tY))$.
We conclude that
$$F_{u,\eta}=F_{u',\eta}$$
and hence (\ref{ibound3}) is valid as claimed.
We conclude the existence of a constant such that
\begin{equation}\label{ibound4}
|R_k(t)|\le Ce^{t (\Lambda'-1)}.
\end{equation}

Based on (\ref{ibound4})
we shall provide the following two estimates for $F_k(t)$:
\begin{equation}\label{estimate}
\sup_{t\ge 0}\, e^{-t\gamma} |F_k(t)|<\infty,\quad
\forall\gamma>\max\{\mu_k,\Lambda-1\}
\end{equation}
for all $k=1,\dots, l$, and
\begin{equation}\label{better estimate}
\sup_{t\ge 0}\, e^{-t\gamma} |F_k(t)|<\infty,\quad \forall\gamma>\Lambda-1
\end{equation}
for those $k=1,\dots$ for which $\mu_k> \Lambda$.

Let $\gamma>\max\{\mu_k,\Lambda-1\}$ and
let $\Lambda'\in (\Lambda,\gamma+1)$.
It is clear that the first term in (\ref{eqk}) is bounded
by a polynomial times $e^{t\mu_k}$.
Applying (\ref{ibound4}) we see that the integrand
$e^{sB_k} R_k (s)$ of the second term
is dominated by a polynomial times
$e^{s(-\mu_k+\Lambda'-1)}$.
It follows that $|F_k(t)|$ is
is dominated by a polynomial times $e^{t\max\{\mu_k,\Lambda'-1\}}$,
and this implies (\ref{estimate}).

Before we prove the second estimate we note the following fact.
If $\mu_k>\Lambda'-1$ then
$e^{sB_k} R_k (s)$ is integrable to infinity. Hence if $\mu_k>\Lambda-1$,
we can replace
the solution formula
(\ref{eqk}) by
\begin{equation} \label{eqk1}
F_k(t) = e^{-tB_k}c_\infty+ e^{-tB_k}\int_t^\infty e^{sB_k} R_k (s) \,ds
\end{equation}
where $c_\infty=\lim_{s\to\infty} e^{sB_k} F_k(s)$.
Note that with (\ref{ibound4}) the equation
(\ref{eqk1}) implies (\ref{better estimate})
for every $k$ for which the limit $c_\infty$ vanishes.

For the proof of (\ref{better estimate}) we assume
$\mu_k > \Lambda$.
Then it follows from (\ref{ibound2}) with $\Lambda<\Lambda'<\mu_k$
that $\lim_{s\to\infty} e^{sB_k} F_k(s)=0$.
Hence $c_\infty=0$ and (\ref{eqk1}) implies
(\ref{better estimate}).

We are now ready to prove our claims (\ref{one}) and (\ref{two}).
Let $\Lambda\in\mathbb{R}$ and let
$\mu_\Lambda=\max\Xi\cap \mathbb{R}_{\le\Lambda}$ (with $\max\emptyset=-\infty$).
Let
$\eta\in\Upsilon_\Lambda$.
We will show
\begin{equation}\label{Upsi}
\eta\in\Upsilon_{\max\{\mu_\Lambda,\Lambda-1\}}.
\end{equation}
By iteration this will then imply that
$\Upsilon_\Lambda=\Upsilon_{\mu_\Lambda}$
as stated in (\ref{one}). Indeed, if $\Lambda-1\le\mu_\Lambda$
then this follows immediately.
Otherwise (\ref{Upsi}) implies
$\Upsilon_\Lambda=\Upsilon_{\Lambda-1}$. Since
$\Lambda$ was arbitrary we reach $\Upsilon_\Lambda=\Upsilon_{\mu_\Lambda}$ after
finitely many steps.

Let $\lambda'>\max\{\mu_\Lambda,\Lambda-1\}$. We need to prove
\begin{equation}\label{claim1}
\sup_{t\ge 0}\, e^{-t\lambda'} |F_{v,\eta}(t)| < \infty
\end{equation}
for all $v\in V$.

Note that for $m\in L\cap K\cap H$ we have $F_{mv,\eta}=F_{v,\eta}$ (see (\ref{XY}))
and hence $F_{v,\eta}=F_{v_0,\eta}$, where $v_0$ is the projection
of $v$ to $V^{L\cap K\cap H}$. Hence for every $v\in V$ we have that $F_{v,\eta}$
is a linear combination the functions $F_{w_j,\eta}$ plus
$F_{u,\eta}$ for some $u\in \bar\mathfrak{q}_0 V$.
Thus it follows from (\ref{ibound3}) that it suffices to establish
(\ref{claim1}) for $v=w_j$ for all $j$, and hence it suffices to show
\begin{equation}\label{claim2}
\sup_{t\ge 0}\, e^{-t\lambda'} |F_k(t)| < \infty
\end{equation}
for $k=1,\dots,l$. If $\mu_k\le\mu_\Lambda$ then this is
clear from (\ref{estimate}).
On the other hand, if $\mu_k>\mu_\Lambda$ then $\mu_k>\Lambda$
and (\ref{claim2}) follows from (\ref{better estimate}).

Next we show (\ref{two}).
Let $1\le k_0\le l$ and $\eta\in \Upsilon_{\Lambda}$
where $\Lambda=\mu_{k_0}$.
Then $\mu_{k_0}>\Lambda-1$ and according to
(\ref{eqk1}) the limit
$\lim_{s\to\infty} e^{sB_{k}}F_{k}(s)$ exists for every $k\le k_0$.
Moreover, for $k< k_0$ we have $\mu_k>\Lambda$ and the limit vanishes
as seen in the proof of (\ref{better estimate}) above.
If we assume that this limit is zero also for $k=k_0$,
then (\ref{better estimate}) is valid for all $k\leq k_0$,
from which we conclude as above that
(\ref{Upsi}) holds, with $\mu_\Lambda$ now replaced by $\mu_{k_0+1}$.
By step (\ref{one}) this implies
that $\eta\in\Upsilon_{\mu_{k_0+1}}$. Hence
$$\Upsilon_{\mu_{k_0+1}}=\{\eta\in \Upsilon_{\mu_{k_0}}\mid \lim_{s\to\infty}
e^{sB_{k_0}}F_{k_0}(s)=0\}$$
and as $\eta\mapsto \lim_{s\to\infty} e^{sB_{k_0}}F_{k_0}(s)$
is linear into $\mathbb{C}^{m_{k_0}}$, (\ref{two}) follows.
\end{proof}

\section{The spherical subrepresentation theorem}\label{sphsub}

\par Let $Z=G/H$ be real spherical with $PH$ open, and recall
the definition (\ref{eq1}) of the subalgebra $\bar\mathfrak{q}_1\subset\bar\mathfrak{q}$.
The algebraic version of the
subrepresentation theorem is:

\begin{cor}\label{C-lem}
Let $V$ be a Harish-Chandra module and assume that $\mathrm{Hom}_H (V^\infty, \mathbb{C})\neq \{0\}$.
Then $(V/\bar\mathfrak{q}_1 V)^{L\cap K\cap H}\neq \{0\}$.
\end{cor}

\begin{proof} Immediate from Theorem \ref{thm1}.
\end{proof}

Let $\bar Q=\theta(Q)$ denote the parabolic subgroup opposite to $Q$, with
Levi decomposition $\bar Q=L\bar U$.
If $\tau$ is a finite dimensional representation of $\bar Q$, then we write
$\mathrm{Ind}_{\bar Q}^G$ for the induced $(\mathfrak{g}, K)$-module of $K$-finite smooth sections
of the $G$-equivariant vector bundle $\tau \times_{\bar Q} G \to \bar Q\backslash G$.

\begin{theorem}
Let $V$ be an irreducible Harish-Chandra module such that
$\mathrm{Hom}_H (V^\infty, \mathbb{C})\neq \{0\}$.

Then there exists a finite
dimensional irreducible representation
$\tau$ of $\bar Q$, trivial on $\bar U$, such that
 $\mathrm{Hom}_{(\mathfrak{l}\cap\mathfrak{h},L\cap H\cap K)}(\tau,\mathbb{C})\neq \{0\}$, and an embedding
$$V \hookrightarrow \mathrm{Ind}_{\bar Q}^G \tau\, .$$
\end{theorem}

\begin{rmk}\label{second remark}
Assume that $L\cap H$ has finitely many connected components, and
recall (see Proposition \ref{lst}) that this assumption
is valid for example when $G$ and $H$ are algebraic. Then
it follows from Lemma \ref{LcapH} that
$\mathrm{Hom}_{L\cap H}(\tau,\mathbb{C})\neq \{0\}$
for
the representation
$\tau$ above.
\end{rmk}

\begin{proof}
Let $\bar\mathfrak{q}_0=\mathfrak{l}_n+\bar\mathfrak{u}$, then $\bar\mathfrak{q}_1=\bar\mathfrak{q}_0+(\mathfrak{l}\cap\mathfrak{h})$.
The space $V/\bar\mathfrak{q}_0 V$ is
finite-dimensional since $\bar\mathfrak{n}\subset\bar\mathfrak{q}_0$,
and it carries compatible actions of
$\bar\mathfrak{q}$ and $L\cap K$ since $\bar\mathfrak{q}_0$ is an
$L\cap K$-invariant ideal in $\bar\mathfrak{q}$.  Since $L$ is reductive
the action of the pair $(\mathfrak{l},L\cap K)$ lifts uniquely
to a representation of $L$ on $V/\bar\mathfrak{q}_0 V$.
We can then extend to an action of $\bar Q=L\bar U$,
which is trivial on $\bar U$ and compatible
with the action of $\bar\mathfrak{q}$.

The quotient map
$V/\bar \mathfrak{q}_0V\to V/\bar\mathfrak{q}_1V$ is clearly a homomorphism for the pair
$(\mathfrak{l}\cap\mathfrak{h},L\cap H\cap K)$. By Corollary \ref{C-lem} the module
$V/\bar \mathfrak{q}_1V$ has a non-zero vector fixed by the compact group
$L\cap K\cap H$. Hence also its dual admits such a vector,
and by composing with the quotient map we obtain that
$\mathrm{Hom}_{(\mathfrak{l}\cap\mathfrak{h},L\cap H\cap K)}(V/\bar\mathfrak{q}_0 V,\mathbb{C})\neq \{0\}$.

Let $\tau$ be an irreducible $\bar Q$-subrepresentation
of $V/\bar\mathfrak{q}_0V$ for which  $\mathrm{Hom}_{(\mathfrak{l}\cap\mathfrak{h},L\cap H\cap K)}(\tau,\mathbb{C})\neq \{0\}$.
Since the quotient map $V/\bar \mathfrak{u} V\to V/\bar \mathfrak{q}_0V$ is $L$-equivariant, we have
$$\operatorname{Hom}_L(V/\bar \mathfrak{u} V,\tau)\neq \{0 \}.$$
As $V$ is irreducible, the desired embedding follows by
Frobenius reciprocity
(see \cite{HS}, Theorem 4.9, and note that our induction is not normalized).
\end{proof}

\section{Regular singularities}\label{appendix}

The goal of this section is to provide a proof for the expansion
(\ref{power}).  To begin with let us first recall that matrix coefficients on $Z$ satisfy
certain systems of differential equations. For that we fix a Harish-Chandra module $V$
and a $K$-type $\tau$ occurring in $V$. We denote by $V[\tau]$ the $\tau$-isotypical part
of $V$, and consider for $\eta\in  (V^{-\infty})^H$ and $v\in V[\tau]$
the matrix coefficient
$$f(a)=m_{v, \eta}(a)= \eta(\pi(a^{-1})v)\, $$
on $A$ (where $A$ originates from the Iwasawa decomposition $G=KAN$
chosen in Section \ref{str}).

For simplicity we assume that $V$ is irreducible and obtain that
the center $\mathcal{Z}(\mathfrak{g})$ of $\mathcal{U}(\mathfrak{g})$
acts by scalars on $V$
(otherwise we replace the annihilating ideal of $V$ in $\mathcal{Z}(\mathfrak{g})$
with an ideal of finite co-dimension, and proceed as before).
The theory of $\tau$-radial parts (see Remark \ref{rmk2} below)
then gives a system of differential equations for $f$ on a
subcone of $A$.

\par
Let $R_1, \ldots, R_n$ be a basis of
root vectors of $\mathfrak{u}$, say $R_j\in\mathfrak{g}^{\alpha_j}$ corresponds to the
root $\alpha_j$. Set $Q_j:=\theta(R_j)\in \mathfrak{g}^{-\alpha_j}$.
For $t\in\mathbb{R}$ we let
\begin{equation*
\mathfrak{a}_t=\{ X\in\mathfrak{a} \mid \alpha_j(X)<-t ,\, j=1,\dots,n\}
\end{equation*}
and $A_t:=\exp(\mathfrak{a}_t)$. Let
$D_\epsilon:=\{|z|<\epsilon\}^n \subset \mathbb{C}^n$
for $\epsilon>0$, and
define
\begin{equation}\label{iota}
\iota: A_t \to \mathbb{C}^n , \ \ a\mapsto \iota(a)\:=(a^{\alpha_j})_{j=1,\dots,n}\,
\end{equation}
where $a^{\alpha_j}=e^{\alpha_j(X)}$ for $a=\exp(X)$,
then $\iota(A_t)\subset D_{e^{-t}}.$

\begin{lemma}\label{app1} There exists $t\ge 0$ such that
for all $a\in A_t$ one has
\begin{equation} \label{decomp}
\mathfrak{g}= \operatorname{Ad}(a^{-1}) \mathfrak{k} +\mathfrak{a}_Z +\mathfrak{h} \end{equation}
\end{lemma}

For symmetric spaces this is obtained with $t=0$ in \cite{vdB2} Lemma 1.5.

\begin{proof} Recall the decomposition
(\ref{ldeco}). As $\mathfrak{m}_Z$ is centralized by $A$, we have
$\mathfrak{m}_Z\subset\operatorname{Ad}(a^{-1})\mathfrak{k}$ for every $a\in A$. Hence
it suffices to show that each $R_j\in\mathfrak{u}$ decomposes according to (\ref{decomp}).
For all $1\leq j \leq n$ we obtain from  (\ref{ldeco})
\begin{equation} \label{decomp5}
Q_j= X_{\mathfrak{m}, j} + X_{\mathfrak{a}, j} + X_{\mathfrak{h}, j} + \sum_{i=1}^n c_{ij} R_i
\end{equation}
with $X_{\mathfrak{m}, j} \in \mathfrak{m}_Z$, $X_{\mathfrak{a}, j} \in \mathfrak{a}_Z$ and $X_{\mathfrak{h}, j} \in \mathfrak{h}$.
For $a\in A$ we have
\begin{equation} \label{decomp3} R_j = a^{\alpha_j} \operatorname{Ad}(a^{-1})(Q_j + R_j)  -
a^{2\alpha_j} Q_j\, ,\end{equation}
and inserting (\ref{decomp5}) for the last term
in (\ref{decomp3}) we obtain
\begin{equation} \label {decomp2}
\begin{aligned}
R_j+  a^{2\alpha_j} &\sum_{i=1}^n c_{ij}R_i
\\&=
a^{\alpha_j} \operatorname{Ad}(a^{-1}) \underbrace{(Q_j +R_j)}_{\in \mathfrak{k}}
-a^{2\alpha_j}(X_{\mathfrak{m}, j} + X_{\mathfrak{a}, j} + X_{\mathfrak{h}, j}) \,
\end{aligned}\end{equation}
Consider the $n\times n$-matrix
\begin{equation*
{\bf1}+ (c_{ij} a^{2\alpha_j})_{ij}.
\end{equation*}
It is invertible for $a\in A_t $ for sufficiently large $t$.
This proves (\ref{decomp}).
\end{proof}

In fact, a more precise version of the decomposition (\ref{decomp}) is obtained as follows.
Let $(X_{j})_j$, $(Y_{k})_k$
and $(Z_{l})_l$ be fixed bases for $\mathfrak{k}$, $\mathfrak{a}_Z$ and $\mathfrak{h}$,
respectively. Then the previous proof shows
(see (\ref{iota}) and (\ref{decomp2}))

\begin{lemma}\label{app2}
There exists $t>0$ such that for every $X\in\mathfrak{g}$ there exist
holomorphic functions
$f_j, g_k, h_l \in \mathcal{O}(D_{e^{-t}})$
such that
for all $a\in A_t$ one has
$$X=
\sum_{j=1}^{\dim\mathfrak{k}} f_j(\iota(a)) \operatorname{Ad}(a^{-1})X_{j} +
\sum_{k=1}^{\dim\mathfrak{a}_Z} g_k(\iota(a))  Y_{k} +
\sum_{l=1}^{\dim \mathfrak{h} } h_l(\iota(a) ) Z_{l}  .$$
Moreover, if $X\in\mathfrak{u}$, then this can be attained with
$$f_j, g_k, h_l \in \mathcal{O}_0(D_{e^{-t}})=\{\varphi\in \mathcal{O}(D_{e^{-t}})\mid \varphi(0)=0\}.$$
\end{lemma}

Let $(U_j)_j$ be a fixed homogeneous basis of $\mathcal{U}(\mathfrak{k})$,
$(V_k)_k$ one of $\mathcal{U}(\mathfrak{a}_Z)$ and $(W_l)_l$ of $\mathcal{U}(\mathfrak{h})$.
We then obtain the following from Lemma \ref{app1} and
Lemma \ref{app2}.

\begin{lemma}\label{app3} There exists $t>0$ such that
for every $u\in \mathcal{U}(\mathfrak{g})$ there exist $f_{j,k,l}\in \mathcal{O}(D_{e^{-t}})$
such that
\begin{equation} \label{deco}
u = \sum_{j,k,l}
f_{j,k,l} (\iota(a)) (\operatorname{Ad}(a^{-1}) U_j) V_k W_l ,\end{equation} for all $a\in A_t$,
where $j,k,l$ extend over all sets of indices with
$$\deg (U_j) + \deg (V_k) +\deg(W_l)\leq \deg(u).$$
\end{lemma}

\begin{proof} Let $t$ be as in Lemma \ref{app2}. We proceed by induction on $\deg(u)$,
the case of degree zero being clear.
Let $u\in \mathcal{U}(\mathfrak{g})$ and assume $u=Xv$ with $X\in\mathfrak{g}$ and $v\in \mathcal{U}(\mathfrak{g})$ of degree one less.
Applying Lemma  \ref{app2} to $X$, we see that
it suffices to treat the case where $X\in\mathfrak{a}_Z+\mathfrak{h}$. We write
$Xv=[X,v]+vX$ and apply the induction hypothesis to $[X,v]$.
For the term $vX$ we
apply the induction hypothesis to $v$ and write it as a
linear combination of elements $(\operatorname{Ad}(a^{-1}) U_j) V_k W_l$
with holomorphic coefficients and $\deg(U_j)+\deg(V_k)+\deg(W_l)\le\deg(v)$. Next we write
$$(\operatorname{Ad}(a^{-1}) U_j) V_k W_l X=(\operatorname{Ad}(a^{-1}) U_j) V_k([W_l,X]+XW_l).$$
The induction hypothesis applies to $V_k[W_l,X]$. The terms with $V_kXW_l$
already have the form asserted in (\ref{deco}).
\end{proof}

\begin{rmk}\label{rmk2} The decomposition (\ref{deco}) allows one to define the $\tau$-radial part of an
element $u\in \mathcal{U}(\mathfrak{g})^{M\cap H}$. This goes as follows:
We let $\mathcal{U}(\mathfrak{g})$ act on $C^\infty(G)$ by right
differentiation:
$$ (X\cdot f)(g) = \frac d{dt}\Big|_{t=0} f(g\exp(tX))\qquad (f\in C^\infty(G), g\in G, X\in \mathfrak{g})\, .$$
Henceforth we regard smooth functions on $Z$ as right $H$-invariant functions on $G$.
For a $K$-type $(\tau,W_\tau)$
we denote by $C^\infty(Z)_\tau$ the space of functions which are of left $K$-type $\tau$, for example the
matrix coefficients $m_{v,\eta}$ with $v$ in $V[\tau]$.
\par It is easily seen from the Peter-Weyl theorem
that for every $f\in C^\infty(Z)_\tau$ there exist a unique
$\Phi\in C^\infty(Z, W_\tau\otimes W_{\tau^*})$
such that
$$f(gH)= \operatorname{Tr} (\Phi(gH)),\quad gH\in Z,$$
where $\operatorname{Tr}$ stands for contraction $W_\tau\otimes W_{\tau^*}\to\mathbb{C}$.
Moreover $\Phi(kgH)=(1\otimes \tau^*(k))\Phi(gH)$ for $k\in K$.
In particular, $\Phi(aH)\in W_\tau\otimes W^{M\cap H}_{\tau^*}$ for $a\in A$.
Hence we have a finite sum
\begin{equation}\label{Fjsum}
f(gH)=\sum_j \langle F_j(gH),w_j\rangle,
\end{equation}
with $F_j\in C^\infty(Z,W_{\tau^*})$ and $w_j\in W_\tau$.
Moreover, $F_j(kgH)=\tau^*(k)(gH)$ for all $k\in K, g\in G$.

\par Let $u\in \mathcal{U}(\mathfrak{g})^{M\cap H}$ and let $t>0$
be as above. According to (\ref{deco}) we can write $u$ as a sum of
$(\operatorname{Ad}(a^{-1})U) V W$ with $U\in \mathcal{U}(\mathfrak{k})^{M\cap H}$, $V\in \mathcal{U}(\mathfrak{a}_Z)$ and $W\in \mathcal{U}(\mathfrak{h})$
and coefficients depending holomorphically on $\iota(a)$ for $a\in A_t$.
In order to compute $ u\cdot f|_{A_t} $ we may assume that $W={\bf1}$ as
$f$ is right $H$-invariant. For $u=(\operatorname{Ad}(a^{-1}))UV$ we then have
$$(u\cdot f)(a) = \sum  \langle  (V\cdot F_j)(a), U^t\cdot w_j\rangle $$
We finally arrive at an action of $\mathcal{U}(\mathfrak{g})^{M\cap H}$ on $C^\infty (A_t,
W_{\tau^*}^{M\cap H})$ given by
$${\rm rad}(u)(F)(a):= \sum \tau^*(U) (V\cdot F)(a)\,,
$$
a differential operator with $\operatorname{End}(W_{\tau^*}^{M\cap H})$-valued coefficients.
With this definition we find
$$uf(a)=\sum_j \langle\operatorname{rad}(u)(F_j)(a),w_j\rangle$$
when $f$ is given by (\ref{Fjsum}).

\par In particular for $f=m_{w,\eta}$
the functions $F_j \in C^\infty(A_t, W_{\tau^*}^{M\cap H})$ obtained as above
satisfy the system of differential equations
\begin{equation*
{\rm rad }(z)(F) = \chi_V(z) F
\qquad (z\in \mathcal{Z}(\mathfrak{g}))\end{equation*}
with $\chi_V: \mathcal{Z}(\mathfrak{g})\to \mathbb{C}$ the infinitesimal character of $V$.
\end{rmk}

In the decomposition (\ref{deco}) of $u\in \mathcal{U}(\mathfrak{g})$ we would like to restrict
the middle parts $V_k$ from $\mathcal{U}(\mathfrak{a}_Z)$ to a fixed finite set of elements, independent of
$u$. This will be done at the cost of enlarging the product with an extra factor from $\mathcal{Z}(\mathfrak{g})$
(which act by scalars on $V$).

For that we need to recall  parts of the  construction of the Harish-Chandra homomorphism.
The decomposition $\mathfrak{g}=\mathfrak{n} +\mathfrak{a} +\mathfrak{m} +\bar \mathfrak{n}$ results in a
direct sum decomposition
$$\mathcal{U}(\mathfrak{g}) = (\mathfrak{n}\, \mathcal{U}(\mathfrak{g}) + \mathcal{U}(\mathfrak{g})\bar\mathfrak{n}) + \mathcal{U}(\mathfrak{a}+ \mathfrak{m})$$
and allows for a linear projection
$$\mu_1: \mathcal{U}(\mathfrak{g}) \to \mathcal{U}(\mathfrak{a}+ \mathfrak{m})\, .$$
Recall (\cite{vdB} Lemma 3.6) that $\mu_1$ restricts to an algebra
homomorphism $\mathcal{Z}(\mathfrak{g})\to\mathcal{Z}(\mathfrak{m}+\mathfrak{a})$ and that for $z\in \mathcal{Z}(\mathfrak{g})$
of degree $d$ one has
$$z-\mu_1(z)\in \mathfrak{n} \,\mathcal{U}(\mathfrak{g})_{d-1},$$
(where $\mathcal{U}(\mathfrak{g})_{d-1}$ signifies elements of degree $\le d-1$).
It follows from the Harish-Chandra isomorphism theorem that
$\mathcal{Z}(\mathfrak{m}+\mathfrak{a})$ is finitely generated over $\mu_1(\mathcal{Z}(\mathfrak{g}))$.
More precisely (see \cite{vdB}, Lemma 3.7), there exist
elements $v_1,\dots,v_r\in\mathcal{Z}(\mathfrak{m}+\mathfrak{a})$ such that every
$v\in\mathcal{Z}(\mathfrak{m}+\mathfrak{a})$ can be written as
\begin{equation}\label{vsum}
v=\sum_{j=1}^r \mu_1(z_j)v_j
\end{equation}
with $z_j\in \mathcal{Z}(\mathfrak{g})$ and
$\deg(z_j)+\deg(v_j)\le\deg(v)$ for each~$j$.

Further, as $\mathfrak{a}$ and $\mathfrak{m}$ commute and as
$\mathfrak{a}=\mathfrak{a}_Z+(\mathfrak{a}\cap\mathfrak{h})$,
we obtain an algebra homomorphism
$$p: \mathcal{U}(\mathfrak{m} +\mathfrak{a}) \to \mathcal{U}(\mathfrak{a}_Z)$$ with
kernel $(\mathfrak{m}+(\mathfrak{a}\cap\mathfrak{h}))\,\mathcal{U}(\mathfrak{m}+\mathfrak{a})$.
We compose with $\mu_1$ and obtain a linear map
$$\mu_2: \mathcal{U}(\mathfrak{g})\to \mathcal{U}(\mathfrak{a}_Z)\, .$$
The restriction
to $\mathcal{Z}(\mathfrak{g})$ is an algebra homomorphism and
will be denoted by $\mu$. It follows from the above that
\begin{equation}\label{zdeco}
z-\mu(z)\in \mathfrak{n}\, \mathcal{U}(\mathfrak{g})_{d-1} + (\mathfrak{m}+(\mathfrak{a}\cap\mathfrak{h}))\,\mathcal{U}(\mathfrak{m} +\mathfrak{a})_{d-1}
\end{equation}
for $z\in\mathcal{Z}(\mathfrak{g})$.
Furthermore by applying $p$ to (\ref{vsum}) we see that
$\mathcal{U}(\mathfrak{a}_Z)$ is finitely generated over $\mu(\mathcal{Z}(\mathfrak{g}))$
with generators $p(v_j)\in \mathcal{U}(\mathfrak{a}_Z)$. Let $\mathcal{Y}$ denote the finite set
$$\mathcal{Y}=\{p(v_1),\dots,p(v_r)\}\subset\mathcal{U}(\mathfrak{a}_Z).$$

\par
\begin{lemma} \label{app4}For all $n\in \mathbb{Z}_{\ge 0}$ there exists $t=t_n>0$ such
that for all $u\in \mathcal{U}(\mathfrak{g})$ with $\deg(u)\le n$
there exist
$$U_j \in \mathcal{U}(\mathfrak{k}),
V_j \in \mathcal{Y}\subset\mathcal{U}(\mathfrak{a}_Z),  W_j \in
\mathcal{U}(\mathfrak{h}), z_j\in\mathcal{Z}(\mathfrak{g})$$
with
$$\deg U_j +\deg V_j + \deg W_j+\deg z_j \leq n,$$
and holomorphic functions $f_j \in \mathcal{O}(D_{e^{-t}})$ such that
\begin{equation}\label{sum with Z}
u=\sum_{j=1}^p f_j(\iota(a))
(\operatorname{Ad}(a^{-1})U_j)V_{j} W_j z_j\end{equation}
for all $a\in A_{t}$.
\end{lemma}

\begin{proof} Note the dependence of $t$ on the degree of $u$, contrary to
what was the case in Lemma \ref{app3}.
The proof will again be by induction on $n$, and the case $n=0$ is again clear.
Furthermore, we see from Lemma \ref{app3} that it suffices to establish
the decomposition for $u\in\mathcal{U}(\mathfrak{a}_Z)_n$, since if $\deg(U_j)>0$ or
$\deg(W_l)>0$ in (\ref{deco}) then $\deg(V_k)<n$ and the induction hypothesis applies.

Using (\ref{vsum})
we write $u=\sum \mu(z_{i}) Y_i$ with $z_{i}\in \mathcal{Z}(\mathfrak{g})$ and $Y_i\in\mathcal{Y}$,
such that
$\deg(z_{i}) + \deg(Y_i)\leq n$ for all $1\leq i \leq r$.
With $\mathfrak{n}=\mathfrak{u}+\mathfrak{n}_L$ and $\mathfrak{n}_L\subset \mathfrak{h}$
we obtain from  (\ref{zdeco}) for each $i$,
$$\mu(z_{i}) =z_{i} +  \sum_m \underbrace{ R_m}_{\in\mathfrak{u}} u_{mi}  +
\sum_{m'} \underbrace{S_{m'}}_{\in\mathfrak{m}} {u_{m'i}'} +
\sum_{m''} \underbrace{T_{m''}}_{\in \mathfrak{h}}
u_{m''i}''$$
with $u_{mi},u_{m'i}',u_{m''i}''\in\mathcal{U}(\mathfrak{g})$  all
of degree $< \deg z_{i}$.
Hence
$$u=\sum_i z_{i}Y_i +  \sum_{i,m}  R_m u_{mi}Y_i  +
\sum_{i,m'} S_{m'} u_{m'i}'Y_i +
\sum_{i,m''} T_{m''} u_{m''i}''Y_i.$$
The terms in the first sum already have the desired form.
The terms with $S_{m'}$ are dealt with directly by the induction hypothesis,
and the terms with $T_{m''}$ are dealt with similarly
after commuting $T_{m''}$ and $u_{m''i}''Y_i$.
Hence only the second sum remains, which is a sum of terms
$Rv$ with $R\in\mathfrak{u}$ and $v\in\mathcal{U}(\mathfrak{g})_{n-1}$.

We use Lemma \ref{app2} to decompose $R$ as a linear combination of
terms $\operatorname{Ad}(a^{-1})X$, with $X\in\mathfrak{k}$, and basis vectors of $\mathfrak{a}_z$ or $\mathfrak{h}$,
and with coefficients from $\mathcal{O}_0(D_{e^{-s}})$ (for some fixed $s>0$).
By the induction hypothesis applied to $v$,
the terms $(\operatorname{Ad}(a^{-1})X)v$ again have the desired form (\ref{sum with Z}).
Here $a\in A_{t_{n-1}}$.
We thus reach the conclusion that for each such $a$, our element
$u\in\mathcal{U}(\mathfrak{a}_Z)_n$ is of form (\ref{sum with Z})
plus elements of the form $g(\iota(a)) w$
with $g\in \mathcal{O}_0(D_{e^{-s}})$ and $w\in\mathcal{U}(\mathfrak{g})_n$ independent of $a$.
The same conclusion then applies to every element $u\in\mathcal{U}(\mathfrak{g})_n$,
as observed in the beginning of the proof.

\par We do the above for a basis $u_1, \ldots , u_N$ of $\mathcal{U}(\mathfrak{g})_n$ and
finally use the fact that an $N\times N$-matrix
of the form ${\bf1}_{N\times N} + F(z)$ with $F$ holomorphic and $F(0)=0$ is invertible for
sufficiently small $z$.
\end{proof}

We can now give the proof of (\ref{power}).
For $v\in V[\tau]$ we consider the matrix coefficient
$f(gH)=m_{v,\eta}(gH)=\eta(\pi(g)^{-1}v)$,
which we recall is a joint eigenfunction for $\mathcal{Z}(\mathfrak{g})$. As explained
in Remark \ref{rmk2}, $f$ is a sum of functions
of the form $\langle F(gH),w\rangle$ where $w\in W_\tau$ and where
$F\in C^\infty(Z, W_{\tau^*})$ is $\tau^*$-spherical, that is,
$F(kgH)=\tau^*(k)F(gH)$ for all $k\in K$, $g\in G$.
Note that this implies
$F(aH)\in W_{\tau^*}^{M\cap H}$ for $a\in A$.
Furthermore, the action of
elements from $\mathcal{U}(\mathfrak{g})^{M\cap H}$ is computed by taking radial parts,
and $F$ is a joint eigenfunction for $\mathcal{Z}(\mathfrak{g})$.

Let $\Pi\subset\mathfrak{a}^*$ denote the set of simple roots. For simplicity
we assume $\mathfrak{g}=[\mathfrak{g},\mathfrak{g}]$ in the remainder of
this section, so that
$\Pi$ is a basis for $\mathfrak{a}^*$.
Extension to the reductive case is
elementary.
Let $d=|\Pi|=\dim\mathfrak{a}$.

\newcommand{\mathcal{D}}{\mathcal{D}}
\begin{proposition}
Let $F\in C^\infty(Z, W_{\tau^*})$ be a $\tau^*$-spherical
joint eigenfunction for $\mathcal{Z}(\mathfrak{g})$.
There exist a neighborhood $\mathcal{D}$ of $0$ in $\mathbb{C}^d$, a number $M\in\mathbb{N}$, a finite set
$S\subset \mathbb{C}^d$,
and for each $s\in S$ and each multiindex $0\le |m|\le M$
a holomorphic $W_{\tau^*}^{M\cap H}$-valued function
$h_{s,m}$ on $\mathcal{D}$ such that
\begin{equation}
F(aH)=\sum_{s\in S}\sum_{0\le |m|\le M} z^s(\log z)^m h_{s,m}(z),
\end{equation}
for all $a\in A$ such that
 $z=(a^{\alpha})_{\alpha\in\Pi}\in \mathcal{D}$.
\end{proposition}

\begin{proof}
Let
$X_1,\dots,X_d$ be the basis
for $\mathfrak{a}$ which is dual to $\Pi$,
and
let $D$ be the maximal degree of the
the finite set of operators from $\mathcal{Y}\subset\mathcal{U}(\mathfrak{a}_Z)$.
For each multi-index $k=(k_1,\dots,k_d)$ with $|k|\le D$
we define
$$F_k(gH)=X_1^{k_1}\cdots X_d^{k_d}F(gH),\quad gH\in Z,$$
a $W_{\tau^*}$-valued function.
For each $j=1,\dots,d$ and each multiindex $l$ with $|l|\le D$, the function $X_j F_l$ can then be determined
by giving the monomial $X_jX_1^{k_1}\cdots X_d^{k_d}\in\mathcal{U}(\mathfrak{a})$ an
expression (\ref{sum with Z}) according to Lemma \ref{app4} with $t=t_{D+1}$,
and applying the theory of radial parts, as explained in
Remark \ref{rmk2}.
It follows that on $A_t$, the function $X_jF_l$ is
a combination of the $F_k$'s, with $|k|\le D$ and with
$\operatorname{End}(W_{\tau^*}^{M\cap H})$-valued coefficients depending on $a\in A_t$.

We thus see that the vector-valued function $\mathbf F(a)=(F_k(aH))_k$
satisfies a first order ordinary differential system
$$X_j\mathbf F (a)= M(a) \mathbf F(a)$$
on $A_t$ with an operator $M(a)$. Furthermore, it follows from
(\ref{sum with Z}) that $M(a)$ extends to a holomorphic function of
 $z=(a^{\alpha_j})_{j=1,\dots,n}\in\mathbb{C}^n$
in a neighborhood of $0$. In particular, $M(a)$ is then a
holomorphic function of the smaller tuple
 $(a^{\alpha})_{\alpha\in\Pi}\in\mathbb{C}^d$
in a neighborhood of $0$, since every positive root is a sum of simple roots.
With coordinates on $A$ given by
 $(a^{\alpha})_{\alpha\in\Pi}\in\mathbb{C}^d$
we thus obtain a system of equations,
which has a simple singularity at $0$ according to
\cite{Kn}, p 702, Example 2. The proposition now follows from
 \cite{Kn}, Thm.~B.16.
\end{proof}

By taking inner products with elements from $W_\tau$ we obtain
a similar expansion of the matrix coefficients $m_{v,\eta}(a)$ for $a\in A_t$,
for all $v\in V[\tau]$.

\begin{cor}
Let $\eta\in (V^{-\infty})^H$ and $v\in V[\tau]$.
There exist a neighborhood $\mathcal{D}$ of $0$ in $\mathbb{C}^d$, a number $M\in\mathbb{N}$, a finite set
$S\subset \mathbb{C}^d$,
and for each $s\in S$ and each multiindex $0\le |m|\le M$
a holomorphic function
$f_{s,m}$ on $\mathcal{D}$ such that
\begin{equation}
m_{v,\eta}(a)=\sum_{s\in S}\sum_{0\le |m|\le M} z^s(\log z)^m f_{s,m}(z),
\end{equation}
for all $a\in A$ such that
 $z=(a^{\alpha})_{\alpha\in\Pi}\in \mathcal{D}$.
\end{cor}

Let now $X\in\mathfrak{a}$ be such that $\alpha(X)\in\mathbb{Z}_-$
for each $\alpha\in\Pi$ and consider $m_{v,\eta}(a_s)$ for $s\to\infty$,
where $a_s=\exp(sX)$. Note that  $(a_s^{\alpha})_{\alpha\in\Pi}\in \mathcal{D}$
for $s\gg 0$ and that
 $a_s^{\alpha}$ is a positive integral
power of $e^{-s}$ for each $\alpha\in\Pi$.
Hence $f_ {s,m}( (a_s^{\alpha})_{\alpha\in\Pi})$ extends to a holomorphic
function of $\zeta=e^{-s}$.
This completes the proof of
(\ref{power}) and hence of
Lemma \ref{lem1}.