The diagonal slice of Schottky space

Abstract: An irreducible representation of the free group on two generators X,Y into SL(2,C) is determined up to conjugation by the traces of X,Y and XY. We study the diagonal slice of representations for which X,Y and XY have equal trace. Using the three-fold symmetry and Keen-Series pleating rays we locate those groups which are free and discrete, in which case the resulting hyperbolic manifold is a genus-2 handlebody.We also compute the Bowditch set, consisting of those representations for which no primitive elements in the group generated by X,Y are parabolic or elliptic, and at most finitely many have trace with absolute value at most 2. In contrast to the quasifuchsian punctured torus groups originally studied by Bowditch, computer graphics show that this set is significantly different from the discreteness locus.

\section{Introduction} \label{intro}

It is well known that an  irreducible representation of the free group $F_2$ on two generators $X,Y $ into $SL(2,\mathbb C)$ is determined up to conjugation by the traces of $X,Y$ and $XY$. More generally, if we take the GIT quotient of   all (not necessarily irreducible) representations, then the resulting $SL(2,\mathbb C)$ character variety of $F_2$ can be identified with $\mathbb C^3$ via these traces, see for example \cite{goldman2} and the references therein.   If the representation is free, discrete, purely loxodromic and geometrically finite, the resulting manifold is a genus-$2$ handlebody.  The collection of all such representations is known as \emph{Schottky space}, denoted  $\mathcal {SCH}$. It is a consequence of Bers' density theorem that $\mathcal {SCH}$ is the interior of the discreteness locus, see for example~\cite{canary}. It is natural to ask, for which values of   $x = \Tr X, y = \Tr Y, z= \Tr XY$  is the corresponding  representation in  $\mathcal {SCH}$?

Let ${\cal P}$ denote the set of primitive elements in $F_2$. For $(x,y,z) \in \mathbb C^3$, let $\rho_{(x,y,z)}$ denote a choice of representation $  F_2 \to SL(2,\mathbb C)$ in the conjugacy class determined by the trace triple.
The \emph{Bowditch set} (or $BQ$-set) $\mathcal B$  is defined in ~\cite{tan_gen}  as the set of $(x,y,z) \in \mathbb C^3$ corresponding to irreducible representations for which
\begin{equation*}
 \begin{split}  & \Tr \rho_{(x,y,z)}(g) \notin [-2,2]  \ \ \forall g \in {\cal P}    \ \ \mbox {\rm and}   \ \ \cr
& \{ g \in {\cal P}: |\Tr \rho_{(x,y,z)}(g)| \leq 2 \} \ \mbox {\rm is finite}.\end{split}\end{equation*}
(The exceptional case in which
$\Tr [X,Y] = 2$ corresponds to reducible representations and is excluded from the discussion, see Remark~\ref{reducible}.)
The Bowditch set is open and ${\rm Out}(F_2)$ acts properly discontinuously on it. Clearly ${\mathcal{SCH}} \subset \mathcal B$.

Bowditch's original work~\cite{bow_mar} was on the case in which  the commutator $[X,Y] = XYX^{-1}Y^{-1}$
is parabolic and $\Tr  [X,Y] =  -2 $. He conjectured that the subsets of  ${\mathcal{SCH}}$ and $ { \mathcal B}$ corresponding to this restriction  coincide.
Although this has not been proven, computer pictures  indicate his conjecture may well be true

In this paper we restrict to the special case in which $x = y = z$, which we call the \emph{diagonal slice} of the character variety, denoted $  {\Delta}$ and parametrized by the single complex variable $x$.  We show that in this slice, the analogue of Bowditch's conjecture is far from being true. This is illustrated in Figure~\ref{Diagonal-and-Riley-mu-3-Ray-BQ} which compares the  intersections of $\Delta$ with ${\mathcal{SCH}}$ and ${ \mathcal B}$.
The discreteness locus is the outer region foliated by rays; these are the Keen-Series pleating rays which relate to the geometry of the convex hull boundary as explained in Section~\ref{sec:pleating} and whose closure is known to be $\overline{ {\Delta} \cap {\mathcal{SCH}}}$, see Theorem~\ref{thm:density}.
The Bowditch set, by contrast, is the complement of the black part.
It  is clear   that ${ \mathcal B} \cap {\Delta}$ contains a large open region not in   $ {\Delta} \cap {\mathcal{SCH}}$, and also has different symmetries.
In particular,   it is not hard to show that  the interval $  (2,3)$ is contained in $   { \mathcal B} \setminus {\mathcal{SCH}}$, see the discussion  in Section~\ref{sec:algorithm1}.

The main content of this paper is an explanation and justification of  how these plots were made, in particular to explain how we enumerated and computed the pleating rays for the symmetric genus $2$ handlebody corresponding to the trace triple $(x,x,x)$.

\begin{figure}[hbt]\includegraphics[width=6cm]{Figs/Diagonal-and-Riley-mu-3-Ray-BQ}
\caption{Superposition of the discreteness locus for $\pi_1({\cal H})$ and the Bowditch set in the $x$-plane. The Bowditch set for the $(x,x,x)$-triple is the complement of the central black region, while the discreteness locus is the closure of the region foliated by rays. The rays  are actually computed as the pleating rays for the quotient orbifold ${\mathcal S}$.}\label{Diagonal-and-Riley-mu-3-Ray-BQ}
\end{figure}

To compute the Bowditch set ${ \mathcal B}$ we use an algorithm based on the ideas in~\cite{bow_mar} and developed further  in~\cite{tan_gen}. This is explained
in Section~\ref{sec:algorithm}.

The discreteness problem is tackled as follows. If $(x,x,x) \in {\mathcal{SCH}}$ then the quotient $3$-manifold $\mathbb H^3/G$ is a handlebody ${\cal H}$ with order $3$ symmetry. We use the symmetry to  reduce  the problem of finding $ {\Delta} \cap {\mathcal{SCH}}$ to a problem very similar to that of determining the so-called \emph{Riley slice  of Schottky space}. This is actually a space of groups on the boundary of $ {\mathcal{SCH}}$, consisting of those free, discrete and geometrically finite groups for which  the two  generators $X,Y$ are parabolic, thus contained in the slice $(2,2,z) \subset \mathbb C^3$. The corresponding  manifold is a handlebody whose conformal boundary is a sphere with four parabolic points.   The problem of finding those $z$-values for which such a group is free, discrete and geometrically finite was  solved using the method of pleating rays in~\cite{ksriley}.  In the present case, the quotient of ${\cal H}$ by the symmetry is an orbifold ${\mathcal S}$ with two order $3$ cone axes, whose conformal boundary is a sphere with four order $3$ cone points. Thus similar methods enable us to find   $ {\Delta} \cap {\mathcal{SCH}}$ here.

Although Figure~\ref{Diagonal-and-Riley-mu-3-Ray-BQ} shows that in $\Delta$, the analogue of Bowditch's conjecture fails since ${ \mathcal B}$  and the interior of the discreteness locus are plainly distinct, in many other  slices, see for example Figure~\ref{Figs/Riley-Ray-BQ},  the (modified) Bowditch set  and the interior of the discreteness locus appear to coincide. This  is connected to the dynamics of the action of a suitable mapping class group on representations and raises many interesting questions which we hope to address elsewhere.

The plan of the paper is as follows.
We begin in Section~\ref{sec:markoff} with a discussion of the Markoff tree and the algorithm used to compute the Bowditch set.
In Section~\ref{sec:basicconfig} we introduce a basic geometrical construction  which conveniently
encapsulates the $3$-fold symmetry. The quotient of the original   handlebody  ${\cal H}$ by the symmetry is a ball with two order $3$ cone axes. This orbifold ${\mathcal S}$ has a further $4$-fold symmetry group whose quotient is again a topological ball.  Our construction allows us to write down specific $SL(2,\mathbb C)$ representations of all the groups involved with ease.
In Section~\ref{sec:discrete} we turn to the discreteness question. After reducing the problem to one on ${\mathcal S}$,
we briefly review material from the Keen-Series theory of pleating rays and
recall  what is needed from~\cite{ksriley}, allowing us to
apply a similar proof in the present context.  Section~\ref{sec:torustree}, not strictly logically necessary for our development, explains how we did our trace computations in practice, by relating the problem to one on a commensurable torus with a single cone point of angle $4\pi/3$.

\section{The Markoff tree and the Bowditch set} \label{sec:markoff}
Let $A =  \begin{pmatrix} a & b \cr c & d\end{pmatrix} \in SL(2,\mathbb C)$ so that $ad-bc= 1$. As usual we define its trace $\Tr A = a+d$.

Let $F_2 = \langle X,Y|  \  \rangle$ be the free group on two generators.
It is well known that a representation $\rho{\colon \thinspace} F_2 \to SL(2,\mathbb C)$ is determined up to conjugation (modulo taking the GIT quotient under the conjugation action, see~\cite{goldman2})
by the three traces $ x = \tr X, y = \tr Y, z = \tr XY$. In fact,
given $x,y,z \in \mathbb C$   we can define a representation $\rho_{x,y,z}{\colon \thinspace} F_2  \to SL(2,\mathbb C) $ by
$\rho(X) = \begin{pmatrix} x & 1 \cr -1 & 0 \end{pmatrix}, \ \rho(Y) = \begin{pmatrix} 0 & \xi \cr -\xi^{-1} & y \end{pmatrix}$ where $ z = -(\xi + \xi^{-1})$.  Clearly with this definition, $ \Tr X = x, \Tr Y = y$ and $\Tr XY = z$.

\subsection{The Markoff Tree}
For matrices $U,V \in SL(2,\mathbb C)$  set $u = \Tr U, v = \Tr V, w = \Tr UV$. Recall the trace relations:
\begin{equation}\label{eqn:inverse} \Tr UV^{-1} = uv-w \end{equation} and
\begin{equation} \label{eqn:commreln}  u^2+v^2+w^2 = uvw + \Tr {[U,V]} +2.
\end{equation}
Setting $\mu =  \Tr {[U,V]} + 2$, this last equation takes the form
$$u^2+v^2+w^2 - uvw = \mu.$$

Let $F_2 = \langle X,Y| \ \  \rangle$ as above. An element $U \in F_2$ is \emph{primitive} if it is a member  of a generating pair; we denote the set of all primitive elements by ${\cal P}$.  The conjugacy classes of primitive elements are enumerated by $\hat \mathbb Q = \mathbb Q \cup \infty$ and are conveniently organised relative to the Farey diagram ${  \mathcal F}$ as shown in Figure~\ref{fig:farey}.
This consists of  the images of the ideal triangle with vertices at $1/0,0/1$ and $1/1$ under the action of $SL(2,\mathbb Z)$ on the upper half plane, suitably conjugated to the position shown in the disk. The label $p/q$ in the disk is just the conjugated image of the actual point $p/q \in \mathbb R$.

\begin{figure}[ht]
\includegraphics[width=5.5cm]{Figs/fareydiagram.pdf}
\hspace{1cm}
\includegraphics[width=5.5cm]{Figs/fareywords.pdf}
\caption{The Farey diagram, showing the arrangement of rational numbers on the left with the corresponding  primitive words on the right.}\label{fig:farey}
\end{figure}

 Since the rational points are precisely the images of $\infty$ under $SL(2,\mathbb Z)$, they correspond bijectively to  the vertices of ${  \mathcal F}$.
A pair  $p/q , r/s \in \hat \mathbb Q$ are the endpoints of an edge  if and only if $pr-qs = \pm 1$; such pairs are called \emph{neighbours}.
A triple of points in $\hat \mathbb Q$ are the vertices of a triangle precisely when they are the images of the vertices of the initial triangle $(1/0,0/1,1/1)$; such triples are always of the form
$(p/q , r/s,(p+r )/( q+s))$ where  $p/q , r/s$ are neighbours.
In other words,  if $p/q , r/s$ are the endpoints of an edge, then the vertex of the triangle on the side away from the centre of the disk is found by `Farey addition' to be $(p+r )/( q+s)$. Starting from $1/0 = -1/0=  \infty$ and $0/1$, all   points in $\hat \mathbb Q$ are obtained recursively  in this way. (Note we need to start with $-1/0=  \infty$  to get the negative fractions on the left side of the left hand diagram in Figure~\ref{fig:farey}.)

The right hand picture in Figure~\ref{fig:farey}
 shows a corresponding arrangement of primitive elements in $F_2$, one in each conjugacy class, starting with initial triple $(A,B, AB)$. Each vertex is labelled by a certain cyclically shortest representative of the corresponding word. Pairs of primitive elements form a generating pair if and only if they are at the endpoints of an edge. Triples at the vertices of a triangle correspond to a generator triple of the form $(U,V, UV)$.  Corresponding to the process of Farey addition, successive words can be found by juxtaposition as indicated on the diagram. Note that for this to work it is important to preserve the order:   if $U, V$ are the endpoints of an edge  with $U$ before $V$ in the anti-clockwise order round the circle, the correct  concatenation is $UV$. Note also that the words on the left side of the diagram involve $B^{-1}A$ corresponding to starting with $\infty = -1/0 $. We denote the particular representative of the conjugacy class corresponding to $p/q \in \hat \mathbb Q$ found by  concatenation by  $W_{p/q}$. Its word length in the generators $A,B$ is a function $F(p/q)$ of $p/q$. A function on $\hat \mathbb Q$ is said to have \emph{Fibonacci growth} if it is comparable with  uniform upper and lower bounds to $F$.

 \begin{figure}[ht] \includegraphics[width=7cm]{Figs/markofftree.pdf}
\vspace{-0.5cm}
\caption{The Markoff tree used to compute traces with an initial triple $(u,v,w)$. }\label{fig:markofftree}
\end{figure}

In this paper we  are largely interested in computing traces of primitive elements.  Following Bowditch~\cite{bow_mar}, these can also be easily computed  by using the  trivalent  tree $\mathbb T$ dual to ${  \mathcal F}$, see the left frame of Figure~\ref{fig:farey}, and Figure~\ref{fig:markofftree}. Let $\Omega$ denote the set of complementary regions of $\mathbb T$, abstractly, a complementary region is the closure of a connected component of the complement of $\mathbb T$. As is apparent from Figure~\ref{fig:farey},   there is a bijection between $\Omega$ and the set of vertices of ${  \mathcal F}$.   Thus the set $\Omega$ can be identified with conjugacy classes of primitive elements and hence with $\hat \mathbb Q$.

Given a representation $\rho {\colon \thinspace}  F_2 \to SL(2,\mathbb C)$, each $U \in \Omega$ is labelled by $u=\Tr \rho(U)$, the trace of the corresponding generator, as shown in
Figure~\ref{fig:markofftree}. Labels on opposite sides of an edge correspond to traces of a generator pair: the three labels round a vertex correspond to a generator triple $(U,V,UV)$.  Crossing an edge adjacent to regions $U,V$  of ${  \mathcal F}$ corresponds to changing the  generator triple from $(U,V,UV)$ to  $(U,V,UV^{-1})$.

Suppose that $(U,V,W)$ are
  the labels of regions round a vertex with $ u = \Tr \rho(U)$, $v = \Tr \rho(V)$, $w = \Tr \rho(W)$. By~\eqref{eqn:commreln}
we have
$u^2+ v^2 + w^2 -uvw = \mu$. By~\eqref{eqn:inverse}, the two vertices  opposite the ends of an edge labelled $(U,V)$ have labels  $w, uv-w$ respectively. More precisely, crossing the $3$ edges of a triangle of $\mathcal F$  gives rise to the three basic moves
 $(u,v, w) \to (u,v, uv-w)$,  $(u,v, w) \to (u,uw-v, w)$,  $(u,v, w) \to (vw-u,v, w)$
which  generates  traces of all possible elements in $\Omega$ (and hence $\mathcal P$).  Note that any of these three moves  leaves $\Tr {[U,V]}$ and hence $\mu$ invariant; in other words, $\mu$ is an invariant of the tree. Bowditch's original paper was mostly confined to the case $\mu = 0$.

In this way, the Markoff tree provides a fast way to compute traces of  elements  in ${\cal P}$ starting from an initial triple $(u,v,w)$.   This is illustrated in Figure~\ref{fig:trace2}  with the initial triple $( \sqrt{x+1}, 0, \sqrt{-x+2})$ which is used in Section~\ref{sec:torustree1}. We denote the tree of traces associated to an initial triple $(u,v,w)$ by $\mathbb T_{(u,v,w)}$.
Later we will use a variant of this construction to compute traces of curves on a four pointed sphere, see Section~\ref{sec:traces}.

\begin{figure}[hbt]    \begin{center}\vspace{-.5cm}
 \includegraphics[height=7cm]{Figs/fareytraces}  \vspace{-.8cm}
 \caption{The Farey tessellation used to  compute  traces. See Section~\ref{sec:torustree1} for a discussion of the choice of sign of the square roots.}\label{fig:trace2}
\end{center}
\end{figure}

 \subsection{The Bowditch set}
\label{sec:tree}

It is convenient to rephrase the above discussion using the terminology introduced in~\cite{bow_mar}.
As above, let $\Omega$ denote the set of complementary regions of the tree $\mathbb T$.
Define a \emph{Markoff map} to be a map $\phi: \Omega \to \mathbb C$ such that $\phi$ satisfies the trace relations~\eqref{eqn:inverse} and~\eqref{eqn:commreln}.   The set of all Markoff maps is denoted $\Phi$.
 Since traces depend only on conjugacy classes,
a representation $\rho {\colon \thinspace} F_2 \to SL(2,\mathbb C)$ defines a Markoff map by setting $\phi(U) = \Tr \rho(U)$ for $U \in \Omega$.  Fixing once and for all an identification of $\Omega$ with $\hat \mathbb Q$ (and recalling that $\Omega$ is identified with conjugacy classes of elements in ${\cal P}$), we have $\phi(p/q) = \Tr \rho(W_{p/q}), p/q \in \hat \mathbb Q$, where $W_{p/q}$ is the special word in the conjugacy class corresponding to $p/q \in \Omega$.

Thus as explained above,  using the trace relations~\eqref{eqn:inverse} and~\eqref{eqn:commreln}, an initial triple $(x,y,z) \in \mathbb C^3$ uniquely determines a Markoff map $\phi = \phi_{x,y,z}$ together with a corresponding labelling of $\mathbb T$.  Conversely a Markoff map $\phi \in \Phi$ determines $(x,y,z)\in \mathbb C^3$ by setting $ x = \phi(0/1),  y = \phi(1/0), z = \phi(1/1)$. In this way, we can identify $\Phi$ with $\mathbb C^3$.
For $\phi \in \Phi$, denote the corresponding tree $\mathbb T_{\phi} = \mathbb T_{(\phi(0/1), \phi(1/0),\phi(1/1))}$.

The Bowditch set $\mathcal B$ is the set of all  $\phi \in \Phi$ with $\mu \neq 4$ which satisfy the following conditions:
\begin{equation}\label{eqn:B1}   \phi (U)  \notin [-2,2]  \ \ \forall  \ U \in \Omega  \ \ \ \mbox {\rm {and}} \end{equation}
\begin{equation} \label{eqn:B2}   \{U \in \Omega:  | \phi (U)| \leq 2 \}   \ \ \mbox {\rm {is finite}}.
\end{equation}
The Bowditch set $\mathcal B$ is open in $\mathbb C^3$ and ${\rm Out}(F_2)$ acts properly discontinuously on $\mathcal B$. Furthermore, if $\phi \in \mathcal B$, then $\log^+|\phi(U)|={\rm max}\{0, \log|\phi(U)|\}$ has Fibonacci growth on $\Omega$ (see \cite{tan_gen}).

\begin{remark}\label{reducible} \rm{The maps $\phi$ for which  $\mu=4$ correspond  to the reducible representations: our definition above automatically excludes them from ${ \mathcal B}$.  For such $\phi$, there are infinitely many $U \in \Omega$ such that $|\phi(U)|<m$  for $m>2$, they can alternatively be excluded from ${ \mathcal B}$ by relaxing condition (\ref{eqn:B2}) to the condition  that $\{U \in \Omega:  | \phi (U)| \leq 2+\epsilon \}   $ be finite for any $\epsilon >0$. }
\end{remark}

\subsubsection{Background to the algorithm}\label{sec:algorithm}
Our algorithm for computing which points lie in ${ \mathcal B}$ is based on results from~\cite{bow_mar, tan_gen} which we summarise here. We consider only $\phi$ for which $\mu \neq 4$.
Following Bowditch~\cite{bow_mar}, we orient the edges of $\mathbb T_{\phi}$ in the following way. Suppose that labels of the regions adjacent to some edge $e$ are $u,v$ and the  labels of the two  remaining regions at the two end vertices are $w,t$, see Figure~\ref{fig:markofftree}. From the trace relations, $t = uv-w$.
Orient $e$ by putting an arrow from $t$ to $w$ whenever $|t| > |w|$ and vice versa. If both moduli are equal,  make either choice; if the inequality is strict, say that the edge is \emph{oriented decisively}.

A  \emph{sink region} of $\mathbb T_{\phi}$ is a connected non-empty subtree $T$ such that  the arrow on any edge not in $T$ points towards $T$ decisively. A sink region may consist of a single \emph{sink vertex} $v$ (the three edges adjacent to $v$ point towards $v$) and no edges. Clearly a sink region is not unique: one can always add further vertices and edges around the boundary of $T_{\phi}$.

For any $m \geq 0$ and $\phi \in \Phi$ define $\Omega_{\phi}(m) = \{ U \in \Omega | |\phi( U)| \leq m\}$.
The following  lemmas from~\cite{tan_gen} show that $\Omega_{\phi}(2)$ is connected, and that from any initial vertex not adjacent to regions in $\Omega_{\phi}(2)$, the arrows determine a descending path
through $\mathbb T$  which either runs into a sink, or meets vertices adjacent to regions in $\Omega_{\phi}(2)$. Furthermore, if $\phi(U)$ takes values away from the exceptional set $E = [-2,2] \cup \{\pm \sqrt{\mu}\} \subset \mathbb C$, then there exists a finite segment of $\partial U$ such that the edges adjacent to $U$ not in this segment are directed towards this segment.

\begin{lemma}[\cite{tan_gen} Lemma 3.7]\label{forkvertex}
Suppose $U,V,W \in \Omega$ meet at a vertex $v$  with the arrows on both the edges adjacent to $U$ pointing away from $v$. Then either $|\phi(U)| \leq 2$ or $\phi(V) = \phi(W) = 0$. \end{lemma}

\begin{corollary}[\cite{tan_gen} Theorem 3.1(2)]\label{connected}
Let $\phi \in \Phi$. Then  $\Omega_{\phi}(2)$ (more generally, $\Omega_{\phi}(m)$ for $m \ge 2$) is connected.
\end{corollary}

\begin{lemma}[\cite{tan_gen} Lemma 3.11 and following comment] \label{infiniteray}
Suppose $\beta$ is an infinite ray consisting of a sequence of edges of $\mathbb T_{\phi}$ all of whose arrows point away from the initial vertex. Then $\beta$ meets at least one region $U \in \Omega$  with $|\phi( U)| < 2$. Furthermore, if the ray does not follow the boundary of a single region, it meets infinitely many  regions with this property.
\end{lemma}

\begin{lemma}[\cite{tan_gen} Lemma 3.20] \label{finiteboundary}
Suppose that $\phi(U) \notin E$ and consider the regions $V_i, i \in \mathbb Z$ adjacent to $U$ in order round $\partial U$. Then away from a finite subset, the values $|\phi(V_i)|$  are increasing and approach  infinity as $ i \to \infty$ in both directions.   Hence there exists a finite segment of $\partial U$ such that the edges adjacent to $U$ not in this segment are directed towards this segment.
 \end{lemma}

We remark that if $\phi(U)=  \pm \sqrt{\mu}$  and $\sqrt{\mu} \not\in [-2,2]$, then the values of $|\phi(V_i)|$ in Lemma \ref{finiteboundary} approach zero in one direction round $\partial U$ (\cite{tan_gen} Lemma 3.10) and hence $\phi \not\in { \mathcal B}$ since condition (\ref{eqn:B2}) will not be satisfied. Hence, for $\phi \in { \mathcal B}$, $\phi(U) \not\in E$ for all $U \in \Omega$.

The set $\Omega_{\phi}(2)$ can be used to construct a sink region $T$ which is finite if and only if $\phi \in \mathcal B$. Essentially, if $\phi \in \mathcal B$, then $T$ consists of finite segments of the boundaries of the (finite number of) elements of $\Omega_{\phi}(2)$. These are the segments alluded to in Lemma \ref{finiteboundary}; they have to be large enough so the conclusion of the lemma holds, and also to contain all edges adjacent to $U,V$ with $U,V \in \Omega_{\phi}(2)$ so that the union is connected. To do this, an explicit function
  $H_{\mu}:\mathbb C \rightarrow \mathbb R^+ \cup \{\infty\}$ is constructed (see~\cite{tan_gen} Lemma 3.20, the following remark and Lemma 3.23)   as follows:
\begin{enumerate}
\item If $x \in E$, define $H_{\mu}(x) =\infty$;
\item For $x \not\in E$, let $x=\lambda +\lambda^{-1}$ with $|\lambda|>1$ (note that $|\lambda| \neq 1$ since $x \not\in [-2,2]$). Define
\begin{equation}\label{eqn:Hx}
H_{\mu}(x) =\max \left\{ 2, \sqrt{\left|\frac{x^2-\mu}{x^2-4}\right|}\frac{2|\lambda|^2}{|\lambda|-1} \right\} .
\end{equation}
\end{enumerate}
Then $H_{\mu}$ is continuous on $\mathbb C \setminus E$. Now we can define  a specific attracting subtree:

\begin{definition} \label{sinkregion} Given $\phi \in \Phi$, let $T$ be the subset of $\mathbb T_{\phi}$ defined as follows:
\begin{enumerate}
\item An edge with adjacent regions $U,V$  is in $T$ if and only if either $|\phi(U)| \leq 2$ and $|\phi(V)| \leq H_{\mu}(\phi(U))$, or vice versa.
\item Any sink vertex is in $T$, as are any vertices which are the end points of two edges in $T$.
\end{enumerate}
\end{definition}

Based on the above lemmas, we have the following theorem (see also the special properties of the function  $H_{\mu}$  and Lemmas 3.21-3.24 in~\cite{tan_gen}).

\begin{theorem} Given $\phi \in \Phi$ (with $\mu \neq 4$), the set  $T$   in Definition \ref{sinkregion} is a non-empty, connected subtree of
$\mathbb T_{\phi}$. Moreover  $T$ is a sink region for $\mathbb T_{\phi}$, that is, all edges not in $T$ are directed decisively towards $T$. Furthermore, $T$ is finite if and only if $\phi \in \mathcal B$.
\end{theorem}

\subsubsection{The algorithm}\label{sec:algorithm1}

Based on the above discussion, our algorithm to decide whether or not $\phi \in { \mathcal B}$ is as follows.
\begin{description}
\item [Step 1] Starting at any vertex, follow the direction of decreasing arrows. On reaching  a  sink vertex,
stop. This vertex is in $T$ by Definition~\ref{sinkregion}. If the input is ${ \mathcal B}$, then this method always finds a sink vertex in finite time because there is a finite sink region.  Otherwise, the process may not terminate in (pre-specified) finite time and the algorithm is indecisive.
\item [Step 2]  Assuming a stopping point is found in Step 1, starting from this point, search outwards by a depth first search using Definition~\ref{sinkregion} to identify whether or not an edge is in $T$. This works because of the connectedness of $T$. If this search terminates in (pre-specified) finite time then  $\phi_{x,y,z} \in { \mathcal B}$. Otherwise, the algorithm is indecisive.
\end{description}

Note that  if the starting point is a sink vertex and the three adjacent edges are not in $T$, then $T$ consists of just the sink vertex  by the connectedness of $T$, hence  $\phi_{x,y,z} \in { \mathcal B}$. This occurs for example for the tree $\mathbb T_{(x,x,x)}$ with $x \in (2,3)$.
\begin{comment}
The algorithm is as follows. It  is justified by Lemmas 3.7, 3.11, 3.21, 3.22 in \cite{tan_gen}.

\begin{description}
\item [Step 1] Starting at any vertex, follow the direction of decreasing arrows. On reaching a vertex with two decreasing arrows (a \emph{fork vertex})
or a vertex with no decreasing arrows (a \emph{sink vertex}),
stop. This vertex is in $T$  since neither a fork nor sink can be contained in the complement of $T$, we are at a vertex in $T$.
\item [Step 2] We now need to check whether or not $T$ is finite.  Starting from the stopping point found in Step 1, search outwards by a depth first search  using Lemma~\ref{test} to identify whether or not an edge is in $T$. If this search terminates then  $\phi_{(x,y,z)} \in { \mathcal B}$.
Otherwise, the process may not terminate in finite time. Note that in the case where the starting point is a sink, if the three edges adjacent to the sink are not in $T$, then $T$ consists of just the sink and  $\phi_{(x,y,z)} \in { \mathcal B}$. This occurs for example for the tree $\mathbb T_{(x,x,x)}$ with $x \in (2,3)$.
\end{description}

\begin{remark}
 \rm{ In the actual program we used, rather than stopping at a fork vertex, we continue to follow any decreasing arrows until we reach a sink vertex, because in practice, it is easier to write a program which does this.
If the input is ${ \mathcal B}$, then this method always finds a sink vertex
in a finite time because there is a finite sink region.
}
\end{remark}
\end{comment}

Figure~\ref{Diagonal-BQ} shows the Bowditch set in the diagonal slice $\Delta$ as determined by this algorithm.

\begin{remark}
 \rm{We do not have an algorithm whose output is  $\phi_{x,y,z} \notin { \mathcal B}$. When $\mu=0$, it was shown in \cite{NT} that if $|\phi(U)|\le 0.5$ for some $U \in \Omega$, then  $\phi_{x,y,z} \notin { \mathcal B}$. Hence in Step 1 above, if $\mu=0$, we can stop when we hit a region satisfying this condition and conclude that  $\phi_{x,y,z} \notin { \mathcal B}$. Using the same methods, a similar upper bound can be found for $\mu$ close to $0$. In particular, there is a neighbourhood of $(0,0,0)$ which is disjoint from ${ \mathcal B}$, as clearly illustrated in Figure ~\ref{Diagonal-BQ}. However, as shown in \cite{GMST}, no such universal positive bound exists for all $\mu$: precisely, for any $\epsilon >0$ and $\mu>4$, there exist  $\phi \in { \mathcal B}_{\mu}$ and $U \in \Omega$ such that $|\phi(U)|<\epsilon$.  Another issue is that the sink region may be extremely large so may not be detected in a program with a given finite number of steps, this occurs when we approach the boundary of ${ \mathcal B}$. Thus the  algorithm is not completely decisive although it appears to give nice results. In particular, there may be false negatives; however points which are determined to be in ${ \mathcal B}$ are correctly marked. }
\end{remark}

  \begin{figure}[ht]
  \includegraphics[width=7cm]{Figs/Diagonal-BQ}
\caption{The Bowditch set  ${ \mathcal B}$ for the Markoff maps $\phi_{(x,x,x)}$, plotted in the $x$-plane. The coloured (grey) points are in  ${ \mathcal B}$ and the black ones are outside. The colours (shades of grey) indicate the size of the sink region $T$.} \label{Diagonal-BQ} \end{figure}

\section{Groups, manifolds, symmetries and quotients} \label{sec:basicconfig}

In this section we detail a construction which allows us conveniently to exploit the three-fold symmetry of groups in the diagonal slice $\Delta$.
As is well known, if the image of a representation $\rho {\colon \thinspace} F_2 \to SL(2,\mathbb C)$ is free and discrete  then $\mathbb H^3/\rho(F_2)$ is a genus two handlebody ${\cal H}$,  see~\cite{hempel} Theorem 5.2. (Note that a hyperbolic $3$-manifold is irreducible, hence prime, and that $\pi_2(M) = 0$.)
Rather than working with ${\cal H}$, however, it is much easier to work with the quotient ${\mathcal S}$ of ${\cal H}$ by the order $3$-symmetry ${\kappa}$ corresponding to cyclic permutation of the parameters. We also introduce  a  commensurable orbifold ${\cal T}$ with a torus boundary  $ \partial {\cal T}$.

Both ${\mathcal S} = {\cal H} /{\kappa}$ and ${\cal T}$  surject to a $3$-orbifold ${\mathcal U}$ with fundamental group  a so-called \emph{$(P,Q,R)$-group}. Its boundary  $\partial {\mathcal U}$ is a sphere with three order $2$ and one order $3$ cone points. A similar construction has been used extensively by Akiyoshi et al, see for example~\cite{akiyoshi}, and is the basis of Wada's program OPTi, hence was convenient for our computations.
In this section we explain these constructions in detail,  using them to find explicit representations of all four groups.

\subsection{The handlebody and related orbifolds} \label{sec:mainhandlebody}

The symmetric handlebody ${\cal H}$ can be thought of as made by gluing two solid pairs of pants each with order $3$-symmetry. More precisely, take a $3$-ball  and remove three open disks from the boundary, placed so as to have order $3$ rotational symmetry. Gluing two such balls along the open disks produces a handlebody ${\cal H}$ with the required order three symmetry ${\kappa}$. Rather than write down a suitably symmetric representation of $\pi_1({\cal H})$ directly, we consider first the quotient orbifold ${\mathcal S} = {\cal H} /{\kappa}$.  As will be justified in retrospect when we have identified the representations explicitly, this is a ball with  two cone axes around each of which the angle is $2\pi/3$. Its boundary $\partial {\mathcal S}$ is a sphere $\Sigma_{0;3,3,3,3}$ with $4$ order $3$ cone points.  We will call ${\mathcal S}$ the \emph{large coned ball}.

The ball ${\mathcal S}$ has a further order $4$ symmetry group.
Consider  the two cone axes  which form the singular locus of ${\mathcal S}$, together with their common perpendicular $C$. This configuration is invariant under the $\pi$-rotation about $C$, and also under $\pi$-rotations about a unique pair of orthogonal
 lines  on the plane orthogonal to $C$ passing through its midpoint $O$, see Section~\ref{sec:largeball}.  Denoting these latter rotations $P,Q$, the $\pi$-rotation about $C$ is $PQ$ and the entire configuration is invariant under
 $\langle P,Q \rangle = \mathbb Z_2 \times \mathbb Z_2$. Thus we obtain a further quotient  orbifold ${\mathcal U} = {\mathcal S}/(\mathbb Z_2 \times \mathbb Z_2)$, also topologically a ball,  which we call  the \emph{small coned ball}.
The singular locus of ${\mathcal U} $ is as follows. Let $\bar O$ and $\bar  E$ be the images in ${\mathcal U}$ of the midpoint $O$ of $C$ and the point where $C$ meets $\Ax K$ respectively, where $K$ is one of the two order three rotations. Let $\bar C$ be the image of $C$, so that $\bar C$ is a line from $\bar O$ to $\bar E$.
 From  $\bar  O$ emanate three mutually orthogonal   lines  corresponding  to the order $2$ axes of $P,Q$ and $PQ$. One of  these is the line  $\bar C$ corresponding to $PQ$ which ends at $\bar  E$. From $\bar  E$ also emanates an order $3$ singular line, the projection of $\Ax K$, perpendicular to $\bar C$. The boundary $\partial {\mathcal U}$ is a sphere $\Sigma_{0;2,2,2,3}$ with $3$ cone points of order $2$  and one of order $3$. The order $3$ cone point is the end point of the order 3 singular line and the order $2$ cone points are the endpoints on $\partial U$ of the axes of $P, Q$ and a third involution $R$ defined below.

Finally, there is a double cover of the small coned ball ${\mathcal U}$ by an orbifold ${\cal T}$ which is topologically a solid  torus. Its
 boundary  is a torus $  \partial \mathcal T$ with a single cone point of angle $4 \pi/3$. Just as the quotient of a once punctured torus $  \Sigma_{1;\infty}$ by the hypelliptic involution is the surface $\Sigma_{0;2,2,2,\infty}$, so the quotient of    $ \partial \mathcal T $  by the hypelliptic involution $\iota$ is the surface $\partial {\mathcal U} = \Sigma_{0;2,2,2,3}$.   The involution $\iota$ extends to an involution, also denoted $\iota$, of $\mathcal T$ such that $\mathcal T/\iota = {\mathcal U}$.

The group $\pi_1({\mathcal U})$  is generated by  $(P,Q, K)$. We can replace $K$
by a further involution $R$ such that $RQP = K$.
 To do this,  let  $R$ be an order $2$ rotation about an axis contained in the plane through $E$ orthogonal to $ \Ax K$, such that the axis makes an angle $\pi/3$ with $C$. (We will fix orientations more precisely below.) Then $R (QP)$ is a $2\pi/3$-rotation about $ \Ax K$, in other words, provided orientations have been chosen correctly,  we can identify $\pi_1({\mathcal U}) $ with a group $  \langle P,Q,R | P^2=Q^2=R^2= (RPQ)^3  = -\rm{id}, PQ=-QP \rangle \subset SL(2,\mathbb C)$. (For a discussion on the choice of signs, see Remark~\ref{signchoices} below.)

In~\cite{akiyoshi} and other papers by the same authors,  groups generated by three involutions $P,Q,R$
with $R QP$ parabolic,
are used as a convenient way of parameterizing representations of once punctured tori, where the
 torus in question is now a two-fold cover  of  the orbifold
 with fundamental group $\langle P,Q,R \rangle$ with quotient induced by the hyperelliptic involution.
  A small modification of their parameterization allows us to write down a convenient  general form for a representation $\pi_1({\mathcal U}) \to SL(2,\mathbb C)$, from which we obtain explicit representations of $\pi_1({\cal H}), \pi_1({\mathcal S})$ and $\pi_1({\cal T})$. This we do in the next section.

 \subsection{The basic configuration and the small coned ball}\label{sec:basicconfig1}

 We start with a general construction for representations $\pi_1({\mathcal U}) \to SL(2,\mathbb C)$, that is, of subgroups $ \langle P,Q,R | P^2=Q^2=R^2= (RPQ)^3  = -\rm{id}, PQ=-QP \rangle \subset SL(2,\mathbb C)$.  For convenience we refer to such a group (or its image in $PSL(2,\mathbb C)$) as a \emph{$(P,Q,R)$-group}.

We will make our calculations using  \emph{line matrices} following \cite{Fenchel}. Note this  will define representations into $SL(2,\mathbb C)$, thus fixing the signs of traces.
Let $u,u' \in \hat {\mathbb C}$, and denote the oriented line from $u$ to $u'$ by $[u,u']$.
The associated line matrix $M([u,u'])\inSL(2,\mathbb C) $ is a matrix which induces an order two rotation about $[u,u']$
and such that $
M([u,u'])^2 = -\rm{id}$, so that in particular
 $$M([0,\infty])  = \begin{pmatrix} i & 0 \\ 0 & -i \end{pmatrix}.$$
 By \cite{Fenchel}, p. 64, equation (1), we have, if $u,u' \in \mathbb C$:
$$
M([u,u']) = \frac{i}{u-u'}
\begin{pmatrix} u+u' & -2uu' \\ 2 & -u-u' \end{pmatrix}.
$$

The representation we require is derived from
a basic configuration shown in Figure~\ref{fig:basicconfig}.
It depends on a single parameter $ \zeta \in \mathbb C$ which we will relate to the original parameter $x$ in \ref{sec:solidtorus} below.

 \begin{figure}[hbt]
\includegraphics[height=4cm]{Figs/PQR}
 \caption{The basic configuration for the $(P,Q,R)$-group $\pi_1({\mathcal U})$.}\label{fig:basicconfig}
\end{figure}

Let $\zeta\in\mathbb C$ and $P, Q, R\inSL(2,\mathbb C) $ be $\pi$-rotations
about the oriented lines $[\zeta, -\zeta]$, $[i\zeta, -i\zeta]$ and $[1,-3]$, respectively. By construction $P^2=Q^2= R^2 = -\rm {id}$. Moreover $ \Ax P$ and $ \Ax Q$ intersect at the point $|\zeta|j \in \mathbb H^3$ on the hemisphere of radius $|\zeta|$ centre $0 \in \mathbb C$, where $ z + tj$ represents the point at height $t>0$ above $z\in \mathbb C$ in the upper half space model of $\mathbb H^3$. Thus $PQ = -QP$ and $PQ$ is an order $2$ rotation about the vertical axis $0 + tj, t>0$.

Let $V$ be the vertical plane  above the real axis  in $\mathbb H^3$. Note that  the oriented axes of the order two rotations $PQ$ and $R$ both lie in  $V$,  intersecting in the point $\sqrt 3 j $   at angle $\pi/3$.
The line $[-\sqrt{3}i, \sqrt{3}i]$ passes through this point and is orthogonal to $V$.
It follows that $RPQ = -RQP$  is anti-clockwise rotation through $2\pi/3$ about the line
$[-\sqrt{3}i, \sqrt{3}i]$.  Using line matrices as above, we can now easily write down the
corresponding representation in $SL(2,\mathbb C)$:

\begin{align*}
P & = M([\zeta, -\zeta])
    = \frac{i}{2\zeta} \begin{pmatrix} 0 & 2\zeta^2 \\ 2 & 0 \end{pmatrix}
    = \begin{pmatrix} 0 & i \zeta \\ i/\zeta & 0\end{pmatrix}\\
Q & = M([i\zeta, -i\zeta])
    = \frac{1}{2\zeta} \begin{pmatrix} 0 & -2\zeta^2 \\ 2 & 0 \end{pmatrix}
    = \begin{pmatrix} 0 & -\zeta \\ 1/\zeta & 0 \end{pmatrix} \\
R & = M([1, -3])
    = \frac{i}{4} \begin{pmatrix} -2 & 6 \\ 2 & 2 \end{pmatrix}
    = \begin{pmatrix} -i/2 & 3i/2 \\ i/2 & i/2 \end{pmatrix}.
\end{align*}
Let $K = RQP$.  Then,
$$
K = \begin{pmatrix} -1/2 & -3/2 \\ 1/2 & -1/2 \end{pmatrix}, \quad
K^3 = \begin{pmatrix}1 & 0 \\ 0 & 1 \end{pmatrix},
$$
 so that as expected, $K$ is a anticlockwise rotation about $[-\sqrt{3}i, \sqrt{3}i]$ by $2\pi/3$.

Note that as matrices in $SL(2,\mathbb C)$, $P^2 = Q^2 = -\rm {id}$ and $PQ = -QP$. As isometries of $\mathbb H^3$, the signs are irrelevant. We could have chosen $K = RPQ$ in which case $K^3 = -\rm {id}$ but see Remark~\ref{signchoices} below.
We denote the group generated by $P,Q,R$ by $G_{{\mathcal U}}(\zeta)$ and the corresponding representation $\pi_1({\mathcal U}) \to SL(2,\mathbb C)$ by $\rho_{{\mathcal U}}(\xi)$.

 \subsubsection{The large coned ball ${\mathcal S}$}\label{sec:largeball}

To relate $\pi_1({\mathcal U})$ to $\pi_1({\mathcal S})$, start with two oriented axes $A_0, A_1$ about each of which we have order $3$ anticlockwise rotations $K_0,K_1$, measured  with respect to the orientation of the axes. Let $C$ denoted the common perpendicular between $A_0$ and $ A_1$, oriented from $A_0$ to $A_1$. We denote this configuration, which is clearly well defined up to isometry, by ${ \mathcal C} {  \mathcal F}$.
As described in
\ref{sec:mainhandlebody},
${ \mathcal C} {  \mathcal F}$ has a further $\mathbb Z_2 \times \mathbb Z_2$ group of symmetries generated by
 the $\pi$-rotations $P,Q$ with axes through the  mid-point of $C$: precisely, let $\Pi$ be the plane through  the mid-point of $C$ and orthogonal to $C$. Then
 the axes of $P,Q$ are the two lines in $\Pi$ which bisect the angles between the projections of   $\Ax K_0, \Ax  K_1 $ onto $\Pi$, chosen so that the angle bisected by  $ \Ax P$ is that between the projection of the lines  $\Ax K_0, \Ax  K_1 $  with the same (say outward) orientation.

This choice of $P$ ensures that $PK_0P^{-1} = K_1$ while $QK_0Q^{-1} = K_1^{-1}$.
Also $PQ$ is the order $2$ rotation about $C$, and $PQK_i Q^{-1}P^{-1} = K_i^{-1}, i = 0,1$.
As in Section~\ref{sec:mainhandlebody},
${\mathcal U} = {\mathcal S}/(\mathbb Z_2 \times \mathbb Z_2)$ and we can take $\pi_1({\mathcal U})$ to be the $(P,Q,R)$-group defined in~Section~\ref{sec:basicconfig1}.
In terms of  $(P,Q,R)$, the
generators of $\pi_1({\mathcal S})$ are $K_0 =  RQP, K_1 = PK_0P^{-1}$.
Thus
\begin{equation*}
K_0 =  -\begin{pmatrix} 1/2 & 3/2 \\ -1/2 & 1/2 \end{pmatrix}, \quad
K_1 =- \begin{pmatrix} 1/2 & -\zeta^2/2 \\ 3/2\zeta^2 & 1/2 \end{pmatrix}.
\end{equation*}

In terms of generators for $ \pi_1(\partial {\mathcal S})$, we have additionally
$ K_2 = QK_0Q^{-1}, K_3 = RK_0R^{-1}$ where
\begin{equation*}
K_2 =  -\begin{pmatrix} 1/2 & \zeta^2/2 \\ -3/2\zeta^2 & 1/2 \end{pmatrix}, \quad
K_3 = -\begin{pmatrix} 1/2 & -3/2 \\ 1/2 & 1/2 \end{pmatrix}
\end{equation*}
so that $K_0 K_3 K_1 K_2 = \rm{id}$.

 We denote  the group generated by $K_0,K_1$ by $G_{{\mathcal S}}(\zeta)$ and the corresponding representation $\pi_1({\mathcal S}) \to SL(2,\mathbb C)$ by $\rho_{{\mathcal S}}(\xi)$.
From now on, we frequently drop the subscript and refer to $K_0$ as $K$.

\subsubsection{\textbf{The handlebody ${\cal H}$}}\label{sec:handlebody}

Observe that the generator $X \in \pi_1({\cal H})$ projects to the loop $K_0K_1$ in ${\cal H} /{\kappa}$. (This latter is  a loop in $\partial {\cal H} /{\kappa}$ which separates one of each pair of the cone points of $K_0,K_1$ from the other pair.)
We arrange that the action of ${\kappa}$ is induced by conjugation by $K_0^{-1} = K^{-1}$, so the generators of $\pi_1({\cal H})$ can be written in terms of the generators of $\pi_1({\mathcal S})$ as $X = K_0 K_1$, $Y = K_0^{-1}XK_0 = K_1 K_0$. Thus we have:
$$
K^{-1} X K = Y, \quad
K^{-1} Y K = (XY)^{-1}, \quad
K^{-1} (XY)^{-1} K = X.
$$
Using the formulae from the previous section, this gives
$$
X = \begin{pmatrix}
    \frac{9}{4\zeta^2} + \frac{1}{4} &
    -\frac{\zeta^2}{4} + \frac{3}{4} \\
    \frac{3}{4\zeta^2} - \frac{1}{4} &
    \frac{\zeta^2}{4} + \frac{1}{4}
    \end{pmatrix},\quad
Y = \begin{pmatrix}
    \frac{\zeta^2}{4} + \frac{1}{4} &
    -\frac{\zeta^2}{4} + \frac{3}{4} \\
    \frac{3}{4\zeta^2} - \frac{1}{4} &
    \frac{9}{4\zeta^2} + \frac{1}{4}
    \end{pmatrix}.
$$
In particular this reveals the relation between the parameter $\zeta$ and $x$:
\begin{equation}\label{eqn:paramreln}
x= \tr X = \tr Y = \tr XY = \frac{\zeta^2}{4}+\frac{9}{4\zeta^2}+\frac{1}{2}.
\end{equation}
We denote the group generated by $X,Y$ by $G_{{\cal H}}(\zeta)$ and the corresponding representation $\pi_1({\cal H}) \to SL(2,\mathbb C)$ by $\rho_{{\cal H}}(\zeta)$; we explain in Section~\ref{sec:whichparameter}  why up to conjugation $\rho_{{\cal H}}(\zeta)$ in fact depends only on $x$.

 \begin{remark} \label{signchoices}   \rm{In the above discussion, we made choices of sign so that $K^3 = \mbox{\rm{id}}, X = K_0K_1$ (where $K=K_0$ as above).   To compute the discreteness locus of a family of representations only  requires looking in $PSL(2,\mathbb C)$, however for computations involving  traces  we need a lift to $SL(2,\mathbb C)$.

 By~\cite{Culler},  any $PSL(2,\mathbb C)$  representation of a Kleinian group can be lifted  to $SL(2,\mathbb C)$ provided there are no elements of order $2$; in particular  this applies to $PSL(2,\mathbb C)$ representations of $\pi_1({\mathcal S})$ and $ \pi_1({\cal H}) $.
 Since the product of the three generating loops corresponding to $X,Y, Z$ is the identity in $\pi_1({\cal H})$,  we should make a choice of lift in which $XYZ = \rm{id}$ in $SL(2,\mathbb C)$.
 We could choose the element $K$ which represents the  $3$-fold symmetry ${\kappa}$ to be  such that either  $K^3 =   \rm{id}$ or
 $K^3 = -\rm{id}$; however since we intend to work with   representations of $\pi_1({\mathcal S}) \to SL(2,\mathbb C)$ we should make the choice $K^3 =   \rm{id}$ because in the quotient orbifold ${\mathcal S}$, $K$ corresponds to a loop round  an order $3$ cone axis.

In the representation we have written down we achieve    $K^3 = \rm{id}$ with  the choice
  $K =  RQP = \begin{pmatrix} -1/2 & -3/2 \\ 1/2 & -1/2 \end{pmatrix}$.
It is easy to check that taking $K^3 = \rm{id}$, if we let $X = K_0K_1$ we get $XYZ = \rm {id}$ as required, but if we choose $X = -K_0K_1$ we get $XYZ = -\rm {id}$ which is wrong.
}
 \end{remark}

 \subsubsection{\textbf{The singular solid torus ${\cal T}$.}}\label{sec:solidtorus}

Finally we discuss the associated singular solid torus  ${\cal T}$, which is constructed in a standard way from the $(P,Q,R)$-group. We do not logically need to use ${\cal T}$ in our further development, however as explained in Section~\ref{sec:torustree},   in practice we used ${\cal T}$ for computations, moreover the interpretation of the problem in the more familiar setting of a torus with a cone point may be helpful.

The boundary $  \partial {\mathcal U}$  is a sphere with $4$ cone points  $x_P, x_Q,x_R$ and $ x_K$ corresponding to $P,Q, R$ and $K = RQP$.  Thus we can take as generators of $\pi_1({\cal T})$ the loop $B=PQ$ separating $x_P, x_Q$ from $x_R, x_K$, and the loop $A = RQ$ separating $x_R, x_Q$ from $x_P, x_K$. Since $P,Q$ have a common fixed point, $B$ is an order $2$ elliptic, while since the axes of $R,Q$ are (generically) disjoint, $A$ is a loxodromic whose axis extends the common perpendicular to $ \Ax R$ and $ \Ax Q$.

Using the formulae above for the $(P,Q,R)$-group we compute:
\begin{equation*}
RQ = A    = \begin{pmatrix} 3i/2\zeta & i\zeta/2 \\ i/2\zeta & -i\zeta/2 \end{pmatrix},  \quad
PQ= B  = \begin{pmatrix} i & 0 \\ 0 & -i\end{pmatrix},
\end{equation*} \quad
so that \begin{equation}\label{eqn:torustraces}
\tr A   = \frac{3i}{2\zeta} - \frac{i\zeta}{2}, \quad
\tr B =  0, \quad
\tr AB = -\frac{\zeta}{2} - \frac{3}{2\zeta}.
\end{equation}

Note that $AB = RP$ and $A^2 = -K_0K_1, B^2=-\rm{id}$.
We also deduce that
\begin{equation*}
ABA^{-1} B^{-1} = [A, B] = \begin{pmatrix} 1/2 & -3/2 \\ 1/2 & 1/2 \end{pmatrix},  \\ \
\mbox{{\rm so that }}
\Tr {[A,B] }= 1.
\end{equation*}

Note that $\Tr A \in [-2,2]$ if and only if $|\zeta| = \sqrt 3$ or $\zeta = it$ with
$ 1 \leq |t| \leq \sqrt 3$,
 justifying the above remark that generically $A$ is loxodromic.
Note also that $A^2 = -K_0K_1 $   is consistent with the direct computation using \eqref{eqn:torustraces}
 that $\Tr (A^2) = -( \frac{\zeta^2}{4}+\frac{9}{4\zeta^2}+\frac{1}{2})$.  Also note that $ [A, B] = -K^2$, so that the commutator is rotation by $ 4\pi/3$ about $ \Ax K$.  Since $\Tr K^2 = (\Tr K)^2 -2$ we find also that $\Tr {[A,B]} =  1$ independently  of the choice of sign for $K$. This is consistent with  $\Tr {[A,B]} = -2 \cos (2\pi/3)$,  the sign being negative by analogy  with the well known fact that for any representation of a once punctured torus group for which the commutator is parabolic, we have  $\Tr  {[A,B]}= -2 $.

 We denote  the group generated by $A,B$ by $G_{{\cal T}}(\zeta)$ and the corresponding representation $\pi_1({\cal T}) \to SL(2,\mathbb C)$ by $\rho_{{\cal T}}(\xi)$.

 \begin{remark} \label{signchoices1}   \rm{Once again there are questions of sign which this time are a little more subtle. If   $\alpha \in PSL(2,\mathbb C)$ corresponds to an element of order $2$ in $\pi_1(M)$, then the corresponding representation cannot be lifted to $SL(2,\mathbb C)$, because for non-trivial $\alpha \in SL(2,\mathbb C)$, necessarily  $\alpha^2 = -\rm{id}$, see~\cite{kra} and \cite{Culler}.
Since in $\pi_1({\cal T})$ the element $B^2$ is trivial,  a $PSL(2,\mathbb C)$ representation of  $\pi_1({\cal T})$ cannot be lifted to $SL(2,\mathbb C)$.
Nevertheless,  we can as above write down a group in  $SL(2,\mathbb C)$ which projects to a $PSL(2,\mathbb C)$ representation for $\pi_1({\cal T})$. See Section~\ref{sec:torustree} for further discussion on this point.
} \end{remark}

 \subsubsection{More on the configuration for the large coned ball ${\mathcal S}$}\label{sec:morelargeball}

 The relation~\eqref{eqn:paramreln} can be given a geometrical interpretation in terms of the perpendicular distance between the axes of $K_0,K_1$ which sheds light on the symmetries of  the configuration  ${ \mathcal C} {  \mathcal F}$ in Section~\ref{sec:largeball}.
To measure complex distance, we use the conventions spelled out in detail in~\cite{ser-wolp} Section 2.1. The signed complex distance $ {\bf d}_{ \alpha} (L_1,L_2)$ between two oriented lines $L_1,L_2$ along their oriented common perpendicular $\alpha$ is defined as follows.
The signed real distance $  d_{  \alpha} (L_1,L_2)$  is the positive real hyperbolic distance between $L_1,L_2$  if $\alpha$ is oriented from $L_1$ to $L_2$ and its negative otherwise. Let ${\bf v_i}, i=1,2$ be unit vectors to $L_i$ at the points $L_i \cap \alpha$ and let $\bf w_1$ be the parallel translate of $\bf v_1$ along $\alpha$ to the point $ \alpha \cap L_2$.
Then $  {\bf d}_{ \alpha} (L_1,L_2)= \delta_{\bf \alpha} (L_1,L_2)  + i\theta$ where $\theta$ is the angle, mod $2\pi i$, from $\bf w_1$ to $\bf v_2$ measured anticlockwise  in the plane spanned by  $\bf w_1$ to $\bf v_2$ and oriented by $\alpha$.

Let
$\sigma$ be the  signed complex distance from the \emph{oriented} axis $\Ax  K_0 $ to
 the oriented axis $ \Ax  K_1 $, measured along the common perpendicular $C$  oriented  from $\Ax  K_0 $ to
$ \Ax  K_1$. Then $\Ax  K_0 ,  \Ax  K_1$ together with  $\Ax  K_0K_1 $ form the alternate sides of a right angled skew hexagon whose other three sides are the common perpendiculars between the three axes taken in pairs. The cosine formula gives $\sigma$  in terms of the complex half translation lengths $\lambda_0, \lambda_1, \lambda_2$  of $K_0, K_1$ and $K_0K_1$ respectively. To get the sides oriented consistently round the hexagon we have to reverse the orientation of $ \Ax  K_0$  so that the complex distance $\sigma$  should be replaced by $ \sigma' = \sigma+ i \pi$ and $\lambda_0$  by $ \lambda'_0  = -\lambda_0 $, see~\cite{ser-wolp}, so the   formula  gives
$$ \cosh \sigma' = \frac{\cosh \lambda_2 - \cosh  \lambda'_0 \cosh  \lambda_1}{ \sinh  \lambda'_0 \sinh  \lambda_1} .$$

As in~\ref{sec:handlebody}, we have $X  = K_0K_1$ so
 $ x = \tr K_0K_1 = 2 \cosh \lambda_2$
while for $i = 0,1$ we have $ \cosh \lambda_i = \cos 2 \pi/3 = -1/2$ and $\sinh \lambda_i = i\sin 2 \pi/3 = i \sqrt 3/2$.  (Note that since $K_0, K_1$ are conjugate we should take $\lambda_0= \lambda_1$ so the possible additive ambiguity of $i\pi$  in the definition of the $\lambda_i$  does not change the resulting equation.) Substituting, we find
\begin{equation}\label{eqn:cxdist}
 -\cosh \sigma = \frac{x/2 - 1/4}{( \sqrt {3}/2)^2} = \frac{2x-1}{3}.\end{equation}

 We can also relate $\sigma$ directly to our parameter $\zeta$.  By construction $ \Ax  K_0$ is the oriented line
 $[-\sqrt 3 i, \sqrt 3i]$, while $K_1 = PK_0P^{-1}$ so that $ \Ax  K_1$ is the oriented line
 $[  i\zeta^2/\sqrt  3,  -i\zeta^2/\sqrt  3]$ and   $C$ is the oriented line from $\infty$ to $0$.
 Thus the real part of the hyperbolic distance from $ \Ax  K_0$ to $ \Ax  K_1$ is $2\log \sqrt 3 /|\zeta|$ and the anticlockwise
 angle, measured in the plane oriented \emph{downwards} along the vertical axis $C$, is $-(\pi + 2\arg \zeta) $.
 Hence $$\sigma = 2\log \sqrt 3 /|\zeta| - 2i\arg \zeta - i\pi= 2 \log \frac{\sqrt 3}{\zeta}.$$
  Comparing  to~\eqref{eqn:cxdist}, we find
   $$ \biggl [  \biggl(\frac{\sqrt 3}{\zeta}\biggr)^2 + \biggl(\frac{\zeta}  {\sqrt 3}\biggr)^2\biggr ] = 2 \cosh (\sigma + i \pi) =\frac{2(2x-1)}{3}$$
     or
   \begin{equation}\label{eqn:paramreln1}  x -1/2 = \frac{3}{4}   \biggl [  \biggl(\frac{\sqrt 3}{\zeta}\biggr)^2 + \biggl(\frac{\zeta}  {\sqrt 3}\biggr)^2\biggr ]   \end{equation}
 recovering and  giving a more satisfactory geometrical meaning to~\eqref{eqn:paramreln}.

 \subsubsection{Dependence on $x$ versus $\xi$.}\label{sec:whichparameter}
 It is not  perhaps  immediately obvious why the groups $G_{{\mathcal S}}(\xi), G_{{\cal H}}(\xi)$ as defined above depend up to conjugation only on our original parameter $x$. This is clarified by the above discussion, because up to conjugation $G_{{\mathcal S}}(\xi)$ depends only on the configuration ${ \mathcal C}{  \mathcal F}$ and hence on $\sigma$ which is related to $x$ as in~\eqref{eqn:cxdist}. An alternative way to see this is the discussion on computing traces in Section~\ref{sec:traces}. Thus from now on, we shall alternatively write
$G_{{\mathcal S}}(x), G_{{\cal H}}(x)$ in place of $G_{{\mathcal S}}(\xi), G_{{\cal H}}(\xi)$.

   \subsubsection{Symmetries}\label{sec:symmetries}
   The discussion in~\ref{sec:morelargeball} gives  insight into various symmetries of the parameters $x$ and $\zeta$.
Equation~\eqref{eqn:paramreln} shows that the map $ \zeta \mapsto x$ is a $4$-fold covering with branch points at $ \zeta = \pm\sqrt{3},  \zeta = \pm i\sqrt{3}$  and $z = 0,\infty$.
Correspondingly, we have a Klein  $4$-group $\mathbb Z_2 \times \mathbb Z_2$ of symmetries which change $\zeta$ but not $x$:
\begin{enumerate}
\item  replacing $\zeta$ by $-\zeta$ leaves the basic construction unchanged but the line matrices defining $P,Q$ change sign.

\item  replacing $\zeta$ by $-3/\zeta$  is an order $2$ rotation about the axis  $[-\sqrt 3i, \sqrt 3 i]$. This fixes $K_0$ and moves $K_1$ into a position on the opposite side of $K_0$ along the vertical line $C$. This changes nothing other than the position we choose for the basic configuration in Section~\ref{sec:basicconfig1}. Note however that
the line matrices defining $P,Q$ change sign.
\end{enumerate}

There is also a symmetry which changes $x$ as well as $\zeta$.
 Say we fix the orientation of one of the two axes $\Ax K_0,  \Ax  K_1 $ while reversing the other. On the level of the  configuration ${ \mathcal C} {  \mathcal F}$ from~\ref{sec:largeball},  this interchanges $P$ and $Q$. Since  $PK_0P^{-1} = K_1$ while $QK_0Q^{-1} = K_1^{-1}$,
 this is equivalent to  fixing the orientation of one of the two axes $ \Ax  K_0 ,  \Ax  K_1 $ while reversing the other.   This symmetry interchanges the marked group $P,Q, R$ with the marked group $Q,P, R$, so that one group is discrete if and only if so is the other.
  In terms of our parameters,  the complex distance $\sigma$ between the axes changes to $ \sigma + i\pi$, so that $\cosh \sigma \mapsto - \cosh \sigma$ giving the symmetry $(x -1/2) \mapsto -(x -1/2)$ of  \eqref{eqn:cxdist}. Note that the diagonal slice of the Bowditch set $\Delta \cap { \mathcal B}$ does not possess this symmetry.
 Interchanging $P$ and $Q$  is  induced by the map $ \zeta \mapsto i\zeta$; more precisely this map sends $P$ to $Q$ and $Q$ to $-P$.  This clearly  induces the same symmetry in equation \eqref{eqn:paramreln}.  Note that by the definition, in this symmetry $R$ remains unchanged.

 On  the level of the torus group $\pi_1({\cal T})$, we have by definition  $RQ= A$, $PQ=B$ so that $AB = RP$. Thus sending $P$ to $Q$ and $Q$ to $-P$ while fixing $R$ sends $B$ to $-B$ and  $A$ to $-AB$. (Recall that on the level of matrices, $PQ = -QP$.) The symmetry should therefore replace the trace triple $(\Tr A, \Tr B, \Tr AB)$  by the triple $(-\Tr AB, -\Tr B, \Tr A)$.
It is easily checked from~\eqref{eqn:torustraces} that this is exactly the change effected by $ \zeta \mapsto i\zeta$.

Finally, we have the  symmetry of complex conjugation induced by $ x \to \bar x$ or equivalently $ \zeta \mapsto \bar \zeta$. This sends $\sigma \mapsto \bar \sigma$ thus replacing  $G_{{\cal H}}(x)$ by a conjugate group in which the distance between $ \Ax  K_0 $ and $ \Ax  K_1 $  is unchanged but the angle measured along their common perpendicular changes sign. Clearly these are different groups but one is discrete if and only if the same is true of the other.

 The diagonal slice of the Bowditch set obviously also enjoys the symmetry by conjugation, however, that is its only symmetry. In particular   $(x,x,x) \to (-x,-x,-x)$  is not  a symmetry and  the corresponding $SL(2,\mathbb C)$ representations project to different representations of $F_2$ into $PSL(2,\mathbb C)$. This is because any two distinct
lifts of a  representation from $PSL(2,\mathbb C)$ to $SL(2,\mathbb C)$  differ by multiplying  exactly two of the parameters $x,y,z$ by $-1$.
The allowed replacement $X \to -X$ and $Y \to -Y$
gives the group $(-x,-x,x)$ with parameters which are not in the diagonal slice $\Delta$.

 The symmetries can be seen in our plots by comparing
 Figure~\ref{Diagonal-BQ}, the Bowditch set for the triple $\phi_{(x,x,x)}$  in the $x$-plane,  with the right hand frame of Figure~\ref{fig:BQ-sets-comparison}, which shows the same set in the
 $\zeta$-plane.  Note the symmetry of complex conjugation in both pictures. In addition, Figure~\ref{fig:BQ-sets-comparison} is invariant under the maps $ \zeta \mapsto -\zeta$ and $ \zeta \mapsto -3/\zeta$, neither of which are seen in  Figure~\ref{Diagonal-BQ}. Thus the upper half plane in
Figure~\ref{fig:BQ-sets-comparison} is a $4$-fold covering of the upper half plane in Figure~\ref{Diagonal-BQ}: as is easily checked from~\eqref{eqn:paramreln},  the imaginary axis in Figure~\ref{fig:BQ-sets-comparison} maps to the negative real axis in
Figure~\ref{Diagonal-BQ} while the real axis in Figure~\ref{fig:BQ-sets-comparison} maps to the positive real axis in
Figure~\ref{Diagonal-BQ}. In particular, note the following branch points and special values:
if $x=3$, then $\zeta=\pm 1, \pm 3$; if $x=2$, then $\zeta=\pm\sqrt{3}$; if $x=-1$, then $\zeta=\pm \sqrt{3}i$; if $x=-2$, then $\zeta=\pm i, \pm 3i$.

Finally, the symmetry $(x -1/2) \mapsto -(x -1/2) $ is not visible in either picture because it does not preserve the property of lying in the Bowditch set. As we shall see later, this symmetry is visible in pictures of the discreteness locus, see the right hand frame of Figure~\ref{Figs/Riley-Ray-BQ} below.

\section{Discreteness} \label{sec:discrete}

 We now turn to the question of finding those values of the parameter $x$ for which the group
$\langle X,Y \rangle$ is free and discrete.
Let $ \mathcal D_{{\mathcal S}},  \mathcal D_{{\cal H}} \subset \mathbb C$ denote the subsets of the complex $x$-plane on which the groups $G_{{\mathcal S}}(x), G_{{\cal H}}(x)$ respectively are discrete and the corresponding representations are faithful, so that in particular,  $G_{{\cal H}}(x)$ is free. (See Section~\ref{sec:whichparameter} for the replacement of $G_{{\mathcal S}}(\xi), G_{{\cal H}}(\xi)$ by
$G_{{\mathcal S}}(x), G_{{\cal H}}(x)$.)
 We first  show  that $  \mathcal D_{{\mathcal S}}= \mathcal D_{{\cal H}} $.

We begin with the easy observation that since all the groups in Section~\ref{sec:basicconfig} are commensurable, they are either all discrete or all non-discrete together:
\begin{lemma} \label{lem:commensurable} Suppose that $G, H$ are subgroups of $PSL(2,\mathbb C)$ with $G \supset H$ and that $[G:H]$ is finite. Then $G$  is discrete if and only if  the same is true of $H$.\end{lemma}
\begin{proof} If $G$ is discrete, clearly so is $H$. Suppose that  $H$ is discrete but $G$ is not. Then   infinitely many distinct  orbit points   in $G \cdot O$ accumulate in some compact set $D \subset \mathbb H^3$. Label the cosets of  $[G:H]$ as $g_1H \ldots  g_kH$. Then  for some $i$ there are infinitely many points
$g_i h_r \cdot O  \in D$, which gives infinitely many distinct points $h_r \in g_i^{-1}D$. This contradicts discreteness of $H$. \end{proof}

\begin{lemma} \label{lem:faithfultogether} The representation  $\rho_{{\mathcal S}}(x) {\colon \thinspace} \pi_1({\mathcal S}) \to G_{{\mathcal S}}(x)$ is faithful if and only if the same is true of  $ \rho_{{\cal H}}(x) {\colon \thinspace} \pi_1({\cal H}) \to G_{{\cal H}}(x)$.
\end{lemma}
\begin{proof}
Note that $\pi_1({\mathcal S})$  is isomorphic to  $\mathbb Z/3\mathbb Z * \mathbb Z/3\mathbb Z = <k_0, k_1 | k_0^3 = k_1^3 = \rm{id}>$, while
$\pi_1({\cal H})$ is the subgroup of $\pi_1({\mathcal S})$ generated by $k_0 k_1$ and $k_1 k_0$
and is isomorphic to a free group of rank $2$. By construction, $\rho_{{\cal H}}(x)$ is the restriction of $\rho_{{\mathcal S}}(x)$  to $\pi_1({\cal H})$.
Thus, if $\rho_{{\mathcal S}}(x)$ is faithful, then so is $\rho_{{\cal H}}(x)$.

Now $\pi_1({\cal H})$ has index three in $\pi_1({\mathcal S})$
 and $ \pi_1({\mathcal S}) = \pi_1({\cal H}) \cup k_0 \pi_1({\cal H}) \cup k_0^{-1} \pi_1({\cal H})$.
Suppose that $\rho_{{\cal H}}(x)$ is faithful but  $\rho_{{\mathcal S}}(x)$ is not.
Then, there exists $g \in \pi_1(S)$ such that $\rho_{{\mathcal S}}(x)(g) = \rm{id}$.
Now $g = k_0^e h$, where $e=\pm  1$ and $h \in \pi_1({\cal H})$.
Thus $ \mbox{\rm{id} }= \rho_{{\mathcal S}}(x)(g) = \rho_{{\mathcal S}}(x)(k_0^e) \rho_{{\cal H}}(x)(h)  $ so that $\rho_{{\cal H}}(x)(h^3) = \rho_{{\mathcal S}}(x)(k_0^{-3e}) = \rm{id}$
 contradicting  the assumption that $\rho_{{\cal H}}(x)$ is faithful.
\end{proof}

\begin{corollary} \label{discretetogether} The  representations $\rho_{{\mathcal S}}(x), \rho_{{\cal H}}(x)$  are faithful and discrete together, that is, $ \mathcal D_{{\mathcal S}}= \mathcal D_{{\cal H}} $.
\end{corollary}
Thus we may write $ \mathcal D =  \mathcal D_{{\mathcal S}}= \mathcal D_{{\cal H}} $. Our aim is this section is to find $ \mathcal D \subset \mathbb C$.

\subsubsection{Fundamental domains}\label{sec:dirichlet}
We can make a rough estimate  for $ \mathcal D $ by exhibiting  a fundamental domain for $G_{{\mathcal S}}(x)$ for sufficiently large $x$.

\begin{prop}  \label{prop:dirichlet}Writing  $x = u + iv$, the region $ \mathcal D$ contains the region  outside the ellipse
$(2u-1)^2/25   +   v^2/4  = 1$ in the $x$-plane.
 \end{prop}
\begin{proof} In view of Corollary~\ref{discretetogether}, we can work with
 the large cone manifold ${\mathcal S}$ with generators $K_0, K_1$ of Section~\ref{sec:largeball}.  By construction, $P$, which conjugates $K_0$ to $K_1$,  is $\pi$-rotation about $0+ j |\zeta|$ in the upper half space $\mathbb H^3$, which is   the mid-point of the common perpendicular between the axes of $K_0$ and $ K_1$.  The axis of $K_0$ is the line through $ j \sqrt 3$ perpendicular to the real axis.  Let $H, H'$ be the hemispheres which meet $\mathbb R$ orthogonally at points $ -3, 1$ and $-1,3$ respectively, and let $E, E'$ be the closed half spaces they cut out which contain infinity. Then $H, H'$ intersect in $ \Ax  K_0 $, moreover  the complement in $\mathbb H^3$ of $E \cup E'$ is a fundamental domain for the group $\langle K_0\rangle$ acting on $\mathbb H^3$. Likewise the images of $H,H'$ under $P$ meet along $ \Ax  K_1 $ and the complement of
$P(E ) \cup P(E')$  is a fundamental domain for  $\langle K_1\rangle$.
 Since $P(z) = \zeta^2/z, z \in \mathbb C$,  $P(H),P( H')$ meet $\mathbb R$ orthogonally in points $-\zeta^2/3, \zeta^2$ and $\zeta^2/3, -\zeta^2$ respectively.
Thus if $|\zeta| < 1$ then the circle of radius $|\zeta|$ separates the regions $ E \cap E'$ and $P(E ) \cap P(E')$. We conclude by Poincar\'e's theorem
 (or a suitable simple version of the Klein-Maskit combination theorem) that in this situation $\langle K_0, K_1\rangle$ is discrete with presentation $\langle K_0, K_1| K_0^3 =  K_1^3 = \rm {id}\rangle$. Thus the representation $\rho_{{\mathcal S}}(x)$ with $ x = x(\zeta)$ as in~\eqref{eqn:paramreln},  is faithful, and hence $x \in  \mathcal D$.

Suppose that $\zeta = e^{i\phi}$. Then from \eqref{eqn:cxdist}, $(2x-1)/3 =  \cosh \sigma = 1/2( e^{2i \phi}/3 + 3 e^{-2i \phi} )$ so that $x = u+iv$   lies on the ellipse
 $(2u-1)^2/25 + v^2/4 = 1$ as claimed. \end{proof}

\begin{figure}[ht]
 \includegraphics[width=6cm]{Figs/dirichletgeneral}
 \caption{The shaded region illustrates the fundamental domain for $\pi_1({\mathcal S})$  acting in its regular set in $\hat {\mathbb C}$ when $|\zeta| >1$, so that $x$ is outside the ellipse of Proposition~\ref{prop:dirichlet}.}\label{fig:dirichletgeneral}
\end{figure}

The configuration when $x \in \mathbb R$ is of particular interest since in this case   $G_{{\cal H}}(x)$ is Fuchsian. The ellipse meets the real axis in points $-2,3$ so that $G_{{\cal H}}(x)$ is  discrete and the representation is faithful on $(\infty, -2]$ (corresponding to $|\zeta| >1, \zeta\in    i\mathbb R$) and
$[3,\infty)$ (corresponding to $|\zeta| >1, \zeta \in    \mathbb R$).  In these two cases the fundamental domains look the same, see Figure~\ref{fig:fuchsianconfigs}.  Note that   the interval $(-2,3)$ is definitely not in $ \mathcal D$: if  $ -2 < x < 2$ then $K_0K_1   $ is elliptic since
 $x= \Tr K_0K_1   $, while if $ -1 < x < 3$ then  $K_0K_1^{-1}   $ is elliptic since
 $\Tr K_0K_1^{-1}  = 1-x   $, see also Section~\ref{sec:fuchsian}.

In the general case, a fundamental domain can be found by
a modification of  Wada's   program OPTi ~\cite{OPTi}.
This program allows one to compute the  limit set and fundamental
 domains for the $PQR$-group $G_{{\mathcal U}}$.  A short python code for doing this is available at \url{http://vivaldi.ics.nara-wu.ac.jp/~yamasita/DiagonalSlice/}

\subsection{The method of pleating rays} \label{sec:pleating}
To determine $ \mathcal D$, we  use  the Keen-Series method of pleating rays applied to the large coned sphere ${\mathcal S}$. This is closely analogous to the problem of computing the Riley slice of Schottky space, that is the parameter space of free discrete groups generated by two parabolics, which was solved in~\cite{ksriley,koms}.

We begin by briefly summarising the elements of pleating ray theory we need.
For more details see various of the first author's papers, for example~\cite{ksriley, chois}.

Suppose that  $G \subset SL(2,\mathbb C)$ is a geometrically finite Kleinian group  with corresponding orbifold  $M = \mathbb H^3/G$  and let ${ \mathcal C}/G$ be its convex core, where ${ \mathcal C}$ is the convex hull in $\mathbb H^3$ of the limit set of $G$, see~\cite{EpM}.
 Then  $\partial { \mathcal C}/G$ is a convex pleated surface (see for example~\cite{EpM}) also homeomorphic to $\partial M$. The
bending of this pleated surface is recorded by means of  a measured geodesic lamination,  the \emph{bending lamination} $\beta = \beta(G)$,
whose support forms the bending lines of the surface   and whose transverse measure records the  total bending angle along short transversals.
 We say  $\beta $ is  \emph{rational} if it is supported on closed curves: note that closed curves in the support of $\beta $ are necessarily simple and pairwise disjoint. If a  bending line is represented by a curve $\gamma \in \pi_1({\mathcal S})$, then by definition it is the projection of a geodesic axis to $\partial { \mathcal C}/G$, so in particular $\beta$ contains no peripheral curves in its support.
Note that any two homotopically distinct non-peripheral simple closed curves on
$\partial  {\mathcal S}$ intersect. Thus in this case,  $\beta $ is  rational only if its support is a single simple essential non-peripheral  closed curve on $\partial { \mathcal C}/G$.

As above, we parameterise representations $ \rho_{{\mathcal S}}(x) {\colon \thinspace} \pi_1({\mathcal S}) \to SL(2,\mathbb C)$ by $x \in \mathbb C$ and denote the image group by $G_{{\mathcal S}}(x)$. From now on, we frequently write $ \rho_{x}$ for $ \rho_{{\mathcal S}}(x)$.

 \begin{definition}
Let $\gamma$ be a homotopy class of simple essential non-peripheral   closed curves on $\partial {\mathcal S}$.
The \emph{pleating ray}  $ {\cal P}_{\gamma}$ of   $\gamma$ is the
set of points $ x \in  \mathcal D$ for which   $ \beta(G_{{\mathcal S}}(x)) = \gamma$.
\end{definition}

Such rays are called \emph{rational pleating rays}; a similar definition can be made for general projective classes of bending lamination, see~\cite{chois}.

The following key lemma is proved in~\cite{chois} Proposition 4.1, see also~\cite{kstop} Lemma 4.6.  The essence is that because the two flat pieces of $\partial { \mathcal C}/G$ on either side of a bending line are invariant under translation along the line, the translation can have no  rotational part.
\begin{lemma}
  \label{lemma:realtrace} If the axis of $g \in G$ is a bending line   of $\partial { \mathcal C}/G_{{\mathcal S}}(x)$, then
$\Tr g(x) \in \mathbb R$.
\end{lemma}
 Notice that the lemma applies even when the bending angle $\theta_{\gamma}$ along $\gamma$ vanishes, so the corresponding surface is flat, or when the angle is $\pi$, in which case either $\gamma$ is  parabolic or $G_{{\mathcal S}}(x)$ is Fuchsian.

If $g \in G$ represents a curve $\gamma$ on $\partial {\mathcal S}$,  define the \emph{real trace locus} $\mathbb R_{\gamma}$ of $\gamma$ to be the locus of points in $\mathbb C$ for which
$\Tr g \in (-\infty, -2] \cup [2,\infty)$. By the above lemma, ${\cal P}_{\gamma} \subset \mathbb R_{\gamma}$.

Our aim is to compute  the locus of faithful discrete representations $ \mathcal D_{{\mathcal S}}=  \mathcal D$. In summary, we do this as follows:
\begin{enumerate}
\item Show that up to homotopy in ${\mathcal S}$, the essential non-peripheral curves on $\partial {\mathcal S}$ are indexed by $ \mathbb Q / \sim $  where $p/q \sim \pm (p+ 2kq)/q, k \in \mathbb Z$. (Proposition~\ref{curveequality}).

\item Given $\gamma \in \pi_1(\partial{\mathcal S})$, give an algorithm for computing $\Tr \rho_x(\gamma)$  as a polynomial in $x$, in particular identifying its  two highest order terms in terms of $p,q$. (Section~\ref{sec:traces} and Proposition~\ref{prop:tracepoly}).

\item Show that ${\cal P}_{0/1} = (-\infty, -3]$ and ${\cal P}_{1/1} = [2, \infty)$ (where
 ${\cal P}_{p/q}$ denotes the pleating ray of the curve $\gamma_{p/q} \in \pi_1(\partial{\mathcal S})$ identified with $p/q$).
(Section~\ref{sec:fuchsian}).

\item Show that  ${\cal P}_{p/q}$ is a union of connected non-singular branches of $\mathbb R_{\gamma}$.
(Theorem~\ref{thm:nonsing}).

\item For   $p,q \neq 0,1$, identify ${\cal P}_{p/q}$ by
showing it has two connected components, namely the branches of $\mathbb R_{\gamma}$ which are asymptotic to the directions $e^{\pm i \pi (p/q + 1)}$ as $|x| \to \infty$.
(Proposition~\ref{prop:connectivity}).

\item Prove that rational rays ${\cal P}_{p/q}$ are dense in $ \mathcal D_{{\mathcal S}}$.  (Theorem~\ref{thm:density}).
\end{enumerate}

One could carry all this out following almost word for word the arguments in~\cite{ksriley}. Rather than do this, we indicate as appropriate how more general results can be put together to provide a  somewhat less ad hoc proof of the results. The claim that ${\cal P}_{p/q}$ has two connected components appears to contradict the results in~\cite{ksriley}, see however  the following remark and Proposition~\ref{prop:connectivity} below.
 The pleating rays are shown on the left in Figure~\ref{Figs/Riley-Ray-BQ}  with the Riley slice rays from~\cite{ksriley} on the right for comparison.

\begin{remark} \label{Komori issue}
\rm{There were two rather subtle errors in~\cite{ksriley} .  The first was  that, in the enumeration of curves on $\partial {\mathcal S}$, we omitted to note that $\gamma_{p/q} $ is homotopic to $ \gamma_{-p/q}$ in ${\mathcal S}$. The second was, that we found only one of the two components of ${\cal P}_{p/q}$.
 Since ${\cal P}_{p/q} =   {\cal P}_{-p/q}$, these two errors in some sense cancelled each other out. They were discussed at length and resolved in ~\cite{koms} and we make corresponding corrections here.}
\end{remark}

 \begin{figure}[ht]
\includegraphics[width=5cm]{Figs/Riley-mu-3-Ray-BQ}\hspace{1cm}
\includegraphics[width=5cm]{Figs/Riley-Ray-BQ}
\caption{Left:  Pleating rays for $G_{{\mathcal S}}(x)$. Right: Pleating rays  for the Riley slice as described in~\cite{ksriley}.
The coloured (grey) regions are the Bowditch sets for the initial triples discussed in~\ref{sec:discussion}; conjecturally these coincide with the closure of the regions filled by the pleating rays. For a discussion of how the rays were actually computed, see Section~\ref{sec:torustree1}
} \label{Figs/Riley-Ray-BQ}
\end{figure}

 \subsection{Step 1: Enumeration of curves on $\partial {\mathcal S}$} \label{sec:enumeration}

We need to enumerate essential non-peripheral unoriented simple curves on $\partial {\mathcal S}$ up to homotopy equivalence in ${\mathcal S}$.
As is well known,  such curves on $\partial {\mathcal S}$ are,  up to homotopy equivalence in $\partial {\mathcal S}$,  in bijective correspondence with lines of rational slope
in the plane, that is, with $\mathbb Q \cup \infty$, see for example~\cite{ksriley, koms}. For $(p,q) $ relatively prime and $q \geq0$,  denote the class corresponding to $p/q$ by $\gamma_{p/q}$.
We have:
\begin{proposition}[\cite{koms} Theorem 1.2] \label{curveequality}The unoriented curves $\gamma_{p/q}, \gamma_{p'/q'}$ are homotopic in ${\mathcal S}$ if and only if
$p'/q' = \pm p/q  + 2k, k \in \mathbb Z$.
\end{proposition}
Missing the identification  $\gamma_{p/q} \sim \gamma_{-p/q}$ was the first of the two errors in~\cite{ksriley} referred to in Remark~\ref{Komori issue}.

Before proving the proposition, we need to explain the identification of curves on $\partial S$ with $\mathbb Q \cup \infty$. In~\cite{ksriley, koms} this was done using the plane punctured at  integer points as an intermediate covering between $\partial S$ and its universal cover.  The idea is sketched in Section~\ref{sec:correspondence}. Here we give a slightly different description of the curve $\gamma_{p/q} $ which leads to a nice proof of the above result.

 Cut ${\mathcal S}$ into two halves along the meridian disk $m$ which is the plane which perpendicularly bisects the common perpendicular $C$ to the two singular axes $\Ax  K_i , i=0,1$.  Each half  is a ball  $\hat B_i$ with a singular axis $\Ax  K_i $. The boundary  $\partial B_i = \partial \hat B_i \cap \partial {\mathcal S}$  is a sphere with two cone points and a hole $\partial m$.  Since the axes of $K_i$ are oriented, we can distinguish one cone point on each $\partial B_i$ as the positive end  of $\Ax  K_i $.
Now $\partial {\mathcal S}$ has a hyperbolic structure inherited from the ordinary set (or from the pleated surface structure on $\partial { \mathcal C}/G_{{\mathcal S}}(x)$), in which $\partial m$ is geodesic. With respect to such a structure,
each $\partial B_i$ has a reflectional symmetry $\iota$ in the plane containing $\Ax  K_i $ and $C$, which maps the cone points to themselves, is an involution on  $\partial m$ with two fixed points which maps the `front' to the `back' as shown in Figure~\ref{fig:fronttoback}.    There is a preferred base point $P_i$ on $\partial m$, namely the foot of the perpendicular from the negative end of  $\Ax  K_i $ to $\partial m$.

 \begin{figure}[ht]
\includegraphics[width=8cm]{Figs/fronttoback}
\caption{The arrangement of arcs on $\partial {\mathcal S}$. The curve shown illustrates the case $p=1, q=3$.}\label{fig:fronttoback}
\end{figure}

Let $\gamma$ be an essential non-peripheral simple curve on $\partial {\mathcal S}$. For each $i= 0,1$,  $\gamma \cap \partial B_i$ consists of $q$ arcs joining $\partial m$ to itself.  Start with the strands of $\gamma$ arranged symmetrically with respect to $\iota$, that is, with front to back symmetry.
 Orient $\partial m$ so that it points `upwards' on the front side of the figure, noticing that $\iota$ reverses the orientation. Lifting $\partial m $ to its cyclic cover $\mathbb R$, enumerate in order the endpoints $X^{k}_i, i=0,1;  k \in \mathbb Z$ of arcs of $\gamma$   starting (say) with the arc meeting
$\partial m$ nearest $P_i$, and so that increasing order is in the direction of the upwards orientation of $\partial m$ viewed from the front side in the figure. Since in fact $X^{k}_i = X^{k+2q}_i $, the enumeration is really $\rm{mod} \  2q\mathbb Z$.

   To reconstruct  $\gamma$ we have to join the endpoints $X^{k}_0$ on $\partial B_0$ to the endpoints $X^{k'}_1$ on $\partial B_1$.
  Since the arcs have to be matched in order round $\partial m$, if  $X^i_0$ is  joined to $X^j_1$ then
$X^{i+k}_0$ is  joined to $X^{j+k}_1$ for all $k \in \mathbb Z$. Set  $ p = j-i$.

It is not hard to see that the resulting curve $\gamma_{p/q}$ is connected if and only if $(p,q)$ are relatively prime. Note that with this description, $\partial m$ is the curve $q=0$, that is $\gamma_{1/0}$. The curve $\gamma_{0/1}$ is the curve $K_0K_1$ and $\gamma_{1/1} = K_0K_1^{-1}$.
We leave it to the reader to see that this description is the same as that obtained from the lattice picture in~\cite{ksriley}, see also Section~\ref{sec:correspondence}.

\medskip
{\sc Proof of Proposition~\ref{curveequality}}
Write $\gamma_{p/q} \sim \gamma_{p'/q'}$ to indicate that $\gamma_{p/q} , \gamma_{p'/q'}$ are homotopic in ${\mathcal S}$.
Since Dehn twisting round $\partial m$ is trivial in ${\mathcal S}$ and sends $X^k_i \to X^{k+2q}_{i}$,  we have $\gamma_{p/q} \sim \gamma_{p/q +2}$. To see why $\gamma_{p/q} \sim \gamma_{-p/q}$, we proceed as follows. Consider the arrangement  shown in Figure~\ref{fig:fronttoback}, in which $B_0$ is on the left  and  $\partial m$ is oriented `upwards'. Each $\partial  B_i$ has  a natural `front' and `back' which are interchanged by the involution $\iota$.
We will prove the result by arranging $\gamma$ on $\partial S$ in two different ways. Start with the strands of $\gamma$ arranged symmetrically front to back, that is, with respect to $\iota$, as on the two sides of Figure~\ref{fig:fronttoback}.

(1) Homotope  all the arcs of $\gamma \cap \partial B_i$ on $\partial B_i$  by dragging the endpoints on the back side around one or other cusp so that they meet $\partial m$ on the front, see Figure~\ref{fig:pullround}. In this case, the choice $p>0$ means that
on an arc approaching $\partial m$ from $\partial B_0$ one turns left (up) for $p$ slots before joining up arcs.

 \begin{figure}[ht]
  \includegraphics[width=4cm]{Figs/pullround}
 \caption{Curves pulled round to the front of $B_0$ ready to be joined as in (1) of the proof of Proposition~\ref{curveequality}.}\label{fig:pullround}
\end{figure}

(2) Now for the second arrangement. Homotope  all the arcs of $\gamma \cap \partial B_i$ on $\partial B_i$, by dragging the endpoints on the front side  so that they meet $\partial m$ on the back of $\partial B_i$. Notice turning left for $p$ slots on the back side, when viewed through ${\mathcal S}$ from the front to the back, looks the same as turning right (down) for $p$ slots on the front side.

Now, starting from situation (2),  homotope $\gamma$ in ${\mathcal S}$ by pulling the arcs meeting $\partial m$ on the back side through $m$ so that they meet $\partial m$ on the front side keeping their `horizontal' level fixed,  so that an endpoint $X$ moves to $\iota(X)$. The endpoints on the front side now have  the same up-down order as they did on the back side. This shows that  the curve obtained connecting the arcs moving up $p$ slots is homotopic in ${\mathcal S}$ to the curve obtained by moving down $p$ slots, as claimed.

We will show that  if $p'/q' \neq \pm p/q  + 2k, k \in \mathbb Z$ then $\gamma_{p/q} \nsim \gamma_{p'/q'}$ after computing traces, see Corollary~\ref{curveinequality}. \qed

\subsection{Step 2: Computation of traces}\label{sec:traces}

Let $V_{p/q} \in SL(2,\mathbb C) = \rho_x(\gamma_{p/q})$, where, since we want to compute $\Tr V_{p/q}$,  we only need consider $V_{p/q}$ up to cyclic permutation and inversion.
Rather than using the associated torus tree, we will work directly with a $4$-holed sphere $\Sigma_{0,4}$ and the associated tree as described in~\cite{MPT}, see also \cite{goldman}.
Let $\alpha, \beta, \gamma,\delta$ denote loops round the four holes, oriented so that $\alpha  \beta \gamma \delta = \rm{id}.$ The fundamental group is identified with the free group $  F_3$ with generators
$\alpha, \beta, \gamma$.
 A representation $\rho{\colon \thinspace} F_3 \to SL(2,\mathbb C) $ is determined up to conjugation by  its values on seven elements as follows (where we use $\hat w$ in place of $w$ in \cite{MPT} etc to distinguish it from a variable $w$ already in other use):
$$ \Tr \rho(\alpha) = a; \Tr \rho(\beta) = b; \Tr \rho(\gamma) = c; \Tr \rho(\delta) = d$$
$$ \Tr \rho(\alpha\beta) = \hat x; \Tr \rho(\beta\gamma) = \hat y; \Tr \rho(\gamma\alpha) = \hat z$$
related by the equation
\begin{equation} \label{conedspheretraces}  \hat x^2 +\hat y^2 +\hat z^2 +\hat x\hat y\hat z = \hat p\hat x + \hat q\hat y+\hat  r\hat z + \hat  s\end{equation}
where
$$ \hat p = ab + cd,  \hat  q = bc+ad,  \hat  r = ac+bd,  \hat  s = 4-a^2 -b^2 - c^2 - d^2 - abcd.$$
We identify our generators $K_i$ as:
$ \alpha = K_0, \beta = K_1, \gamma = K_2, \delta = K_3$.
Thus we find:
$$ a = b = c = d = 1, \hat x = \Tr K_0K_1 = x, \hat y =  \Tr K_1K_2 = 2, \hat z =  \Tr K_2K_0 = -x+1.$$
As a check, it is easy to verify that the trace identity \eqref{conedspheretraces} holds. Notice that none of the expressions  $\hat p, \ldots, \hat z$ depend on the sign choices made in Section~\ref{sec:handlebody}.

The traces can be arranged in a trivalent tree in the usual way. As explained above, we have $\gamma_{0/1} = K_0K_1, \gamma_{1/0} =   \rm{id}, \gamma_{1/1} =K_0K_1^{-1} $.  As explained in~\cite{MPT} Section 2.10, there are now $3$ moves, depending on the values of $\hat p,\hat q, \hat r$. In our case $\hat p=\hat q= \hat r = 2$ so the three moves described there coincide.
Following~\cite{MPT}, if $u,v,w$ are labels round a vertex, with $v,w$ labels adjacent along a common edge $e$, then the label at the vertex at the opposite end of $e$ is $u' = 2-vw-u$, compare Figure~\ref{fig:markofftree} in which $u ' = vw-u$.

Clearly this procedure gives an algorithm for arranging curves and computing  traces on a trivalent tree by analogy with that described in Section~\ref{sec:markoff}.
Curves generated  in this way inherit  a natural labelling from  the usual procedure of Farey addition as described in Section~\ref{sec:markoff}. Denote the curve  which inherits the label $p/q$   by $\delta_{p/q}$; we say this curve is in \emph{Farey position $p/q$} on the tree. We shall refer to this tree together with its new rule for computing traces  as the ${\mathcal S}$-tree, to distinguish it from the Markoff tree of Section~\ref{sec:markoff}.

 We need to show that $\delta_{p/q}$ is the same as the curve  $\gamma_{p/q}$ described in the previous section, namely  the
curve with $2q$ intersections with the meridian $\partial m$ and a twist by $p$.
 \begin{lemma} \label{lem:treeisok}  With the above notation,  $\delta_{p/q} = \gamma_{p/q}$.\end{lemma}
\begin{proof} (Sketch) By definition we have $\delta_{p/q} = \gamma_{p/q}$ for $p,q \in \{0,  1 \}$. With the notation above,  these are the curves $\alpha \beta, \beta \gamma, \gamma \alpha$, each of which separates the punctures in pairs.

Call two essential simple non-peripheral curves on $\partial {\mathcal S}$ \emph{neighbours} if they intersect exactly twice.
Note that of the initial triple,  each pair adjacent along an initial edge are neighbours, so that the triple round the initial vertex are neighbours in pairs.

Now we check inductively that this is always the case.  Give a pair of neighbours $\delta, \delta'$ along an edge, the remaining curves at the vertices at the opposite ends of this edge are obtained by Dehn twisting $\delta$ about $\delta'$ (or vice versa) in each of the two possible directions, see~\cite{MPT}. Thus for example $\delta_{1/1}$ is obtained by  Dehn twisting $\delta_{0/1}$ about $\delta_{1/0}$ while $\delta_{-1/1}$ is obtained by  Dehn twisting $\delta_{0/1}$ about $\delta_{1/0}$ in the opposite direction. Moreover if $\delta, \delta'$ are neighbours, then so are the pairs
$\delta, D^{\pm}_{\delta}( \delta')$ and $\delta', D^{\pm}_{\delta}( \delta')$, where $D_{\delta}( \delta')$ denotes a Dehn twist of $\delta'$ about $\delta$ and by abuse of notation we write $D^{+}, D^{-}$ for $D, D^{-1}$ respectively.

Now we show inductively that $\delta_{p/q} = \gamma_{p/q}$. Suppose that it is true for neighbours $p/q,r/s$ where $|ps-rq| = 1$. By induction we may assume that $\delta_{p/q} , \delta_{r/s}$ are adjacent along an edge $e$ of the tree. The two remaining curves $D^{\pm}_{\delta_{p/q}}(\delta_{r/s})$ at the vertices at the ends of $e$   result from Dehn twisting $\delta_{p/q}$ about $\delta_{r/s}$ in each of the two possible directions, hence are also neighbours. From the inductive hypothesis, we can assume that one of these two curves is $\delta_{p-r/q-s} = \gamma_{p-r/q-s}$ or $\delta_{r-p/s-q} = \gamma_{r-p/s-q}$  depending on whether  $ ps-rq  = \pm 1$. In accordance with the Farey labelling system, the curve at the other vertex is $\delta_{p+r/q+s} $.

 \begin{figure}[ht]
  \includegraphics[width=8cm]{Figs/surgery}
 \caption{Here  $\gamma_{1/3}$ and $\gamma_{1/2}$ are surgered to give $\gamma_{2/5}$, see the proof of Lemma~\ref{lem:treeisok}. The inset circles show the direction of surgery. }\label{fig:surgery}
\end{figure}

Now  the curves
 $D^{\pm}_{\delta_{p/q}}(\delta_{r/s})$ can be found by surgery. On each $B_i$, arrange both curves symmetrically with respect to the front and back of ${\mathcal S}$ as described above, then join the strands so that they have minimal intersection as in Figure~\ref{fig:surgery}. With $\gamma_{p/q}$ in this position, its twist $p$  is its intersection number with the geodesic joining the two positive cone points of the axes $K_i$,  and likewise for
$\gamma_{r/s}$.  To perform the Dehn twist we have to cut the curves at their intersection points and then make a consistent choice of which direction to rejoin the resulting arcs.  One of the two choices will give a curve with $2(q+s)$ intersection  points with the meridian $\partial m$. Clearly the curve with the `positive' surgery (see the inset circles in Figure~\ref{fig:surgery}) will have intersection number $p+r$ with this line, and hence is the curve $\gamma_{p+r/q+s}$. Since we already know the curve at the other vertex is  $\delta_{r-p/s-q} =\gamma_{p-r/q-s}$ (or $\gamma_{r-p/s-q}$), this shows that $\delta_{p+r/q+s} = \gamma_{p+r/q+s}$. This completes the proof.
\end{proof}

\begin{prop} \label{prop:tracepoly} Let $V_{p/q}(x) = \rho_{{\mathcal S}}(x) (\gamma_{p/q})$ as above. Then:
\begin{enumerate}
\item $\Tr V_{p/q}$ is a  polynomial  in $x$ whose top two terms are $(-1)^{p-q-1} (x^{q} - px^{q-1} )$.

\item $\Tr V_{p/q} = \Tr V_{(p/q) + 2} = \Tr V_{-p/q}$.

\item $\Tr V_{p/q} (x) = \Tr V_{p+q/q } (1-x)$.

\end{enumerate}

\end{prop}

\begin{remark}\rm{(1) should be compared to~\cite{ksriley} Corollary 4.3 in which we showed   that the leading term is of the form $(-1)^{p-q-1}c x^q$ for some $c>0$, see also the remark following the corollary in that paper. }\end{remark}
\begin{proof}  (1)
Note that (1) holds for the three initial traces of $  \gamma_{0/1},  \gamma_{1/0}
,  \gamma_{1/1} $. If curves $  \gamma_{p/q},  \gamma_{r/s}$ are adjacent along an edge, then the two curves at the remaining vertices at the ends of the edge are $  \gamma_{p\pm r/q \pm s}$.
The result then follows easily by induction on the tree.

(2) This follows immediately from Proposition~\ref{curveequality} and can also be proved easily  by looking at the symmetries of the ${\mathcal S}$-tree.

(3) This results from the symmetry $ x \mapsto 1-x$ which interchanges $  \gamma_{0/1},  \gamma_{1/1}$.

\end{proof}

Now we can prove the `only if' assertion of Proposition~\ref{curveequality}:
\begin{corollary} \label{curveinequality}  If $p'/q' \neq \pm p/q  + 2k, k \in \mathbb Z$ then $\gamma_{p/q} \not\sim \gamma_{p'/q'}$.

\end{corollary}
\begin{proof} This follows immediately by comparing the top two terms of $\Tr W_{p/q}, \Tr W_{p'/q'}$.
\end{proof}

\subsection{Step 3. The exceptional Fuchsian case: computation of ${\cal P}_{0/1}, {\cal P}_{1/1}$} \label{sec:fuchsian}

 As above, let ${\cal P}_{p/q}$ denote the pleating ray of $\gamma_{p/q}$.
The rays ${\cal P}_{0/1}, {\cal P}_{1/1}$ are exceptional. We have $\Tr \gamma_{0/1}= x, \Tr \gamma_{1/1} = 1-x$. Thus
the real locus for both trace polynomials is exactly the real axis, and on this locus, the group $G_{{\mathcal S}}(x)$, if discrete, is Fuchsian.
This is exactly the situation discussed in~\cite{ksriley} p. 84.

In the ball model of $\mathbb H^3$,  identify the extended real axis with the equatorial circle. Since the limit set is contained in $\hat \mathbb R$, the convex core (the Nielsen region) of $G_{{\mathcal S}}(x)$
 is contained in the equatorial plane. We can think that the convex core has been squashed flat and the bending lines are just the boundary of the Nielsen region, that is, the boundary of the surface $\mathbb H^2/G_{{\mathcal S}}(x)$. Thus to find the bending lamination we just have to determine  the boundary of $\mathbb H^2/G_{{\mathcal S}}(x)$.

 Now if  $x \in \mathbb R$ then either $\zeta \in \mathbb R$ and  $x>0$, or $\zeta \in i\mathbb R$ and $x<0$.
 In both cases, we find a fundamental domain for $G_{{\mathcal S}}(x)$ as  described in~\ref{sec:dirichlet}, see Figure~\ref{fig:fuchsianconfigs}. Thus regarded as a Fuchsian group acting on the upper half plane $\mathbb H$, $G_{{\mathcal S}}(x)$ represents a sphere with two order $3$ cone points and one hole.  However the cases $ x<0, x>0$ are slightly different, because of the relative directions of rotation of $K_0$ and $K_1$.

 In both cases, the axis $K_0$ has fixed points $\pm i \sqrt 3$ and its axis is oriented so that it is anticlockwise rotation about $  i \sqrt 3$. Thus $K_1 = PK_0P^{-1}$ rotates anticlockwise  about
 $P(  i \sqrt 3) = -i \zeta ^2/\sqrt 3 $. If $x<0$ then  $P(  i \sqrt 3)$ is in the upper half plane $\mathbb H$ while  if $x>0$ then  $P(  i \sqrt 3)$ is in the lower  half plane. Hence if $x<0$ then $K_0, K_1$ rotate in the same sense about their fixed points in $\mathbb H$ while if $x>0$ their rotation  directions are opposite. This leads to the two different configurations shown in Figure~\ref{fig:fuchsianconfigs}.

 As is easily checked, if $x>0$ the boundary of the hole is thus $K_0K_1^{-1}$ while if $x<0$ the boundary of the hole is $K_0K_1$.  Since $K_0K_1 = \gamma_{0/1}$ and $K_0K_1^{-1} = \gamma_{1/1}$, combining this with information about the discreteness locus in the Fuchsian case from~\ref{sec:dirichlet},  we conclude that ${\cal P}_{0/1} = (-\infty, -2]$ and $ {\cal P}_{1/1}= [3, \infty)$.

 \begin{figure}[ht]
 \includegraphics[width=6cm]{Figs/dirichletxleq0}
  \includegraphics[width=6cm]{Figs/dirichletxgeq0}
 \caption{Configurations  for $x \in \mathbb R$. Left: $\zeta \in i\mathbb R, x\leq -2$. $K_0$ and $ K_1$ rotate in the same  directions   $\mathbb H$ and the hole is $K_0K_1$.
Right: $\zeta \in  \mathbb R, x\geq 3$.  $K_0$ and $ K_1$ rotate in opposite directions in $\mathbb H$ and the hole is $K_0K_1^{-1}$.}\label{fig:fuchsianconfigs}
\end{figure}

\subsection {Step 4. Non-singularity of pleating rays}\label{sec:nonsingular}

 This is the part of the argument which contains the deepest mathematics. Fortunately the results needed have been proved elsewhere.

 \begin{theorem}[\cite{ksqf, ksriley, chois}] \label{thm:nonsing}
 Suppose that $\gamma \in \pi_1({\mathcal S})$. Then ${\cal P}_{\gamma}$ is open and closed in the real trace locus $\mathbb R_{\gamma}$.
 Moreover $\Tr \rho_x(\gamma)$ is a local coordinate for $\mathbb C$  in a neighbourhood of ${\cal P}_{\gamma}$, and  is a global coordinate for ${\cal P}_{\gamma}$ on any non-empty connected component of ${\cal P}_{\gamma}$.
\end{theorem}
\begin{proof} The  statement  that   ${\cal P}_{\gamma}$ is open in $\mathbb R_{\gamma}$ is essentially \cite{ksriley} Proposition 3.1, see also   \cite{ksqf} Theorems 15  and 26.  The fact that $\Tr \rho_x(\gamma)$ is a local parameter is equivalent to the fact, also proved in both \cite{ksriley}  and   \cite{kstop}, that ${\cal P}_{\gamma}$ is a non-singular $1$-manifold. The open-ness and the final statement are actually  a
special case of  Theorems B and C  of \cite{chois} which state that for general hyperbolic manifolds, if the support of the bending lamination is a union of closed curves  is rational, then the traces of  these  curves are local parameters for the deformation space in a neighbourhood of  the corresponding pleating variety.

That   ${\cal P}_{\gamma}$ is closed in $\mathbb R_{\gamma}$ can be  proved as in~\cite{ksriley} Theorem 3.7. Here is a slightly more sophisticated version of the same idea. Suppose $x_n \to x_{\infty}$ with $x_n \in {\cal P}_{\gamma}$. The limit group
$G_{{\mathcal S}}(x_{\infty})$ is an algebraic limit of groups $G_{{\mathcal S}}(x_n)$ and hence the corresponding representation is discrete and faithful. Each of the two components of $(\partial { \mathcal C}/G_{{\mathcal S}}(x_n)) \setminus \gamma $ is a flat surface corresponding to a conjugacy class of Fuchsian subgroup $F_j(x_n), j=1,2$ (the $F$-peripheral subgroups of~\cite{ksriley}).  Since the limit is algebraic, $F_j(x_n)$ limits on a Fuchsian subgroup $F_j(x_{\infty})$, and similarly for all its conjugates in $G_{{\mathcal S}}(x_{\infty})$.

The  limit sets $\Lambda_{\alpha}$ of each of these subgroups $F_{\alpha}$ is spanned by a hyperbolic  plane $H_{\alpha}$ in $\mathbb H^3$. The Nielsen regions  of  $F_{\alpha}$ in $H_{\alpha}$ fit together along the
 lifts of the bending line $\gamma$ to $\mathbb H^3$, forming a pleated surface $\Pi$ in $\mathbb H^3$. We claim that $\Pi =  \partial { \mathcal C} (x_{\infty})$. This follows since the closure of the union of the $\Lambda_{\alpha}$ is the limit set of $G(x_{\infty})$, see also
 Proposition 7.2 in~\cite{kshowtobend}. The  result follows.   \end{proof}

\begin{remark} \rm{The closure of ${\cal P}_{\gamma}$ in $\mathbb R_{\gamma}$ is a simple case of both the `local limit theorem', Theorem 15 in~\cite{ksqf} and the  `lemme de fermeture' of~\cite{BonO}. These  much more sophisticated results allow that the bending lines may be part of an irrational lamination. Our argument above, in which the bending lamination is supported on closed curves,  is very close to that in the first part of the proof of Theor{\`e}me 6 in~\cite{otal}. }   \end{remark}

\begin{corollary}[\cite{kstop, ksriley, chois}]  If ${\cal P}_{\gamma} \neq \emptyset$, then it is a union of  connected non-singular branches of the real trace locus $\mathbb R_{ \gamma}$.
\end{corollary}
\begin{proof}  Suppose that ${\cal P}_{\gamma} \neq \emptyset$ and let $x \in {\cal P}_{\gamma}$, so that by Lemma~\ref{lemma:realtrace}, $x \in \mathbb R_{\gamma}$.  By Theorem~\ref{thm:nonsing},  ${\cal P}_{\gamma}$ is open and closed in $\mathbb R_{\gamma}$. Since $\Tr \gamma$ is a local coordinate, in a neighbourhood of $x$ the locus $\mathbb R_{\gamma}$
is a $1$-manifold.
\end{proof}

Notice that the theorem says that $\Tr \rho_x(\gamma)$ is a local parameter even in a neighbourhood of a cusp where $ \rho_x(\gamma)$ is parabolic,~\cite{ chois} Theorem C.
Thus we have
\begin{corollary} Suppose that $ x \in {\cal P}_{\gamma}$. Then there is a neighbourhood of $x$ in $\mathbb C$ on which $x \in \mathbb R_{\gamma}$ implies that $ x \in  \mathcal D$.
\end{corollary}

 \begin{corollary}\label{cor:unbounded} If ${\cal P}_{\gamma} \neq \emptyset$, then $\Tr \rho_x(\gamma)$ is unbounded on ${\cal P}_{\gamma}$.
\end{corollary}
\begin{proof}
Since $\Tr \rho_x(\gamma)$  is a local coordinate on connected components of ${\cal P}_{\gamma}$, this follows from the maximum principle on the branch, see~\cite{ksriley} Theorem 4.1.
\end{proof}

\subsection{Step 5. Finding the non-empty pleating rays} \label{sec:general rays}

Now we determine the pleating rays. As above, let ${\cal P}_{p/q}$ denote the ray corresponding to
the curve $\gamma_{p/q}$ and write $\mathbb R_{p/q}$ for the real locus of $\Tr V_{p/q}$. From Proposition~\ref{curveequality} we have ${\cal P}_{p/q} = {\cal P}_{(p+2q)/q}= {\cal P}_{-p/q}$.

By~\ref{sec:nonsingular}, ${\cal P}_{p/q}  $ is a union of non-singular branches of $\mathbb R_{p/q}$.
We  now find those $p/q \neq \{0,1\}$ for which ${\cal P}_{p/q} \neq \emptyset$,  at the same time resolving the connectivity issue.   We follow the method of \cite{ksriley}, using  an inductive argument on position of the pleating rays and their asymptotic directions as $|x| \to \infty$, and at the same time correcting the second of the two errors referred to in Remark~\ref{Komori issue}.   We have:
 \begin{proposition}[c.f.~\cite{ksriley} Theorem 4.1]\label{prop:asympdirn}
The set ${\cal P}_{p/q}$ is the union of the two branches  of $\mathbb R_{p/q} $ which are asymptotic to the half lines $\rho e^{\pm i \pi (p-q)/q}$ as $\rho \to \infty$.
 \end{proposition}
\begin{proof}
Denote by $R(\theta)$ the ray  $te^{i \theta}, t >0$, in the $x$-plane. By Proposition~\ref{prop:tracepoly}, $\Tr V_{p/q}$ is a polynomial in $x$ whose top term is $(-1)^{p-q-1} x^q$. Now $\Tr V_{p/q}$ takes real values on ${\cal P}_{p/q}$,  moreover by Corollary~\ref{cor:unbounded} it is unbounded on ${\cal P}_{p/q}$. It follows that
as $|x| \to \infty$, ${\cal P}_{p/q}$ must be asymptotic to one of the rays $R(k\pi/q), k \in \mathbb Z$.

We have already identified ${\cal P}_{0/1}$ and ${\cal P}_{1/1}$ as the real intervals $(-\infty, -3]$ and $[2,\infty)$   respectively. It follows from Section~\ref{sec:dirichlet} that the  semicircular arc from $-4$ to $4$ (say) in $\mathbb H$ is a continuous path in $ \mathcal D$ from
${\cal P}_{0/1}$ to ${\cal P}_{1/1}$. Hence by the continuity theorem of~\cite{kscont},  if  $0 < p/q <1$ there is a point on  ${\cal P}_{p/q}$ in the upper half plane $\mathbb H$. Likewise there is a point
on  ${\cal P}_{p/q}$ in the lower half plane. (This was missed in~\cite{ksriley}.) Since ${\cal P}_{0/1} \cup {\cal P}_{1/1}$ separates $ \mathcal D$ into two connected components, this shows in particular that  ${\cal P}_{p/q}$ must have at least two connected components.

 Now we proceed by induction on the Farey tree. Suppose we have shown the result for two Farey neighbours $p/q, r/s$. Consider the locus ${\cal P}_{p+r/q+s}$. By the inductive hypothesis,  $\mathbb H$ contains exactly one component  of  each of  ${\cal P}_{p/q}, {\cal P}_{r/s}$, asymptotic to the rays  $R(\pi (p-q)/q), R( \pi (r-s)/s)$ respectively. Exactly as in~\cite{ksriley}  it is easy to check that there is exactly one integer $k \in \{0, 1, \ldots, 2(q+s)-1 \}$ for which $R(k\pi/(q+s))$ lies between $R(\pi (p-q)/q)$ and $ R( \pi (r-s)/s)$, namely $k= (p+r)/(q+s)$. By the same continuity theorem as before, a path in this sector joining suitable points on ${\cal P}_{p/q}, {\cal P}_{r/s}$ must meet ${\cal P}_{p+r/q+s}$. Thus ${\cal P}_{p+r/q+s}$ has at least one connected component asymptotic to $ R( \pi (p + r -q-s)/(q+s))$.
 A similar argument in the lower half plane gives another connected component asymptotic to $ R( \pi (p + r +q+s)/(q+s))$. Since ${\cal P}_{p+r/q+s}$ has exactly two components by Proposition~\ref{prop:connectivity} below, the result follows.
\end{proof}

The issue of connectivity of ${\cal P}_{\gamma}$ is a bit subtle. In the general theory, see~\cite{BonO, chois}, one shows that ${\cal P}_{\gamma}$ has one connected component.   However this result holds in a space of manifolds which are consistently oriented throughout the space and
all of whose  convex cores have non-zero volume. In our case we have:

\begin{proposition} \label{prop:connectivity} If $\gamma  \neq 0/1, 1/1$ and ${\cal P}_{\gamma} \neq \emptyset$,  then  ${\cal P}_{\gamma}$ has exactly two connected components in $ \mathcal D$.
\end{proposition}
\begin{proof}
The usual argument that the pleating ray of a rational lamination  has one connected component goes as follows.  Given a point on ${\cal P}_{\gamma}$, double the convex core along its boundary to obtain a cone manifold with a singular axis of angle $2(\pi - \theta)$ along $\gamma$, where $\theta$ is the bending angle along $\gamma$. (Notice that the convention on defining bending angles differs between papers by the first author and~\cite{BonO}. In our convention, a bending line contained in  flat subsurface has bending angle $0$ but cone angle $2 \pi$, whereas in~\cite{BonO},  the bending angle along a line in a flat  surface is defined to be $\pi$.) By \cite{HK},  such a  hyperbolic cone manifold is  parametrized by its cone angle. One shows that one can continuously deform the cone angle to $0$, at which point the curve whose axis is the bending line has to become parabolic.  The doubled manifold is an oriented hyperbolic manifold with a rank two cusp and finite volume. As long as we are working in a space in which all manifolds have consistent orientation, such a manifold  is  unique up to orientation preserving isometry, from which one deduces that ${\cal P}_{\gamma}$ is connected.

In our case, the parameter space $ \mathcal D$ is separated by two lines along which $G$ is Fuchsian so that ${ \mathcal C}(G)/G$ has zero volume and  the above argument fails.
Note however that, provided that $G$ is not Fuchsian,   ${\mathcal S}$ can be oriented by the triple consisting of the \emph{oriented} axes of $P,Q$ and the oriented line $C$ from $ \Ax  K_1 $ to $ \Ax  K_1 $.
The map $ \zeta \to \bar \zeta$ reverses the relative orientations of $ \Ax P,  \Ax Q$ while fixing that of $C$.   Thus $ \mathcal D \setminus \mathbb R$ has two connected components in which ${\mathcal S}$ has naturally opposite orientations. The above argument shows that  ${\cal P}_{\gamma} $ has at
most  one component  in each component of $ \mathcal D$.
Since we have already shown in Proposition~\ref{prop:asympdirn}
 that ${\cal P}_{\gamma} $ has at least one component in each of the upper and lower half planes, this completes the proof.

This proposition can alternatively be proved by the more \emph{ad hoc} methods used in~\cite{ksriley}.
\end{proof}

\begin{remark}
\rm {Proposition~\ref{prop:asympdirn} shows that ${\cal P}_{p/q} \neq \emptyset$ for all $p/q \in \mathbb Q$. This can be viewed as a special case of the general result of~\cite{BonO} Theorem 1,
see also~\cite{chois} Theorem 2.4.  We have to be careful to include the case, excluded in \cite{BonO}, that the group $G_{{\mathcal S}}(x)$ is Fuchsian so that ${ \mathcal C}/G$ has zero volume.
The conclusion is the following:
\begin{proposition} Let $\gamma$  be an essential simple non-peripheral closed curve on $\partial {\mathcal S}$. Then
${\cal P}_{\gamma} \neq \emptyset$ if and only if $\gamma$ is non-trivial in $\pi_1({\mathcal S})$ and intersects the meridian disk $\gamma_{1/0}$ at least twice.
If $\gamma$ meets  $\gamma_{1/0}$ exactly twice then the bending angle is identically $\pi$ and $G_{{\mathcal S}}(x)$ is Fuchsian.
\end{proposition}
}
\end{remark}

\subsection{ Step 6. Density of rational pleating rays}

Finally, we justify the claim that  the rational pleating rays are dense in $ \mathcal D$:
\begin{theorem}[\cite{kstop} Corollary 6.2, \cite{ksriley} Theorem 5.2] \label{thm:density} Rational pleating rays  are dense in $ \mathcal D_{{\mathcal S}}$.
\end{theorem}
\begin{proof} The proof of this result in any one complex dimensional parameter space   is the same. Here is a quick sketch.
Suppose that $\nu$ is an irrational lamination  with corresponding pleating variety ${\cal P}_{\nu}$, and that $x \in {\cal P}_{\nu} \cap  \mathcal D$. Pick a sequence of rational measured laminations $\nu_n = c_n \delta_{\gamma_n}$ where $c_n \in \mathbb R^+$ so that $\nu_n  \to \nu $ in the space of projective measured laminations on $\partial {\mathcal S}$, where $\delta_{\gamma_n}$ is the unit point mass on $\gamma_n$. Replace the traces of $\gamma_n$ by complex length functions $\lambda_n$ and  scale to get complex analytic functions $c_n \lambda_n$. One shows that
in a neighbourhood of $x \in {\cal P}_{\nu}$ these functions form a normal family  which converges to a non-constant analytic function (\cite {kstop} Theorem 20),   whose real locus
contains the pleating ray ${\cal P}_{\nu}$ (\cite {kstop} Theorem 23).
By   Hurwitz' theorem,  there are nearby points at which the approximating functions $c_n \lambda_n$ must take on real values. In a small enough neighbourhood of $x$, this is enough to force $ y \in {\cal P}_{\gamma_n}$ (\cite {kstop} Theorem 31). This gives density in $\mathop{{Int}}  \mathcal D$. By the result quoted in the introduction that $ \mathcal D = \overline{\mathop{{Int}}  \mathcal D}$ we are done. \end{proof}

  \subsection{The pleating rays for ${\cal H}$}
\label{sec:rays}

By Corollary~\ref{discretetogether}, $ \mathcal D_{{\cal H}}=  \mathcal D_{{\mathcal S}}$. Thus  the rational rays for $ \mathcal D_{{\mathcal S}}$ are also dense in $ \mathcal D_{{\cal H}}$. However it is easy to see that
 a rational pleating laminations on  $\partial {\cal H}(x)$ correspond exactly to those on ${\mathcal S}(x)$,  and that although the actual bending curves differ, their traces are related by a simple formula.

\begin{lemma} \label{bendlinesagree} Suppose that the bending lamination  $\beta_{{\cal H}}(x)$ of ${\cal H}(x)$ is rational so that its support $\lambda$ is a union of disjoint simple closed curves on $\partial {\cal H}$.
Let $\gamma$ be a connected component of
$\lambda$. Then either ${\kappa}(\gamma) = \gamma$ or the three curves $\gamma, {\kappa}(\gamma) , {\kappa}^2(\gamma)$
are disjoint. The support of the bending lamination   $\beta_{{\mathcal S}}(x)$ is exactly the projection of $\gamma$ to ${\mathcal S}$ and all rational bending laminations of ${\mathcal S} $ arise in this way.
\end{lemma}
\begin{proof} The limit set of $G_{{\cal H}}(x)$ and hence its convex core are invariant under the symmetry  ${\kappa}$.
Hence  the support $\lambda$ of $\beta_{{\cal H}}(x)$   is also ${\kappa}$-invariant.
Let $\gamma$ be a connected component of
$\lambda$. Since connected components of $\lambda$ are pairwise disjoint, either ${\kappa}(\gamma) = \gamma$ or the three curves $\gamma, {\kappa}(\gamma), {\kappa}^2(\gamma)$
are disjoint.
In either case,  $\gamma$ cannot pass through a fixed point of ${\kappa}$: at the fixed point $P$ the images of $\gamma$ would meet at angles $2 \pi/3$ so that $\gamma, {\kappa}(\gamma) , {\kappa}^2(\gamma)$ would intersect at $P$, which is impossible.

Let $\pi_{{\kappa}}$ be the projection ${\cal H} \to {\mathcal S}$. In a neighbourhood of a bending line $\pi_{{\kappa}}$ is a covering map hence a local isometry. Since being a bending line can be characterised locally, $\beta_{{\mathcal S}}(x)$ is the projection of $\gamma$ to ${\mathcal S}$.

Let $\gamma$ be a simple closed curve on $\partial {\mathcal S}$.
Clearly, by the same observation about local characterisation of bending lines,  if $\gamma$ is a bending line then so is any connected component of its lift  to $\partial {\cal H}$. This proves the converse.
 \end{proof}

 We remark that if $p/q$ is congruent to $1/0$ or $0/1$ mod $\mathbb Z_2$ then the lift of $\gamma_{p/q}$ has three connected components which are permuted among themselves by ${\kappa}$, while if $p/q$ is congruent to $1/1$ then its lift  has one ${\kappa}$-invariant connected component.
 To see this, check by hand for the curves $\gamma_{1/0}, \gamma_{0/1}, \gamma_{1/1}$ and then note that the lifting property is invariant under the mapping class group of $\partial {\mathcal S}$ which at the same time acts transitively on $p/q$   congruence classes mod $\mathbb Z_2$.

To actually compute the pleating rays for $ \mathcal D_{{\mathcal S}}$, we  computed the traces $\Tr V_{p/q}(x)$ corresponding to the  curves $\gamma_{p/q} \in \pi_1({\mathcal S})$. The above discussion shows that it is unnecessary to actually compute traces of lifted curves   in $\pi_1({\cal H})$. If for some reason one wanted to do this, either one could start again enumerating the curves on ${\cal H}$, or one could note that the complex length of a lift of $\gamma_{p/q}$ in ${\cal H}$ would be either the same as or three times that of the curve $\gamma_{p/q}$ in $  {\mathcal S}$, depending on the $\mathbb Z_2$-parity of $p/q$.

 \section{Computing traces} \label{sec:torustree}

  To compute traces of the elements $V_{p/q}$, rather than use the ${\mathcal S}$-tree as in Section~\ref{sec:traces},  we actually  performed computations on the  associated Markoff tree corresponding to the associated torus ${\cal T}$ of Section~\ref{sec:solidtorus}, referred to in this section as the ${\cal T}$-tree. To justify this, we need to compare the curves in Farey position $p/q$ on the two trees to ensure that they do indeed correspond geometrically as expected. We also need to address the issue about lifting representations  to $SL(2,\mathbb C)$ raised in Remark~\ref{signchoices1}.

  \subsubsection{Correspondence of curves}\label{sec:correspondence}

Homotopy classes of essential simple non-peripheral  loops on $\partial {\cal T}$ are  well known to be in bijective correspondence to unoriented lines of rational slope  in the plane, see for example~\cite{serint,kstop}.
 In fact the word $W_{p/q}$ generated by the concatenation process following the ${\cal T}$-tree described in Section~\ref{sec:markoff} is the cutting sequence of a line of slope $p/q \in \hat \mathbb Q$ across the lattice, see~\cite{serint}.
 \begin{figure}[hbt] \begin{center}
 \includegraphics[height=6cm]{Figs/slope}
 \caption{Lattice representation of a cover of $\partial {\mathcal S}$. The integer  vertices (white circles) correspond to the end points of the order $3$ axes on $\partial {\mathcal S}$; the endpoints of the order $2$ elliptics $P,Q, R$ are coloured blue, green and red respectively.} \label{fig:lattice}
\end{center}\end{figure}

The key point here is that the plane with a cone singularity of angle $2\pi/3$ at integer lattice points, see Figure~\ref{fig:lattice}, is an intermediate covering between the universal cover $\mathbb H$ of $\partial {\cal T}$ and $\partial {\cal T}$ itself.   As described in for example~\cite{ksriley}, the same lattice can also be viewed as an intermediate covering between
$\mathbb H$ and $\partial {\mathcal S}$: the rectangle with vertices at $0, 1, 2i, 2i+1$  projects bijectively to $\partial {\mathcal S}$ while the rectangle with vertices   $0, 1,   i/2, 1+  i/2$ projects bijectively to $\partial {\mathcal U}$ and the unit square projects  bijectively to   the torus $\partial {\cal T}$. The lattice points correspond to the cone points belonging to $K_i, i=1, \ldots, 4$ arranged as shown. Thus there  is also a bijective correspondence between lines of rational slope in the punctured plane and simple essential non-peripheral curves on  $\partial {\mathcal S}$. In this way, one can easily relate the words $W_{p/q}$ (on $\partial {\cal T}$) and $ V_{p/q}$ (on $\partial {\mathcal S}$); this is explained in detail in~\cite{ksriley}.

In this picture,  the meridian loop $\partial m$ of Section~\ref{sec:enumeration} is identified as the `vertical' line of slope $1/0$.
One sees easily that the line of slope $p/q$ in the plane projects to a curve on $\partial {\mathcal S}$ which has exactly $2q$ intersections with $\partial m$ and a twist of $p$ as described in Section~\ref{sec:enumeration}.   It follows from Lemma~\ref{lem:treeisok}  that  the labelling of curves by lines of rational slope $p/q$  exactly corresponds  to  the Farey labelling of curves  by their position on the ${\mathcal S}$-tree.

As above, the curve in Farey position $p/q$ on the ${\mathcal S}$-tree is denoted $\gamma_{p/q}$, corresponding to a word $V_{p/q}$; while the curve in Farey position $p/q$ on the ${\cal T}$-tree is denoted $\omega_{p/q}$, corresponding to a word $W_{p/q}$. Now $  {\mathcal S}$ projects to ${\mathcal U}$ by a four-fold cover and  ${\cal T}$ projects to ${\mathcal U}$ by a two-fold cover. Hence we have:

\begin{proposition}\label{coverings}
The complex length of $\gamma_{p/q} $ is twice that of  $\omega_{p/q}$, hence $   \Tr V_{p/q}(\zeta) = \pm (\Tr W_{p/q} (\zeta))^2 -2 )$. \end{proposition}
Note that this allows for an ambiguity in the signs of the traces since the two lifts of $\pi_1({\cal T})$ and $\pi_1({\mathcal S})$ to $SL(2,\mathbb C)$ are not (indeed cannot be) chosen consistently.
\begin{corollary} Up to sign, the trace of  $\gamma_{p/q} \in \pi_1({\mathcal U})$   may be computed using the formula of Proposition~\ref{coverings} and  the ${\cal T}$-tree.  \end{corollary}

Since we are aiming to compute pleating rays which are a geometrical construct and hence only depend on a $PSL(2,\mathbb C)$ representation, this would be sufficient for our purposes. However it is more  satisfying    to prove the following more precise result which  shows that working with the $SL(2,\mathbb C)$ lift of the representation of $\pi_1({\cal T})$ described in Section~\ref{sec:solidtorus}, we can fix the choice of sign.

 \begin{prop} \label{relatingtraces} With $W_{p/q}, V_{p/q}$  as above, let $f_{p/q} (z) = \Tr V_{p/q}(\zeta)$  and  $g_{p/q}  (\zeta)  = \Tr W_{p/q} (\zeta) $.  Then $-f_{p/q} (\zeta) =(g_{p/q}  (\zeta))^2 -2$ for all $p/q \in \hat \mathbb Q$.   \end{prop}
 \begin{proof} It is easy to check that this is correct for $p/q = 0/1,1/0, 1/1$.
 In detail:

 $\omega_{0/1} = A, \gamma_{0/1}  = K_0K_1$ and we have shown that $  A^2 = -K_0K_1$.
Thus $ f_{0/1} (\zeta) =x$, $(g_{p/q}  (\zeta))^2 -2 = (-x+2)-2 = -x$.

 $\omega_{1/0} =B, \gamma_{1/0}  = \rm{id}$ and $B^2 = - \rm{id}$.
So $ f_{1/0} (\zeta) =2$, $(g_{1/0}  (\zeta))^2 -2 = -2$.

  $\omega_{1/1} =AB, \gamma_{1/1}  = \rm{K_0K_1^{-1}}$.
So $ f_{1/1} (\zeta) =1-x$, $(g_{1/1}  (\zeta))^2 -2 = x-1$.

Now we do an inductive proof.
Suppose that in the ${\mathcal S}$-tree  labels $u,v$ are adjacent along an edge $e$ with $w$ the remaining label at one of the two vertices at the ends of $e$. By the formula in Section~\ref{sec:traces}
 the  label at the other vertex is $2 - uv -w$.

 Suppose that the corresponding labels on the ${\cal T}$-tree  are $u',v',w'$. Then the remaining label at the vertex at the other end of $e$ is  $ u'v'-w'$. Replace these labels by the negatives of the traces of the doubled curves to get labels
 $2-u'^2, 2-v'^2, 2-w'^2 ,2- (u'v'-w')^2 $ around the same $4$ vertices. If we can show that
 $$2 -(2-u'^2)(2-v'^2) - (2-w'^2 )  = -   (2- (u'v'-w')^2 )    $$ we will be done.
This is easily checked by multiplying out, noting that the trace identity~\eqref{conedspheretraces} round a vertex of the ${\cal T}$-tree gives
 $$ u'^2 +   v'^2   + w'^2 = u'v'w' + \Tr [A,B] +2 =  u'v'w' + 3.$$  \end{proof}

  \subsubsection{The actual computations}\label{sec:torustree1}

The above discussion justifies the method we actually used to perform  computations involving traces on ${\mathcal S}$.
Instead of computing on the ${\mathcal S}$-tree with initial traces   $\Tr \gamma_{0/1} =  \Tr  K_0K_1 = x, \Tr \gamma_{1/0}= \Tr   \mbox{\rm{id}}= 2, \Tr \gamma_{1/1} = \Tr K_0K_1^{-1}  = 1-x$, we used the
 ${\cal T}$-tree  with initial triple
$\Tr A = \pm i(3/\zeta - \zeta),  \Tr B =0$ and $\Tr AB = \pm i(3/\zeta + \zeta)$ corresponding to the generators $A,B$ of $G_{{\cal T}}$.  As in Section~\ref{sec:solidtorus}
 Observe that $A^2 = -K_0K_1$, so that $\Tr A^2 = -x$. Since $\Tr B = 0$, we can find $\Tr AB$ from the  identity
$(\Tr A)^2 +  (\Tr AB)^2 = \Tr {[A,B]} +2 = 3$.
Thus setting $(a,b,c) = (\Tr A, \Tr B, \Tr AB)$ we have
$$
a^2 - 2 = -x, \quad c^2  =   1 + x.
$$
It is easily checked that this is in accord with~\eqref{eqn:torustraces}.
Thus associated to $G_{{\cal T}}(x)$  we have the torus tree
 $(a,0,c) = (\sqrt{-x+2}, 0, \sqrt{x+1})$. This is the method we actually used to compute the
 pleating rays shown in Figure~\ref{Diagonal-and-Riley-mu-3-Ray-BQ}.

\begin{remark} \label{rmk:signs} \rm{The sign of the square roots in the above can be  uniquely determined by the formulae for traces in terms of $\xi$. What we actually did was to make an arbitrary choice and plot rays corresponding to curves in the range $0 \leq p/q \leq 1$, thus making a picture in the upper half plane which we could then reflect. As can be seen from Figure~\ref{fig:trace2}, the  signs of the square roots in fact  alternate periodically with period $4$ rather than period $2$, so that, for example, $\Tr \gamma_{3} = -\Tr \gamma_{1}$.}
\end{remark}

\subsubsection{Computations for the Riley slice} \label{sec:discussion}  The traces needed to find the pleating  rays for the Riley slice on the right  in Figure~\ref{Figs/Riley-Ray-BQ} were computed by a method similar to that described above.   Our parameter $x$ can be related to the parameter $\rho$ of~\cite{ksriley}
 by comparing the traces of the word in Farey position $0/1$: these are $K_0K_1$ in our case and $XY$ in the notation of~\cite{ksriley}.   Thus we find   the correct correspondence is $ x \leftrightarrow \rho + 2 $.
For the Riley group  a similar computation to the one above with $\Tr [A,B] = -2$ gives immediately
 $(\Tr A)^2 = -(\rho + 2)$ and  $(\Tr AB)^2 =  \rho + 2 $. Thus writing in terms of the $x$-coordinate we find the initial triple $(\sqrt {-x}, 0, \sqrt {x})$.

\subsubsection{Comparison of Bowditch sets}  It is interesting to compare the Bowditch sets associated to  the two  initial triples $(x,x,x)$ and $(\sqrt{-x+2}, 0, \sqrt{x+1})$.  In the latter case, one needs to modify the definition of the Bowditch set: since $\phi(U)=0$ for some $U \in \Omega$, there is a trace preserving $\mathbb Z$-action on the associated tree $\mathbb T_{(\sqrt{-x+2}, 0, \sqrt{x+1})}$ corresponding to the action of a subgroup of ${\rm Aut}(F_2)$ generated by a parabolic, see  for example \cite{tan_gd} Theorem 1.9.   The Bowditch condition should actually be specified  on $\Omega \setminus \{U\} / \sim$, where $\sim$ is   the equivalence coming from this symmetry.

The results, plotted in the $\zeta$-plane, are shown in Figure~\ref{fig:BQ-sets-comparison}. On the right  the initial triple is  $(x,x,x)$  (with $x$ related to $\xi$ as in \eqref{eqn:paramreln}) corresponding to the handlebody group $G_{{\cal H}}(x)$.
On the left,  the initial triple is $(\sqrt{-x+2}, 0, \sqrt{x+1})$  corresponding to the torus group $G_{{\cal T}}(x)$.
The two regions are clearly distinct: the grey region on the right contains that on the left.
Conjecturally, the left hand grey region is also the discreteness locus for the groups $G_{{\cal T}}(x)$, see Figure \ref{Figs/Riley-Ray-BQ} for the parametrization in terms of $x$.

Note the various symmetries as discussed in Section~\ref{sec:symmetries}, in particular note how  Figure~\ref{Diagonal-BQ} loses the left-right reflectional symmetry seen in Figure~\ref{fig:BQ-sets-comparison}.
The coloured region in Figure~\ref{Diagonal-BQ} is the same region as the right frame of Figure~\ref{fig:BQ-sets-comparison},  drawn in the $x$-plane.

\begin{figure}[ht]\label{xxx}
\includegraphics[width=4cm]{Figs/z.pdf}
\hspace{1cm}
\includegraphics[width=4cm]{Figs/xxx.pdf}
\caption{Bowditch sets (grey) plotted in the $\zeta$-plane with range $[-4,4]\times[-4i,4i]$.  Left:  Initial triple  $(\sqrt{-x+2}, 0, \sqrt{x+1})$
corresponding to the torus group $G_{{\cal T}}(x)$. Right: Initial triple  $(x,x,x)$ corresponding to the handlebody group $G_{{\cal H}}(x)$.
The two regions are clearly distinct: the grey region on the right contains that on the left. }
\label{fig:BQ-sets-comparison}
\end{figure}

\bibliographystyle{alpha}

The diagonal slice of Schottky space