by David Dumas and François Guéritaud

Thurston's embedding was introduced in the 1986 paper *Minimal stretch maps between hyperbolic surfaces*. Thurston showed that for any hyperbolic structure *X* on the surface, the map which takes a class [λ] in PML(S) to the covector d_{X}(log L_{λ}) embeds PML(S) as the boundary of a convex neighborhood of the origin in the cotangent space T_{X}^{*}Teich(S). Here L_{λ} : Teich(S) → **R** is the hyperbolic length function.

For visualization purposes it is natural to consider cases in which the real dimension of the cotangent space of Teichmüller space is as small as possible. However, for closed orientable surfaces, this real dimension is at least six. We therefore consider *non-orientable surfaces* and *surfaces with punctures*.

There are three types of simple closed curves on N_{3}:

- The special one-sided curve fixed by the mapping class group (shown in green)
- Other one-sided curves (shown in blue)
- Two-sided curves (shown in red)

A multi-curve on N_{3} is either the union of the two curves in a dual pair, a two-sided curve and the special one-sided curve, or the union of two or three disjoint one-sided curves. Examples are shown below.

To make a picture of Thurston's embedding, we mark the point in T_{X}^{*}Teich(N_{3})=**R**^{3} corresponding to a simple closed curve with a sphere whose radius is inversely proportional to the hyperbolic length of the corresponding geodesic. The spheres are colored according to the topological type of the associated curve as indicated above.

In the resulting pictures (see below), we see that the two-sided simple curves on N_{3} give a dense subset of a Jordan curve in two-dimensional sphere PML(N_{3}). This *equator* curve is the projective measured lamination space PML(S_{1,1}), which naturally embeds in PML(N_{3}) through the torus picture described above. The one-sided curves in N_{3} dual to these two-sided simple curves accumulate on the equator from one side (i.e. they lie in one "hemisphere"). The other open hemisphere of PML(N_{3})) contains only one simple curve: the special one-sided curve.

**Rotating PML(N**_{3})**Earthquake path of PML(N**_{3}) embeddings

As before, the Teichmüller space has real dimension three, and there are one-sided simple curves, two-sided simple curves, and multi-curves. While PML(N_{1,3}) is again homeomorphic a 2-sphere, it was shown by M. Scharlemann (see *The complex of curves on non-orientable surfaces.* J. London Math. Soc. 25 (1982) 171-184) that in this case the closure of the set of two-sided simple curves is an Apollonian gasket. In contrast, each one-sided curve is an isolated point in the set of simple curves; in fact, there is a single one-sided simple curve in each connected component of the complement of the gasket.

**Rotating PML(N**_{1,3}) in the octahedral Thurston embedding

We use the "pentagonal" hyperbolic structure as the base point X for the Thurston embedding. This is the hyperbolic structure obtained by doubling a regular ideal pentagon across its boundary. We identify the tangent space T_{X}Teich(S_{0,5}) with **R**^{4} using a basis obtained from infinitesimal earthquake deformations along simple closed geodesics. In order to produce a symmetric picture, some care must be taken to choose a basis respecting the symmetry of the hyperbolic structure X. Specifically, we start with the blue and red geodesics shown below:

We therefore choose a basis {*v*_{1},*v*_{2},*v*_{3}, *v*_{4}} for the tangent space in which the blue earthquake vectors lie in the span of {*v*_{1},*v*_{2}} and the red earthquake vectors lie in the span of {*v*_{3},*v*_{4}}. Moreover the basis is chosen so that either collection forms the vertex set of a regular pentagon inscribed in the unit circle of its corresponding plane (considering the vectors *v _{i}* to be orthonormal for the purposes of this geometric description).

After stereographic projection, we proceed as before, marking each of the resulting points with a sphere whose radius is inversely proportional to the hyperbolic length of the corresponding geodesic. Note that the sizes of the spheres are computed *after* stereographic projection: The length of the curve determines the Euclidean radius of the sphere drawn **R**^{3}, rather than a spherical radius in **S**^{3}. (While using spherical radii would seem more natural, the resulting pictures are cluttered and confusing. The few large spheres near infinity dominate the view and make it difficult to see the local structure near the origin.)

Despite the inherent distortion when using stereographic projection to visualize the boundary of this four-dimensional convex set, the intricate structure of PML and of the complex of curves is once again evident. Families of simple closed curves that lie in various four-punctured-sphere subsurfaces appear as a network of thin "filaments". Dehn twist families are seen approaching short curves from many directions.

**Rotating stereographic projection**

For better quality: Download the video file – 51MB, 8000kbit/s h.264 MP4**Rotating stereographic projection II (a closer view)**

For better quality: Download the video file – 78MB, 8000kbit/s h.264 MP4

- Full-size PDF for poster printing – 15MB
- Preview image – 3.2MB, 1728x2592 PNG
- Source image for top panel of the poster – 13MB, 6600x4125 PNG
- Source image for bottom panel of the poster – 10MB, 6600x4125 PNG

The poster is intended for printing on "ARCH D" size paper (24in x 36in or 610mm x 914mm).

- Find explicit generators for the Fuchsian group corresponding to the base hyperbolic structure
- Compute a basis for the tangent space to Teichmüller space (as cocycles for a Fuchsian representation of the surface group into PSL(2,
**R**)) - Generate a list of words in π
_{1}S corresponding to simple closed curves - Compute lengths and length derivatives for these words
- Convert the length derivative data to a list of geometric primitives (centers and radii of spheres)
- Render animation frames and encode video files

Lists of words representing simple closed curves were generated by programs written in Python, by applying mapping class group elements to a set of "seed" words representing the different orbits. The visualizations shown above include 50,000 to 2,000,000 simple closed curves. The hyperbolic lengths and derivatives were then calculated by a finite difference method, resulting in data files representing Thurston's embedding.

A final set of Python programs handled stereographic projection, rotations, computing radii, and generating geometric primitives in the input format expected by the ray tracer POV-Ray. A separate template POV-Ray scene file set lighting, background, and texture parameters for the images and animations.

To create the videos, individual frame images were rendered with POV-Ray and then combined and encoded as h.264 video using FFmpeg. Rendering was the most time-consuming step in the process; each frame took between 5 and 30 minutes to render with a 3.7Ghz Intel Xeon E3-1290 V2 CPU.

**Combien de courbes sur une surface**— article by François Guéritaud in the series Images des mathématiques from CNRS.**Visualizing PML**— NSF Mathematical Sciences Institutes research highlight associated with a 2016 summer workshop at ICERM.- This
**3D graphics program built with Unity 3D**provides a real-time "first person" view of these point-cloud data sets in the 3-sphere. It was created by undergraduate researchers in Fall 2015 and Summer 2016 projects in UIC's Mathematical Computing Laboratory.

**Visualizing PML**— lecture series by David Dumas from the ICERM workshop Illustrating Mathematics (June 30 and July 1, 2016).**Python for mathematical visualization: a four-dimensional case study**— Talk at PyCon 2017 by David Dumas about the use of the Python programming language in this project (May 20, 2017).

- The
**pmls05-demo**repository on github has code samples from talks about the project and full code to generate five-punctured sphere images and animations. - An older source release (Nov 2012) has code for generating the length derivative data for the
five-punctured sphere, but lacks some documentation and corrections present in the github repository.
- normball-0.0.1.tar.gz (108KB)

- (
*Planned: release source code for computing simple curves and length derivatives on the other surfaces, including N*.)_{1,3}and N_{3} **ddencode**is a simple front-end to FFmpeg to simplify the process of encoding a sequence of image files to a high-quality h.264 MP4 video file. It was used to generate the videos linked above (and the source videos for uploading to YouTube).

- 2017-05-21 – PyCon video and slides
- 2017-05-17 – CNRS and math institutes article links; corrections to description of methods
- 2017-05-16 – Lecture video and github links
- 2016-06-23 – Added PML(S
_{0,5}) poster - 2015-09-30 – Major update to include N
_{3}and N_{1,3}visualizations, collaboration with F. Guéritaud - 2012-11-29 – PML(S
_{0,5}) Source code released. - 2012-11-19 – Re-encoded videos to address playback problems in Quicktime Player for Mac OS X; minor text edits.
- 2012-11-09 – Added second rotation stereographic projection video.
- 2012-11-08 – Initial release: Rotating stereographic projection video.

This material is based upon work supported by the National Science Foundation. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author and do not necessarily reflect the views of the National Science Foundation.