# Category: geometry

A Belyi-extender (or dessinflateur) is a rational function $q(t) = \frac{f(t)}{g(t)} \in \mathbb{Q}(t)$ that defines a map
$q : \mathbb{P}^1_{\mathbb{C}} \rightarrow \mathbb{P}^1_{\mathbb{C}}$
unramified outside $\{ 0,1,\infty \}$, and has the property that $q(\{ 0,1,\infty \}) \subseteq \{ 0,1,\infty \}$.

An example of such a Belyi-extender is the power map $q(t)=t^n$, which is totally ramified in $0$ and $\infty$ and we clearly have that $q(0)=0,~q(1)=1$ and $q(\infty)=\infty$.

The composition of two Belyi-extenders is again an extender, and we get a rather mysterious monoid $\mathcal{E}$ of all Belyi-extenders.

Very little seems to be known about this monoid. Its units form the symmetric group $S_3$ which is the automrphism group of $\mathbb{P}^1_{\mathbb{C}} – \{ 0,1,\infty \}$, and mapping an extender $q$ to its degree gives a monoid map $\mathcal{E} \rightarrow \mathbb{N}_+^{\times}$ to the multiplicative monoid of positive natural numbers.

If one relaxes the condition of $q(t) \in \mathbb{Q}(t)$ to being defined over its algebraic closure $\overline{\mathbb{Q}}$, then such maps/functions have been known for some time under the name of dynamical Belyi-functions, for example in Zvonkin’s Belyi Functions: Examples, Properties, and Applications (section 6).

Here, one is interested in the complex dynamical system of iterations of $q$, that is, the limit-behaviour of the orbits
$\{ z,q(z),q^2(z),q^3(z),… \}$
for all complex numbers $z \in \mathbb{C}$.

In general, the 2-sphere $\mathbb{P}^1_{\mathbb{C}} = S^2$ has a finite number of open sets (the Fatou domains) where the limit behaviour of the series is similar, and the union of these open sets is dense in $S^2$. The complement of the Fatou domains is the Julia set of the function, of which we might expect a nice fractal picture.

Let’s take again the power map $q(t)=t^n$. For a complex number $z$ lying outside the unit disc, the series $\{ z,z^n,z^{2n},… \}$ has limit point $\infty$ and for those lying inside the unit circle, this limit is $0$. So, here we have two Fatou domains (interior and exterior of the unit circle) and the Julia set of the power map is the (boring?) unit circle.

Fortunately, there are indeed dynamical Belyi-maps having a more pleasant looking Julia set, such as this one

But then, many dynamical Belyi-maps (and Belyi-extenders) are systems of an entirely different nature, they are completely chaotic, meaning that their Julia set is the whole $2$-sphere! Nowhere do we find an open region where points share the same limit behaviour… (the butterfly effect).

There’s a nice sufficient condition for chaotic behaviour, due to Dennis Sullivan, which is pretty easy to check for dynamical Belyi-maps.

A periodic point for $q(t)$ is a point $p \in S^2 = \mathbb{P}^1_{\mathbb{C}}$ such that $p = q^m(p)$ for some $m > 1$. A critical point is one such that either $q(p) = \infty$ or $q'(p)=0$.

Sullivan’s result is that $q(t)$ is completely chaotic when all its critical points $p$ become eventually periodic, that is some $q^k(p)$ is periodic, but $p$ itself is not periodic.

For a Belyi-map $q(t)$ the critical points are either comlex numbers mapping to $\infty$ or the inverse images of $0$ or $1$ (that is, the black or white dots in the dessin of $q(t)$) which are not leaf-vertices of the dessin.

Let’s do an example, already used by Sullivan himself:
$q(t) = (\frac{t-2}{t})^2$
This is a Belyi-function, and in fact a Belyi-extender as it is defined over $\mathbb{Q}$ and we have that $q(0)=\infty$, $q(1)=1$ and $q(\infty)=1$. The corresponding dessin is (inverse images of $\infty$ are marked with an $\ast$)

The critical points $0$ and $2$ are not periodic, but they become eventually periodic:

$2 \rightarrow^q 0 \rightarrow^q \infty \rightarrow^q 1 \rightarrow^q 1$
and $1$ is periodic.

For a general Belyi-extender $q$, we have that the image under $q$ of any critical point is among $\{ 0,1,\infty \}$ and because we demand that $q(\{ 0,1,\infty \}) \subseteq \{ 0,1,\infty \}$, every critical point of $q$ eventually becomes periodic.

If we want to avoid the corresponding dynamical system to be completely chaotic, we have to ensure that one of the periodic points among $\{ 0,1,\infty \}$ (and there is at least one of those) must be critical.

Let’s consider the very special Belyi-extenders $q$ having the additional property that $q(0)=0$, $q(1)=1$ and $q(\infty)=\infty$, then all three of them are periodic.

So, the system is always completely chaotic unless the black dot at $0$ is not a leaf-vertex of the dessin, or the white dot at $1$ is not a leaf-vertex, or the degree of the region determined by the starred $\infty$ is at least two.

Going back to the mystery Manin-Marcolli sub-monoid of $\mathcal{E}$, it might explain why it is a good idea to restrict to very special Belyi-extenders having associated dessin a $2$-coloured tree, for then the periodic point $\infty$ is critical (the degree of the outside region is at least two), and therefore the conditions of Sullivan’s theorem are not satisfied. So, these Belyi-extenders do not necessarily have to be completely chaotic. (tbc)

A Belyi-extender (or dessinflateur) $\beta$ of degree $d$ is a quotient of two polynomials with rational coefficients
$\beta(t) = \frac{f(t)}{g(t)}$
with the special properties that for each complex number $c$ the polynomial equation of degree $d$ in $t$
$f(t)-c g(t)=0$
has $d$ distinct solutions, except perhaps for $c=0$ or $c=1$, and, in addition, we have that
$\beta(0),\beta(1),\beta(\infty) \in \{ 0,1,\infty \}$

Let’s take for instance the power maps $\beta_n(t)=t^n$.

For every $c$ the degree $n$ polynomial $t^n – c = 0$ has exactly $n$ distinct solutions, except for $c=0$, when there is just one. And, clearly we have that $0^n=0$, $1^n=1$ and $\infty^n=\infty$. So, $\beta_n$ is a Belyi-extender of degree $n$.

A cute observation being that if $\beta$ is a Belyi-extender of degree $d$, and $\beta’$ is an extender of degree $d’$, then $\beta \circ \beta’$ is again a Belyi-extender, this time of degree $d.d’$.

That is, Belyi-extenders form a monoid under composition!

In our example, $\beta_n \circ \beta_m = \beta_{n.m}$. So, the power-maps are a sub-monoid of the Belyi-extenders, isomorphic to the multiplicative monoid $\mathbb{N}_{\times}$ of strictly positive natural numbers.

In their paper Quantum statistical mechanics of the absolute Galois group, Yuri I. Manin and Matilde Marcolli say they use the full monoid of Belyi-extenders to act on all Grothendieck’s dessins d’enfant.

But, they attach properties to these Belyi-extenders which they don’t have, in general. That’s fine, as they foresee in Remark 2.21 of their paper that the construction works equally well for any suitable sub-monoid, as long as this sub-monoid contains all power-map exenders.

I’m trying to figure out what the maximal mystery sub-monoid of extenders is satisfying all the properties they need for their proofs.

But first, let us see what Belyi-extenders have to do with dessins d’enfant.

In his user-friendlier period, Grothendieck told us how to draw a picture, which he called a dessin d’enfant, of an extender $\beta(t) = \frac{f(t)}{g(t)}$ of degree $d$:

Look at all complex solutions of $f(t)=0$ and label them with a black dot (and add a black dot at $\infty$ if $\beta(\infty)=0$). Now, look at all complex solutions of $f(t)-g(t)=0$ and label them with a white dot (and add a white dot at $\infty$ if $\beta(\infty)=1$).

Now comes the fun part.

Because $\beta$ has exactly $d$ pre-images for all real numbers $\lambda$ in the open interval $(0,1)$ (and $\beta$ is continuous), we can connect the black dots with the white dots by $d$ edges (the pre-images of the open interval $(0,1)$), giving us a $2$-coloured graph.

For the power-maps $\beta_n(t)=t^n$, we have just one black dot at $0$ (being the only solution of $t^n=0$), and $n$ white dots at the $n$-th roots of unity (the solutions of $x^n-1=0$). Any $\lambda \in (0,1)$ has as its $n$ pre-images the numbers $\zeta_i.\sqrt[n]{\lambda}$ with $\zeta_i$ an $n$-th root of unity, so we get here as picture an $n$-star. Here for $n=5$:

This dessin should be viewed on the 2-sphere, with the antipodal point of $0$ being $\infty$, so projecting from $\infty$ gives a homeomorphism between the 2-sphere and $\mathbb{C} \cup \{ \infty \}$.

To get all information of the dessin (including possible dots at infinity) it is best to slice the sphere open along the real segments $(\infty,0)$ and $(1,\infty)$ and flatten it to form a ‘diamond’ with the upper triangle corresponding to the closed upper semisphere and the lower triangle to the open lower semisphere.

In the picture above, the right hand side is the dessin drawn in the diamond, and this representation will be important when we come to the action of extenders on more general Grothendieck dessins d’enfant.

Okay, let’s try to get some information about the monoid $\mathcal{E}$ of all Belyi-extenders.

What are its invertible elements?

Well, we’ve seen that the degree of a composition of two extenders is the product of their degrees, so invertible elements must have degree $1$, so are automorphisms of $\mathbb{P}^1_{\mathbb{C}} – \{ 0,1,\infty \} = S^2-\{ 0,1,\infty \}$ permuting the set $\{ 0,1,\infty \}$.

They form the symmetric group $S_3$ on $3$-letters and correspond to the Belyi-extenders
$t,~1-t,~\frac{1}{t},~\frac{1}{1-t},~\frac{t-1}{t},~\frac{t}{t-1}$
You can compose these units with an extender to get anther extender of the same degree where the roles of $0,1$ and $\infty$ are changed.

For example, if you want to colour all your white dots black and the black dots white, you compose with the unit $1-t$.

Manin and Marcolli use this and claim that you can transform any extender $\eta$ to an extender $\gamma$ by composing with a unit, such that $\gamma(0)=0, \gamma(1)=1$ and $\gamma(\infty)=\infty$.

That’s fine as long as your original extender $\eta$ maps $\{ 0,1,\infty \}$ onto $\{ 0,1,\infty \}$, but usually a Belyi-extender only maps into $\{ 0,1,\infty \}$.

Here are some extenders of degree three (taken from Melanie Wood’s paper Belyi-extending maps and the Galois action on dessins d’enfants):

with dessin $5$ corresponding to the Belyi-extender
$\beta(t) = \frac{t^2(t-1)}{(t-\frac{4}{3})^3}$
with $\beta(0)=0=\beta(1)$ and $\beta(\infty) = 1$.

So, a first property of the mystery Manin-Marcolli monoid $\mathcal{E}_{MMM}$ must surely be that all its elements $\gamma(t)$ map $\{ 0,1,\infty \}$ onto $\{ 0,1,\infty \}$, for they use this property a number of times, for instance to construct a monoid map
$\mathcal{E}_{MMM} \rightarrow M_2(\mathbb{Z})^+ \qquad \gamma \mapsto \begin{bmatrix} d & m-1 \\ 0 & 1 \end{bmatrix}$
where $d$ is the degree of $\gamma$ and $m$ is the number of black dots in the dessin (or white dots for that matter).

Further, they seem to believe that the dessin of any Belyi-extender must be a 2-coloured tree.

Already last time we’ve encountered a Belyi-extender $\zeta(t) = \frac{27 t^2(t-1)^2}{4(t^2-t+1)^3}$ with dessin

But then, you may argue, this extender sends all of $0,1$ and $\infty$ to $0$, so it cannot belong to $\mathcal{E}_{MMM}$.

Here’s a trick to construct Belyi-extenders from Belyi-maps $\beta : \mathbb{P}^1 \rightarrow \mathbb{P}^1$, defined over $\mathbb{Q}$ and having the property that there are rational points in the fibers over $0,1$ and $\infty$.

Let’s take an example, the ‘monstrous dessin’ corresponding to the congruence subgroup $\Gamma_0(2)$

with map $\beta(t) = \frac{(t+256)^3}{1728 t^2}$.

As it stands, $\beta$ is not a Belyi-extender because it does not map $1$ into $\{ 0,1,\infty \}$. But we have that
$-256 \in \beta^{-1}(0),~\infty \in \beta^{-1}(\infty),~\text{and}~512,-64 \in \beta^{-1}(1)$
(the last one follows from $(t+256)^2-1728 t^3=(t-512)^2(t+64)$).

We can now pre-compose $\beta$ with the automorphism (defined over $\mathbb{Q}$) sending $0$ to $-256$, $1$ to $-64$ and fixing $\infty$ to get a Belyi-extender
$\gamma(t) = \frac{(192t)^3}{1728(192t-256)^2}$
which maps $\gamma(0)=0,~\gamma(1)=1$ and $\gamma(\infty)=\infty$ (so belongs to $\mathcal{E}_{MMM}$) with the same dessin, which is not a tree,

That is, $\mathcal{E}_{MMM}$ can at best consist only of those Belyi-extenders $\gamma(t)$ that map $\{ 0,1,\infty \}$ onto $\{ 0,1,\infty \}$ and such that their dessin is a tree.

Let me stop, for now, by asking for a reference (or counterexample) to perhaps the most startling claim in the Manin-Marcolli paper, namely that any 2-coloured tree can be realised as the dessin of a Belyi-extender!

I’m trying to get into the latest Manin-Marcolli paper Quantum Statistical Mechanics of the Absolute Galois Group on how to create from Grothendieck’s dessins d’enfant a quantum system, generalising the Bost-Connes system to the non-Abelian part of the absolute Galois group $Gal(\overline{\mathbb{Q}}/\mathbb{Q})$.

In doing so they want to extend the action of the multiplicative monoid $\mathbb{N}_{\times}$ by power maps on the roots of unity to the action of a larger monoid on all dessins d’enfants.

Here they use an idea, originally due to Jordan Ellenberg, worked out by Melanie Wood in her paper Belyi-extending maps and the Galois action on dessins d’enfants.

To grasp this, it’s best to remember what dessins have to do with Belyi maps, which are maps defined over $\overline{\mathbb{Q}}$
$\pi : \Sigma \rightarrow \mathbb{P}^1$
from a Riemann surface $\Sigma$ to the complex projective line (aka the 2-sphere), ramified only in $0,1$ and $\infty$. The dessin determining $\pi$ is the 2-coloured graph on the surface $\Sigma$ with as black vertices the pre-images of $0$, white vertices the pre-images of $1$ and these vertices are joined by the lifts of the closed interval $[0,1]$, so the number of edges is equal to the degree $d$ of the map.

Wood considers a very special subclass of these maps, which she calls Belyi-extender maps, of the form
$\gamma : \mathbb{P}^1 \rightarrow \mathbb{P}^1$
defined over $\mathbb{Q}$ with the additional property that $\gamma$ maps $\{ 0,1,\infty \}$ into $\{ 0,1,\infty \}$.

The upshot being that post-compositions of Belyi’s with Belyi-extenders $\gamma \circ \pi$ are again Belyi maps, and if two Belyi’s $\pi$ and $\pi’$ lie in the same Galois orbit, then so must all $\gamma \circ \pi$ and $\gamma \circ \pi’$.

The crucial Ellenberg-Wood idea is then to construct “new Galois invariants” of dessins by checking existing and easily computable Galois invariants on the dessins of the Belyi’s $\gamma \circ \pi$.

For this we need to know how to draw the dessin of $\gamma \circ \pi$ on $\Sigma$ if we know the dessins of $\pi$ and of the Belyi-extender $\gamma$. Here’s the procedure

Here, the middle dessin is that of the Belyi-extender $\gamma$ (which in this case is the power map $t \rightarrow t^4$) and the upper graph is the unmarked dessin of $\pi$.

One has to replace each of the black-white edges in the dessin of $\pi$ by the dessin of the expander $\gamma$, but one must be very careful in respecting the orientations on the two dessins. In the upper picture just one edge is replaced and one has to do this for all edges in a compatible manner.

Thus, a Belyi-expander $\gamma$ inflates the dessin $\pi$ with factor the degree of $\gamma$. For this reason i prefer to call them dessinflateurs, a contraction of dessin+inflator.

In her paper, Melanie Wood says she can separate dessins for which all known Galois invariants were the same, such as these two dessins,

by inflating them with a suitable Belyi-extender and computing the monodromy group of the inflated dessin.

This monodromy group is the permutation group generated by two elements, the first one gives the permutation on the edges given by walking counter-clockwise around all black vertices, the second by walking around all white vertices.

For example, by labelling the edges of $\Delta$, its monodromy is generated by the permutations $(2,3,5,4)(1,6)(8,10,9)$ and $(1,3,2)(4,7,5,8)(9,10)$ and GAP tells us that the order of this group is $1814400$. For $\Omega$ the generating permutations are $(1,2)(3,6,4,7)(8,9,10)$ and $(1,2,4,3)(5,6)(7,9,8)$, giving an isomorphic group.

Let’s inflate these dessins using the Belyi-extender $\gamma(t) = -\frac{27}{4}(t^3-t^2)$ with corresponding dessin

It took me a couple of attempts before I got the inflated dessins correct (as i knew from Wood that this simple extender would not separate the dessins). Inflated $\Omega$ on top:

Both dessins give a monodromy group of order $35838544379904000000$.

Now we’re ready to do serious work.

Melanie Wood uses in her paper the extender $\zeta(t)=\frac{27 t^2(t-1)^2}{4(t^2-t+1)^3}$ with associated dessin

and says she can now separate the inflated dessins by the order of their monodromy groups. She gets for the inflated $\Delta$ the order $19752284160000$ and for inflated $\Omega$ the order $214066877211724763979841536000000000000$.

It’s very easy to make mistakes in these computations, so probably I did something horribly wrong but I get for both $\Delta$ and $\Omega$ that the order of the monodromy group of the inflated dessin is $214066877211724763979841536000000000000$.

I’d be very happy when someone would be able to spot the error!

Let’s try to identify the $\Psi(n) = n \prod_{p|n}(1+\frac{1}{p})$ points of $\mathbb{P}^1(\mathbb{Z}/n \mathbb{Z})$ with the lattices $L_{M \frac{g}{h}}$ at hyperdistance $n$ from the standard lattice $L_1$ in Conway’s big picture.

Here are all $24=\Psi(12)$ lattices at hyperdistance $12$ from $L_1$ (the boundary lattices):

You can also see the $4 = \Psi(3)$ lattices at hyperdistance $3$ (those connected to $1$ with a red arrow) as well as the intermediate $12 = \Psi(6)$ lattices at hyperdistance $6$.

The vertices of Conway’s Big Picture are the projective classes of integral sublattices of the standard lattice $\mathbb{Z}^2=\mathbb{Z} e_1 \oplus \mathbb{Z} e_2$.

Let’s say our sublattice is generated by the integral vectors $v=(v_1,v_2)$ and $w=(w_1.w_2)$. How do we determine its class $L_{M,\frac{g}{h}}$ where $M \in \mathbb{Q}_+$ is a strictly positive rational number and $0 \leq \frac{g}{h} < 1$?

Here’s an example: the sublattice (the thick dots) is spanned by the vectors $v=(2,1)$ and $w=(1,4)$

Well, we try to find a basechange matrix in $SL_2(\mathbb{Z})$ such that the new 2nd base vector is of the form $(0,z)$. To do this take coprime $(c,d) \in \mathbb{Z}^2$ such that $cv_1+dw_1=0$ and complete with $(a,b)$ satisfying $ad-bc=1$ via Bezout to a matrix in $SL_2(\mathbb{Z})$ such that
$\begin{bmatrix} a & b \\ c & d \end{bmatrix} \begin{bmatrix} v_1 & v_2 \\ w_1 & w_2 \end{bmatrix} = \begin{bmatrix} x & y \\ 0 & z \end{bmatrix}$
then the sublattice is of class $L_{\frac{x}{z},\frac{y}{z}~mod~1}$.

In the example, we have
$\begin{bmatrix} 0 & 1 \\ -1 & 2 \end{bmatrix} \begin{bmatrix} 2 & 1 \\ 1 & 4 \end{bmatrix} = \begin{bmatrix} 1 & 4 \\ 0 & 7 \end{bmatrix}$
so this sublattice is of class $L_{\frac{1}{7},\frac{4}{7}}$.

Starting from a class $L_{M,\frac{g}{h}}$ it is easy to work out its hyperdistance from $L_1$: let $d$ be the smallest natural number making the corresponding matrix integral
$d. \begin{bmatrix} M & \frac{g}{h} \\ 0 & 1 \end{bmatrix} = \begin{bmatrix} u & v \\ 0 & w \end{bmatrix} \in M_2(\mathbb{Z})$
then $L_{M,\frac{g}{h}}$ is at hyperdistance $u . w$ from $L_1$.

Now that we know how to find the lattice class of any sublattice of $\mathbb{Z}^2$, let us assign a class to any point $[c:d]$ of $\mathbb{P}^1(\mathbb{Z}/n\mathbb{Z})$.

As $gcd(c,d)=1$, by Bezout we can find a integral matrix with determinant $1$
$S_{[c:d]} = \begin{bmatrix} a & b \\ c & d \end{bmatrix}$
But then the matrix
$\begin{bmatrix} a.n & b.n \\ c & d \end{bmatrix}$
has determinant $n$.

Working backwards we see that the class $L_{[c:d]}$ of the sublattice of $\mathbb{Z}^2$ spanned by the vectors $(a.n,b.n)$ and $(c,d)$ is of hyperdistance $n$ from $L_1$.

This is how the correspondence between points of $\mathbb{P}^1(\mathbb{Z}/n\mathbb{Z})$ and classes in Conway’s big picture at hyperdistance $n$ from $L_1$ works.

Let’s do an example. Take the point $[7:3] \in \mathbb{P}^1(\mathbb{Z}/12\mathbb{Z})$ (see last time), then
$\begin{bmatrix} -2 & -1 \\ 7 & 3 \end{bmatrix} \in SL_2(\mathbb{Z})$
so we have to determine the class of the sublattice spanned by $(-24,-12)$ and $(7,3)$. As before we have to compute
$\begin{bmatrix} -2 & -7 \\ 7 & 24 \end{bmatrix} \begin{bmatrix} -24 & -12 \\ 7 & 3 \end{bmatrix} = \begin{bmatrix} -1 & 3 \\ 0 & -12 \end{bmatrix}$
giving us that the class $L_{[7:3]} = L_{\frac{1}{12}\frac{3}{4}}$ (remember that the second term must be taken $mod~1$).

If you do this for all points in $\mathbb{P}^1(\mathbb{Z}/12\mathbb{Z})$ (and $\mathbb{P}^1(\mathbb{Z}/6\mathbb{Z})$ and $\mathbb{P}^1(\mathbb{Z}/3 \mathbb{Z})$) you get this version of the picture we started with

You’ll spot that the preimages of a canonical coordinate of $\mathbb{P}^1(\mathbb{Z}/m\mathbb{Z})$ for $m | n$ are the very same coordinate together with ‘new’ canonical coordinates in $\mathbb{P}^1(\mathbb{Z}/n\mathbb{Z})$.

To see that this correspondence is one-to-one and that the index of the congruence subgroup
$\Gamma_0(n) = \{ \begin{bmatrix} p & q \\ r & s \end{bmatrix}~|~n|r~\text{and}~ps-qr=1 \}$
in the full modular group $\Gamma = PSL_2(\mathbb{Z})$ is equal to $\Psi(n)$ it is useful to consider the action of $PGL_2(\mathbb{Q})^+$ on the right on the classes of lattices.

The stabilizer of $L_1$ is the full modular group $\Gamma$ and the stabilizer of any class is a suitable conjugate of $\Gamma$. For example, for the class $L_n$ (that is, of the sublattice spanned by $(n,0)$ and $(0,1)$, which is of hyperdistance $n$ from $L_1$) this stabilizer is
$Stab(L_n) = \{ \begin{bmatrix} a & \frac{b}{n} \\ c.n & d \end{bmatrix}~|~ad-bc = 1 \}$
and a very useful observation is that
$Stab(L_1) \cap Stab(L_n) = \Gamma_0(n)$
This is the way Conway likes us to think about the congruence subgroup $\Gamma_0(n)$: it is the joint stabilizer of the classes $L_1$ and $L_n$ (as well as all classes in the ‘thread’ $L_m$ with $m | n$).

On the other hand, $\Gamma$ acts by rotations on the big picture: it only fixes $L_1$ and maps a class to another one of the same hyperdistance from $L_1$.The index of $\Gamma_0(n)$ in $\Gamma$ is then the number of classes at hyperdistance $n$.

To see that this number is $\Psi(n)$, first check that the classes at hyperdistance $p^k$ for $p$ a prime number and for all $k$ for the $p+1$ free valent tree with root $L_1$, so there are exactly $p^{k-1}(p+1)$ classes as hyperdistance $p^k$.

To get from this that the number of hyperdistance $n$ classes is indeed $\Psi(n) = \prod_{p|n}p^{v_p(n)-1}(p+1)$ we have to use the prime- factorisation of the hyperdistance (see this post).

The fundamental domain for the action of $\Gamma_0(12)$ by Moebius tranfos on the upper half plane must then consist of $48=2 \Psi(12)$ black or white hyperbolic triangles

Next time we’ll see how to deduce the ‘monstrous’ Grothendieck dessin d’enfant for $\Gamma_0(12)$ from it

Dedekind’s Psi-function $\Psi(n)= n \prod_{p |n}(1 + \frac{1}{p})$ pops up in a number of topics:

• $\Psi(n)$ is the index of the congruence subgroup $\Gamma_0(n)$ in the modular group $\Gamma=PSL_2(\mathbb{Z})$,
• $\Psi(n)$ is the number of points in the projective line $\mathbb{P}^1(\mathbb{Z}/n\mathbb{Z})$,
• $\Psi(n)$ is the number of classes of $2$-dimensional lattices $L_{M \frac{g}{h}}$ at hyperdistance $n$ in Conway’s big picture from the standard lattice $L_1$,
• $\Psi(n)$ is the number of admissible maximal commuting sets of operators in the Pauli group of a single qudit.

The first and third interpretation have obvious connections with Monstrous Moonshine.

Conway’s big picture originated from the desire to better understand the Moonshine groups, and Ogg’s Jack Daniels problem
asks for a conceptual interpretation of the fact that the prime numbers such that $\Gamma_0(p)^+$ is a genus zero group are exactly the prime divisors of the order of the Monster simple group.

Here’s a nice talk by Ken Ono : Can’t you just feel the Moonshine?

For this reason it might be worthwhile to make the connection between these two concepts and the number of points of $\mathbb{P}^1(\mathbb{Z}/n\mathbb{Z})$ as explicit as possible.

Surely all of this is classical, but it is nicely summarised in the paper by Tatitscheff, He and McKay “Cusps, congruence groups and monstrous dessins”.

The ‘monstrous dessins’ from their title refers to the fact that the lattices $L_{M \frac{g}{h}}$ at hyperdistance $n$ from $L_1$ are permuted by the action of the modular groups and so determine a Grothendieck’s dessin d’enfant. In this paper they describe the dessins corresponding to the $15$ genus zero congruence subgroups $\Gamma_0(n)$, that is when $n=1,2,3,4,5,6,7,8,9,10,12,13,16,18$ or $25$.

Here’s the ‘monstrous dessin’ for $\Gamma_0(6)$

But, one can compute these dessins for arbitrary $n$, describing the ripples in Conway’s big picture, and try to figure out whether they are consistent with the Riemann hypothesis.

We will get there eventually, but let’s start at an easy pace and try to describe the points of the projective line $\mathbb{P}^1(\mathbb{Z}/n \mathbb{Z})$.

Over a field $k$ the points of $\mathbb{P}^1(k)$ correspond to the lines through the origin in the affine plane $\mathbb{A}^2(k)$ and they can represented by projective coordinates $[a:b]$ which are equivalence classes of couples $(a,b) \in k^2- \{ (0,0) \}$ under scalar multiplication with non-zero elements in $k$, so with points $[a:1]$ for all $a \in k$ together with the point at infinity $[1:0]$. When $n=p$ is a prime number we have $\# \mathbb{P}^1(\mathbb{Z}/p\mathbb{Z}) = p+1$. Here are the $8$ lines through the origin in $\mathbb{A}^2(\mathbb{Z}/7\mathbb{Z})$

Over an arbitrary (commutative) ring $R$ the points of $\mathbb{P}^1(R)$ again represent equivalence classes, this time of pairs
$(a,b) \in R^2~:~aR+bR=R$
with respect to scalar multiplication by units in $R$, that is
$(a,b) \sim (c,d)~\quad~\text{iff}~\qquad \exists \lambda \in R^*~:~a=\lambda c, b = \lambda d$
For $\mathbb{P}^1(\mathbb{Z}/n \mathbb{Z})$ we have to find all pairs of integers $(a,b) \in \mathbb{Z}^2$ with $0 \leq a,b < n$ with $gcd(a,b)=1$ and use Cremona’s trick to test for equivalence:
$(a,b) = (c,d) \in \mathbb{P}^1(\mathbb{Z}/n \mathbb{Z})~\quad \text{iff}~\quad ad-bc \equiv 0~mod~n$
The problem is to find a canonical representative in each class in an efficient way because this is used a huge number of times in working with modular symbols.

Perhaps the best algorithm, for large $n$, is sketched in pages 145-146 of Bill Stein’s Modular forms: a computational approach.

For small $n$ the algorithm in $\S 1.3$ in the Tatitscheff, He and McKay paper suffices:

• Consider the action of $(\mathbb{Z}/n\mathbb{Z})^*$ on $\{ 0,1,…,n-1 \}=\mathbb{Z}/n\mathbb{Z}$ and let $D$ be the set of the smallest elements in each orbit,
• For each $d \in D$ compute the stabilizer subgroup $G_d$ for this action and let $C_d$ be the set of smallest elements in each $G_d$-orbit on the set of all elements in $\mathbb{Z}/n \mathbb{Z}$ coprime with $d$,
• Then $\mathbb{P}^1(\mathbb{Z}/n\mathbb{Z})= \{ [c:d]~|~d \in D, c \in C_d \}$.

Let’s work this out for $n=12$ which will be our running example (the smallest non-squarefree non-primepower):

• $(\mathbb{Z}/12\mathbb{Z})^* = \{ 1,5,7,11 \} \simeq C_2 \times C_2$,
• The orbits on $\{ 0,1,…,11 \}$ are
$\{ 0 \}, \{ 1,5,7,11 \}, \{ 2,10 \}, \{ 3,9 \}, \{ 4,8 \}, \{ 6 \}$
and $D=\{ 0,1,2,3,4,6 \}$,
• $G_0 = C_2 \times C_2$, $G_1 = \{ 1 \}$, $G_2 = \{ 1,7 \}$, $G_3 = \{ 1,5 \}$, $G_4=\{ 1,7 \}$ and $G_6=C_2 \times C_2$,
• $1$ is the only number coprime with $0$, giving us $[1:0]$,
• $\{ 0,1,…,11 \}$ are all coprime with $1$, and we have trivial stabilizer, giving us the points $[0:1],[1:1],…,[11:1]$,
• $\{ 1,3,5,7,9,11 \}$ are coprime with $2$ and under the action of $\{ 1,7 \}$ they split into the orbits
$\{ 1,7 \},~\{ 3,9 \},~\{ 5,11 \}$
giving us the points $[1:2],[3:2]$ and $[5:2]$,
• $\{ 1,2,4,5,7,8,10,11 \}$ are coprime with $3$, the action of $\{ 1,5 \}$ gives us the orbits
$\{ 1,5 \},~\{ 2,10 \},~\{ 4,8 \},~\{ 7,11 \}$
and additional points $[1:3],[2:3],[4:3]$ and $[7:3]$,
• $\{ 1,3,5,7,9,11 \}$ are coprime with $4$ and under the action of $\{ 1,7 \}$ we get orbits
$\{ 1,7 \},~\{ 3,9 \},~\{ 5,11 \}$
and points $[1:4],[3:4]$ and $[5,4]$,
• Finally, $\{ 1,5,7,11 \}$ are the only coprimes with $6$ and they form a single orbit under $C_2 \times C_2$ giving us just one additional point $[1:6]$.

This gives us all $24= \Psi(12)$ points of $\mathbb{P}^1(\mathbb{Z}/12 \mathbb{Z})$ (strangely, op page 43 of the T-H-M paper they use different representants).

One way to see that $\# \mathbb{P}^1(\mathbb{Z}/n \mathbb{Z}) = \Psi(n)$ comes from a consequence of the Chinese Remainder Theorem that for the prime factorization $n = p_1^{e_1} … p_k^{e_k}$ we have
$\mathbb{P}^1(\mathbb{Z}/n \mathbb{Z}) = \mathbb{P}^1(\mathbb{Z}/p_1^{e_1} \mathbb{Z}) \times … \times \mathbb{P}^1(\mathbb{Z}/p_k^{e_k} \mathbb{Z})$
and for a prime power $p^k$ we have canonical representants for $\mathbb{P}^1(\mathbb{Z}/p^k \mathbb{Z})$
$[a:1]~\text{for}~a=0,1,…,p^k-1~\quad \text{and} \quad [1:b]~\text{for}~b=0,p,2p,3p,…,p^k-p$
which shows that $\# \mathbb{P}^1(\mathbb{Z}/p^k \mathbb{Z}) = (p+1)p^{k-1}= \Psi(p^k)$.

Next time, we’ll connect $\mathbb{P}^1(\mathbb{Z}/n \mathbb{Z})$ to Conway’s big picture and the congruence subgroup $\Gamma_0(n)$.