Conserved quantities in the theory of discrete surfaces
Wayne Rossman
June 13, 2007

1

BackgroundSuppose you are given a simple first order smooth ordinary differential equation with a given
initial condition. If you cannot write down its solution explicitly, you might find a discrete
approximate solution by using the Euler or Runga-Kutta algorithm, just to have some initialidea how the smooth solution behaves. In this case, your interest in the approximate solution
is only as a stepping stone for understanding the smooth true solution. We can think of theequation (i.e. the algorithm) for the discrete approximate solution as a finite dimensional
problem because the full space of objects (a vector space of discrete functions) that can beinserted to test for validity in the equation is finite dimensional. Likewise, we can call the smoothdifferential equation an infinite dimensional problem (this might be somewhat unconventional),
because the objects insertable into the equation form an infinite dimensional vector space.
Or you might instead look at a related ordinary difference equation, with little concern thatthe resulting discrete solution approximates the smooth solution, and rather be more concernedthat the difference equation maintains some property found in the smooth differential equationthat you deem important. In this case, as your primary interest is the ”finite dimensional”
difference equation situation itself, you might discard the smooth equation altogether, or you
might acknowledge the existence of the smooth equation but regard it only as an incidental
limiting case of the difference equation you care much more about.
Both approaches are of interest, and are now common in surface theory, though usually involving partial differential equations, not ordinary ones. Discrete analogs of smooth minimal andconstant mean curvature surfaces are being studied. But there is no single definitive way todefine these analogs; the definition one chooses depends on which properties of smooth minimal
and constant mean curvature surfaces one wishes to emulate in the discrete case.
A good example of the first approach is the Costa surface, with the following story: The simplest
example of a minimal surface in R3 is the plane, and two other rather simple examples are: 1)
the catenoid, a surface of revolution produced by revolving a catenary and parametrized by{(cosh u cos v, cosh u sin v, u) ∈ R3 | u ∈ R, v ∈ [0, 2π)} ,Figure 1: The Costa surface.
1

2

RossmanFigure 2: Cut-aways of three constant mean curvature surfaces, a Delaunay unduloid, a Delaunay
nodoid and a Wente torus, in R3 . The first two are surfaces of revolution.

where R denotes the real numbers, and 2) the helicoid, foliated by straight lines and parametrizedby{(sinh u cos v, sinh u sin v, v) ∈ R3 | u, v ∈ R} .
The famous Weierstrass representation says that all minimal surfaces can be locally parametrizedby pairs of meromorphic functions f , g defined on Riemann surfaces Σ with local complex coordinates z, by using path integrals:
Z z√
Re(1 − g 2 , i + ig 2 , 2g)f dz , i = −1 .
z0

The Costa surface, found in 1984 by Costa [5] is a complete minimal surface homeomorphic toa torus minus three points. It has the Weierstrass data
Σ = {(z, w) ∈ (C ∪ {∞})2 | w2 = z(z 2 − 1)} \ {(−1, 0), (1, 0), (∞, ∞)} ,andg = B/w ,f = w/(z 2 − 1) ,

Z 1µ

¶1/2

where B is the constant
s
B=

2
0

t
1 − t2

ÁZ 1
dt
0

dt
.
t(1 − t2 )1/2

There had been a long standing conjecture that the only complete embedded minimal surfaces
with finite topology in R3 are the plane and catenoid and helicoid. In 1985, Hoffman and Meeks[9] confirmed that the Costa surface is a counterexample by proving it is embedded, but only
after numerics led them to see that the surface possessed certain lines and planes of symmetrythat were useful for the proof. Though their final proof used no numerics, the numerics helped
them to find it.
There are also examples of this first approach amongst the constant mean curvature (CMC)
surfaces in R3 . Simple examples of these surfaces are the round sphere and round cylinder.
Less trivial examples are Delaunay surfaces of revolution, parametrizable explicitly in terms ofthe nonconstant periodic Jacobi elliptic function v(x) satisfying(v 0 )2 = −(v 2 − 4s2 )(v 2 − 4t2 ) , v(0) = 2|t| ,with s, t ∈ R \ {0}, s 6= t and s + t = 1/2, and the elliptic integral of the third kind
Z x
4st
dρ .
2
0 4st + v (ρ)
Hopf [11] asked whether any compact CMC surface without boundary in R3 must be a roundsphere. He proved it is true when the surface is simply connected. Alexandrov proved it when
the surface is embedded, using the maximum principle for second-order elliptic differential
equations. However, Wente [17] showed it is false in general, by finding compact nonembedded

Symposium Valenciennes

3

CMC surfaces in R3 without boundary and of genus 1. These tori can be described in terms of
Jacobi elliptic functions and integrals, as the Delaunay surfaces are.
As an example of the first approach, Abresch saw from numerical experimentation that one
family of curvature lines on the Wente tori appeared to be planar curves. He then proceeded tomathematically prove the existence of such CMC tori [1]. After that, Spruck [16] showed thatthe tori found by Abresch are exactly the same collection of surfaces as found by Wente.
To apply the second approach, on the other hand, one must decide what properties of smooth
CMC surfaces one would like to see preserved in the discrete CMC surfaces, and different choicesof those properties result in genuinely different theories.
One choice, following the definitions given above, is to demand that the discrete CMC surfaces
also locally minimize area with respect to variations of the surface that preserve volume toeach side. For this, one would choose the discrete surfaces to be triangulated, i.e. as surfaces
made by gluing triangles together along edges, and then consider variations that continuously
move the vertices while preserving the simplicial structure. If any variation that moves just one
interior vertex and preserves the volume to each side of the surface will never decrease area, we
can say that this discrete surface is of constant mean curvature. See works of K. Polthier [14],[15], and the Surface Evolver by K. Brakke.
But another choice of the property to be preserved, now often used, is as follows: The governing
equation, i.e. the Gauss equation, when using conformal coordinates, for the local existence ofa CMC surface is the sinh-Gordon equation∂z̄ ∂z u + sinh u = 0 ,where z is a local complex coordinate on a Riemann surface. (Let us ignore umbilic points here.
Anyway, these umbilics are isolated on any CMC surface other than the round sphere.)
So we can now obtain discrete CMC surfaces via a discretization of that integrable system (thesinh-Gordon equation). In this case, the area-minimizing property has been lost, so there is no
reason to think about variations of the surface, and so one can consider the surfaces as a meshof planar quadrilaterals (for which one cannot freely move the vertices, as planarity will thenbe lost) rather than of triangles. Perhaps the first breakthrough in this field was by Wunderlich[18], with further developments by the TU-Berlin geometry group.
At first, preserving the area-minimizing property for the discrete surfaces might seem preferable
to preserving a relationship with integrable systems, as the former property is fundamental to
how the smooth surfaces are defined, while the latter property is something that one later
discovers in their mathematical structure. The first way is clearly important, but there aregood reasons for considering the second way too. The second way, by preserving relations tointegrable systems, preserves much of the interesting underlying mathematical structure [2], [8].
(For example, only the second way gives discrete versions of the Bianchi permutability theorem,
and discrete transformations of Backlund, Darboux or Ribaucour type.)
Interestingly, it seems that one cannot simultaneously preserve the area-minimizing property
and the relationships to integrable systems. This can already be seen in the discrete minimal
catenoids coming from each way, as these two types of catenoids really do not coincide.
The integrable systems viewpoint for smooth CMC surfaces itself leads to both approaches
for discretizing. There is a method, called the DPW method after its founders Dorfmeister,Pedit and Wu [7], that is based on integrable systems methods and produces smooth CMC
surfaces (or more generally harmonic maps into symmetric spaces). Central to the method is
an introduction of a spectral parameter λ lying in the unit circle S 1 in the complex plane. In
this method one uses techniques in integrable systems theory, dating back at least to Kričever[12], to construct an object called an extended frame depending on λ, from which a CMC surface
can be constructed. In fact, any CMC surface can be constructed this way. One is implicitly
finding a solution to the sinh-Gordon equation, but the beauty of the method is that the sinhGordon equation itself is essentially bypassed and the surface is constructed without needing
to know anything specific about that solution.

4

RossmanFigure 3: Graphics of a smooth trinoid made with N. Schmitt’s CMCLab software.

Figure 4: Discrete surfaces made from a variational viewpoint. Left: a discrete CMC Delaunay surface.
Middle: a discrete triply-periodic Schwarz P minimal surface. Right: a discrete minimal surface similar
to the smooth Costa minimal surface.
Now one can take either of the two approaches described here for discretizing the DPW method.
The first approach is to find discrete approximations to the desired smooth surface. In thiscase the smooth surface is doubly an infinite dimensional problem, because you have both thesurface parameter and the spectral parameter. The surface parameter can be discretized in theusual way, and the spectral parameter λ can be discretized by chopping away all but a finite
number of terms in the Fourier series for the extended frame. What is left is a finite dimensional
problem, and this is exactly what is solved numerically in the CMCLab program by N. Schmitt.
The second approach is to create a discrete formulation of the DPW method that preserves
integrable systems properties. This approach has been taken by T. Hoffman [10].
By the above avenues, and others as well, computers and discrete methods and integrable
systems methods have come to play a central role in surface theory.
We now describe the DPW recipe for constructing any nonminimal CMC surface in R3 . On
a Riemann surface Σ with local complex coordinates z, we define a holomorphic potential as a
trace-free matrix-valued λ-dependent 1-form, λ ∈ S 1 ,
P∞
µP∞
¶
j
cj (z)λjj=−1 aj (z)λP
ξ = Pj=0
dz ,∞∞j− j=0 cj (z)λj
j=0 bj (z)λwhere the aj dz, bj dz, cj dz are all holomorphic 1-forms defined on Σ, and a−1 is never zero.
Choose an SL2 C-valued solution φ ofdφ = φξ ,analytic in λ, and write φ = F B (this is Iwasawa splitting) so that F is SU2 -valued for all
λ ∈ S 1 and B extends holomorphically to {λ ∈ C | |λ| ≤ 1} and B|λ=0 is upper-triangular.
Although φ is holomorphic in z, F and B are only real-analytic in z. We call F an extended
frame because it is represents a framing for a CMC surface at any λ ∈ S 1 . Then we insert FSymposium Valenciennes

5

Figure 5: Smooth and discrete versions of the triply-periodic Fischer-Koch type minimal surface. The
discrete version is made from a variational viewpoint.

into the Sym-Bobenko formula, choosing λ = 1,£¤
f = −2iH −1 λ∂λ F · F −1 λ=1 ,which is of the form

−i
f=
2

µ

−x3
x1 − ix2

x1 + ix2
x3¶,

for real-valued functions x1 = x1 (z, z̄), x2 = x2 (z, z̄), x3 = x3 (z, z̄), and then one can prove that
Σ 3 z 7→ (x1 , x2 , x3 ) ∈ R3
becomes a conformal parametrization of a CMC H surface.
Among the simple examples of holomorphic potentials, for the sphere, cylinder and Delaunay
surfaces, respectively, areµ¶−1 0 1ξ=λdz , Σ = C ,
0 0
µ
¶
1
0
λ−1 + 1ξ=dz , Σ = C \ {0} ,
0
4 1+λ
¶
µ
0
sλ−1 + t dz
ξ=ξ=, Σ = C \ {0} ,sλ + t0z
for s, t ∈ R \ {0, 1/4} and s + t = 1/2.

2

A conserved quantities approach to smooth CMC surfaces

In the next section we introduce an approach to discrete CMC surfaces coming from joint
work with F. Burstall, U. Hertrich-Jeromin, and S. Santos. But to motivate that discussion,in this section we first explain a result of Burstall and Calderbank [4] for the case of smooth
CMC surfaces. We begin by describing the 3-dimensional space forms using the 5-dimensional
Minkowski space R4,1 .
Minkowski 5-space. We give a 2×2 matrix formulation for Minkowski 5-space. Let H denotethe quaternions and Im H the imaginary quaternions.
½µ¶¯¾
x x∞ ¯¯
4,1
R = X=
x
∈
Im
H,
x
,
x
∈
R
0
∞
x0 −x ¯

6Rossmanwith Minkowski metric hX, Y i such that hX, Y i · I = − 21 (XY + Y X), I = identity matrix. This
metric has signature (+, +, +, +, −) with respect to the (orthonormal) basis
µ

¶ µ
0
j
,
−i
0

i
0

0
−j

¶

µ
,

k
0

0
−k

¶

µ
,

0 1
−1 0

¶

µ
,

0 1
1 0

¶
.

If we set x4 = 12 (x∞ − x0 ), x5 = 21 (x∞ + x0 ), we can write X asµ
X = x1

¶
µ
0
j
+ x2
−i
0

i
0

0
−j

¶

µ
k
+ x3
0

¶
µ
¶
µ
¶
0
0 1
0 1
+ x4
+ x5
,
−k
−1 0
1 0

where x = x1 i + x2 j + x3 k, and then we have the correspondence X ↔ (x1 , x2 , x3 , x4 , x5 ) to
the more usual way{(x1 , x2 , x3 , x4 , x5 ) ∈ R5 | ||(x1 , x2 , x3 , x4 , x5 )|| = x21 + x22 + x23 + x24 − x25 }of denoting R4,1 . The 4-dimensional light cone is
L4 = {X ∈ R4,1 | ||X||2 = 0} .
We can make the 3-dimensional space forms as follows: A space form M is M = L4 ∩
{X | hX, Qi = −I} for any nonzero Q ∈ R4,1 . It will turn out that M has curvature κ, where
Q2 = κ · I, so without loss of generality we can obtain any space form by choosing
µ
¶
0 1
Q=
,
(1)
κ 0
and then

µ
¶¾
2
x −x2
·
,
X=
1 −x
1 − κx2½
M=

which is equivalent to {(x1 , x2 , x3 ) ∈ R3 ∪{∞} |x21 +x22 +x23 6= −κ−1 }, where x = x1 i+x2 j+x3 k ∈Im H. Note that when κ < 0, M becomes two copies of hyperbolic 3-space with sectional
curvature κ.
The tangent space of M at X is½
TX M = Ta =µ¶¾2
a + κxax−xa − ax·,κ(xa + ax) −a − κxax
(1 − κx2 )2

for a ∈ Im H. When X = X(t) ∈ M is a smooth function of a real variable t, and when 0
denotes differentiation with respect to t, we have
X 0 = Tx0 .
A computation gives
hTa , Tb i =−4Re(ab) ,(1 − κx2 )2(2)||Ta || = 1 ⇔ |a| = 21 |1 − κx2 | .
Also,

µ
¶
4(x0 )2
κx −1
X = T 2κ(xx +x x)
·
+
.
κ −κx
(1 − κx2 )2·x0 +x00
1−κx2
00

0

0

Note that generally X 00 is not contained in TX M .
The following lemma follows from (2).
Lemma 1. The M determined by the Q in (1) has constant sectional curvature κ.

(3)

Symposium Valenciennes

7

We see from (2) that the collection of M given by the above choice (1) for Q, for various κ, are
all conformally equivalent (or Moebius equivalent).
Smooth surfaces in space forms. Letx = x1 (u, v)i + x2 (u, v)j + x3 (u, v)k ≈ X ∈ M
be a surface in M . Assume (u, v) is a conformal curvature-line coordinate system (every CMC
surface can be parametrized this way). We call such coordinates isothermic coordinates.
Note that x1 , x2 and x3 can be chosen before the space form M is chosen, and only once M(and hence κ) is chosen do we know the form of X. In particular, the surface can be defined
before the space form is chosen. Because we will always choose Q as in (1), we will indicate
this by denoting M as Mκ , with the subscript κ.
We let n denote the unit normal vector for x, once Mκ is chosen. n0 denotes the unit normal withrespect to Euclidean 3-space M0 , where κ = 0. We sometimes write Hκ for the mean curvatureof the surface x with respect to the space form Mκ , to denote that the mean curvature depends
on the choice of space form. H0 is the mean curvature in the case of Euclidean 3-space M0 .
Lemma 2. The mean curvature H = Hκ of x with respect to the space form Mκ given by Q
as in (1), with 4x = ∂u ∂u x + ∂v ∂v x, is−2
H = −1Re{4x · n} −2 |xu |

κ(xn + nx) =1 − κx22−2
Re{4x · n0 } − κ(xn0 + n0 x) == −1
2 (1 − κx )|xu |(1 − κx2 )H0 − κ(xn0 + n0 x) .
Then H is constant exactly when ∂u H = ∂v H = 0, which is equivalent to2
1(∂u H0 ) · (1 − κx2 ) = κ k2 −k
2 ∂u (x ) ,2
2(∂v H0 ) · (1 − κx2 ) = κ k1 −k
2 ∂v (x ) ,(4)

where the kj ∈ R are the principal curvatures with respect to the Euclidean space form M0 , i.e.
∂u n0 = −k1 ∂u x and ∂v n0 = −k2 ∂v x.
d
Proof. Letting x1u denote du(x1 ), and similarly taking other notations, the unit normal vector
to the surface is Tn , where n = (1 − κx2 )n0 and

n0 =1 (x2u x3v − x3u x2v )i + (x3u x1v − x1u x3v )j + (x1u x2v − x2u x1v )k·p.
2(x2u x3v − x3u x2v )2 + (x3u x1v − x1u x3v )2 + (x1u x2v − x2u x1v )2The first fundamental form (gij ) satisfies hTxu , Txv i = 0 = g12 = g21 , andg11 = hTxu , Txu i =4|xv |2
4|xu |2== hTxv , Txv i = g22 .
(1 − κx2 )2(1 − κx2 )2Then using (3), with the symbol 0 denoting either ∂u or ∂v , we have (where the superscript ”T ”
denotes the part of a vector tangent to TX M )T
b11 = hXuu
, Tn i = hXuu , Tn i =

4κx2u
−4
Re{xuu · n} +(xn + nx) ,2
2(1 − κx )(1 − κx2 )3T
b12 = b21 = hXuv
, Tn i = hXuv , Tn i = 0 ,T
b22 = hXvv
, Tn i = hXvv , Tn i =

−4
4κx2v
Re{x
·
n}
+
(xn + nx) .
vv(1 − κx2 )2(1 − κx2 )3The result follows, using H0 = (k1 + k2 )/2.
We now define the Christoffel transformation x∗ , which for a CMC surface in R3 gives theparallel CMC surface. Let x give a surface in R3 with mean curvature H0 and unit normal n0 .
The Christoffel transformation x∗ satisfies that

8

Rossman• x∗ is defined on the same domain as x,• x∗ has the same conformal structure as x,• x and x∗ have opposite orientations,• and x and x∗ have parallel tangent planes at corresponding points.

This definition above turns out to be equivalent to the following definition, and the existenceof the integrating factor ρ below is equivalent to the existence of isothermic coordinates. Then,−1
once we have x∗ , we will see that we can take x∗ so that dx∗ = x−1
u du − xv dv.
Definition 1. A Christoffel transformation x∗ of x in R3 is such that dx∗ = ρ(dn0 + H0 dx)
for some real-valued function ρ on the surface x (here x∗ is determined only up to translations
and homotheties).
Lemma 3. x∗ exists if and only if x is isothermic.
Proof. We prove only one direction here. Assume x is isothermic, and take isothermic coordinates u, v for x, so xuv = Axu + Bxv for some A, B. Then
−2
−1
d(x−1
u du − xv dv) = 16g11 (xu xuv xu + xv xuv xv )du ∧ dv = 0 .

This implies that there exists an x∗ such that−1dx∗ = x−1
u du − xv dv .

Also,

−1
dn0 + H0 dx = 18 (b11 − b22 )(x−1
u du − xv dv) ,implying that x∗ is a Christoffel transform.
−1
Corollary 1. Christoffel transformations can be taken as solutions of dx∗ = x−1
u du − xv dv.

As a result of the above corollary, we can now simply take the definition of x∗ as follows:
Definition 2. The Christoffel transformation of x is any x∗ (defined in R3 up to translation)−1such that dx∗ = x−1
u du − xv dv.
Lemma 4.
dx∗ =Proof.

µ=

2(dn0 + H0 dx) .
(k1 − k2 )|xu |2

¶
2
−1
−1
(dn0 + H0 dx) − xu du + xv dv |xu |2 =(k1 − k2 )|xu |22
2(−k1 xu du − k2 xv dv + k1 +k
2 (xu du + xv dv)) + xu du − xv dv = 0 .
k1 − k2

The following corollary shows that the Christoffel transform is what it should be when theambient space is R3 , and this surface is CMC, i.e. it gives the parallel CMC surface.
Corollary 2. If H0 is constant for the surface x in R3 , then x∗ is equal to a real constanttimes x + H0−1 n0 .
Proof. Because H0 is constant and we have isothermic coordinates, the Q̂ in the Hopf differentialQ̂(d(u + iv))2 is real. Because Q̂ is holomorphic with respect to u + iv, it is therefore a realconstant. So
Q̂ = hn0 , xuu − xvv i = (k1 − k2 )|xu |2
is constant. We then apply Lemma 4.

Symposium Valenciennes

9

In the next definition, the nonzero real constant c can be chosen freely, and we are once again
considering general space forms M .
µ¶xdx∗ −xdx∗ x
Definition 3. For some c ∈ R \ {0}, we set τ = c. If there exist smooth Q
dx∗−dx∗ xand Z in R4,1 depending on (u, v) such thatd(Q + λZ) = (Q + λZ)λτ − λτ (Q + λZ)(5)holds for all λ ∈ R, then we call Q + λZ a linear conserved quantity of x.
Some properties of linear conserved quantities are immediate. For example, Q and Z 2 areconstant, Xτ = τ X = 0, X ⊥ Z and X ⊥ dZ. We now show these properties:
Lemma 5. Q is constant.
Proof. Set λ = 0 in the conserved quantity equation (5).
Lemma 6. Xτ = τ X = 0.
Proof.
Xτ =2c
1 − κx2

since

µ ¶
x ¡
1
1
¡

µ ¶
¢ x
¡
−x
dx∗ 1
1

¢
−x = 0 ,

µ ¶
¢ x
1 −x
=0.
1

Similarly one can show τ F = 0.
Lemma 7. If Q + λZ is a linear conserved quantity, then Z 2 is constant.
Proof. First note that d(Z 2 ) = Z · dZ + dZ · Z = Z(Qτ − τ Q) + (Qτ − τ Q)Z = (QZ + ZQ)τ −τ (QZ + ZQ), since Zτ = τ Z. Because QZ + ZQ is real, we have d(Z 2 ) = 0.
Lemma 8. X is perpendicular to both Z and dZ.
Proof. XZ + ZX is a real multiple of the identity, and is zero because τ (XZ + ZX) = τ ZX =Zτ X = 0. Thus, X ⊥ Z. Next, X · dZ + dZ · X = X(Qτ − τ Q) + (Qτ − τ Q)X = XQτ − τ QX =(−QX − I)τ − τ (−XQ − I) = −τ + τ = 0. Thus X ⊥ dZ.
Properties like these will be utilized to prove the following theorem:
Theorem 1. [4] The surface x is constant mean curvature in a space form M (produced by
Q 6= 0) if and only if there exists (for that Q) a linear conserved quantity Q + λZ.
Furthermore, when x is not totally umbilic, then Z is unique with positive norm. In fact, in thefollowing proof we can see that Z is uniquely determined by the mean curvature and normal
vector, so Z represents the central sphere congruence.
Proof. In the case that x is part of a sphere, then dx∗ = 0, so τ = 0, and so we can take Z = 0.
So we now assume x is not totally umbilic.
Assume that x has a linear conserved quantity. We can take Q as in (1), and denote thecomponents of Z by
µ
¶
z z∞
Z=
.
z0 −z
The above lemmas tell us that XZ + ZX = 0 and X ⊥ dZ, which, respectively, imply
xz − x2 z0 + zx + z∞ = 0 and x dz − x2 dz0 + dz x + dz∞ = 0 .

10RossmanSymposium Valenciennes

11

Differentiating the first of these two equations, and then applying the second one, we havedx · z − (xdx + dx · x)z0 + zdx = 0 ,which impliesz = z0 · x + h · n0
for some real-valued function h. Thenx(z0 x + hn0 ) − x2 z0 + (z0 x + hn0 )x + z∞ = hxn0 + z0 x2 + hn0 x + z∞ = 0 ,soz∞ = −h(xn0 + n0 x) − z0 x2 .
Thus
Z = z0

µ
¶
µ
x −x2
n0
+h
1 −x
0

−n0 x − xn0
−n0

¶
.

Because Z 2 is constant,(z0 x + hn0 )2 − z0 h(xn0 + n0 x) − z02 x2 = −h2
is constant, and so h is constant, and then also |Z| is constant and nonnegative. A directcomputation, using n0 dx∗ + dx∗ n0 = 0, now shows that Zτ = τ Z, so the condition Zτ = τ Z
coming from Equation (5) provides no extra information. The relation dZ = Qτ − τ Q from (5)gives thatdz0 = cκ(x · dx∗ + dx∗ · x) and dz0 · x + z0 dx + hdn0 = c(dx∗ + κxdx∗ x) .
These two equations tell us thath=2c(1 − κx2 ),
x2u (k2 − k1 )(6)which we know to be constant, and thatz0 = 21 h(k2 + k1 ) = h · H0 .

(7)

Equations (6) and (7) tell us that (4) holds, and so Hκ is constant. One direction of the theorem
now follows.
To prove the other direction, assume that x is a CMC surface with isothermic coordinatez = u + iv, then the Hopf differential is a constant multiple of dz 2 , sob11 − b22 =

4x2u (k2 − k1 )1 − κx2

is constant, and soh=2c(1 − κx2 ), c ∈ R \ {0} ,x2u (k2 − k1 )is also constant. Take Q as in (1), and then take
µ
¶
µ
x −x2
n0
Z = hH0
+h
1 −x
0

¶
−xn0 − n0 x.
−n0Then set the candidate for the conserved quantity to be P = Q + λZ, where dx∗ = x−1
u du −
x−1
dv,
and
τ
is
as
in
Definition
3.
Then
a
computation
gives
dP
+
λτ
P
−
P
λτ
=
0,
by
Equationv(4).
We now explain the conserved quantity equation in terms of the Calapso transformation, inorder to motivate a definition used in the discrete setting.

12

RossmanDefinition 4. Let x be a surface with isothermic coordinates. A Calapso transformation T ∈M ob(3) is a solution ofT −1 dT = λτ .
Then the transformation x → T x is a Calapso transformation. (We can also call it a Ttransform or ”conformal deformation”.)
The Calapso transformation is classical, and was studied by Calapso, Bianchi and Cartan. It
preserves the conformal structure and is thus of interest in Moebius geometry. However, in thecase that the starting surface is CMC, it is the same as the Lawson correspondence, which isquite important in the differential geometry of CMC surfaces.
Lemma 9. If x is isothermic, then the Calapso transform exists.
Proof. The systemµ ¶¢x −1 ¡T Tu = λU , U =xu 1 −x ,1µ ¶¢x −1 ¡T −1 Tv = λV , V = −xv 1 −x1−1

has a solution if and only if U V − V U + Vu − Uv = 0, and this equation holds precisely becauseof the conditions for isothermicness, that isx2u = x2v , xu xv + xv xu = 0 , xuv = Axu + Bxv
for some functions A, B.
Now if we set P = Q + λZ, thendP + λτ P − P λτ = 0
if and only if dP + T −1 dT · P − P · T −1 dT = 0 if and only ifd(T P T −1 ) = 0
if and only if T P T −1 is constant. It is this last condition of T P T −1 being constant that we will
use to define discrete CMC surfaces, just as it defines smooth CMC surfaces, by Theorem 1.
Darboux transformations. For smooth surfaces, a Darboux transform is one such that• there exists an envelope of spheres between the original surface and the transform,• the envelope (i.e. the transform) preserves curvature lines, and• the transform preserves conformality.

3

A conserved quantities approach to discrete CMC surfaces

Our purpose in this section is to present a definition for discrete constant mean curvature(CMC) H surfaces in any of the three space forms Euclidean 3-space R3 , spherical 3-space
S 3 and hyperbolic 3-space H 3 . This new definition is equivalent to the previously known
definitions [2] in the case of R3 . It also satisfies a Calapso transformation relation (the Lawson
correspondence), suggesting the definition is also natural for the space form S 3 , and for CMC
surfaces with H ≥ 1 in H 3 . The definition is the first one for CMC surfaces with −1 < H < 1
in H 3 .
Isothermic discrete surfaces and their Christoffel transforms. Consider a discretesurface fp ∈ ImH ≈ R3 , where p is any point in a discrete lattice domain (locally alwaysSymposium Valenciennes

13

Figure 7: Discrete profile curves for discrete CMC surfaces of revolution. The meanings of thesegraphics are explained in Example 1.
a subdomain of Z 2 ). Consider one quadrilateral in the lattice with vertices p, q, r, s (forexample the points (m, n), (m + 1, n), (m + 1, n + 1), (m, n + 1) for some m, n ∈ Z) orderedcounterclockwise about the quadrilateral. We define the cross ratio of this quadrilateral as
qpqrs = (fq − fp )(fr − fq )−1 (fs − fr )(fp − fs )−1 .
When, for every quadrilateral, we can write the cross ratio as
qpqrs = apq /aps ∈ R
so that the function apq defined on the edges of f satisfies
apq = asr ∈ R and aps = aqr ∈ R ,then we say that f is isothermic. Note that the apq are symmetric, i.e. apq = aqp for any
adjacent p and q.
Remark 1. When considering a smooth surface x(u, v) in R3 and defining
Q = (x(u + ², v − ²) − x(u − ², v − ²))(x(u + ², v + ²) − x(u + ², v − ²))−1 ×(x(u − ², v + ²) − x(u + ², v + ²))(x(u − ², v − ²) − x(u − ², v + ²))−1 ,Bobenko and Pinkall [3] proved that
Q = −I + O(²)
if and only if x is conformal, and
Q = −I + O(²2 )if and only if x is isothermic. This leads to the following definition for discrete isothermic
surfaces in the narrow sense: f is discrete isothermic if(fq − fp )(fr − fq )−1 (fs − fr )(fp − fs )−1 = −I
for all quadrilaterals. However, with this definition, transformations, such as the Calapso transform, of isothermic surfaces will not remain isothermic. Hence the broader definition given
above has been taken up.
When f is isothermic, we can define the Christoffel transform f∗ of f bydf∗pq dfpq = apq .
We can then prove the following:
Lemma 10. [2] f is isothermic if and only if there exists a discrete surface f∗ satisfying theabove equation for df∗ , and f∗ is isothermic with the same cross ratios as f.

14

RossmanCalapso transformations. Like in the smooth case, we can define Calapso transformationsin the discrete case. For adjacent vertices p, q, we define T by
µ ¶
¡
¢
f
Tq = Tp (1 + λτpq ) , τpq = c p (f∗q − f∗p ) 1 −fq .
1Lemma 11. If f is discrete isothermic, then T exists.
Linear conserved quantities. We can now discretize (5) to obtain(1 + λτpq )(Q + λZ)q = (Q + λZ)p (1 + λτpq ) ,(8)where λ ∈ R and Q, Z ∈ R4,1 are functions on the lattice domain.
We can derive this discretization as follows: We say that f is CMC if if there exists a P = Q+λZ
so that T P T −1 is constant, just like in the smooth case. This is equivalent toTq Pq Tq−1 = Tp Pp Tp−1
for all adjacent vertices p, q, which is equivalent to(1 + λτpq )Pq = Pp (1 + λτpq ) ,which becomes Equation (8) above.
We now come to our goal:
Definition 5. If a linear conserved quantity Q + λZ, Q 6= 0, exists for an isothermic discretesurface f, we say that f is of constant mean curvature in the space form M determined by Q.
Remark 2. One can see, in the case that M = R3 , that the above definition is equivalent
to the definition found by Bobenko and Pinkall [2]: f is CMC if |fp − f∗p |2 is constant, andthen it is the constant H0−2 . Also, the property of being discrete CMC is preserved by Calapso
transformations, so the definition here is the right one for the space form M1 = S 3 , and also
for the space form M−1 = H 3 when the mean curvature H−1 has absolute value at least 1.
Furthermore, we can define the constant mean curvature of f to be the lower left entry of Z,
and the normal vector of f to be the upper left entry of Z. Note that here the mean curvatureand the normal vector will be unique if Z is.
Looking at the coefficients in front of the λk in Equation (8) for k = 0, 1, 2, we immediately
have the following lemma:
Lemma 12. Equation (8) is equivalent to dQ = 0 and dZpq = Qp τpq −τpq Qq and τpq Zq = Zp τpq .
Example 1. In the last figure, we show discrete CMC surfaces of revolution. The first two
curves are profile curves for discrete nonminimal CMC surfaces of revolution in R3 , the firstbeing unduloidal and the second nodoidal. (For each of these two curves, the axis of rotation
producing the surface is a vertical line drawn to the left of the curve, and is not shown in the
figure.) The third picture shows the profile curve for a discrete CMC surface of revolution in
S 3 , where S 3 is stereographically projected to R3 , and the shown circle is a geodesic of S 3 thatis also the axis of the surface – and furthermore, this example has a periodicity that causes
it to close on itself and form a torus. The final three pictures show discrete CMC surfaces ofrevolution in H 3 . The first two, with H > 1 and H = 1 respectively, are shown in the Poincare
model, and the first is unduloidal while the second looks similar to a smooth embedded catenoid
cousin. (For these two curves, the corresponding axis of revolution is the vertical line between
the uppermost and lowermost points of the circle shown, and this circle lies in the boundary
sphere at infinity of H 3 .) The last picture is a minimal surface that lies in both copies of
M−1 = H 3 ∪ H 3 , and the horizontal plane shown here is the virtual boundary at infinity of two
copies of the halfspace model for H 3 . This example was not known before, because the notionof discrete CMC was not defined before in this case.

Symposium Valenciennes

15

Darboux transformations. The Darboux transforms of discrete surfaces have similar enveloping properties to the case of smooth surfaces. In the discrete case, the eight vertices of
two corresponding quadrilaterals (one on the original surface and the corresponding one on the
Darboux transform) all lie in one sphere.
Polynomial conserved quantities. Equation (8) can be extended to define surfaces withpolynomial conserved quantities, as follows:(1 + λτpq )(Q + λP1 + λ2 P2 + ... + λn−1 Pn−1 + λn Z)q == (Q + λP1 + λ2 P2 + ... + λn−1 Pn−1 + λn Z)p (1 + λτpq ) .

(9)

16

Rossman

References[14] U. Pinkall and K. Polthier, Computing
discrete minimal surfaces and their conjugates, J. Exp. Math., 2 (1993), 15–36.
[1] U. Abresch, Constant mean curvature tori
in terms of elliptic functions, J. Reine [15] K. Polthier and W. Rossman, Discrete
Angew. Math., 374 (1987), 169–192.
Constant Mean Curvature Surfaces and[2] A. Bobenko and U. Pinkall, Discretizationof surfaces and integrable systems, Oxford Lecture Ser. Math. Appl., 16, Oxford
Univ. Press (1999), 3–58.
[3] A. Bobenko and U. Pinkall, Discrete
isothermic surfaces, J. reine angew.
Math., 475 (1996), 187–208.

their Index, J. Reine. U. Angew. Math.,
549 (2002), 47–77.
[16] J. Spruck, The elliptic sinh-Gordon equation and the construction of toroidal soap
bubbles, Lect. Notes in Math., 1340,
Springer, Berlin (1988), 275-301.

[17] H. C. Wente, Counterexample to a conjecture of H. Hopf, Pac. J. Math., 121[4] F. Burstall, D. Calderbank, Conformal(1986), 193–243.
submanifold geometry, in preparation.
[5] C. Costa, Example of a complete minimal [18] W. Wunderlich, Zur Differenzengeometrie
der Flachen konstanter negativer Krumimmersion of R3 of genus one and threemung, Osterreich. Akad. Wiss. Math.embedded ends, Bull. Soc. Bras. Mat., 15
Nat. Kl. S.-B. IIa., 160 (1951), 39–77.
(1984), 47–54.
[6] C. Delaunay, Sur la surface de révolution
dont la courbure moyenne est constante,
J. Math. pures et appl. Sér. 1 6 (1841),
309–320.
[7] J. Dorfmeister, F. Pedit and H. Wu,
Weierstrass type representation of harmonic maps into symmetric spaces,
Comm. Anal. Geom., 6 (1998), 633–668.
[8] U. Hertrich-Jeromin, Introduction toMöbius Differential Geometry, Cambridge
University Press, London Math. Society
Lect. Note Series 300 (2003).
[9] D. Hoffman and W. H. Meeks III, A
complete embedded minimal surface withgenus one, three ends and finite total curvature, J. Diff. Geom., 21 (1985), 109–127.
[10] T. Hoffmann, Discrete curves and surfaces, Ph.D. thesis, TU-Berlin (2000).
[11] H. Hopf, Differential geometry in thelarge, Lect. Notes in Math., 1000,
Springer, Berlin (1983).
[12] I. M. Kričever, An analogue of thed’Alembert formula for the equations of aprincipal chiral field and the sine-Gordon
equation, Soviet Math. Dokl., 22 (1980),
79–84.
[13] M. Melko and I. Sterling, Integrable systems, harmonic maps and the classicaltheory of surfaces, Aspects Math., 21,Vieweg (1994), 129–144.