---
abstract: |
    We study multidimensional quadrilateral lattices satisfying simultaneously two integrable constraints: a quadratic constraint and the projective Moutard constraint. When the lattice is two dimensional and the quadric under consideration is the Möbius sphere one obtains, after the stereographic projection, the discrete isothermic surfaces defined by Bobenko and Pinkall by an algebraic constraint imposed on the (complex) cross-ratio of the circular lattice. We derive the analogous condition for our generalized isthermic lattices using Steiner’s projective structure of conics and we present basic geometric constructions which encode integrability of the lattice. In particular we introduce the Darboux transformation of the generalized isothermic lattice and we derive the corresponding Bianchi permutability principle. Finally, we study two dimensional generalized isothermic lattices, in particular geometry of their initial boundary value problem.\
    [*Keywords:*]{} discrete geometry; integrable systems; multidimensional quadrilateral lattices; isothermic surfaces; Darboux transformation\
    \
address: 'Wydzia[ł]{} Matematyki i Informatyki, Uniwersytet Warmińsko-Mazurski w Olsztynie, ul. Żo[ł]{}nierska 14, 10-561 Olsztyn, Poland'
author:
- 'Adam Doliwa$^\ddagger$'
title: Generalized Isothermic Lattices
---

[^1]

Introduction
============

Isothermic surfaces
-------------------

In the year 1837 Gabriel Lamé presented [@Lame-isoth] results of his studies on distribution of temperature in a homogeneous solid body in thermal equilibrium. He was interested, in particular, in description of the *isothermic surfaces*, i.e. surfaces of constant temperature within the body; notice that his definition makes sense only for families of surfaces, and not for a single surface. Then he found a condition under which one parameter family of surfaces in (a subset of) ${{\mathbb E}}^3$ consists of isothermic surfaces, and showed (for details see [@Lame] or [@DarbouxOS]) that the three families of confocal quadrics, which provide elliptic coordinates in ${{\mathbb E}}^3$, meet that criterion. Subsequently, he proposed to determine all triply orthogonal systems composed by three isothermic families (triply isothermic systems). Such a program was fulfiled by Gaston Darboux [@DarbouxOS] (see also [@Eisenhart-SV]).

Another path of research was initiated by Joseph Bertrand [@Bertrand] who showed that the surfaces of triply isothermic systems are divided by their lines of curvature into “infinitesimal squares”, or in exact terms, they allow for conformal curvature parametrization. This definition of isothermic surfaces (or surfaces of isothermic curvature lines), which can be applied to a single surface, was commonly accepted in the second half of the XIX-th century (see [@Bianchi; @DarbouxIV]). We mention that the minimal surfaces and the constant mean curvature surfaces are particular examples of the isothermic surfaces. The theory of isothermic surfaces was one of the most favorite subjects of study among prominent geometers of that period. Such surfaces exhibit particular properties, for example there exists a transformation, described by Gaston Darboux in [@Darboux-transf], which produces from a given isothermic surface a family of new surfaces of the same type.

The Gauss-Mainardi-Codazzi equations for isothermic surfaces constitute a nonlinear system generalizing the $\mathrm{sinh}$-Gordon equation (the latter governs the constant mean curvature surfaces), and the Darboux transformation can be interpreted as Bäcklund-type transformation of the system. Soon after that Luigi Bianchi showed [@Bianchi-isoth] that two Darboux transforms of a given isothermic surface determine *in algebraic terms* new isothermic surface being their simultaneous Darboux transform. The Bianchi permutability principle can be considered as a hallmark of integrability (in the sense of soliton theory) of the above-mentioned system. Indeed, the isothermic surfaces were reinterpreted by Cieśliński, Goldstein and Sym [@CGS] within the theory of soliton surfaces [@Sym]. More information on isothermic surfaces and their history the Reader can find in the paper of Klimczewski, Nieszporski and Sym [@KNS], where also a more detailed description of the relation between the “ancient” differential geometry and the soliton theory is given, and in books by Rogers and Schief [@RogersSchief] and by Hertrich-Jeromin [@UHJ].

Discrete isothermic surfaces and discrete integrable geometry
-------------------------------------------------------------

In the recent studies of the relation between geometry and the integrable systems theory a particular attention is payed to discrete (difference) integrable equations and the corresponding discrete surfaces or lattice submanifolds. Also here the discrete analogs of isothermic surfaces played a prominent role in the development of the subject. Bobenko and Pinkall [@BP2] introduced the integrable discrete analogoue of isothermic surfaces as mappings built of “conformal squares”, i.e., maps ${{\boldsymbol x}}:{{\mathbb Z}}^2\to{{\mathbb E}}^3$ with all elementary quadrilaterals circular, and such that the complex cross-ratios (with the plane of a quadrilateral identified with the complex plane ${{\mathbb C}}$) $$q(m,n) = {\mathrm{cr}}\big({{\boldsymbol x}}(m,n),{{\boldsymbol x}}(m+1,n+1);{{\boldsymbol x}}(m+1,n),{{\boldsymbol x}}(m,n+1) \big)_{{\mathbb C}}$$ are equal to $-1$. Soon after that it turned out [@BP-V] that it is more convenient to allow for the cross-rations to satisfy the constraint $$\label{eq:isoth-cr}
q(m,n) q(m+1,n+1) = q(m+1,n)q(m,n+1).$$ Then the cross-ratio is a ratio of functions of single variables, which corresponds to allowed reparametrization of the curvature coordinates on isothermic surfaces.

After the pioneering work of Bobenko and Pinkall, which was an important step in building the geometric approach to integrable discrete equations (see also [@DS-AL; @BP1; @DCN] and older results of the difference geometry (*Differenzengeometrie*) summarized in Robert Sauer’s books [@Sauer2; @Sauer]), the discrete isothermic surfaces and their Darboux transformations were studied in a number of papers [@HeHP; @Ciesl; @Schief-C]. Distinguished integrable reductions of isothermic lattices are the discrete constant mean curvature surfaces or the discrete minimal surfaces [@BP-V; @UHJ]. It should be mentioned that the complex cross-ratio condition was extended to circular lattices of dimension three [@BP-V; @Ciesl] placing the Darboux transformations of the discrete isothermic surfaces on equal footing with the lattice itself.

In the present-day approach to the relation between discrete integrable systems and geometry [@DS-EMP; @BobSur] the key role is played by the integrable discrete analogue of conjugate nets – multidimensional lattices of planar quadrilaterals (the quadrilateral lattices) [@MQL]. These are maps $x:{{\mathbb Z}}^N\to{{\mathbb P}}^M$ of $N$-dimensional integer lattice in $M\geq N$-dimensional projective space with all elementary quadrilaterals planar. Integrability of such lattices (for $N>2$) is based on the following elementary geometry fact (see Figure \[fig:TiTjTkx\]).

![The geometric integrability scheme[]{data-label="fig:TiTjTkx"}](TiTjTkx.eps)

\[lem:gen-hex\] Consider points $x_0$, $x_1$, $x_2$ and $x_3$ in general position in ${{\mathbb P}}^M$, $M\geq 3$. On the plane $\langle x_0, x_i, x_j \rangle$, $1\leq i < j \leq 3$ choose a point $x_{ij}$ not on the lines $\langle x_0, x_i \rangle$, $\langle x_0,x_j
\rangle$ and $\langle x_i, x_j \rangle$. Then there exists the unique point $x_{123}$ which belongs simultaneously to the three planes $\langle x_3, x_{13}, x_{23} \rangle$, $\langle x_2, x_{12}, x_{23} \rangle$ and $\langle x_1, x_{12}, x_{13} \rangle$.

Various constraints compatible with the geometric integrability scheme define integrable reductions of the quadrilateral lattice. It turns out that such geometric notion of integrability very often associates integrable reductions of the quadrilateral lattice with classical theorems of incidence geometry. We advocate this point of view in the present paper.

Among basic reductions of the quadrilateral lattice the so called quadratic reductions [@q-red] play a distinguished role. The lattice vertices are then contained in a hyperquadric (or in the intersection of several of them). Such reductions of the quadrilateral lattice can be often associated with various subgeometries of the projective geometry, when the quadric plays the role of the absolute of the geometry (see also corresponding remarks in [@q-red; @W-cong]). In particular, when the hyperquadric is the Möbius hypersphere one obtains, after the stereographic projection, the circular lattices [@Bobenko-O; @CDS; @DMS; @KoSchief2], which are the integrable discrete analogue of submanifolds of ${{\mathbb E}}^M$ in curvature line parametrization. Because of the Möbius invariance of the complex cross-ratio it is also more convenient to consider (discrete) isothermic surfaces in the Möbius sphere (both dimensions of the lattice and of the sphere can be enlarged) keeping the “cross-ratio definition”.

For a person trained in the projective geometry it is more or less natural to generalize the Möbius geometry approach to discrete isothermic surfaces (lattices) in quadrics replacing the Möbius sphere by a quadric, and correspondingly, the complex cross-ratio by the Steiner cross-ratio of four points of a conic being intersection of the quadric by the plane of elementary quadrilateral of the quadrilateral lattice. However the “cross-ratio point of view” doesn’t answer the crucial question about integrability (understood as compatibility of the constraint with the geometric integrablity scheme) of such discrete isothermic surfaces in quadrics. Our general methodological principle in the integrable discrete geometry, applied successfuly earlier, for example in [@DS-sym; @DNS-Bianchi], which we would like to follow here is (i) to isolate basic reductions of the quadrilateral lattice and then (ii) to incorporate other geometric systems into the theory considering them as superpositions of the basic reductions.

In this context we would like to recall another equivalent characterization of the classical isothermic surfaces which can be found in the classical monograph of Darboux [@DarbouxIV vol. 2, p. 267]: [*Les cinq coordonnées pentasphériques d’un point de toute surface isothermique considérées comme fonctions des paramètres $\rho$ et $\rho_1$ des lignes de courbure satisfont à une équation linéaire du second ordre dont les invariants sont égaux. Inversement, si une équation de la forme $$\label{eq:Moutard-diff}
\frac{\partial^2\theta}{\partial \rho \partial \rho_1} = \lambda \theta$$ ou, plus généralement, une équation à invariants égaux, admet cinq solutions particulières $x_1$, $x_2$, …, $x_5$ liées par l’équation $$\sum_{1}^{5} x_i^2 = 0,$$ les quantités $x_i$ sont les coordonnées pantasphériques qui définissent une surface isothermique rapportée à ses lignes de courbure.* ]{}

In literature there are known two (closely related) discrete integrable versions (of Nimmo and Schief [@NiSchief] and of Nieszporski [@Nieszporski-dK]) of the Moutard equation . It turns out that for our purposes it suits the discrete Moutard equation proposed in [@NiSchief]. Its projectively invariant geometric characterization has been discovered [@DGNS] only recently (for geometric meaning of the adjoint Moutard equation of Nieszporski in terms of the so called Koenigs lattice see [@Dol-Koe]). Indeed, it turns out that the generalized isothermic lattices can be obtained by adding to the quadratic constraint the projective Moutard constraint. Finally, the quadratic reduction and the Moutard reduction, when applied simultaneously, give *a posteriori* the cross-ratio condition .

In fact, the direct algebraic discrete counterpart of the above description of the isothermic surfaces, i.e. existence of the light-cone lift which satisfies the (discrete) Moutard equation, appeared first in a preprint by Bobenko and Suris [@BobSur]. However, the pure geometric characterization of the discrete isothermic surfaces was not given there.

Because integrability of the discrete Moutard equation can be seen better when one considers a system of such equations for multidimensional lattices, there was a need to find the projective geometric characterization of the system. The corresponding reduction of the quadrilateral lattice was called in [@BQL], because of its connection with the discrete BKP equation, the B-quadrilateral lattice (BQL). In fact, research in this direction prevented me from publication of the above mentioned generalization of discrete isothermic surfaces, announced however in my talk during the Workshop “Geometry and Integrable Systems” (Berlin, 3-7 June 2005). I suggested also there that the discrete S-isothermic surfaces of Tim Hoffmann [@Hoffmann] (see also [@BobSur]) should be considered as an example of the generalized isothermic lattices where the quadric under consideration is the Lie quadric. The final results of my research on generalized isothermic lattices were presented on the Conference “Symmetries and Integrability of Difference Equations VII” (Melbourne, 10-14 July 2006).

When my paper was almost ready there appeared the preprint of Bobenko and Suris [@BobSur-gen-isoth] where similar ideas were presented in application to the sphere (Möbius, Laguerre and Lie) geometries. I would like also to point out a recent paper by Wallner and Pottman [@WP] devoted, among others, to discrete isothermic surfaces in the Laguerre geometry.

Plan of the paper
-----------------

As it often happens, the logical presentation of results of a research goes in opposite direction to their chronological derivation. In Section \[sec:BQL+QQL\] we collect some geometric results from the theory of the B-quadrilateral lattices (BQLs) and of quadrilateral lattices in quadrics (QQLs). Some new results concerning the relation between (Steiner’s) cross-ratios of vertices of elementary guadrilaterals of elementary hexahedrons of the QQLs are given there as well. Then in Section \[sec:gen-isoth-latt\] we define generalized isothermic lattices and discuss their basic properties. In particular, we give the synthetic-geometry proof of a basic lemma (the half-hexahedron lemma) which immediately gives the cross-ratio characterization of the lattices. We also present some algebraic consequences (some of them known already [@BobSur]) of the system of Moutard equations supplemented by a quadratic constraint. In Section \[sec:Darboux\] we study in more detail the Darboux transformation of the generalized isothermic lattices and the corresponding Bianchi permutability principle. Finally, in Section \[sec:isothermic\] we consider two dimensional generalized isothermic lattices. In two Appendices we recall necessary information concerning the cross-ratio of four points on a conic curve and we perform some auxilliary calculations.

The B-quadrilateral lattices and the quadrilateral lattices in quadrics {#sec:BQL+QQL}
=======================================================================

It turns out that compatibility of both BQLs and QQLs with the geometric integrablity scheme follows from certain classical geometric facts. We start each section, devoted to a particular lattice, from the corresponding geometric statement.

The B-quadrilateral lattice [@BQL]
----------------------------------

\[lem:BKP-hex\] Under hypothesess of Lemma \[lem:gen-hex\], assume that the points $x_0$, $x_{12}$, $x_{13}$, $x_{23}$ are coplanar, then the points $x_1$, $x_2$, $x_3$, and $x_{123}$ are coplanar as well (see Figure \[fig:moutard\]).

![Elementary hexahedron of the B-quadrilateral lattice[]{data-label="fig:moutard"}](moutard)

As it was discussed in [@BQL] the above fact is equivalent to the Möbius theorem (see, for example [@Coxeter]) on mutually inscribed tetrahedra. Another equivalent, but more symmetric, formulation of Lemma \[lem:BKP-hex\] is provided by the Cox theorem (see [@Coxeter]): *Let $\sigma_1$, $\sigma_2$, $\sigma_3$, $\sigma_4$ be four planes of general position through a point $S$. Let $S_{ij}$ be an arbitrary point on the line $\langle \sigma_i,
\sigma_j \rangle$. Let $\sigma_{ijk}$ denote the plane $\langle S_{ij}, S_{ik},
S_{jk} \rangle$. Then the four planes $\sigma_{234}$, $\sigma_{134}$, $\sigma_{124}$, $\sigma_{123}$ all pass through one point $S_{1234}$.*

\[def:BQL\] A quadrilateral lattice $x:{{\mathbb Z}}^N\to{{\mathbb P}}^M$ is called the *B-quadrilateral lattice* if for any triple of different indices $i,j,k$ the points $x$, $x_{(ij)}$, $x_{(jk)}$ and $x_{(ik)}$ are coplanar.

Here and in all the paper, given a fuction $F$ on ${{\mathbb Z}}^N$, we denote its shift in the $i$th direction in a standard manner: $F_{(i)}(n_1,\dots, n_i, \dots , n_N) = F(n_1,\dots, n_i + 1, \dots , n_N)$. One can show that a quadrilateral lattice $x:{{\mathbb Z}}^N\to{{\mathbb P}}^M$ is a B-quadrilateral lattice if and only if it allows for a homogoneous representation ${{\boldsymbol x}}:{{\mathbb Z}}^N\to{{\mathbb R}}^{M+1}_{*}$ satisfying the system of discrete Moutard equations (the discrete BKP linear problem) $$\label{eq:BKP-lin}
{{\boldsymbol x}}_{(ij)} - {{\boldsymbol x}}= f^{ij} ({{\boldsymbol x}}_{(i)} - {{\boldsymbol x}}_{(j)}) , \quad 1\leq i< j\leq N,$$ for suitable functions $f^{ij}:{{\mathbb Z}}^N\to{{\mathbb R}}$.

The compatibility condition of the system implies that the functions $f^{ij}$ can be written in terms of the potential $\tau:{{\mathbb Z}}^N\to{{\mathbb R}}$, $$\label{eq:tau}
f^{ij} = \frac{\tau_{(i)}\tau_{(j)}}{\tau \, \tau_{(ij)}}, \qquad i\ne j ,$$ which satisfies Miwa’s discrete BKP equations [@Miwa] $$\label{eq:BKP-nlin}
\tau\, \tau_{(ijk)} = \tau_{(ij)}\tau_{(k)} - \tau_{(ik)}\tau_{(j)} + 
\tau_{(jk)}\tau_{(i)}, \quad 1\leq i< j < k \leq N.$$

The trapezoidal lattice [@BobSur] is another reduction of the quadrilateral lattice being algebraically described by the discrete Moutard equations . Geometrically, the trapezoidal lattices are characterized by parallelity of diagonals of the elementary quadrilaterals, thus they belong to the affine geometry. Moreover, because the trapezoidal constraint is imposed on the level of elementary quadrilaterals then, from the point of view of the geometric integrability scheme, one has to check its three dimensional consistency. In contrary, the BQL constraint is imposed on the level of elementary hexahedrons, and to prove geometrically its integrability one has to check four dimensional consistency.

The quadrilateral lattices in quadrics
--------------------------------------

\[lem:qred-hex\] Under hypotheses of Lemma \[lem:gen-hex\], assume that the points $x_0$, $x_1$, $x_2$, $x_3$, $x_{12}$, $x_{13}$, $x_{23}$ belong to a quadric $\mathcal{Q}$. Then the point $x_{123}$ belongs to the quadric $\mathcal{Q}$ as well.

The above fact is a consequence of the classical *eight points theorem* (see, for example [@Coxeter]) which says that *seven points in general position determine a unique eighth point, such that every quadric through the seven passes also through the eighth*. In our case the point $x_{123}$ is contained in the three (degenerate) quadrics being pairs of opposite facets of the hexahedron.

Lemma \[lem:BKP-hex\] can be considered as a “reduced” resion of Lemma \[lem:qred-hex\] when the quadric $\mathcal{Q}$ degenerates to a pair of planes.

\[def:QQL\] A quadrilateral lattice $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ contained in a hyperquadric $\mathcal{Q}$ is called the *$\mathcal{Q}$-reduced quadrilateral lattice* (QQL).

Integrability of the QQLs was pointed out in [@q-red], where also the corresponding Darboux-type transformation (called in this context the Ribaucour transformation) was constructed in the vectorial form. When the quadric is irreducible then generically it cuts the planes of the hexahedron along conics.

A quadrilateral lattice $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ in a hyperquadric $\mathcal{Q}$ such that the intersection of the planes of elementary quadrilaterals of the latice with the quadric are irreducible conic curves is called *locally irreducible*.

The following result, which will not be used in the sequel and whose proof can be found in Appendix \[sec:aux-calc\], generalizes the relation between complex cross-ratios of the opposite quadrilaterals of elementary hexahedrons of the circular lattices [@Bobenko-O].

\[prop:qred-cr\] Given locally irreducible quadrilateral lattice $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ in a hyperquadric $\mathcal{Q}$, denote by $$\lambda^{ij} = {\mathrm{cr}}(x_{(i)},x_{(j)};x,x_{(ij)}), \qquad 1\leq i<j\leq N,$$ the cross-ratios (defined with respect to the corresponding conic curves) of the four vertices of the quadrilaterals. Then the cross-ratios are related by the following system of equations $$\label{eq:qred-cross-r}
\lambda^{ij}\lambda^{ij}_{(k)}\lambda^{jk}\lambda^{jk}_{(i)} =
\lambda^{ik}\lambda^{ik}_{(j)}, \quad 1\leq i<j<k\leq N.$$

The system can be considered as the gauge invariant integrable difference equation governing QQLs.

Generalized isothermic lattices {#sec:gen-isoth-latt}
===============================

Because simultaneous application of integrable constraints preserves integrability we know *a priori* that the following reduction of the quadrilateral lattice is integrable.

A B-quadrilateral lattice in a hyperquadric $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ satisfying the local irreducibility condition is called a *generalized isothermic lattice*.

The half hexahedron lemma and its consequences
----------------------------------------------

We start again from a geometric result, which leads to the cross-ratio characterization of the generalized isothermic lattice.

\[lem:half-hex\] Under hypotheses of Lemmas \[lem:BKP-hex\] and \[lem:qred-hex\] and assuming irreducibility of the conics of the intersection of the corresponding planes with the quadric we have $$\label{eq:product-cr}
{\mathrm{cr}}(x_1,x_3;x_0,x_{13}) = {\mathrm{cr}}(x_1,x_2;x_0,x_{12}) \,
{\mathrm{cr}}(x_2,x_3;x_0,x_{23}),$$ where the cross-ratios are defined with respect to the corresponding conics.

Denote (see Figure \[fig:hex-lem\])

![Geometric proof of Lemma \[lem:half-hex\][]{data-label="fig:hex-lem"}](hex-lem)

the plane $\sigma=\langle x_1, x_2, x_3 \rangle$, and represent points of the conic $\mathcal{C}_{ij}={{\mathcal Q}}\cap\langle x_0, x_i, x_j \rangle$, $1\leq i < j \leq 3$ by points of the line $\langle x_i, x_j \rangle$, via the corresponding planar pencil with the base at $x_0$. In this way the projective structure of the conics conicides with that of the corresponding lines.

By $\ell_{ij}$ denote the intersection of the tangent line to $\mathcal{C}_{ij}$ at $x_0$ with $\sigma$. Notice that the points $\ell_{ij}$ belong to the intersection line of the tangent plane to the quadric at $x_0$ with $\sigma$, and all they represent $x_0$ but from the point of view of different conics. Denote by $d_{ij}$ the intersection point of the line $\langle x_0, x_{ij} \rangle$ with the plane $\sigma$. Notice that coplanarity of the points $x_0$, $x_{12}$, $x_{23}$ and $x_{23}$ is equivalent to collinearity of $d_{12}$, $d_{23}$ and $d_{13}$. Moreover, by definition of the cross-ratio on conics, we have $$\label{eq:cr-c-l}
{\mathrm{cr}}(x_i,x_j;x_0,x_{ij})  = {\mathrm{cr}}(x_i,x_j;\ell_{ij},d_{ij}).$$

To find a relation between the cross-ratios consider perspectivity between the lines $\langle x_1, x_2 \rangle$ and $\langle x_1, x_3 \rangle$ with the center $\ell_{23}$. It transforms $x_2$ into $x_3$, $x_1$ into $x_1$, $d_{12}$ into $\tilde{d}_{12}$ (this is just definition of $\tilde{d}_{12}$) and $\ell_{12}$ into $\ell_{13}$, therefore $${\mathrm{cr}}(x_1,x_2;\ell_{12},d_{12}) = 
{\mathrm{cr}}(x_1,x_3;\ell_{13},\tilde{d}_{12}).$$ Smilarly, considering perspectivity between the lines $\langle x_2, x_3 \rangle$ and $\langle x_1, x_3 \rangle$ with the center $\ell_{12}$ we obtain $${\mathrm{cr}}(x_2,x_3;\ell_{23},d_{23}) = 
{\mathrm{cr}}(x_1,x_3;\ell_{13},\tilde{d}_{23}),$$ where again $\tilde{d}_{23}$ is the projection of $d_{23}$. The comparison of Figures \[fig:hex-lem\] and \[fig:multiplication\] gives $$\label{eq:cr-l-pr}
{\mathrm{cr}}(x_1,x_3;\ell_{13},\tilde{d}_{12})
{\mathrm{cr}}(x_1,x_3;\ell_{13},\tilde{d}_{23})=
{\mathrm{cr}}(x_1,x_3;\ell_{13},d_{13}),$$ which because of equations - implies the statement.

For those who do not like synthetic geometry proofs we give the algebraic proof of the above Lemma in Appendix \[sec:aux-calc\].

Equation can be written in a more symmetric form $$\label{eq:product-cr-symm}
{\mathrm{cr}}(x_1,x_2;x_0,x_{12}) \,
{\mathrm{cr}}(x_2,x_3;x_0,x_{23}) \, {\mathrm{cr}}(x_3,x_1;x_0,x_{13})  = 1.$$

The cross-ratio of the four (coplanar) points $x_0$, $x_{12}$, $x_{13}$ and $x_{23}$ can be expressed by the other cross-ratios as $${\mathrm{cr}}(x_0,x_{12};x_{13},x_{23})={\mathrm{cr}}(x_0,x_{2};x_{3},x_{23})
{\mathrm{cr}}(x_1,x_{0};x_{3},x_{13}).$$

Consider the line $\langle d_{13}, d_{23} \rangle$, which is the section of the planar pencil containing lines $\langle x_0, x_{13} \rangle$ and $\langle x_0, x_{23} \rangle$ with the plane $\sigma$. Denote by $\ell$ the intersection point of the line with the line $\langle \ell_{13}, \ell_{23} \rangle$, then $${\mathrm{cr}}(x_0,x_{12};x_{13},x_{23}) = {\mathrm{cr}}(\ell,d_{12};d_{13},d_{23}).$$ After projection from $\ell_{12}$ we have, in notation of Figure \[fig:hex-lem\], $${\mathrm{cr}}(\ell,d_{12};d_{13},d_{23}) = {\mathrm{cr}}(\ell_{13},x_1;d_{13},\tilde{d}_{23})=
{\mathrm{cr}}(1,\infty;\lambda\nu,\nu).$$ Then the standard permutation properties of the cross-ratio give the statement.

\[cor:hex\] Under assumption of Lemma \[lem:half-hex\] the cross-ratios on opposite quadrilaterals of the hexahedron are equal, i.e. $$\begin{aligned}
\nonumber
{\mathrm{cr}}(x_1,x_2;x_0,x_{12}) & = {\mathrm{cr}}(x_{13},x_{23};x_3,x_{123}) , \\
\label{eq:equalities-cr}
{\mathrm{cr}}(x_2,x_3;x_0,x_{23}) & = {\mathrm{cr}}(x_{12},x_{13};x_1,x_{123}) , \\
\nonumber
{\mathrm{cr}}(x_1,x_3;x_0,x_{13}) & = {\mathrm{cr}}(x_{12},x_{23};x_2,x_{123}) . \end{aligned}$$

Equation written for the three quadrilaterals meeting in $x_3$ reads $${\mathrm{cr}}(x_0,x_{12};x_3,x_{1}) \,
{\mathrm{cr}}(x_{13},x_{23};x_3,x_{123}) \, {\mathrm{cr}}(x_{23},x_0;x_3,x_{2})  = 1,$$ which compared with gives, after using elementary properties of the cross-ratio, the first equation of . Others can be obtained similarly.

By symmetry we have also $${\mathrm{cr}}(x_0,x_{12};x_{13},x_{23}) =  {\mathrm{cr}}(x_{1},x_{2};x_{3},x_{123}).$$

Notice that two neighbouring facets of the above hexahedron determine the whole hexahedron via construction visualized on Fig. \[fig:twotothree\].

![Construction of the point $x_{13}$ from points $x_0$, $x_{1}$, $x_{2}$, $x_{3}$, $x_{12}$ and $x_{23}$. It belongs to the intersection line of two planes $\langle x_0, x_{1}, x_{3} \rangle$ with the plane $\langle x_0, x_{12}, x_{23} \rangle$. Because the line intersects the quadric at $x_0$, it must have also the second intersection point.[]{data-label="fig:twotothree"}](twotothree)

It is easy to see that, unlike in the case of isothermic lattices, three vertices of a quadrilateral of trapezoidal lattice in a quadric [@BobSur] determine the forth vertex.

\[prop:cross-ratio-gen\] A quadrilateral lattice in a quadric $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ satisfying the local irreducibility condition is a generalized isothermic lattice if and only if there exist functions $\alpha^{i}:{{\mathbb Z}}\to{{\mathbb R}}$ of single arguments $n_i$ such that the cross-ratios $\lambda^{ij} = {\mathrm{cr}}(x_{(i)},x_{(j)};x,x_{(ij)})$ can be factorized as follows $$\label{eq:isot-cross-res}
\lambda^{ij} = \frac{\alpha^i}{\alpha^j}, \qquad 1\leq i<j\leq N.$$

Equations and can be rewritten as $$\label{eq:izot-cross-r-1}
\lambda^{ij}\lambda^{jk} = \lambda^{ik}, \quad 1\leq i<j<k\leq N,$$ and $$\label{eq:izot-cross-r-op}
\lambda^{ij}=\lambda^{ij}_{(k)}, \quad \lambda^{jk} = \lambda^{jk}_{(i)}, \quad
\lambda^{ik}=\lambda^{ik}_{(j)}, \quad 1\leq i<j<k\leq N;$$ notice their consistency with the general system .

Equations - imply that cross-ratios of two dimensional sub-lattices of the generalized isothermic lattice satisfy condition of the form , i.e., $$\lambda^{ij}_{(ij)}\lambda^{ij} = \lambda^{ij}_{(i)}\lambda^{ij}_{(j)},
\quad 1\leq i<j\leq N.$$ For a fixed pair $i,j$, the above relation can be resolved as in (the first equation in asserts that the functions $\alpha^i$ and $\alpha^j$ are the same for all $i,j$ sublattices). Finally, equations imply the the functions $\alpha$ can be defined consistently on the whole lattice.

For convenience of the Reader we present also the algebraic proof of the above properties of the generalized isothermic lattice (see also [@BobSur] for analogous results concerning T-nets in a quadric).

Assume that solutions of the system of the discrete Moutard equations satisfy the quadratic constraint $$\label{eq:q-constr}
({{\boldsymbol x}}| {{\boldsymbol x}}) = 0,$$ where $(\cdot | \cdot ) $ is a symmetric nondegenerate bilinear form. Then the coefficients of the Moutard equations should be of the form $$\label{eq:f-isoth}
f^{ij} = \frac{({{\boldsymbol x}}| {{\boldsymbol x}}_{(i)} - {{\boldsymbol x}}_{(j)})}{({{\boldsymbol x}}_{(i)} | {{\boldsymbol x}}_{(j)})},
\quad 1\leq i<j\leq N.$$ Moreover by direct calculations one shows that $$({{\boldsymbol x}}_{(i)} | {{\boldsymbol x}})_{(j)} = ({{\boldsymbol x}}_{(i)} | {{\boldsymbol x}}),
\quad 1\leq i<j\leq N,$$ which implies that the products $({{\boldsymbol x}}_{(i)} | {{\boldsymbol x}}) $, which we denote by $\alpha_i$, are functions of single variables $n_i$.

Consider the points $x, x_{(i)}, x_{(j)}$ as the (projective) basis of the plane $\langle x, x_{(i)}, x_{(j)} \rangle$. Then the homogeneous coordinates of points of the plane can be written as $${{\boldsymbol y}}= t {{\boldsymbol x}}+ t_i {{\boldsymbol x}}_{(i)} + t_j {{\boldsymbol x}}_{(j)}, \qquad (t,t_i,t_j) \in {{\mathbb R}}^3_*,$$ modulo the standard common proportionality factor. In particular, the line $\langle x, x_{(i)} \rangle$ is given by equation $t_j = 0$, and the line $\langle x, x_{(j)} \rangle$ is given by equation $t_i = 0$. Due to the discrete Moutard equation the line $\langle x, x_{(ij)} \rangle$ is given by equation $t_i + t_j = 0$.

To find the cross-ratio $ {\mathrm{cr}}(x_{(i)},x_{(j)};x,x_{(ij)})=\lambda^{ij} $ via lines of the planar pencil with the base point $x$ we need equation of the tangent to the conic $({{\boldsymbol y}}| {{\boldsymbol y}}) =0$ at that point. It is easy to check that the conic is given by $$t t_i \alpha_{i} + t t_j \alpha_{j} + t_i t_j  ({{\boldsymbol x}}_{(i)} | {{\boldsymbol x}}_{(j)}) =0.$$ The tangent to the conic at $x$ is then given by $$t_i \alpha_{i} + t_j \alpha_{j} =0,$$ which implies equation .

It should be mentioned the a similar quadratic reduction of the discrete Moutard equation (for $N=2$) appeared in a paper of Wolfgang Schief [@Schief-C] under the name of discrete vectorial Calapso equation, as an integrable discrete vectorial analogue of the Calapso equation [@Calapso], which is of the fourth order and describes isothermic surfaces. It turns out that the discrete Calapso equation describes also the so called Bianchi reduction of discrete asymptotic surfaces [@DNS-Bianchi-ass].

Isothermic lattices in the Möbius sphere and the so called Clifford configuration
---------------------------------------------------------------------------------

In [@KoSch-Clifford] Konopelchenko and Schief observed that the “complex cross-ratio definition” of the discrete isothermic surfaces, when extended to three dimensional lattices, is related with the so called Clifford configuration of circles (see Figure \[fig:Clifford\]).

![The so called Clifford configuration of circles (the second Miquel configuration)[]{data-label="fig:Clifford"}](Clifford.eps){width="8cm"}

In this Section we would like to explain that fact geometrically.

Our point of view on the Clifford configuration is closely related to the geometric definition of the isothermic lattice in the Möbius sphere. Therefore we start with another, less restrictive, configuration of circles on the plane, the Miquel configuration (see Figure \[fig:Miquel\]), which provides geometric explanation of integrability of the circular lattice [@CDS]. When the quadric in Lemma \[lem:qred-hex\] is the standard sphere then the intersection curves of the planes of the quadrilaterals with the sphere are circles. After the stereographic projection from a generic point of the sphere we obtain the classical Miquel theorem [@Miquel], which can be stated as follows (given three distinct points $a$, $b$ and $c$, by $C(a,b,c)$ we denote the unique circle-line passing through them).

\[th:Miquel\] Given four coplanar points $s_0$, $s_i$, $i=1,2,3$. On each circle $C(s_0,s_i, s_j)$, $1\leq i < j \leq 3$ choose a point, denoted correspondingly by $s_{ij}$. Then there exists the unique point $s_{123}$ which belongs simultaneously to the three circles $C(s_1, s_{12}, s_{13})$, $C(s_2, s_{12}, s_{23})$ and $C(s_3, s_{13}, s_{23})$.

![The (first) Miquel configuration of circles[]{data-label="fig:Miquel"}](Miquel.eps){width="8cm"}

The additional assumption about coplanarity of the points $x_0$, $x_{12}$, $x_{13}$, $x_{23}$ on the sphere is then equivalent to the additional assumption about concircularity of the corresponding points $s_0$, $s_{12}$, $s_{13}$, $s_{23}$. In view of Lemma \[lem:BKP-hex\] we obtain therefore another configuration of circles, the so called Clifford configuration, which can be described as follows.

\[th:Clifford\] Under hypotheses of Theorem \[th:Miquel\] assume that the points $s_0$, $s_{12}$, $s_{13}$, $s_{23}$ are concircular, then the points $s_1$, $s_{2}$, $s_{3}$ and $s_{123}$ are concircular as well.

The original Clifford’s formulation of the above result was more symmetric. Its relation to Theorem \[th:Clifford\] is analogous to relation of the Cox theorem to Lemma \[lem:BKP-hex\]. We remark that although the above theorem is usually attributed (see for example [@Coxeter]) to Wiliam Clifford [@Clifford] it appeared in much earlier paper [@Miquel] of Auguste Miquel, where we read as *Thèorème II* the following statement: *Lorsqu’un quadrilatère complet curviligne ABCDEF est formé par quatre arcs de cercle AB, BC, CD, DA, qui se coupent tous quatre en un même point P, si l’on circonscrit des circonférences de cercle à chacun des quatre triangles curvilignes que forment les côtés de ce quadrilatère, les circonférences de cercle AFB, EBC, DCF, DAE ainsi obtenues se couperont toutes quatre en un même point G*. Points on Figure \[fig:Clifford\] are labelled in double way to visualize simultaneously the configuration in formulation of Theorem \[th:Clifford\] and in Miquel’s formulation.

The Darboux transformation of the generalized isotermic lattice {#sec:Darboux}
===============================================================

The fundamental, Moutard and Ribaucour transformations
------------------------------------------------------

Usually, on the discrete level there is no essential difference between integrable lattices and their transformations. The analogue of the fundamental transformation of Jonas for quadrilateral lattices is defined as construction of a new level of the lattice [@TQL] keeping the basic property of planarity of elementary quadrilaterals. Below we recall the relevant definitions of the fundamental transformation and its important reductions – the BQL reduction [@BQL] (algebraically equivalent to the Moutard transformation [@NiSchief]), and the QQL reduction [@q-red] called the Ribaucour transformation.

*The fundamental transform* of a quadrilateral lattice $x:{{\mathbb Z}}^N\to{{\mathbb P}}^M$ is a new quadrilateral lattice $\hat{x}:{{\mathbb Z}}^N\to{{\mathbb P}}^M$ constructed under assumption that for any point $x$ of the lattice and any direction $i$, the four points $x$, $x_{(i)}$, $\hat{x}$ and $\hat{x}_{(i)}$ are coplanar.

The fundamental transformation of a B-quadrilateral lattice $x:{{\mathbb Z}}^N\to{{\mathbb P}}^M$ constructed under additional assumption that for any point $x$ of the lattice and any pair $i,j$ of different directions, the four points $x$, $x_{(ij)}$, $\hat{x}_{(i)}$ and $\hat{x}_{(j)}$ are coplanar is called *the BQL (Moutard) reduction* of the fundamental transformation of $x$.

Algebraic description of the above transformation is given as follows [@NiSchief]. Given solution ${{\boldsymbol x}}$ of the system of discrete Moutard equations and given its scalar solution $\theta$, then the solution $\hat{{{\boldsymbol x}}}$ of the system $$\label{eq:Moutard-transf}
\hat{{\boldsymbol x}}_{(i)} - {{\boldsymbol x}}= \frac{\theta}{\theta_{(i)}}\left( \hat{{\boldsymbol x}}-
{{\boldsymbol x}}_{(i)}\right),$$ satisfies equations with the new potential $$\hat{f}^{ij} = f^{ij}\frac{\theta_{(i)}\theta_{(j)}}{\theta\theta_{(ij)}},
\qquad i<j,$$ and new $\tau$-function $$\label{eq:transf-Mout-tau}
\hat\tau = \theta\tau.$$

The fundamental transformation of a quadrilateral lattice $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ in a quadric, constructed under additional assumption that also $\hat{x}$ satifies the same quadratic constraint is called *the Riboucour transformation* of $x$.

The Darboux transformation
--------------------------

The fundamental transformation of a generalized isothermic lattice which is simultaneously the Ribaucour and the Moutard transformation is called *the Darboux transformation*.

Notice that there is essentially no difference between the Darboux transformation and construction of a new level of the generalized isothermic lattice. Therefore, given points $x$, $x_{(i)}$, $x_{(j)}$, $x_{(ij)}$, $i\ne j$, of the initial lattice, and given points $\hat{x}$, $\hat{x}_{(j)}$ of its Darboux transform then the point $\hat{x}_{(i)}$ is determined by the “half-hexahedron construction” visualized on Figure \[fig:twotothree\], i.e., $\hat{x}_{(i)}$ is the intersection point of the line $\langle x, x_{(i)}, \hat{x} \rangle \cap 
\langle x, x_{(ij)}, \hat{x}_{(j)} \rangle$ with the quadric. Moreover, Lemma \[lem:half-hex\] implies $$\label{eq:cr-prod-transf}
{\mathrm{cr}}(x_{(i)},\hat{x};x,\hat{x}_{(i)}) = {\mathrm{cr}}(x_{(i)},x_{(j)};x,x_{(ij)})
{\mathrm{cr}}(x_{(j)},\hat{x};x,\hat{x}_{(j)}),$$ while Corollary \[cor:hex\] gives $${\mathrm{cr}}(\hat{x}_{(i)},\hat{x}_{(j)};\hat{x},\hat{x}_{(ij)})=
{\mathrm{cr}}(x_{(i)},x_{(j)};x,x_{(ij)}).$$

The algebraic derivation of the above results is given below. The Darboux transformations of the discrete isothermic surfaces in the light-cone description were discussed in a similar spirit in [@BobSur].

\[prop:Darboux-prod\] If $\hat{x}:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ is a Darboux transform of the generalized isothermic lattice $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$ then the product $(\hat{{{\boldsymbol x}}}|{{\boldsymbol x}})$ of their homogeneous coordinates in the gauge of the linear problem and the corresponding Moutard transformation , with respect to the bilinear form $(\cdot|\cdot)$ defining the quadric $\mathcal{Q}$, is constant.

The homogeneous coordinates ${{\boldsymbol x}}$ and $\hat{{{\boldsymbol x}}}$ in the Moutard transformation satisfy equation of the form $$\hat{{{\boldsymbol x}}}_{(i)} - {{\boldsymbol x}}= f^i(\hat{{{\boldsymbol x}}} - {{\boldsymbol x}}_{(i)}), \qquad 1\leq i \leq N,$$ with appropriate functions $f^i:{{\mathbb Z}}^N\to{{\mathbb R}}$. The quadratic condition $(\hat{{{\boldsymbol x}}}_{(i)}|\hat{{{\boldsymbol x}}}_{(i)})=0$ together with other quadratic conditions give $$\label{eq:fi-isoth}
f^i = \frac{({{\boldsymbol x}}| \hat{{{\boldsymbol x}}} -{{\boldsymbol x}}_{(i)})}{(\hat{{{\boldsymbol x}}}|{{\boldsymbol x}}_{(i)})},$$ which implies $(\hat{{{\boldsymbol x}}}_{(i)}|{{{\boldsymbol x}}}_{(i)}) = (\hat{{{\boldsymbol x}}}|{{\boldsymbol x}})$.

Notice the above proof goes along the corresponding reasoning in the first part of the algebraic proof of Proposition \[prop:cross-ratio-gen\]. The analogous reasoning as in its second part gives the following statement.

Under hypothesis of Proposition \[prop:Darboux-prod\] denote $\zeta=(\hat{{{\boldsymbol x}}}|{{\boldsymbol x}})$ and $\lambda^i = {\mathrm{cr}}(\hat{x},x_{(i)};x,\hat{x}_{(i)})$ then equations and should be replaced by $$\label{eq:cr-lambda-prod}
\lambda^i = \frac{\zeta}{\alpha_i}, \qquad \lambda^i \lambda^{ij} = \lambda^j.$$

The above reasoning can be reversed giving the algebraic way to find the Darboux transform of a given generalized isothermic lattice.

\[th:Darboux\] Given a solution ${{\boldsymbol x}}:{{\mathbb Z}}^N\to {{\mathbb R}}^{M+1}_*$ of the system of Moutard equations satisfying the constraint $({{\boldsymbol x}}|{{\boldsymbol x}})=0$, considered as homogeneous coordinates of generalized isothermic lattice $x:{{\mathbb Z}}^N\to\mathcal{Q}\subset{{\mathbb P}}^M$, denote $\alpha_i = ({{\boldsymbol x}}_{(i)}|{{\boldsymbol x}})$. Given a point $[\hat{{{\boldsymbol x}}}_0]=\hat{x}_0\in\mathcal{Q}$, denote $\zeta =
(\hat{{{\boldsymbol x}}}_0|{{\boldsymbol x}}(0))$. Then there exists unique solution of the linear system $$\label{eq:Darboux}
\hat{{{\boldsymbol x}}}_{(i)} = {{\boldsymbol x}}+ \frac{\zeta - \alpha_i}{(\hat{{{\boldsymbol x}}}|{{\boldsymbol x}}_{(i)})}
(\hat{{{\boldsymbol x}}} - {{\boldsymbol x}}_{(i)}), \qquad, \quad 1\leq i \leq N,$$ with initial condition $\hat{{{\boldsymbol x}}}(0) = \hat{{{\boldsymbol x}}}_0$ which gives the Darboux transform of the lattice $x$. In particular $$\begin{aligned}
 \label{eq:cr-zeta-alpha}
{\mathrm{cr}}(\hat{x},x_{(i)};x,\hat{x}_{(i)}) & = \frac{\zeta}{\alpha_i} \\
{\mathrm{cr}}(\hat{x}_{(i)},\hat{x}_{(j)};\hat{x},\hat{x}_{(ij)}) & = 
{\mathrm{cr}}(x_{(i)},x_{(j)};x,x_{(ij)}) =
\frac{\alpha^i}{\alpha^j}.\end{aligned}$$

Before proving the Theorem let us state a Lemma relating the parameter $\zeta$ of the Darboux transformation with the functional parameter $\theta$ of the Moutard transformation .

\[eq:zeta-theta\] Under hypotheses of Theorem \[th:Darboux\] the solution $\theta$ of the system $$\theta_{(i)}=\theta\frac{(\hat{{{\boldsymbol x}}}|{{\boldsymbol x}}_{(i)})}{\zeta - \alpha_i}, \qquad
\quad 1\leq i \leq N,$$ satisfies the system of discrete Moutard equations with the coefficients given by .

By direct verification. Notice that both ways to calculate $\theta_{(ij)}$, $i\ne j$, from $\theta$ give the same result, and to do that we do not use compatibility of the system .

By direct calculation one can check that the system preserves the constraints $(\hat{{{\boldsymbol x}}}|\hat{{{\boldsymbol x}}}) =0$ and $(\hat{{{\boldsymbol x}}}|{{\boldsymbol x}})=\zeta$, moreover $(\hat{{{\boldsymbol x}}}_{(i)}|\hat{{{\boldsymbol x}}}) =({{\boldsymbol x}}_{(i)}|{{\boldsymbol x}})$. Compatibility of the system can be checked by direct calculation, but in fact it is the consequence of Lemma \[eq:zeta-theta\] and properties of the Moutard transformation .

Notice that because there is essentially no diffrence between the lattice directions and the transformation directions, the tranformation equations can be guessed by keeping the Moutard-like form supplementing it by calculation of the coefficient $f^i$ from the quadratic constraint. We will use this observation in the next Section where we consider the permutability principle for the Darboux transformations of generalized isothermic lattices.

The Bianchi permutability principle
-----------------------------------

The original Bianchi superposition principle for the Darboux transformations of the isothermic surfaces reads as follows [@Bianchi-isoth]: *Se dalla superficie isoterma $S$ si ottengono due nuove superficie isoterme $S_1$, $S_2$ mediante le trasformazioni di Darboux $D_{m_1}$, $D_{m_2}$ a costanti $m_1$, $m_2$ differenti, esiste una quarta superficie isoterma $\overline{S}$, pienamente determinata e costruibile in termini finiti, che è legata alla sua volta alle medesime superficie $S_1$, $S_2$ da due trasformazioni di Darboux $\overline{D}_{m_2}$, $\overline{D}_{m_1}$ colle costanti invertite $m_2$, $m_1$.*

![Geometric construction of the superposition of two Darboux transformations[]{data-label="fig:bianchi"}](bianchi)

Its version for generalized isothermic lattices can be formulated analogously.

\[prop:superp\] When from given generalized isothermic lattice $x$ there were constructed two new isothermic lattices $\hat{x}^1$ and $\hat{x}^2$ via the Darboux transformations with different parameters $\zeta_1$ and $\zeta_2$, then there exists the unique forth generalized isothermic lattice $\hat{x}^{12}$, determined in algebraic terms from the three previous ones, which is connected with two intermediate lattices $\hat{x}^1$ and $\hat{x}^2$ via two Darboux transformations with reversed parameters $\zeta_2$, $\zeta_1$.

The algebraic properties of the B-reduction of the fundamental transformation (the discrete Moutard transformation) imply that in the gauge of the linear problem and of the transformation equations the superposition of two such transformations reads $$\label{eq:-transf-sup}
\hat{{{\boldsymbol x}}}^{12} - {{\boldsymbol x}}= f (\hat{{{\boldsymbol x}}}^1 - \hat{{{\boldsymbol x}}}^2),$$ where $f$ is an appropriate function [@NiSchief; @BQL]. Because of the additional quadratic constraints the function is given by (compare also equations and ) $$\label{eq:f-isoth-sup}
f = \frac{({{\boldsymbol x}}| \hat{{{\boldsymbol x}}}^1 - \hat{{{\boldsymbol x}}}^2)}{(\hat{{{\boldsymbol x}}}^1 | \hat{{{\boldsymbol x}}}^2)}.$$ The lattice $\hat{x}^{12}$ with homogeneous coordinates given by and is superposition of two Darboux transforms. Finally, direct calculation shows that $$(\hat{{{\boldsymbol x}}}^{12} | \hat{{{\boldsymbol x}}}^2) =
(\hat{{{\boldsymbol x}}}^1 | {{\boldsymbol x}}) = \zeta_1, \qquad
(\hat{{{\boldsymbol x}}}^{12} | \hat{{{\boldsymbol x}}}^1)=
(\hat{{{\boldsymbol x}}}^2 | {{\boldsymbol x}})= \zeta_2 .$$

The final algebraic superposition formula reads $$\hat{{{\boldsymbol x}}}^{12} - {{\boldsymbol x}}= \frac{\zeta_1 - \zeta_2}{(\hat{{{\boldsymbol x}}}^1 | \hat{{{\boldsymbol x}}}^2)}
(\hat{{{\boldsymbol x}}}^1 - \hat{{{\boldsymbol x}}}^2),$$ while the cross-ratio of the four corresponding points calculated with respect to the conic intersection of the plane $\langle x,\hat{x}^1,\hat{x}^2 \rangle$ and the quadric is given by $$\label{eq:cr-zeta-sup}
 {\mathrm{cr}}(\hat{x}^1,\hat{x}^2;x,\hat{x}^{12}) = \frac{\zeta_1}{\zeta_2}.$$

To find the lattice $\hat{x}^{12}$ geometrically we can use again the “half-hexahedron construction” in the new context visualized on Figure \[fig:bianchi\] (compare with Figure \[fig:twotothree\]). Moreover, Lemma \[lem:half-hex\] gives $$\label{eq:cr-prod-sup}
{\mathrm{cr}}(\hat{x}^1,\hat{x}^2;x,\hat{x}^{12}) =
{\mathrm{cr}}(\hat{x}^1, x_{(i)};x,\hat{x}^1_{(i)}) \,
{\mathrm{cr}}(x_{(i)}, \hat{x}^2;x,\hat{x}^2_{(i)}),$$ which, due to equation , is in agreement with .

Two dimensional generalized isothermic lattice {#sec:isothermic}
==============================================

In the previous Sections we were mainly interested in generalized isothermic lattices of dimension greater then two. However, simultaneous application of the B-constraint and the quadratic constraint lowers dimensionality of the lattice (in the sense of the initial boundary value problem). One can see it from Figure \[fig:twotothree\], which implies that two intersecting strips made of planar quadrialterals with vertices in a quadric (see Figure \[fig:isoth-init\]) can be extended to a two dimensional quadrilateral lattice in the quadric. Because of Lemma \[lem:half-hex\] such lattice satisfies Steiner’s version of the cross-ratio constraint .

![Two intersecting initial strips of a two dimensional generalized isothermic lattice allow to build all the lattice. The additional transverse quadrilateral allows to build the lattice (together wih its Darboux transform) in a three dimensional fashion[]{data-label="fig:isoth-init"}](isoth-init)

One can define however geometrically two dimensional generalized isothermic lattices (generalized discrete isothermic surfaces) without using the three dimensional construction. An important tool here is the projective interpretation of the discrete Moutard equation [@DGNS] as representing quadrilateral lattice with additional linear relation between any of its points $x$ and its four second-order neighbours $x_{(\pm 1 \pm 2)}$. Geometrically, such five points of a *two dimensional B-quadrilateral lattice* are contained in a subspace of dimension three; for generic two dimensional quadrilateral lattice such points are contained in a subspace of dimension four. To exclude further degenerations we assume that no of the four points $x_{\pm 1}$, $x_{\pm 2}$ belongs to that three dimensional subspace.

A two dimensional B-quadrilateral lattice in a hyperquadric $x:{{\mathbb Z}}^2\to\mathcal{Q}\subset{{\mathbb P}}^M$ satisfying the local irreducibility condition is called a *generalized discrete isothermic surface*.

Notice that the above Definition gives (the conclusion was drawn by Alexander Bobenko) a geometric characterization of the classical discrete isothermic surfaces of Bobenko and Pinkall. Mainly, the non-trivial intersection of a three dimensional subspace with the Möbius sphere is a two dimensional sphere. After the stereographic projection, which preserves co-sphericity of points, a discrete isothermic surface in the Möbius (hyper)sphere gives *circular two dimensional lattice $s:{{\mathbb Z}}^2\to{{\mathbb E}}^M$ such that for any of its points $s$ there exists a sphere containig the point and its four second-order neighbours $s_{(\pm 1 \pm 2)}$*.

Notice that, actually, all calculations where we used simultaneously both the discrete Moutard equation and the quadratic constraint (algebraic proofs of Proposition \[prop:cross-ratio-gen\], Theorem \[th:Darboux\] and Proposition \[prop:superp\]) remain true for $N=2$. Therefore the corresponding results on the cross-ratio characterization of generalized discrete isothermic surfaces, their Darboux transformation and the Bianchi superposition principle are still valid.

To complete this Section let us present the geometric construction of a generalized discrete isothermic surface (see Figure \[fig:geom-constr-surf\]).

![Geometric construction of generalized discrete isothermic surfaces[]{data-label="fig:geom-constr-surf"}](geom-constr-surf)

The basic step of the construction, which allows to build the generalized discrete isothermic surface from two initial quadrilateral strips in a quadric in a two dimensional fashion can be discribed as follows. Consider the four dimensional subspace $V_4 = \langle x, x_{(1)}, x_{(2)}, x_{(-1)}, x_{(-2)} \rangle$, where the basic step takes place. Denote by $V_3 = \langle x, x_{(-1-2)}, x_{(-12)}, x_{(1-2)} \rangle$ its three dimensional subspace passing through the points $x$, $x_{(-1-2)}$, $x_{(-12)}$ and $x_{(1-2)}$, and by $V_2 = \langle x, x_{(1)}, x_{(2)} \rangle$ the plane of the elementary quatrilateral whose fourth vertex $x_{12}$ we are going to find. In the construction of the two dimensional B-quadrilateral lattice the vertex must belong to the line $V_1 = V_3 \cap V_2$. In our case it should also belong to the conic $\mathcal{C}=V_2 \cap \mathcal{Q}$. Because the conic contains already one point $x$ of the line $V_1$, the second point is unique. Notice that although the points $x_{(-1)}$ and $x_{(-2)}$ do not play any role in the construction, they can be easily recovered in a similar way as above.

Conclusions and discussion
==========================

In the paper we defined new integrable reduction of the lattice of planar quadrilaterals, which contains as a particular example the discrete isothermic surfaces. We studied, by using geometric and algebraic means, various aspects of such generalized isothermic lattices. In particular, we defined the (analogs of the) Darboux transformations for the lattices and we showed the corresponding permutablity principle.

The theory of integrable systems is deeply connected with results of geometers of the turn of XIX and XX centuries. The relation of integrability and geometry is even more visible on the discrete level, where into the game there enter basic results of the projective geometry. In our presentation of the generalized isothermic lattices the basic geometric results were a variant of the Möbius theorem and the generalization of the Miquel theorem to arbitrary quadric, which combined together gave the corresponding generalization of the Clifford theorem (known already to Miquel). An important tool in our research was also Steiner’s description of conics and the geometric properties of von Staudt’s algebra (see Appendix \[sec:cr\]).

Acknowledgements {#acknowledgements .unnumbered}
================

I would like to thank to Jaros[ł]{}aw Kosiorek and Andrzej Matraś for discusions cencerning incidence geometry and related algebraic questions. The main part of the paper was prepared during my work at DFG Research Center MATHEON in Institut für Mathematik of the Technische Universität Berlin. The paper was supportet also in part by the Polish Ministry of Science and Higher Education research grant 1 P03B 017 28. Finally, it is my pleasure to thank the organizers of the SIDE VII Conference for support.

The cross-ratio and the projective structure of a conic {#sec:cr}
=======================================================

For convenience of the Reader we have collected some facts from projective geometry (see, for example [@Samuel; @VeblenYoung]) used in the main text of the paper. Let $a,b,c$ be distinct points of the projective line $D$ over the field ${{\mathbb K}}$. Given $d\in D$, the cross-ratio ${\mathrm{cr}}(a,b;c,d)$ of the four points $a,b,c,d$ is defined as $h(d)\in
\hat{{\mathbb K}}={{\mathbb K}}\cup\{\infty\}$ where $h$ is the unique projective transformation $D\to\hat{{\mathbb K}}$ that takes $a$, $b$ and $c$ to $\infty$, $0$ and $1$, respectively. For $D=\hat{{\mathbb K}}$, with the usual conventions about operations with $0$ and $\infty$, the cross-ratio is given by $${\mathrm{cr}}(a,b;c,d) = \frac{(c-a)(d-b)}{(c-b)(d-a)}.$$ Denote by $\boldsymbol{a}$ and $\boldsymbol{b}$ homogenous coordinates of the points $a$ and $b$. If the homogenous coordinates of the points $c$ and $d$ collinear with $a$, $b$ are, respectively, $\boldsymbol{c}=\alpha\boldsymbol{a}+ \beta\boldsymbol{b}$ and $\boldsymbol{d}=\gamma\boldsymbol{a}+ \delta\boldsymbol{b}$, then $$\label{eq:cr-proj}
{\mathrm{cr}}(a,b;c,d) = \frac{\beta\gamma}{\alpha\delta}.$$ Let $D$ and $D^\prime$ be projective lines, $a,b,c,d$ distinct points on $D$, and $a^\prime,b^\prime,c^\prime,d^\prime$ distinct points on $D^\prime$. There exists a projective transformation $u:D\to D^\prime$ taking $a,b,c,d$ into $a^\prime,b^\prime,c^\prime,d^\prime$, respectively, if and only if the cross-ratios ${\mathrm{cr}}(a,b;c,d)$ and ${\mathrm{cr}}(a^\prime,b^\prime;c^\prime,d^\prime)$ are equal.

![Multiplication of cross-ratios on the line[]{data-label="fig:multiplication"}](multiplication)

Following von Staudt one can perform geometrically algebraic operations on cross-ratios (see [@VeblenYoung] for details). We will be concerned with geometric multiplication, which can be considered as “projectivization” of the Thales theorem – see the self-explanatory Figure \[fig:multiplication\].

The planar pencil of lines has the natural projective structure inherited from any line not intersecting its base. Let $\mathcal{C}$ be an irreducible conic in a projective plane, and $a\in\mathcal{C}$ a point. To each line $D$ of the pencil $F_a$ of the base $a$, we associate the second point where $D$ intersects $\mathcal{C}$ (see Figure \[fig:conic\]); we denote this point by $j_a(D)$. When $D$ is the tangent to $\mathcal{C}$ at $a$, let $j_a(D)$ be the point $a$. Thus $j_a$ is a bijection from $F_a$ to $\mathcal{C}$. If $b$ is another point on $\mathcal{C}$, the composition $j_b^{-1}\circ j_a$ is a projective transformation from $F_a$ to $F_b$.

![The projective structure of a conic[]{data-label="fig:conic"}](conic)

Thus the bijection from $F_a$ to $\mathcal{C}$ allows us to transport to $\mathcal{C}$ the projective structure of $F_a$. This structure does not depend on the point $a$. Conversely, any projective transformation between two pencils defines a conic.

Finally we present the relation between the complex cross-ratio of the Möbius geometry, and Steiner’s conic cross-ratio.

\[prop:Moebius-and-Steiner\] Four points $a,b,c,d\in\hat{{\mathbb C}}$ are cocircular or collinear if and only if their cross-ratio ${\mathrm{cr}}(a,b;c,d)_{{\mathbb C}}$ computed in $\hat{{\mathbb C}}$, is real. The cross-ratio ${\mathrm{cr}}(a,b;c,d)_{{\mathbb C}}$ is equal to the cross-ratio ${\mathrm{cr}}(a,b;c,d)$ computed using the real projective line structure of the line or the circle considered as a conic.

Auxiliary calculations {#sec:aux-calc}
======================

In this Appendix we would like to prove the “cross-ratio” characterization (Proposition \[prop:qred-cr\]) of locally irreducible quadrilateral lattices in quadrics, and then to give algebraic proof of the basic Lemma \[lem:half-hex\]. It turns out that in the course of our calculations will give also algebraic “down to earth” proofs of the basic Lemmas \[lem:gen-hex\], \[lem:BKP-hex\] and \[lem:qred-hex\].

Proposition \[prop:qred-cr\] is immediate consequence of the following result.

\[lem:qred-hex-cr\] Under hypotheses of Lemma \[lem:qred-hex\] and irreducibility of the conics, being intersections of the planes of the quadrilaterals with the quadric, the cross-ratios (defined with respect to the conics) of points on opposite sides of the hexahedron are connected by the following equation $$\label{eq:qred-cross-ratios} \begin{split}
{\mathrm{cr}}(x_1,x_2;x_0,x_{12})& 
\, {\mathrm{cr}}(x_{13},x_{23};x_3,x_{123}) {\mathrm{cr}}(x_2,x_3; x_0,x_{23}) \, {\mathrm{cr}}(x_{12},x_{13};x_1,x_{123}) = \\ 
&={\mathrm{cr}}(x_1,x_2;x_0,x_{13}) \, {\mathrm{cr}}(x_{13},x_{23};x_2,x_{123}).
\end{split}$$

Let us choose the points $x_0$, $x_1$, $x_2$ and $x_3$ as the basis of projective coordinate system in the corresponding three dimensional subspace, i.e., $${{\boldsymbol x}}_0 = [1:0:0:0], \quad {{\boldsymbol x}}_1 = [0:1:0:0],
\quad {{\boldsymbol x}}_2 = [0:0:1:0],\quad {{\boldsymbol x}}_3 = [0:0:0:1],$$ and denote by $[t_0 : t_1 : t_2 : t_3 ]$ the corresponding homogeneous coordinates. For generic points $x_{ij}\in\langle x_0, x_{i}, x_{j} \rangle$, $1\leq i < j
\leq 3$, with homogeneous coordinates $${{\boldsymbol x}}_{12} =  [a_0:a_1:a_2:0], \quad
{{\boldsymbol x}}_{13} =  [b_0:b_1:0:b_3], \quad
{{\boldsymbol x}}_{23} =  [c_0:0:c_2:c_3],$$ one obtains, via the standard linear algebra, equations of the planes $\langle x_1, x_{12}, x_{13} \rangle$, $\langle x_2, x_{12}, x_{23} \rangle$, $\langle x_3, x_{13}, x_{23} \rangle$ respectively, $$\begin{aligned}
a_2 b_3 t_0 = & a_0 b_3 t_2 + b_0 a_2 t_3,\\
a_1 c_3 t_0 = & a_0 c_3 t_1 + c_0 a_1 t_3,\\
\label{eq:til_pi_12}
b_1 c_2 t_0 = & b_0 c_2 t_1 + c_0 b_1 t_2.\end{aligned}$$ The intersection point $x_{123}$ of the planes has the following coordinates ${{\boldsymbol x}}_{123}=[y_0:y_1:y_2:y_3]$ $$\label{eq:x_123-gen} \begin{split}
y_0 = & a_0 b_0 c_0 
\left( \frac{1}{a_2 b_1 c_3} + \frac{1}{a_1 b_3 c_2} \right), \\
y_1 = &  \frac{b_0 c_0}{b_3 c_2} +  \frac{a_0 c_0}{a_2 c_3} -
\frac{c_0^2}{c_2 c_3}  , \\
y_2 = &  \frac{a_0 b_0}{a_1 b_3} +  \frac{b_0 c_0}{b_1 c_3} -
\frac{b_0^2}{b_1 b_3}  , \\
y_3 = &  \frac{a_0 c_0}{a_1 c_2} +  \frac{a_0 b_0}{a_2 b_1} -
\frac{a_0^2}{a_1 a_2}  .
\end{split}$$ Up to now we have not used the additional quadratic restriction, and what we have done was just the algebraic proof of Lemma \[lem:gen-hex\].

Any quadric ${{\mathcal Q}}$ passing through $x_0$, $x_1$, $x_2$ and $x_3$ must have equation of the form $$\label{eq:quadric}
a_{01}t_0t_1 + a_{02}t_0t_2 + a_{03}t_0t_3 +
a_{12}t_1t_2 + a_{13}t_1t_3 + a_{23}t_2t_3 = 0.$$ The homogeneous coordinates of the points $x_{12}$, $x_{23}$ and $x_{13}$ can be parametrized in terms of the corresponding cross-ratios $\lambda = {\mathrm{cr}}(x_1,x_2;x_0,x_{12})$, $\nu = {\mathrm{cr}}(x_2,x_3;x_0,x_{23})$ and $\mu = {\mathrm{cr}}(x_1,x_3;x_0,x_{13})$ as $$\begin{aligned}
 \label{eq:x12}
{{\boldsymbol x}}_{12} &= 
\Big[\, \frac{\lambda a_{12}}{1-\lambda} \,:\, -\lambda a_{02} \,:\, 
a_{01}\,: \,0 \, \Big], \\
\label{eq:x23}
{{\boldsymbol x}}_{23} &=  
\Big[\, \frac{\nu a_{23}}{1-\nu}\,: \, 0\,: \, -\nu  a_{03} \,: \,a_{02}
\,\Big],\\
\label{eq:x_13-quadr} 
{{\boldsymbol x}}_{13} &= 
\Big[\, \frac{\mu a_{13}}{1-\mu} \,: \, -\mu  a_{03}\,: \, 0 \,: \, a_{01}
\,\Big].\end{aligned}$$ We will only show how to find the homogeneous coordinates of $x_{12}$ in terms of $\lambda$. Let us parametrize points of the conic $$\label{eq:conic}
a_{01}t_0t_1 + a_{02}t_0t_2 +  a_{12}t_1t_2 = 0,$$ being intersection of the quadric with the plane $t_3=0$, by the planar pencil with base at $x_0$. The point $x_0$ corresponds to the tangent to the conic at $x_0$ $$a_{01}t_1 + a_{02}t_2 = 0,$$ while the points $x_1$ and $x_2$ correspond to lines $t_2=0$ and $t_1 = 0$, respectively. The line $\langle x_0, x_{12} \rangle$ must have equation (see equation in Appendix \[sec:cr\]) of the form $$a_{01}t_1 + \lambda a_{02}t_2 = 0,$$ which inserted into equation of the conic gives, after exluding the point $x_0$, the homogeneous coordinates of the point $x_{12}$.

Inserting expressions , and into formulas we obtain homogeneous coordinates $[y_0:y_1:y_2:y_3]$ of the point $x_{123}$ parametrized in terms of the cross-ratios $\lambda$, $\nu$, $\mu$ $$\label{eq:x_123-quadr} 
\begin{split}
y_0 = & \frac{a_{12}a_{23}a_{13}}{a_{01} a_{02} a_{03}} 
\frac{\lambda \nu - \mu}{(1-\lambda)(1-\mu)(1-\nu)}, \\
y_1 = &  \frac{a_{23}}{1-\nu} \left( 
\frac{a_{13}}{a_{01}a_{03}}\frac{\mu}{1-\mu} -
\frac{a_{23}}{a_{02}a_{03}}\frac{\nu}{1-\nu} -
\frac{a_{12}}{a_{01}a_{02}}\frac{\lambda\nu}{1-\lambda} \right)  , \\
y_2 = & \frac{a_{13}}{1-\mu} \left( 
-\frac{a_{13}}{a_{01}a_{03}}\frac{\mu}{1-\mu} +
\frac{a_{23}}{a_{02}a_{03}}\frac{\nu}{1-\nu} +
\frac{a_{12}}{a_{01}a_{02}}\frac{\mu}{1-\lambda} \right) , \\
y_3 = &  \frac{a_{12}}{1-\lambda} \left( 
\frac{a_{13}}{a_{01}a_{03}}\frac{\lambda}{1-\mu} -
\frac{a_{23}}{a_{02}a_{03}}\frac{1}{1-\nu} -
\frac{a_{12}}{a_{01}a_{02}}\frac{\lambda}{1-\lambda} \right) .
\end{split}$$ One can check that such expressions do satisfy the quadric equation , i. e., we have obtained the direct proof of Lemma \[lem:qred-hex\] under additional assumption of ireducibility of the conics.

Further calculations give the cross-ratios on remaining three sides of the hexahedron $$\begin{aligned}
 \label{eq:qred-cross-ratios-opposite}
{\mathrm{cr}}(x_{13}, x_{23}; x_3, x_{123}) & = 
 - \frac{\mu (1-\nu) a_{13} y_1}{\nu (1-\mu) a_{23} y_2} , \nonumber \\
{\mathrm{cr}}(x_{12}, x_{13}; x_1, x_{123}) & =
- \frac{(1-\mu) a_{12} y_2}{(1-\lambda) a_{13} y_3}, \\ 
{\mathrm{cr}}(x_{12}, x_{23}; x_2, x_{123}) & = 
\lambda \frac{(1-\nu) a_{12} y_1}{(1-\lambda) a_{23} y_3}, \nonumber\end{aligned}$$ with $y_i$ given by , which implies equation .

We can express coplanarity of four points in ${{\mathbb P}}^3$ as vanishing of the determinant of the matrix formed by their homogeneous coordinates. In notation of the proof above (up to equation ) we have $$y_0 = \det \left(
{{\boldsymbol x}}_{123} \: {{\boldsymbol x}}_{1} \: {{\boldsymbol x}}_{2} \: {{\boldsymbol x}}_{3} \right) = 
 \frac{a_0b_0c_0}{a_1 a_2 b_1  b_3 c_2 c_3} 
\det \left( 
{{\boldsymbol x}}_0 \: {{\boldsymbol x}}_{12} \: {{\boldsymbol x}}_{23} \: {{\boldsymbol x}}_{13}  \right),$$ which implies the statement of Lemma \[lem:BKP-hex\].

Finally, when vertices of the hexahedron are contained in the quadric, condition $y_0=0$ inserted into equation gives $\mu = \lambda \nu$.

[10]{}

J. Bertrand, *Mémoire sur les surfaces isothermes orthogonales*, J. Math. Pur. Appl. (Liouville J.) **9** (1844), 117-130, `http://portail.mathdoc.fr/JMPA/`.

L. Bianchi, *Il teorema di permutabilità per le trasformazioni di Darboux delle superficie isoterme*, Rend. Acc. Naz. Lincei **13** (1904), 359–367.

L. Bianchi, *Lezioni di geometria differenziale*, Zanichelli, Bologna, 1924.

A. Bobenko, *Discrete conformal maps and surfaces*, Symmetries and Integrability of Difference Equations (P. Clarkson and F. Nijhoff, eds.), Cambridge University Press, 1999, pp. 97–108.

A. Bobenko and U. Pinkall, *Discrete surfaces with constant negative Gaussian curvature and the Hirota equation*, J. Diff. Geom. **43** (1996), 527–611.

A. Bobenko and U. Pinkall, *Discrete isothermic surfaces*, J. Reine Angew. Math. **475** (1996), 187–208.

A. Bobenko and U. Pinkall, *Discretization of surfaces and integrable systems*, Discrete integrable geometry and physics, (A. Bobenko and R. Seiler, eds.), Clarendon Press, Oxford, 1999, pp. 3–58.

A. Bobenko and Yu. Suris, *Discrete differential geometry. Consistency as integrability*, [arXiv:math.DG/0504358]{}.

A. Bobenko and Yu. Suris, *Isothermic surfaces in sphere geometries as Moutard nets*, `arXiv:math.DG/0610434`.

P. Calapso, *Sulla superficie a linee di curvature isoterme*, Rend. Circ. Mat. Palermo **17** (1903), 275–286.

J. Cieśliński, *The Bäcklund transformations for discrete isothermic surfaces*, Symmetries and Integrability of Difference Equations (P. A. Clarkson and F. W. Nijhoff, eds.), University Press, Cambridge, 1999, pp. 109–121.

J. Cieśliński, P. Goldstein and A. Sym, *Isothermic surfaces in $\mathbb{E}^3$ as soliton surfaces*, Phys. Lett. A **205** (1995), 37–43.

J. Cieśliński, A. Doliwa and P. M. Santini, *The integrable discrete analogues of orthogonal coordinate systems are multidimensional circular lattices*, Phys. Lett. A **235** (1997), 480–488.

W. K. Clifford, *A synthetic proof of Miquels’s theorem*, Oxford, Cambridge and Dublin Messenger of Math. **5** (1871), 124–141.

H. S. M. Coxeter, *Introduction to geometry*, Wiley, New York, 1969.

G. Darboux, *Sur les surfaces isothermiques*, Ann. Sc. É. N. S., **16** (1899), 491–508,\
`http://www.numdam.org/item?id=ASENS_1899_3_16__491_0`.

G. Darboux, *Leçons sur la théorie générale des surfaces. [I–IV]{}*, Gauthier – Villars, Paris, 1887–1896.

G. Darboux, *Leçons sur les systémes orthogonaux et les coordonnées curvilignes*, Gauthier-Villars, Paris, 1910.

A. Doliwa, [*Geometric discretisation of the Toda system*]{}, Phys. Lett. A [**234**]{} (1997), 187–192.

A. Doliwa, *Quadratic reductions of quadrilateral lattices*, J. Geom. Phys. **30** (1999), 169–186.

A. Doliwa, *Discrete asymptotic nets and W-congruences in Plücker line geometry*, J. Geom. Phys. **39** (2001), 9–29.

A. Doliwa, *Geometric discretization of the [K]{}oenigs nets*, J. Math. Phys. **44** (2003), 2234–2249.

A. Doliwa, *The B-quadrilateral lattice, its transformations and the algebro-geometric construction*, J. Geom. Phys. **57** (2007), 1171–1192.

A. Doliwa, *Generalized isothermic lattices*, talk given at the Workshop “Geometry and Integrable Systems” (Berlin, 3-7 June 2005), `http://www.math.tu-berlin.de/geometrie/MISGAM/workshop2005/abstracts.html`.

A. Doliwa and P. M. Santini, *Integrable dynamics of a discrete curve and the Ablowitz-Ladik hierarchy*, J. Math. Phys. **36** (1995), 1259–1273.

A. Doliwa and P. M. Santini, *Multidimensional quadrilateral lattices are integrable*, Phys. Lett. A **233** (1997), 365–372.

A. Doliwa and P. M. Santini, *The symmetric, [D]{}-invariant and [E]{}gorov reductions of the quadrilateral lattice*, J. Geom. Phys. **36** (2000), 60–102.

A. Doliwa and P. M. Santini, *Integrable systems and discrete geometry*, Encyclopedia of Mathematical Physics, (J. P. Françoise, G. Naber and T. S. Tsun, eds.) Vol. III, Elsevier, 2006, pp. 78-87.

A. Doliwa, S. V. Manakov and P. M. Santini, *$\bar\partial$-reductions of the multidimensional quadrilateral lattice: the multidimensional circular lattice*, Comm. Math. Phys. **196** (1998), 1–18.

A. Doliwa, P. M. Santini and M. Mañas, *Transformations of quadrilateral lattices*, J. Math. Phys. **41** (2000), 944–990.

A. Doliwa, M. Nieszporski and P. M. Santini, *Asymptotic lattices and their integrable reductions I: the Bianchi and the Fubini-Ragazzi lattices*, J. Phys. A **34** (2001), 10423–10439.

A. Doliwa, M. Nieszporski and P. M. Santini, *Geometric discretization of the Bianchi system*, J. Geom. Phys. **52** (2004), 217–240.

A. Doliwa, P. Grinevich, M. Nieszporski, and P. M. Santini, *Integrable lattices and their sub-lattices: from the discrete Moutard (discrete Cauchy–Riemann) 4-point equation to the self-adjoint 5-point scheme*, J. Math. Phys. (in press) `arXiv:nlin.SI/0410046`.

L. P. Eisenhart, *Separable systems of Stäckel*, Ann. of Math. **35** (1934), 284–305.

U. Hertrich-Jeromin, *Introduction to Möbius differential geometry*, Cambridge University Press, 2003.

U. Hertrich-Jeromin, T. Hoffmann and U. Pinkall, *A discrete version of the Darboux transform for isothermic surfaces*, Discrete integrable geometry and physics, (A. Bobenko and R. Seiler, eds.), Clarendon Press, Oxford, 1999, pp. 59–81.

T. Hoffmann, *Discrete S-isothermic and S-cmc surfaces*, talk given at the Workshop “Geometry and Integrable Systems” (Berlin, 3-7 June 2005),\
`http://www.math.tu-berlin.de/geometrie/MISGAM/workshop2005/abstracts.html`.

P. Klimczewski, M. Nieszporski and A. Sym, *Luigi Bianchi, Pasquale Calapso and solitons*, Rend. Sem. Mat. Messina, Atti del Congresso Internazionale in onore di Pasquale Calapso, Messina 1998 (2000), 223–240.

B. G. Konopelchenko and W. K. Schief, *Three-dimensional integrable lattices in [Euclidean]{} spaces: Conjugacy and orthogonality*, Proc. Roy. Soc. London A **454** (1998), 3075–3104.

B. G. Konopelchenko and W. K. Schief, *Conformal geometry of the (discrete) Schwarzian Davey-Stewartson II hierarchy*, Glasgow Math. J. **47A** (2005), 121–131.

G. Lamé, *Mémoire sur les surfaces isothermes dans les corps solides homogènes en équilibre de température*. J. Math. Pur. Appl. (Liouville J.) **2** (1837), 147-183, `http://portail.mathdoc.fr/JMPA/`.

G. Lamé, [*Leçons sur les coordonnées curvilignes et leurs diverses applications*]{}, Mallet–Bachalier, Paris, 1859, `http://gallica.bnf.fr/ark:/12148/bpt6k99670t`.

A. Miquel, *Théorèmes sur les intersections des cercles et des sphères*, J. Math. Pur. Appl. (Liouville J.) **3** (1838), 517-522, `http://portail.mathdoc.fr/JMPA/`.

T. Miwa, *On Hirota’s difference equations*, Proc. Japan Acad. **58** (1982), 9–12.

M. Nieszporski, *A Laplace ladder of discrete Laplace equations*, Theor. Math. Phys. **133** (2002), 1576–1584.

J. J. C. Nimmo and W. K. Schief, *Superposition principles associated with the [Moutard]{} transformation. [An]{} integrable discretisation of a (2+1)-dimensional sine-[Gordon]{} system*, Proc. R. Soc. London A **453** (1997), 255–279.

D. Pedoe, *Geometry, a comprehensive course*, Dover, New York, 1988.

C. Rogers and W. K. Schief, *Bäcklund and Darboux transformations. Geometry and modern applications in soliton theory*, Cambridge University Press, Cambridge, 2002.

P. Samuel, *Projective geometry*, Springer, New York–Berlin–Heidelberg, 1988.

R. Sauer, *Projective Liniengeometrie*, de Gruyter, Berlin–Leipzig, 1937.

R. Sauer, *Differenzengeometrie*, Springer, Berlin, 1970.

W. K. Schief, *Isothermic surfaces in spaces of arbitrary dimension: integrability, discretization and [B]{}äcklund transformations – [A]{} discrete [C]{}alapso equation*, Stud. Appl. Math. **106** (2001), 85–137.

A. Sym, *Soliton sufaces and their applications*, Geometric aspects of the Einstein equations and integrable systems, Lecture Notes in Physics **239**, (R. Martini, ed.), Springer, 1985, pp. 154–231.

O. Veblen and J. W. Young, *Projective geometry*, Ginn and Co., Boston, 1910–1918.

J. Wallner and H. Pottmann, *Infinitesimally flexible meshes and discrete minimal surfaces*, Geometry Preprint 162, TU Wien, 2006, `http://www.geometrie.tuwien.ac.at/wallner/cmin.pdf`.

[^1]: $\ddagger$ Supported in part by the DFG Research Center MATHEON
