SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS MATTEO PETRERA AND YURI B SURIS Institut für Mathemati, MA 7- Technische Universität Berlin, Str des 7 Juni 36, 063 Berlin, Germany arxiv:08365v [math-ph] 7 Aug 0 Abstract We prove that the cosine law for spherical triangles and spherical tetrahedra defines integrable systems, both in the sense of multidimensional consistency and in the sense of dynamical systems Introduction Nowadays, it is a well accepted fact that many fundamental notions and results of differential geometry are closely related with the theory of integrable systems Moreover, the recent development of discrete differential geometry [] led to uncovering and identifying the true roots of such relations in the most fundamental incidence theorems of projective geometry and its metric subgeometries On the other hand, this development led also to a clarification of the notion of integrability, which is conceptually much more clear at the level of discrete systems In particular, the notion of multidimensional consistency as integrability crystallized in the course of this development (see []) In the present paper, we offer a new evidence in favor of the thesis that the very basic structures of geometry are integrable in their nature We show that the cosine laws for spherical simplices define integrable systems in several senses Our main results are as follows: We show that the map sending planar angles of a spherical triangle to its side lengths is integrable in the sense of multidimensional consistency (Sect 3) Moreover, we identify this map with the symmetric reduction of the well nown discrete Darboux system, which appeared for the first time in the context of discrete conjugate nets We also discuss the relation of this result with a geometric realization of tetrahedral angle structures, due to Feng Luo [] We show that the same map, interpreted as a dynamical system, is completely integrable in the Liouville-Arnold sense (Sect 4) This map can be considered as a non-rational discretization of the famous Euler top, which is a basic example of integrable system Moreover, we identify the second iterate of this map as a birational discretization of the Euler top, introduced by Hirota and Kimura [5] and studied in detail in our previous publication [5] This second iterate is shown to admit a nice geometric interpretation as a particularly simple and nice map between spherical triangles, which we call switch (the angles of the switched triangle are set to be equal to the sides of the original one) A closely related involutive map between spherical triangles was studied in a paper by Jonas [7], which can be thus considered as an early instance of appearance of integrable systems in elementary geometry E-mail: petrera@mathtu-berlinde, suris@mathtu-berlinde
MATTEO PETRERA AND YURI B SURIS We show that the map sending dihedral angles of a spherical tetrahedron to its side lengths is a further example of an integrable dynamical system (Sect 5) For this map, we establish conserved quantities, which turn out to be a re-interpretation of the sine laws, and an invariant volume form The latter is shown to be closely related, via the famous Schläfli formula, to the determinant of the Hesse matrix of the volume of the spherical tetrahedron as a function of dihedral angles For convenience of the reader, we recall (in Sect ) the standard derivation of the spherical cosine and sine laws via Gram matrices Cosine law for spherical simplices The maps we are dealing with in the present paper have their origin in the cosine law for spherical simplices We understand a spherical (n )-simplex as an intersection of n hemispheres in S n in general position Let v,, v n S n be the vertices of a spherical (n )-simplex, and let v,, v n S n be the vertices of the polar simplex, defined by the conditions v i, v j = 0 if i j, v i, v i > 0 The edge lengths l ij (0, π) and the dihedral angles α ij (0, π) are given by cos l ij = v i, v j, cos α ij = v i, v j Each set of data {l ij } or {α ij } defines a spherical simplex uniquely, and the spherical cosine law gives a transition from one set of data to another The most direct way to derive the cosine law is through the so called Gram matrices G = ( cos α ij ) n i,j=, G = (cos l ij ) n i,j= According to Lemma in [0], Gram matrices of spherical simplices are characterized by the following properties: they are positive definite symmetric matrices with diagonal entries equal to Introducing the matrices V and W whose columns are v i, resp v i, we can write: W T W = G, V T V = G, V T W = D, where D is a diagonal matrix with positive entries There follows: G = DG D, () from which the cosine laws follow immediately They can be put as cos l ij = g ij, cos α ij = g ij, () gii g jj g ii g jj where g ij, g ij are the cofactors of the Gram matrices G, G, respectively, see, for instance, [8] We now consider in more detail the particular cases n = 3, 4 of this construction Spherical triangles For a spherical triangle, α ij is just the plane (inner) angle opposite to the edge l ij, see Fig The angles of a spherical triangle belong to the tetrahedron-shaped domain T = {α = (α, α 3, α 3 ) (0, π) 3 : α + α 3 + α 3 > π, α ij + α i + α j < π} Here and below we use (i, j, ) to denote any permutation of (,,3) The sides belong to the domain T = {l = (l, l 3, l 3 ) (0, π) 3 : l + l 3 + l 3 < π, l ij + l i + l j > 0},
Spherical triangles Let (α i ) S be a spherical triangle with interior angles α i (0, π) andedges i, i =,, 3 The edge i, whose opposite angle is α i, lins the vertices v j v, see Figure below Hereafter (i, j, ) denotesanycyclicpermutationof(,, 3) We refer to [5] for elementary properties of spherical triangles SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 3 α 3 3 α α 3 3 3 Figure A spherical triangle The following conditions are necessary and sufficient for the existence of : α + α + α 3 > π, α i + α j + α < π which can be obtained from T via the reflection I : T T, I(α) = π α The space The quantities i and α i are are related by the law of cosines, which can be derived of spherical triangles (up to orthogonal transformations of R starting from the Gram matrix associated with : 3 ) is parametrized by either of T or T : a spherical triangle exists if and only if α T G := cos α, resp if and only if l T The cosine laws for spherical triangles, as given 3 in cos (), α read: cos α 3 cos α cos α cos α cos l ij = cos α ij + cos α i cos α j, (3) We have: sin α i sin α j δ := det G = cos α α α 3 cosα cos α cos α 3 α ij = cos l ij cos l i cos l j (4) s s α sin l i s sin α l j s α3 = 4cos cos cos cos We will consider formulas (3 4) as defining the mutually inverse maps where s := α + α + α 3 > π Thus we have δ > 0 We construct the matrix of cofactors (a ij ) i,j 3, with a ij := Φ : ( ) T i+j G ij, T where, G ij Φ is the : T (i, j)-th minor T, of G The law of cosines is given by α l, cos i = a l α j, () ajj a We note that Φ = I Φ I The cosine law can be supplemented by the sine law for spherical triangles: sin l = sin l 3 = sin l 3 d = sin α sin α 3 sin α 3 γ = γ, (5) d where d = det G = cos α cos α 3 cos α 3 cos α cos α 3 cos α 3, d = det G = cos l cos l 3 cos l 3 + cos l cos l 3 cos l 3, and γ = sin α sin α 3 sin α 3, γ = sin l sin l 3 sin l 3 The sine law is derived by considering principal minors of the relation () Upon the change of variables x ij = cos α ij, y ij = cos l ij,
4 MATTEO PETRERA AND YURI B SURIS the cosine law for spherical triangles taes the form We denote these maps by y ij = x ij = x ij + x i x j, (6) x i x j y ij y i y j (7) y i yj φ : τ τ, x y, φ : τ τ, y x, and call them the algebraic forms of Φ, Φ, respectively They are defined on the corresponding subsets τ, τ of [, ] 3, and are related by φ = i φ i, i(x) = x The sine law in these new variables reads: where y x = y 3 x 3 = y 3 x 3 = d γ = (γ ) d, (8) d = x x 3 x 3 x x 3 x 3, d = y y 3 y 3 + y y 3 y 3, (9) and γ = ( x )( x 3)( x 3), (γ ) = ( y )( y 3)( y 3) For further reference, we mention that the sine law can be represented in either of the following forms: d = ( y ij)( x i)( x j), d = ( x ij)( y i)( y j) (0) Spherical tetrahedra For a spherical tetrahedron, see Fig, the cosine law is expressed by the formula () with g ij = x ij + x i x j + x im x jm + x m (x i x jm + x im x j x ij x m ), () g ii = (x j + x jm + x m) x j x jm x m, () where we adopt the notation x ij = cos α ij and (i, j,, m) stands for any permutation of (,, 3, 4) These formulae are supplemented (see [, 9]) by the sine law, which can be put as and by the formulas sin l sin l 34 sin α sin α 34 = sin l 3 sin l 4 sin α 3 sin α 4 = sin l 4 sin l 3 sin α 4 sin α 3 = d γ = γ d, (3) cos l cos l 34 cos l 3 cos l 4 cos α cos α 34 cos α 3 cos α 4 = cos l 3 cos l 4 cos l 4 cos l 3 cos α 3 cos α 4 cos α 4 cos α 3 = cos l 4 cos l 3 cos l cos l 34 cos α 4 cos α 3 cos α cos α 34 = d γ = γ d, (4) where d = det G, d = det G, and γ = g g g 33 g 44, γ = g g g 33g 44
SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 5 MATTEO PETRERA AND YURI B SURIS SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 5 4 4 α (4) 3 α (4) 3 α (4) (4) 3 α (4) 3 < α < α (4) < α 3 α 3 < 3 α < α 3 < α 3 Figure Spherical triangle Edges i and angles α i are related by (-) Figure A spherical tetrahedron cos i = cos α i +cosα j cos α We now consider a somewhat fancy combinatorial interpretation, of formulas (6)() for the cosine law for spherical triangles In this sin interpretation, α j sin α we combinatorially assign the three inner angles of a spherical costriangle α i = cos and i its cos three j cos side, lengths to the six -faces () of an elementary cube, in such a manner thatsin each j sin angle and an opposite side are assigned to two see, for opposite instance, faces [] ofhere the and cube below More we use precisely, (i, j, ) we todenoteacyclicpermutationof assign the quantities x j to the three (,,3) facesthe of asine 3-dimensional law for spherical cubetriangles parallel reads: to the coordinate plane j, and the quantities y j = T i x j to the three opposite faces, sin = sin see Fig = sin 3 Here, 3 D T i stands for the unit shift in the i-th coordinate direction = sin α sin α sin α 3 γ, (3) where D = cos α cos α cos α 3 cosα cos α cos α 3, (4) x 3 y γ = sinα sin α sin α 3 x 3 We will consider formulas (), () as defining the mutually inverse maps Φ : T T, Φ : T y 3 T, (5) α, α We note that Φ = I Φ I Upon the change of variables 3 3 3 Figure A spherical tetrahedron 3 Spherical triangles and symmetric discrete Darboux system We nowthe consider inner angles a somewhat α i determine fancy the sides combinatorial i uniquely (and interpretation vice versa), according of formulas to the() for the cosine so called 3 lawspherical cosine for spherical law: triangles triangles andinsymmetric this interpretation, discrete Darboux we combinatorially system assign the three inner angles of a spherical triangle and its three side lengths to the six -faces of an elementary cube, in such a manner that each angle and an opposite side are assigned to two opposite faces of the cube More precisely, we assign the quantities x j to the three faces of a 3-dimensional cube parallel to the coordinate plane j, and the quantities y j = T i x j to the three opposite faces, cf Fig 3 Here, T i stands for the unit shift in the i-th coordinate direction x x i =cosα i, y i =cos i, Figure 3 A map on an elementary 3D cube with fields assigned to -faces Figure 3 A map on an elementary 3D cube with fields assigned to -faces We interpret y ij as T x ij, the shift of x ij in the coordinate direction Thus, what we consider now is, in the terminology of [], a 3D system with fields assigned to elementary squares, given by the formulas Thus, what we consider(tnow 3 x, is, T xin 3, the T x 3 terminology ) = φ (x, x 3 of, x 3 [], ), a 3D system with fields assigned to elementary squares, given by the formulas (T 3 x,t x 3,T x 3 )=φ(x,x 3,x 3 ), (3) y 3 3 3
6 MATTEO PETRERA AND YURI B SURIS or, explicitly, T 3 x = x + x 3 x 3 x 3 x 3, T x 3 = x 3 + x x 3 x x 3, (5) T x 3 = x 3 + x x 3 x x 3 This 3D system will be called the symmetric discrete Darboux system It can be seen as the special case (symmetric reduction) x ij = x ji of the general discrete Darboux system, given by x ij + x i x j T x ij = (6) xi x i xj x j System (6) is well nown in the theory of discrete integrable systems of geometric origin It seems to have appeared for the first time in [8, eq (70)] as a gauge version of a simpler but less symmetric form of the discrete Darboux system, T x ij = x ij + x i x j x j x j, (7) derived earlier in [3] Either form of the discrete Darboux system describes the so called rotation coefficients of discrete conjugate nets f : Z 3 RP n, ie, nets with planar elementary quadrilaterals, see [, p 4] The version of rotation coefficients satisfying system (6) has an advantage of being defined in more local terms, see [, p 95 96] Besides, system (6) admits, unlie system (7), the symmetric reduction x ij = x ji, and is therefore better suited for the description of the so called symmetric conjugate nets and discrete Egorov nets, see [4] The symmetric discrete Darboux system is a close relative of the discrete CKP system, see [7] Thus, we have proved that the symmetric discrete Darboux system (5) admits an alternative interpretation in terms of spherical triangles This interpretation leads also to a new proof of integrability of the symmetric discrete Darboux system Let us recall that integrability of discrete 3D systems is synonymous with their 4D consistency [] The meaning of this notion can be described as follows Consider the initial value problem with the data x ij = x ji, i, j =,, 3, 4, prescribed at six squares adjacent to one common vertex of the 4D cube Then the application of a 3D map lie (5) to the four 3D cubes adjacent to this vertex allows us to determine all T x ij At the second stage, the map is applied to the other four 3-faces of the 4D cube, with the result being all T m (T x ij ) computed in two different ways (since each of the corresponding squares is shared by two 3-faces), see Fig 4 Now, 4D consistency of the map means that T m (T x ij ) = T (T m x ij ) for any permutation i, j,, m of,, 3, 4 and for arbitrary initial data Theorem Symmetric discrete Darboux system (5) is 4D consistent Proof We will show that this is just a corollary of the geometry of spherical tetrahedra Consider a spherical tetrahedron in S 3 with the edge lengths l ij (connecting the vertices i and j) and dihedral angles α ij (at the edges (m), where, m is the complement of i, j in,, 3, 4) One can tae either α ij or l ij as independent data parametrizing the spherical tetrahedron The two sets of data are related by the cosine law for spherical tetrahedra, see section What seems to be less nown is that this cosine law for
SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 7 SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 77 T 3 T 4 x T 3 T 4 x x 34 x 34 x 3 x 3 x3 x3 T 3 x 4 T 3 x 4 T TT 3 x 34 x 4 T T xt 34 x 34 xx x 4 x 4 T TT 4 x 43 x 3 x x 4 4 T T 4 x T 3 4 x 3 Figure 4 Initial data and results of two-stage application of a 4D consis- Figure tent map 4 Initial with fields data onand -faces results of two-stage application of a 4D consistent map with fields on -faces sphericaltetrahedra triangles, ifcan thebe combinatorics derived via of a two-fold this two-fold application application of the is incosine agreement law for with the spherical spherical Fig 4 Intriangles, other words, if the if the a spherical combinatorics tetrahedron of of thiscombinatorially two-fold application corresponds is is in in agreement to agreement a 4D cube with with Fig Fig 4 In 4 other In other words, words, a a spherical tetrahedron combinatorially corresponds to to a 4D a 4D cube cube α (m) i α j α ij α (m) ij α α (m) i j α (m) j i α () ij α () ij α (i) j ij α (j) i ij α (j) i Figure 5 Left: Sides and angles in the spherical triangle L(m) Right: Sides Figure and 5 angles Left: in sides the spherical and angles triangle in the (ij) spherical triangle L(m); right: sides and angles in the spherical triangle (ij) To show this, we start with the dihedral angles α ij as initial data From these, we canto Figure determine show this, 5 allleft: we the start planar sides withand angles the angles dihedral α () ij (the inangles the latter spherical α ij isas the initial notation triangle data for From L(m); thethese, planar right: weangle can opposite determine sidesto and all the the angles side planar (ij) in in angles the the spherical α () ij (the latter triangle is the (ij) notation withfor the the vertices planari, angle j,, see opposite Fig 5 (right)) to the side Indeed, (ij) indraw the spherical a great -sphere triangle with S m the S 3 vertices with thei, pole j,, at seethe Fig vertex ) Indeed, m It draw will cut a To the great show tetrahedron -sphere this, we Salong m start Sa 3 with spherical the pole dihedral triangle at the which angles vertex is αcalled m ij It aswill L(m), initial cut data the tetrahedron lin From of the these, vertex alongwe am can (see determine spherical []) all triangle This triangle the planar which angles is has, called as α () one L(m), easily ij (the latter thesees, linthe is the of the angles notation vertex α ij for m, α(see i and the [3]) α planar j This and angle triangle the side opposite lengths to the has, side asα (ij) one (m), in easily α (m) i and the sees, α spherical the (m) j, angles see Fig triangle α ij 5, (left) with α i and Therefore, the vertices α j and the we find: i, j, side, see lengths Fig ) α (m) ij Indeed,, α (m) i and draw a greatα -sphere (m) j, see Fig S m 3 (left) S 3 with Therefore, ( the (m) pole we at find: α the vertex m It will cut the tetrahedron along a ij, α (m) i spherical triangle which is called (α (m) L(m),, ) α(m) j = Φ(αij, α i, α j ) the lin of the vertex m (see [3]) This triangle In the above combinatorial interpretation, ij, α i has, as one easily sees, the angles α ij,, α(m) j α )=Φ(α the dihedral ij, α i, angles α j ) α i and α j and the side ij are assigned lengths α (m) to six ij, α (m) - faces of the 4D cube, adjacent to one of its vertices The four 3-faces (3D cubes) adjacent i and, see Fig 3 (left) Therefore, we find: α (m) j α (m) i α j α ij (α (m) ij α (m) ij α i i α (i) j, α (m) i, α(m) j )=Φ(α ij, α i, α j ) j j
8 MATTEO PETRERA AND YURI B SURIS to this vertex, correspond combinatorially to the four spherical triangles L(m) The result of application of our 3D system to these four 3D cubes can be written as α (m) ij =T α ij At the second stage, we apply the cosine law to the four faces (i, j, ) of the tetrahedron, see Fig 5 (right) In the above combinatorial interpretation, these four spherical triangles correspond to the four 3-faces of the cube adjacent to its opposite vertex In particular, the side lengths l ij are assigned to the six -faces of the 4D cube, adjacent to this opposite vertex We find: (l ij, l i, l j ) = Φ ( α () ij, α(j) i, α(i) j ) This can be written as l ij = T α () ij = T (T m α ij ) Clearly, each edge is shared by two -faces, so we get two answers for each l ij Of course, they must coincide as definite geometric quantities in a given spherical tetrahedron The proof of Theorem given above is completely free of computations and it is based on the existence of a spherical tetrahedron with given dihedral angles We mention that very similar arguments were used in [] in the opposite direction: the coincidence of the two answers for each l ij was demonstrated there by a direct computation and used to provide a geometric realization for a given angle structure Remar In the combinatorial interpretation adopted in the present section, the map φ : α l is interpreted as solving a spherical triangle by its three planar angles, and combinatorially corresponds to one (directed) long diagonal of the cube The inverse map φ : l α would be interpreted as solving a spherical triangle by its three sides and would combinatorially correspond to the opposite orientation of this same long diagonal It is easy to realize that other long diagonals of the cube would correspond to solving a spherical triangle by two sides and a planar angle enclosed by them, or by a side and two adjacent planar angles The corresponding maps could be called companion maps to φ It is well nown that companion maps for a multidimensionally consistent map with fields on -faces satisfy the so called functional tetrahedron equation, see, eg, [, Sect 66] We come to a conclusion that the cosine law for spherical triangles delivers an example of solutions of the functional tetrahedron equation According to a private communication by R Kashaev, he was aware of this example since many years 4 Spherical triangles and discrete time Euler top In the present section, it will be more convenient to change notation for the edges and angles of a spherical triangle to (l, l, l 3 ) for (l 3, l 3, l ), resp to (α, α, α 3 ) for (α 3, α 3, α ) We study a dynamical system generated by iterations of the map φ from (6), resp its inverse map φ from (7), which now tae the form φ : y i = x i + x j x, φ : x i = x j x y i y j y (8) yj y Due to homogeneity, one can insert a small parameter ε into the map φ by a scaling transformation x i εx i, y i εy i leading to y i = x i + εx j x (9) ε xj ε x
SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 9 Geometrically, small values of ε in (9) correspond to small values of x i = cos α i in (8), ie, to spherical triangles with all three angles (and all three sides) being close to π/ The form (9) of the map φ maes it obvious that it can be considered as a time discretization of the famous Euler top, described by the following system of differential equations: ẋ i = x j x (0) The latter system can be brought to the form ẋ i = c i x j x, involving arbitrary real parameters c i, by a simple scaling transformation x i b i x i, where, however, some of b i may be imaginary, depending on the signs of c i It is well-nown that system (0) is integrable in the Liouville-Arnold sense and can be solved in terms of elliptic functions, see for instance [6] In particular, it is bi-hamiltonian and admits integrals of motion I ij = x i x j, two of which are functionally independent We will show that discretization (9) of the Euler top inherits integrability Theorem The map φ admits integrals of motion E ij = x i, () x j two of which are functionally independent It has an invariant volume form given by ω = dx dx dx 3, ϕ(x) where ϕ(x) is any of the functions ( x j)( x ) or ( x i ) (the quotient of any two such functions is an integral of motion) Proof The first statement (about integrals of motion) follows immediately from the sine law (8) As for the second statement (about an invariant volume form), a direct computation leads to the following expressions: and y i x j = y i x i = x + x i x j = ( x j )3 x, x j x y x j x x i x j The second factor in the last formula can be interpreted as a conjugation of the Jacobi matrix y/ x by the diagonal matrix with entries x i, and is therefore irrelevant for the computation of its determinant, which is thus equal to ( ) y y 3 y det = x ( x )( x )( x y 3 y 3) y y = d ( x )( x )( x 3) According to the sine law (8), the latter expression is equal to ( ) y det = ( y j )( y ) x ( x j )( x ), which finishes the proof
0 MATTEO PETRERA AND YURI B SURIS Corollary The map φ is Poisson with respect to Poisson bracets (bi-vector fields) obtained by contraction of the tri-vector field ϕ(x) 3 with any of the one-forms de ij (x) Explicitly, a family of compatible invariant Poisson bracets is given by {x i, x j } = C i x ( x j) C j x ( x i ), where C, C, C 3 are arbitrary constants Thus, the map φ is bi-hamiltonian and is completely integrable in the Liouville-Arnold sense A somewhat unpleasant feature of the map φ is its non-rational character, due to the square roots in the denominators A striing and unexpected fact is that the second iterate of φ is a rational (actually, a birational) map, which appeared previously in the literature on a different occasion Theorem 3 On the set τ φ (τ), where the second iterate of φ is defined, this second iterate φ : x x is a birational map given by x i = x i + x j x + x i ( x i + x j + x ) (x + x + x 3) x x x 3 () These functions are the (unique) solutions of the following system of equations: x i x i = x j x + x j x, (3) which constitute the Hirota-Kimura discretization of the Euler top [5] Proof We let φ(x) = y, φ(y) = x According to formulas (8), we have x i = y i y j y, x i = yj y y i + y j y (4) yj y Using (4), we find, on one hand, x i x i = y j y, yj y and, on the other hand, x j x + x j x = (y j + y y i )(y y i y j ) + (y j y y i )(y + y i y j ) ( yi ) = yj y y j y yj y This proves equations (3) for φ These equations are linear with respect to x, x, x 3, and thus they can be solved to give the explicit birational map () Actually, the integrability attributes of the map φ given in Theorem and Corollary were found previously for the rational map φ : integrals in [5] and the invariant volume form along with the bi-hamiltonian structure in [5] We refer also to [6, 3, 4] for our recent study on Hirota-Kimura-type discretizations While the map Φ relates the data within one spherical triangle, a natural interpretation of the map Φ is given by the following map between spherical triangles Definition Let be a spherical triangle with the angles α i and the sides l i switched triangle = S( ) is the spherical triangle with inner angles The α i = l i (5)
SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS (whenever such a triangle exists) The switched triangle is defined up to a transformation from O(3) Clearly, S( ) exists if and only if l T Φ(T ), ie, if and only if α Φ (T ) T In this case, we have l = Φ (α) Thus, the quantities x i = cos l i are related to x i = cos α i by the Hirota-Kimura discrete time Euler top equations () or (3) By abusing notations slightly, we write ( α, l) = S(α, l) By induction, we have: S n (α, l) = (Φ n (α), Φ n+ (α)) Therefore, S n ( ) is defined iff α T Φ (T ) Φ n (T ) For any n, the latter domain contains an open neighborhood of the point (π/, π/, π/) When taling about the switched triangles to our colleagues, we have been often ased: But isn t this just the polar triangle? The answer is: Of course not, since the switch generates a non-trivial dynamical system, while the polarity is an involution Nevertheless, it turns out that these two transformations of spherical triangles are closely related We denote the angles and the sides of the polar triangle = P ( ) by (α, l ) = P (α, l) They are given by α i = π l i, l i = π α i (6) Clearly, in the latter interpretation, P maps T T to itself Comparing formulas (5) and (6), we are led to the following definition Definition Let be a spherical triangle with the angles α i and the sides l i The side flip = F ( ) is the spherical triangle with sides l i = π l i (7) (whenever such a triangle exists) The side flip F ( ) is defined up to a transformation from O(3) The map F is an involution (when considered on equivalence classes modulo O(3)) An immediate consequence of the above definition is the following claim Theorem 4 The switch of spherical triangles, considered as a map on equivalence classes modulo O(3), can be represented as a composition of two involutions, namely as the side flip followed by the polarity transformation: S = P F The angles α i of the side flipped triangle = F ( ) are given by α = π Φ (α) The quantities x i = cos α i and x i = cos α i satisfy or, explicitly, x i + x i = x j x + x j x, (8) x i = x i + x j x + x i ( x i + x j + x ) (x + x + x 3) x x x 3 (9) The map f : x x defined by (9) is a birational involution It is very easy to find integrals of this map Indeed, from the defining relation (7) yielding cos l i = cos l i and the cosine law (3) there follows: cos α i + cos α j cos α sin α j sin α = cos α i + cos α j cos α sin α j sin α
MATTEO PETRERA AND YURI B SURIS As a consequence, one gets three functionally independent conserved quantities of the map (9): cos l i = (x i + x j x ) ( x j )( x ) sin l i = d ( x j )( x ), with the quantity d defined by (9) Note that the quotients E ij = sin l i / sin l j depend only on x i rather than on x i, and are therefore integrals of motion for the map φ = i f (these are nothing but functions ()) Map (9) appeared in an old paper [7] by H Jonas Actually, he considered a closely related angle flip of spherical triangles, defined by the relation α i = π α i, and proved that the cosines of the side lengths of the spherical triangles, are related by x i + x i + x j x + x j x = 0 (30) The latter map is conjugated to (8) via i : x x Jonas results include integrals of the map (30) and its solution in terms of elliptic functions 5 Spherical tetrahedra and discretization of two linearly coupled Euler tops We now consider the algebraic form of the cosine law for spherical tetrahedra as a dynamical system Recall that it is given by y ij = g ij, (3) gii g jj with the expressions for the polynomials g ij given in ( ) This map, which we denote by ψ, can be considered as a time discretization of the following system of ordinary differential equations: ẋ ij = x i x j + x im x jm (3) Here, recall, (i, j,, m) stands for an arbitrary permutation of (,, 3, 4) This sixdimensional system turns out to consist of two linearly coupled copies of the Euler top Namely, each group of three variables p ij = x ij + x m and q ij = x ij x m with ij =, 3, 3, satisfies its own copy of the Euler top equations: ṗ ij = p i p j, q ij = q i q j (33) As a consequence, equations (3) are integrable in the Liouville-Arnold sense, with four independent integrals of motion which can be chosen as p p 3 = (x + x 34 ) (x 3 + x 4 ), p 3 p 3 = (x 3 + x 4 ) (x 3 + x 4 ), q q 3 = (x x 34 ) (x 3 x 4 ), q 3 q 3 = (x 3 x 4 ) (x 3 x 4 ), or, alternatively, as x + x 34 x 3 x 4, x 3 + x 4 x 3 x 4, x x 34 x 3 x 4, x 3 x 4 x 3 x 4 Theorem 5 The map ψ admits four independent integrals of motion ( x 3)( x 4) ( x )( x 34), ( x 3)( x 4) ( x )( x 34), x x 34 x 3 x 4 ( x )( x 34), x 3 x 4 x 3 x 4 ( x )( x 34)
SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 3 It has an invariant volume form given by ω = dx dx 3 dx 3 dx 4 dx 4 dx 34, ϕ(x) where ϕ(x) is any of the functions ( x ij) 5/ ( x m )5/ Proof The claim about integrals of motion follows directly from the sine laws (3 4) Let us turn to the proof of the claim about the invariant volume form Our aim is to compute the determinant of the Jacobi matrix ( y ij / x i j ) We use the following ordering of the index pairs:, 3, 3, 4, 4, 34 The entries of the Jacobi matrix are given by the following expressions: y ij x ij = x m gii g jj, (34) y ij = x m g y j x i gii g jj y ij = x m p ij m x m gii g jj, (35) g jj g g mm, (36) g ii g jj where p ij m = y iy jm + y im y j y ij y i y im y ij y j y jm yij Indeed, the first two expressions, (34 35), are easily obtained by a direct computation, while the third one, (36), is a result of [] The last factors in the above expressions can be interpreted as a conjugation of the Jacobi matrix ( y ij / x i j ) by a diagonal matrix ( ) g33 g33 g44 g44 g33 g 44 diag,,,,,, g g g g g g and are therefore inessential for the computation of the determinant We come to the following result: y 3 y 3 y 4 y 4 p 34 ( ) y y <m det = ( x m ) 3 y y 34 p 4 3 y 4 y 3 y p 4 3 y 34 y 4 x (g g g 33 g 44 ) 3/ y 4 y 34 p 3 4 y y 3 y 4 p 3 4 y 34 y y 3 p 34 y 4 y 4 y 3 y 3 A direct computation with a symbolic manipulation system lie Maple leads to the following factorized expression for the determinant on the right-hand side: y 3 y 3 y 4 y 4 p 34 y 3 y y 34 p 4 3 y 4 y 3 y p 4 3 y 34 y 4 y 4 y 34 p 3 = g g g 33g 44 4 y y 3 y 4 p 3 4 y 34 y y i<j ( y ij ) d 3 p 34 y 4 y 4 y 3 y 3 Therefore, ( ) y <m det = ( x m ) (g g g 33g 44) 3/ d x i<j ( y ij ) (g g g 33 g 44 ) 3/ γ
4 MATTEO PETRERA AND YURI B SURIS At the next step, we use the following result (proved below): (g g g 33g 44) 3/ i<j (g g g 33 g 44 ) = ( y ij) 3/ <m ( (37) x m ) As a consequence, we find with the help of (3): ( ) y i<j det = ( ( ) y ij) d γ 6 ( ) d γ 5 x <m ( x m ) = =, γ d γ d and one more reference to relations (3) proves the claim about the invariant volume form To finish the proof we need to prove formula (37) Comparing expression (9) with (), we see that g ii is the analogue of the quantity d for the spherical triangle with the angles α j, α jm, α m, ie, for the triangle L(i) Then, from (0), we obtain: g / ii = sin α (i) j sin α jm sin α m Multiplying all three such expressions for L(i), we find: g 3/ ii = sin α (i) j sin α(i) jm sin α(i) m sin α j sin α jm sin α m Similarly, g ii is the analogue of the quantity d for the spherical triangle with the sides l j, l jm, l m, ie, for the triangle (j,, m) From (0), we obtain: (g ii) / = sin α (m) j sin l jm sin l m Multiplying all three such expressions for (j,, m), we find: As a consequence, (g ii) 3/ = sin α (m) j sin α() jm (g g g 33 g 44 ) 3/ = (g g g 33g 44) 3/ = sin α(j) m sin l j sin l jm sin l m ( 4 ) sin α (i) j sin 4 α j, i= j< ( 4 m= j< sin α (m) j j< ) sin 4 l j Indeed, each dihedral angle α j participates as a planar angle of two triangles L(i) and L(m), and each side l j is shared by two triangles (j,, m) and (i, j, ) This proves (37) Remar According to the famous Schläfli formula (see, eg, [, 8]), if V = V (α) is the volume of the spherical tetrahedron with the dihedral angles α ij, then dv = ij l mdα ij Thus, our proof of Theorem 5 actually leads to an evaluation of the determinant of the Hesse matrix of the volume function, det( V/ α ij α i j ) We are not aware of such an evaluation in the literature Remar The reduction of the system of ordinary differential equations (3) to two copies of the Euler top (33) yields that the original system admits two invariant 3-forms, dp dp 3 dp 3 and dq dq 3 dq 3 We conjecture that the discretization ψ also admits an invariant 3-form This property is stronger than existence of an invariant volume form established in Theorem 5 If true, this 3-form (or, better, the corresponding 3-vector field) j<
SPHERICAL GEOMETRY AND INTEGRABLE SYSTEMS 5 would yield, by contraction with gradients of integrals of motion, a family of compatible Poisson structures for ψ Acnowledgment The authors are partly supported by DFG (Deutsche Forschungsgemeinschaft) in the frame of Sonderforschungsbereich/Transregio 09 Discretization in Geometry and Dynamics References [] AI Bobeno, YuB Suris, Discrete differential geometry Integrable structure, Graduate Studies in Mathematics, Vol 98 AMS, 008 [] DA Derevnin, AD Mednyh, MG Pashevich, On the volume of symmetric tetrahedron, Sib Math Jour, 45/5, 840 848, 004 [3] LV Bogdanov, BG Konopelcheno, Lattice and q-difference Darboux-Zaharov-Manaov systems via -dressing method, Jour Phys A, 8, L73 Ll78, 995 [4] A Doliwa, The C-(symmetric) quadrilateral lattice, its transformations and the algebro-geometric construction, Jour Geom Phys, 60, 690 707, 00 [5] R Hirota, K Kimura, Discretization of the Euler top, Jour Phys Soc Japan, 69, 67 630, 000 [6] ANW Hone, M Petrera, Three-dimensional discrete systems of Hirota-Kimura type and deformed Lie-Poisson algebras, Jour Geom Mech, /, 55 85, 009 [7] H Jonas, Deutung einer birationalen Raumtransformation im Bereiche der sphärischen Trigonometrie, Math Nach, 6, 303 34, 95 [8] B G Konopelcheno, W K Schief, Three-dimensional integrable lattices in Euclidean spaces: conjugacy and orthogonality, Proc R Soc Lond A, 454, 3075 304 (998) [9] AA Kolpaov, AD Mednyh, MG Pashevich, A volume formula for Z -symmetric spherical tetrahedra, Sib Math Jour, 5/3, 456 470, 0 [0] F Luo, On a problem of Fenchel, Geom Dedicata, 64, 77 8, 997 [] F Luo, Volume and angle structures on 3-manifolds, Asian Jour Math, /4, 555 566, 007 [] F Luo, 3-Dimensional Schläfli formula and its generalization, Comm Contemp Math, 0, 835 84, 008 [3] M Petrera, A Pfadler, YuB Suris, On integrability of Hirota-Kimura type discretizations Experimental study of the discrete Clebsch system, Exp Math, 8/, 3 47, 009 [4] M Petrera, A Pfadler, YuB Suris, On integrability of Hirota-Kimura type discretizations, Reg Chaot Dyn, 6/3-4, 45 89, 0 [5] M Petrera, YuB Suris, On the Hamiltonian structure of Hirota-Kimura discretization of the Euler top, Math Nachr, 83/, 654 663, 0 [6] AG Reyman, MA Semenov-Tian-Shansy, Group theoretical methods in the theory of finitedimensional integrable systems, in: Dynamical Systems VII, Springer, Berlin, 994 [7] W K Schief, Lattice geometry of the giscrete Darboux, KP, BKP and CKP equations Menelaus and Carnot Theorems, Jour Nonlin Math Phys, 0/, 94 08, 003 [8] EB Vinberg (Ed), Geometry II, Encyclopaedia of mathematical sciences, Vol 9, Springer-Verlag, 993