Schläfli – for lush, voluminous polyhedra

Because of some turbulence in real life (™) over the last few months, it’s been somewhat hard to concentrate on research. Fortunately the obligation of teaching sometimes exerts the right sort of psychological pressure to keep my mind on mathematics in short bursts. Concerning teaching I believe it was Bott who said (roughly): the trouble with teaching is that when you’ve done it you feel like you’ve accomplished something. But I think this is exactly wrong, and especially absurd coming from a gifted teacher like Bott.

This quarter I’m teaching an introductory graduate class on Kleinian groups. It’s something I could teach standing on my head, and during a couple of the classes I half suspected that I was. But every time I teach something, no matter how “elementary” or “familiar”, I find that I get something new out of it. This time around I have been thinking about the Schläfli formula for the variation of volume in a smooth family of hyperbolic polyhedra, and the way in which it relates to some other well-known and important volume formulae relating to hyperbolic manifolds and geometry, especially in 3 dimensions. It turns out that there are some elegant and easy ways to derive many otherwise quite complicated statements directly from Schläfli; probably this is well known to experts, but it wasn’t to me, and I think it might make an interesting blog post.

0. Volumes.

The subject of volumes is one where mathematics makes contact with daily life in many different ways. For instance, taking a shower today, I was confronted with this:


No doubt other readers with lush, voluminous hair like mine have had a similar experience.

In two dimensions volume is area. For a hyperbolic polyhedron P there is a beautiful relation between the area of P and the angles at the vertices. This relation (due to Gauss) is cleanest to state when P is a triangle, with angles \alpha, \beta, \gamma. For an “ordinary” triangle, the angles should all be strictly positive numbers of course. But in hyperbolic space it makes sense to consider triangles with some (or all) vertices at infinity — such vertices are called ideal, and a triangle with ideal vertices is called semi-ideal (if there are also some ordinary vertices) or ideal (if all three vertices are ideal). The group of isometries of hyperbolic space acts transitively on the set of distinct triples of points at infinity, leading to the remarkable conclusion that all ideal triangles are isometric. In particular, they all have the same area, which turns out to be \pi.

The following diagram shows how to decompose an ideal triangle into four smaller triangles — one with angles \alpha,\beta,\gamma, and the other three semi-ideal with a single non-ideal vertex with angles \pi-\alpha, \pi-\beta and \pi-\gamma respectively.


A semi-ideal triangle with one regular vertex can be moved by an isometry so that the regular vertex is at the origin (in the unit disk model); thus one sees that such triangles are determined up to isometry by the angle at the regular vertex. Let’s use the notation A(\theta) for the area of such a semi-ideal triangle with regular angle \theta. The decomposition above gives the formula

\text{area}(P ) = \pi - A(\pi-\alpha) - A(\pi-\beta) - A(\pi-\gamma)

So to determine \text{area}(P ) we just need to understand the function A.

The next diagram shows how to decompose two semi-ideal triangles with angles \lambda and \mu into a semi-ideal triangle with angle \mu+\lambda, and an ideal triangle.



A(\lambda) + A(\mu) = \pi + A(\lambda+\mu)

This gives a functional equation for A. Since A(\pi)=0 we may inductively compute A(\theta)=\pi-\theta for all \theta of the form \pi/2^k and then for all \theta of the form \pi n/2^k for integers n and k. By continuity we obtain A(\theta)=\pi-\theta for all \theta, and therefore we deduce the angle defect formula:

\text{area}(P ) = \pi - \alpha - \beta -\gamma

for a hyperbolic triangle P with angles \alpha, \beta, \gamma as above.

1. The Schläfli formula.

The Schläfli formula is a variational formula that applies to the volumes of a 1-parameter family of geodesic polyhedra P(t). There is a version of the formula that holds in any dimension. We suppose that all the polyhedra P(t) are combinatorially equivalent to some fixed P. For each codimension 2 face e of P, let \ell_e(t) denote the (n-2)-dimensional “volume” of the corresponding face e(t) of P(t). For example, if P is 3-dimensional, e(t) is an edge, and \ell_e(t) is its length. Let \theta_e(t) denote the dihedral angle along the face e(t). Lastly, denote the n-dimensional volume of P(t) by V(t). Then the formula is:

{dV(t)}/{dt} = -{1}/{(n-1)} \sum_e \ell_e(t){d\theta_e(t)}/{dt}

 If we agree that the “0-dimensional volume” of a point is 1, then the angle defect formula follows by integration, using the fact that an infinitesimally small hyperbolic triangle looks nearly Euclidean (and therefore has area close to 0 and angle sum close to \pi).

There is a short and slick proof of this formula due to Milnor, contained in the first volume of his collected works; but I am not sure how much insight it really gives. Anyway, here is the argument.

Proof of Schläfli formula: First, since both sides are additive under decomposition into pieces, it suffices to prove the formula for a simplex. Second, since both sides are linear in first derivatives, it suffices to prove the formula for a set of variations which linearly span the space of deformations of a simplex. An n-simplex is cut out by (n+1) totally geodesic hyperplanes, and the variations of each hyperplane are spanned by parabolic motions perpendicular to the hyperplane based at a point at infinity in the plane.

We now choose simple coordinates; let’s concentrate on the case n=3 for concreteness. In the upper half-space model, arrange so that our simplex P has four vertices, three of which are contained in a vertical plane \pi with constant x coordinate. Let \Delta be the face of P contained in \pi. Cyclically label the oriented edges of \Delta as e_1,e_2,e_3 so that e_1 is contained in the intersection of \pi with the unit hemisphere centered at the origin in the x-y plane. We deform P by adjusting the x-coordinate on \pi by dx. We compute

dV/dx = \int_\Delta dydz/z^3 = 1/2 \int_{\partial \Delta} dy/z^2

where the second equality is Stokes’ theorem.

If \theta denotes the dihedral angle along e_1, then x=\cos(\theta) and z_0=\sin(\theta) where z_0 is the maximum of the z coordinate on the geodesic containing e_1. Parameterize e_1 by angle \phi so that y=z_0\cos(\phi) and z=z_0\sin(\phi) along e_1, and then we obtain

\ell_{e_1} = \int_{e_1} dy/z\sin(\phi) = \int_{e_1} dy z_0/z^2 = -dx/d\theta \int_{e_1} dy/z^2

and similarly for e_2 and e_3. Putting this together the formula follows. qed.

2. Schläfli from Crofton via Hodgson.

When I recently corresponded with Martin Bridgeman, I opined that Milnor’s derivation did not seem to really “explain” the formula, and I wondered aloud whether there was a more insightful derivation. Martin helpfully pointed me to Craig Hodgson’s thesis (which I read a very long time ago) where there is a derivation using some version of the Crofton formula.

For those who don’t know, the Crofton formula is one of the most beautiful and useful in all of geometry. The simplest version expresses the length of a (rectifiable) finite plane curve \gamma in terms of the “average” number of times it intersects a random straight line. This should be interpreted in the sense that there is a natural isometry-invariant (infinite) measure on the space of straight lines, and intersection number with a finite curve \gamma defines an integrable function on this space.

If P is a compact hyperbolic polyhedron in n-dimensional hyperbolic space, there is a natural isometry-invariant measure on the space of hyperplanes, and if \pi is a hyperplane, the (n-1) dimensional area of P \cap \pi defines an integrable function with respect to this measure. The integral, up to a constant, is the volume of P. Using this formula, one sees that Schläfli for n-dimensional polyhedra follows from the same formula for (n-1)-dimensional polyhedra, at least up to some constant (which can be determined by looking at an example, for instance). Turning this around, the angle defect formula inductively proves Schläfli in every dimension (again up to a determination of the constant).

3. Infinitesimal volume rigidity.

I would now like to explain some corollaries and easy derivations of the formula. The first is a weak version of Gromov Proportionality, a key step in modern proofs of the Mostow Rigidity Theorem.

Gromov Proportionality is the statement that in each dimension n>1 there is a positive constant v_n so that if M is a closed, oriented hyperbolic n-manifold, there is an equality

\|[M]\|_1 = \text{vol}(M)/v_n

where [M] denotes the fundamental class of M in n-dimensional homology, and where \|\cdot\|_1 denotes the infimum of the sum \sum |t_i| over all singular n-cycles \sum t_i \sigma_i representing the class [M]. Since this quantity manifestly depends only on the topology of M (in fact, only on its fundamental group!) this equality shows that the volume of any hyperbolic structure on M is a topological invariant. For even dimensions, this fact is much more elementary, since it is a special case of Chern-Gauss-Bonnet, which says that the Euler characteristic of M can be obtained by integrating the Pfaffian of the curvature form; for a manifold of constant curvature, this says that volume and Euler characteristic are proportional (and in even dimensions, neither are zero).

A weaker statement is the observation that if M(t) is a 1-parameter family of hyperbolic metrics on M then the volume is constant along the family. To see why this follows from Schläfli, just cut M(t) up into a family of geodesic polyhedra P_i(t) and apply Schläfli to each polyhedron. For every edge (codimension 2 face), the various polyhedra fitting around it have dihedral angles which sum to 2\pi, so the contributions to variation of volume all cancel.

4. Neumann-Zagier formula.

The next application is a quick derivation of a famous theorem of Neumann-Zagier from this paper (with its famous math review by Jorgenson), for the leading order change in volume under hyperbolic Dehn filling of the cusp of a complete finite volume hyperbolic 3-manifold.

Let M denote a hyperbolic 3-manifold which is complete and finite volume, and has a single torus cusp. Let m,l denote the meridian and longitude on the torus. Dehn filling is the operation of gluing in a new solid torus along the boundary of M to obtain a closed manifold. The only relevant parameter is the slope on \partial M to which the meridian of the solid torus is attached; in terms of the coordinates m,l this is a curve pm+ql for coprime integers p,q. One denotes the filled manifold by M_{p/q} and says that it is the result of p/q Dehn filling on M.

At the complete structure, the holonomy representation \rho:\pi_1(M) \to \text{PSL}(2,\mathbb{C}) takes m and l to nontrivial parabolic elements. At a nearby representation \rho' the eigenvalues of \rho'(m) and \rho'(l) become real and distinct; call these eigenvalues \mu^\pm and \lambda^\pm respectively where \mu,\lambda are eigenvalues for the same eigenvector. A deformation gives rise to an incomplete hyperbolic structure which can be completed by adding a geodesic to give a structure on M_{p/q}, if and only if there is a formula

p\log(\mu) + q\log(\lambda) = 2\pi i

where we take the branch of the logarithm which is zero at the complete structure. Thurston showed (by an elementary computation) that for small deformations of the complete structure, the ratio \log{\lambda}/\log{\mu} \sim c where c:=a+bi is the “shape” of the Euclidean structure on the torus at the complete structure. That is, up to conjugacy we may suppose that \rho(m):z \to z+1 and \rho(l):z \to z+c, where we are writing isometries as (complex) fractional linear transformations in the usual way. If instead we have

p\log(\mu)+q\log(\lambda) = t2\pi i

then we obtain (by taking metric completion) the structure of a cone manifold on M_{p/q} with cone angle t2\pi along the core geodesic.

Where does Schläfli come in? Let’s let M(t) denote the 1-parameter family of cone manifolds with angle t2\pi as above for t \in (0,1), interpolating between M and M_{p/q}. Decompose M(t) into polyhedra in such a way that the cone geodesic becomes an edge. From the formulae above and using the approximation \log{\lambda} \sim c\log{\mu} we obtain formulae

\log(\mu(t)) \sim t2\pi i(p+qa-qbi)/((p+qa)^2+(qb)^2)

\log(\lambda(t)) \sim t2\pi i(p+qa-qbi)(a+bi)/((p+qa)^2+(qb)^2)

Thus the length of the core geodesic \ell(t) is the greatest common “divisor” of the real parts of these two quantities, which is approximately

\ell(t) \sim t2\pi b/((p+qa)^2+(qb)^2)

Using Schläfli and integrating immediately gives the following

Theorem (Neumann-Zagier): with notation as above there is an estimate

\text{vol}(M) -\text{vol}(M_{p.q}) = \pi^2 b/((p+qa)^2+(qb)^2) + O(p^{-4}+q^{-4})

The error term comes from the fact that the volume is even in p,q since p/q = (-p)/(-q). The quadratic form Q(p,q):=((p+qa)^2+(qb)^2)/b has an intrinsic definition as the length squared of the curve pm+ql on the cusp torus, divided by the area of the torus.

5. Bloch-Wigner dilogarithm.

One last application is to give an integral formula for the volume of an ideal simplex. An ideal simplex has 4 vertices, and by an isometry we can put these vertices (in the upper half-space model) at 0,1,\infty,z for some complex number z called the simplex parameter. Different (orientation-preserving) orderings of the vertices replace z by 1/(1-z) or (z-1)/z. We can therefore define a function D(z) to be the (oriented) volume of the ideal simplex with parameter z.

This function satisfies D(z)=D(1/(1-z))=D((z-1)/z) by what we just said, and D(z)=-D(1-z)=-D(z^{-1}), by thinking about orientations.

Five distinct points 0,1,\infty,z,w span five different ideal simplices, and with the natural orientation, their algebraic volumes sum to zero. Thus there is a 5-term relation


If you have a book of special functions handy, this will give a strong clue as to the identity of the function D(z). It turns out to be equal to the so-called Bloch-Wigner dilogarithm, defined by the integral formula

D(z):=\text{arg}(1-z)\log{|z|}-\text{Im}\left(\int_0^z \log(1-z)dz/z\right)

This is not at all easy to derive directly from the definition. But it falls out effortlessly from Schläfli, by the following trick.

The simplex \Delta is non-compact, but we can truncate it by cutting off four neighborhoods of the vertices, given by their intersections with suitable horoballs. Fix some big real constant T. Let H_\infty be the horoball centered at infinity with boundary the Euclidean plane at height T, and let H_0,H_1,H_z be horoballs centered at 0,1,z with Euclidean height 1/T. The distances between horoballs on the edges \infty 0,\infty 1, \infty z, 01 are all 2\log{T} and the distances on the edges 0z, z1 are (respectively) 2\log{T}+2\log{|z|}, 2\log{T}+2\log{|1-z|}. If T is very big, cutting off the horoballs doesn’t change the volume very much, and we can estimate the variation of volume as a function of z by looking only at the contribution to Schläfli from these six edges. This estimate gets better and better as T \to \infty. Now, the great thing about an ideal simplex is that the sum of its dihedral angles always adds up to 2\pi. This means that to calculate the contribution to Schläfli, we may subtract the same constant 2\log{T} from the length of each of the six edges. But now the dependence on T goes away altogether, and we obtain

dD(z)/dz = \log{|z|}d\text{arg}(1-z)/dz -\log{|1-z|}d\text{arg}(z)/dz

(this uses the easy calculation that the dihedral angles along the edges 0z and z1 are \text{arg}(1/(1-z)) and \text{arg}(z) respectively).

Integrate this expression from 0 to z, and use the fact that D(0)=0. Then integrate the first term by parts to get


The two terms under the integral together sum to \text{Im}(\log(1-z)d(\log z)/dz) and we are done.

This entry was posted in 3-manifolds, Hyperbolic geometry, Special functions and tagged , , , , , . Bookmark the permalink.

6 Responses to Schläfli – for lush, voluminous polyhedra

    • Thanks Ian! This sort of stuff is increasingly hard to find these days, when everyone has chucked out the contents of their filing cabinets (and usually the filing cabinets too).

      • ianagol says:

        I discussed Schläfli’s formula during my seminar this semester, so I scanned Hodgson’s thesis (don’t tell Craig, since I didn’t ask his permission ;). Hodgson does it a bit differently, looking at the intersection with 2-planes instead of hyperplanes (so this doesn’t require induction, but does require a careful analysis of error terms).

        There is also a curious non-differential formula for volume one may glean from this approach. The area of a 2-dimensional cross-section is ∑(π-\theta) – 2π, so one gets a formula which is C ∑(π-\theta_e)l_e -2π Vol(P*), where P* is the dual polyhedron of all 2-planes that meet P, and C is some constant.

  1. (this is a reply to ianagol’s second comment; wordpress doesn’t allow replies nested too deeply)

    OK, I’ll keep this secret between you, me and the internet. :)

    I like the dual formula; probably the analog is easier to interpret in spherical geometry.

  2. melissa says:

    Danny, this is beautiful. Thanks for sharing.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s