arxiv: v2 [math-ph] 13 Feb 2013

Transcription

1 FACTORIZATIONS OF RATIONAL MATRIX FUNCTIONS WITH APPLICATION TO DISCRETE ISOMONODROMIC TRANSFORMATIONS AND DIFFERENCE PAINLEVÉ EQUATIONS ANTON DZHAMAY SCHOOL OF MATHEMATICAL SCIENCES UNIVERSITY OF NORTHERN COLORADO GREELEY, CO arxiv: v2 [math-ph] 3 Feb 203 Abstract. We study factorizations of rational matrix functions with simple poles on the Riemann sphere. For the quadratic case two poles) we show, using multiplicative representations of such matrix functions, that a good coordinate system on this space is given by a mix of residue eigenvectors of the matrix and its inverse. Our approach is motivated by the theory of discrete isomonodromic transformations and their relationship with difference Painlevé equations. In particular, in these coordinates, basic isomonodromic transformations take the form of the discrete Euler-Lagrange equations. Secondly we show that dpv equations, previously obtained in this context by D. Arinkin and A. Borodin, can be understood as simple relationships between the residues of such matrices and their inverses.. Introduction In this paper we are concerned with various ways of introducing coordinates on the space of rational matrix functions Lz) with simple poles) on the Riemann sphere. Examples like this play an important role in various applications, like Yang-Baxter maps and matrix solitons [, 2], Lax equations and isomonodromy transformations on algebraic curves, [3, 4], discrete integrable systems, [5, 6], and others. Our main interest is related to the study of discrete isomonodromic transformations and their relationship with discrete Painlevé equations. The theory of discrete isomonodromic transformations was developed by A. Borodin for polynomial matrices in [7] and adapted to the notion of local monodromy using rational matrix functions by I. Krichever in [8]. In [9] D. Arinkin and A. Borodin used the theory of d-connections on vector bundles to explained the appearance of difference Painlevé equations, considered from the geometric point of view of H. Sakai [0], in the theory of discrete isomonodromic transformations and later, in [], they introduced the notion of the τ-function of such transformation. These τ-functions appear as the gap probabilities in the discrete probabilistic models of random matrix type. I. Krichever conjectured that discrete isomonodromic transformations can be written in the Lagrangian form and that they should be related to the universal symplectic form of Krichever-Phong, [2]. In [3], using the methods of [5], we verified this conjecture for the quadratic two-pole) case by using the multiplicative coordinates on the space of these matrix functions and finding the explicit formula for the Lagrangian. Recently F. Soloviev showed that our Lagrangian symplectic form coincides with the reduction of the quadratic symplectic form of Krichever-Phong to certain symplectic leaves [4]. Unfortunately, the multiplicative coordinates are not easily obtained from other characteristic properties of Lz), such as its residue matrices. This is an obstacle to the generalization of our results to the higher-order case. In this paper we argue that one way around this difficulty is to consider, in addition to the residues of Lz), the residues of the inverse matrix L z). Then good coordinates on the space of Lz), again in the quadratic case, are given by half of the residues data of Lz) and half of the residue data of L z). Such residues also explain the almost symmetric form of the expressions for the multiplicative coordinates of Lz) and allow us to recognize dpv equations as simple relations between the residues of Lz) ± and the residues of the transformed matrices Lz) ±. The paper is organized as follows. In Section 2 we give a short overview of the additive representation of Lz). In Section 3 we study the multiplicative representation of Lz) and establish the relationship between 2000 Mathematics Subject Classification. 39A0, 37K20.

2 the eigenvectors of the residues of Lz) ± and the eigenvectors of the left and right-divisors of Lz). We also significantly simplify many of the arguments and formulas of [3]. Finally, in Section 4 we restrict our attention to the rank-two case. It this case it is possible to introduce the so-called spectral coordinates on the space of Lz) in such a way that the equations relating the spectral coordinates of Lz) and Lz) are precisely the difference Painlevé equations. Our main observation here is the following. In addition to spectral coordinates, divisor of zeroes and poles, and some asymptotic behavior at infinity, the entries of the matrix Lz) also depend on a choice of a gauge with respect the action of the group of constant nondegenerate diagonal matrices. Understanding the change of this gauge from Lz) to L z) and Lz) allows us to easily obtain the expressions of the multiplicative representation of Lz) in the spectral coordinates and the dpv equations. 2. Additive form of Rational Matrix Functions Let Lz) be a rational matrix function on the Riemann sphere, rank Lz) = m, satisfying the following general conditions. First, we require that there exists a normalization point z 0 at which Lz) is regular, lim z z0 Lz) = L 0, det L 0 0, and all eigenvalues of L 0 are distinct. We can then do a gauge transformation to make L 0 diagonal, thus reducing the global gauge group to the group of diagonal matrices. Without any loss of generality we can assume that z 0 =, and so 2.) lim z Lz) = L 0 = diag{ρ,..., ρ m }. Second, we impose the following conditions on the pole structure of Lz) and its inverse Mz) = Lz). We require that Lz) is holomorphic except for simple poles at the points z,..., z k, Mz) is holomorphic except for simple poles at the points ζ,..., ζ k, all z i and ζ j are distinct, and the determinant det Lz) has also only simple poles at z i and simple zeroes at ζ j. These conditions mean that the residues L i = res zi Lz) and M j = res ζj Mz) the negative sign here is for future convenience) are matrices of rank one. Using the symbol to indicate a row vector, we have: 2.2) 2.3) 2.4) Lz) = L 0 + k i= det Lz) = ρ ρ m k i= z ζ i) k j= z z j), Lz) = Mz) = M 0 k i=j L i z z i, where L 0 = diag{ρ,..., ρ m } and L i = a i b i, M j z ζ j, where M 0 = L 0, M j = c j d j, We call the above representations of Lz) and Mz) additive representations and the vectors a i, b i resp. c i, d i ) additive eigenvectors of Lz) resp. Mz)). Note that in terms of Lz) the vectors c j resp. d j ) can be characterized as the left resp. right) null-vectors of Lζ j ), and the similar statement is true for Mz i ) and a i, b i. By the divisor of Lz) we mean the divisor D of its determinant, D = i z i i ζ i. We denote the space of matrices Lz) satisfying the conditions 2.2) 2.4) by M D r. An important question is how to choose a good coordinate system on the space M D r. For example, given L i, we can determine a i and b i only up to a common scaling factor. This factor has to be adjusted to ensure that det Lz) has zeros at ζ i, which a complicated condition on trl i ). Same problem is present for the collection {c i, d i }. Some insight for a good choice of coordinates is provided by the study of the isospectral and isomonodromic transformations, which suggests that half of the coordinates should be taken from the residues of Lz) and half from the residues of Mz). In fact, in the quadratic case, we have the following result. Theorem 2.. When Lz) has k = 2 poles, the vectors c 2, d ; a 2, b ), considered up to rescaling i.e., as points in P r ), are coordinates on the space M D r. To recover L ± z), consider the function 2.5) Lx 2, x ), y 2, y )) = z 2 z ) logx L 0x 2 ) + z ζ 2 ) logy x 2) + ζ 2 ζ ) logy L 0 y 2) + ζ z 2 ) logx y 2). 2

3 Then 2.6) 2.7) a = L c = x c 2, d ), a 2, b )); b 2 = L x 2 c 2, d ), a 2, b )); L c 2, d ), a 2, b )); y d 2 = L y 2 c 2, d ), a 2, b )). The proof of this theorem is based on the multiplicative representations of Lz), which we consider next. 3. Multiplicative form of Rational Matrix Functions 3.. Elementary Divisors. To change from the additive to a multiplicative representation, we first need to define the building blocks for it. We call such building blocks the elementary divisors. Definition 3.. An elementary divisor is a rational m m-matrix function Bz) on the Riemann sphere of the form 3.) Bz) = I + G, where G = pq is a matrix of rank one. A direct calculation establishes the following elementary facts. Lemma 3.2. Let Bz) be an elementary divisor. Then i) det Bz) = z ζ 0 )/ ), where ζ 0 = z 0 q p; ii) Bz) = I G/z ζ 0 ). In fact, for us it will be more convenient to fix the points ζ 0 and z 0 on CP. Thus, we say that a pair ζ 0, z 0 ) corresponds to an elementary divisor Bz) with the determinant det Bz) = z ζ 0 )/ ). Any generic matrix Bz) with such determinant and normalized by the condition Bz) I as z is of the form Bz) = I + G/ ) with G of rank one and det Bζ 0 ) = 0 more carefully, Bz) should be regular at z 0 and Bz) should be regular at ζ 0 ). Thus, Bζ 0 ) has a left null-vector q and a right null-vector p. We then immediately get that G = pq, where we need to normalize the vectors p and q so that q p = z 0 ζ 0. In a more invariant form this can be written as 3.2) Bz) = I + z 0 ζ 0 pq q p. From that point of view, the formula for the inverse matrix follows from the vanishing of the residue of the identity Bz)Bz) = I = Bz) Bz) at z 0, 3.3) Bz) = I + ζ 0 z 0 z ζ 0 pq q p. The following easy properties of elementary divisors are very useful for what follows. Lemma 3.3. Let 3.4) Bz) = I + z 0 ζ 0 pq q p. Then i) ) z ζ0 3.5) Bz)p = p and q Bz) = z ζ0 ) q. ii) Suppose that at some point z we have Bz )w = v for some vectors v and w. Then 3.6) Bz) = I + ) z 0 z ) wq q w + z ζ 0 ) vq q v i.e., we can determine p from v, w, and q ). Similarly, if w Bz ) = v, 3.7) Bz) = I + ) z 0 z ) pw w p + z ζ 0 ) pv v. p 3

4 Proof. Part i) is immediate. To establish the first formula in part ii), we solve Bz )w = v for p/q p) and then use i): p 3.8) q p = z 0 z w z 0 ζ 0 q w + z z 0 z 0 ζ 0 The second formula is proved in the similar way. v q w = z 0 z z 0 ζ 0 w q w + z ζ 0 z 0 ζ 0 v q v. Finally, we need the following notation. Notation: We define a twisting of an elementary divisor Bz) by a constant non-degenerate matrix A to be a new elementary divisor A Bz) such that A Bz)A = ABz), i.e., 3.9) A Bz) = ABz)A = I + z 0 ζ 0 Ap)q A ) q A )Ap) Factors and Divisors. We begin with the following important remark. For the additive representation of Lz), the ordering of zeroes and poles of det Lz) is not important, but for any multiplicative representation choosing such an ordering is crucial. Thus, from now on our labeling will reflect the fact that ζ s, z s )-pair corresponds to some elementary divisor in Lz). There are two ways to look at the multiplicative structure of Lz) we can look at factors or at divisors. Definition 3.4. We say that elementary divisors C s z) = I + H s /z z s ) corresponding to pairs ζ s, z s ) with H s = m s n s are the factors of Lz) if 3.0) Lz) = L 0 C z) C k z). Definition 3.5. We say that elementary divisors B r sz) resp. B l sz)) corresponding to pairs ζ s, z s ) are right resp. left) divisors of Lz) if Lz) = L r sz)b r sz) resp. Lz) = B l sz)l l sz)) where L r sz) resp. L l sz)) is regular at z s. The main advantage of the divisors as opposed to the factors is that they can be written explicitly in terms of the residues of Lz) and Mz). Note that C k z) = B r k z) and L0 C z) = B l z). In particular, in the quadratic case k = 2, there is no essential difference between divisors and factors. Lemma 3.6. Let L s = res zs Lz) = a s b s and M s = res ζs Mz) = c s d s. Then 3.) 3.2) B r sz) = I + z s ζ s z z s c s b s b sc s B l sz) = I + z s ζ s z z s a s d s d sa s. Proof. To obtain the formula for the right divisor, we take the residue of Lz) = L r sz)b r sz) at z s to get b s q r s). Then we take the residue of Mz) = B r sz)) L r sz)) at ζ s to get c s p r s. The formula then follows. The expression for the left divisors is obtained in the same way Refactorization Transformations. As shown in [7], [8], isomonodromic transformations on the space of rational matrix functions Lz) have the form 3.3) Lz) Lz) = Rz + )Lz)Rz), where Rz) is a certain rational matrix function. Similarly, we can consider isospectral transformations 3.4) Lz) Lz) = Rz)Lz)Rz). In the isospectral case we have to choose Rz) in such a way that the singularity structure of Lz) is the same as Lz), and in the isomonodromic case we require that Rz) induces integral shifts of certain parameters corresponding to the linear difference system given by Lz) see [7], [9] for details). In particular, these requirements are satisfied if we take Rz) to be one of the divisors of Lz), which corresponds to changing the order of the factors in the multiplicative representation of Lz); in the isospectral case the divisor D is fixed and in the isomonodromic case D is shifted by an integer vector. In what follows we only consider 4

5 transformations of this particular type. We focus on the isospectral case, since the isomonodromic case is very similar. Let us now restrict our attention to the quadratic case and take Rz) = B r z): 3.5) Lz) = B l 2z)L 0 B r z) Lz) = B r z)b l 2z)L 0 = B l 2z)L 0 Br z). To study these transformation from the point of view of the discrete Euler-Lagrange equations it is necessary to construct a configuration space Q, a Lagrangian function L FQ Q), and a map η : Q Q M D r such that η maps the Lagrangian dynamics to the refactorization dynamics. That is, denoting by Q the previous and by Q the next point of the discrete dynamics, we want ηq, Q) = Lz) and ηq, Q) = Lz), where the shift map Φ : Q, Q) Q, Q) should satisfy the discrete Euler-Lagrange equations L 3.6) Y Q, Q) + L X Q, Q) = 0, see [5], [5] for the general description of this approach, and [3] for our specific situation. Writing 3.7) Lz) = B r z)b l 2z)L 0 = B l z)l 0 Br 2 z) and using the uniqueness of divisors, we see that B r z) = B l z) and B l 2z)L 0 = L 0 Br 2 z). From Lemma 3.6 it then follows that c = ã, b = d, a 2 = L 0c 2, and d 2 L 0 = b 2. Since we want to parametrize Lz) by Q, Q), we see that if we take Q = a 2, b ) as one half of our coordinates, the second half should be Q = a 2, b ) = L 0 c 2, d ). Now the proof of Theorem 2. is the same as the proof of Theorem 3. of [3], and itis sketched below. Proof. Theorem 2.) Taking the residue of 3.8) Lz) = B l z)l 0 B r 2z) = B l 2z)L 0 B r z) at the point z and comparing the row spaces of the resulting rank-one matrices gives d L 0B r 2z ) = b. Using Lemma 3.3ii) we can find the formula for b 2, 3.9) b 2 = z 2 z ) d L 0 d L + z ζ 2 ) b 0c 2 b c 2 = L x 2 c 2, d ), a 2, b )), where Lx 2, x ), y 2, y )) is given by 2.5). Performing similar calculations at the point z 2 for Lz) and the points ζ, ζ 2 for Mz) gives the rest of the formulas ) and completes the proof. In these coordinates Theorem 3. of [3] takes the following form here we improve the formulas from [3] by using L 0 instead of its root). Theorem 3.7. i) The map η : Q Q M D r is given by ηq, Q) = Lz) = I + z 2 ζ ) a 2b z z 2 b 3.20) I + z z a 2 z ζ 2 ) L 0 ã 2b b L 0 ã 2 )) + ζ ζ 2 ) a 2b L 0 b L 0 a L ζ 2 ζ ) L 0 a 2b 2 b 2 L 0 a 2 ii) The equations of motion Q, Q) = ΦQ, Q) of both the isospectral and isomonodromic dynamics in these coordinates are given by the discrete Euler-Lagrange equations with the Lagrangian function 3.2) LX, Y, t) = z 2 z t)) logx x 2) + z t) ζ 2 ) logy L 0 x 2) + ζ 2 ζ t)) logy L 0 y 2) + ζ t) z 2 ) logx y 2), where in the isomonodromic case z t) = z t, ζ t) = ζ t, and in the isospectral case z t) = z, ζ t) = ζ and LX, Y) is time-independent. 5 )) ;

6 4. Isomonodromic Transformations and dpv In [9], Arinkin and Borodin showed that for rank-two matrices isomonodromic transformations above, when written in a special coordinate system, are given by the difference Painlevé equations. In this section, choosing the dpv case as an example, we show that these equations appear explicitly as relations between the residues of L ± z) and L ± z). Similar computation for the q-pvi case was done earlier by Jimbo and Sakai, [6]. 4.. Spectral Coordinates. In the quadratic two-pole) case the space of the rank-two matrices Lz) satisfying the requirements ) can be described using different parameters. These parameters come in two groups. The first group, that we call the type of Lz), consists of the zeroes and poles of the determinant of Lz) and some asymptotic data at z =. The space of Lz) of the fixed type is two-dimensional, and the second group of parameters is a special coordinate system on this space, called the spectral coordinates. Expressing isomonodromy transformation in those coordinates gives rise to the difference Painlevé equations. Definition 4.. Let Lz) be a rational 2 2 matrix function on the Riemann sphere satisfying the following conditions: 4.) In addition, put Lz) = diag{ρ, ρ 2 } + L z z + L 2 z z 2, L i = a i b i, det Lz) = ρ ρ 2 z ζ )z ζ 2 ) z z )z z 2 ). 4.2) ρ k = L ), ρ 2 k 2 = L ) 22, µ = L ) 2, where L = res Lz) dz = L + L 2. We call ρ, ρ 2, ζ, ζ 2, z, z 2, k, k 2, µ) the type of Lz). These parameters are not independent, since 4.3) k + k 2 = tr L 0 L = z ζ ) + z 2 ζ 2 ). The conditions on ρ i k i correspond to fixing the formal type of the solution of the difference equation at infinity, and the choice of µ corresponds to fixing the gauge under the global action by constant nondegenerate diagonal matrices. Definition 4.2. The spectral coordinates γ, π) are defined by the conditions that Lγ) 2 = 0 and therefore Lz) 2 = µz γ) z z )z z 2) ); π = γ z) γ ζ 2) Lγ). The normalization conditions in this definition are chosen to match the formulas in [9] and []. Notation: To find the expression of Lz) and Mz) in the spectral coordinates, it is convenient to introduce the notation ϕa, b) = πγ a) ρ γ b). Lemma 4.3. Let Lz) be a rational 2 2 matrix function on the Riemann sphere that has the type ρ, ρ 2, ζ, ζ 2, z, z 2, k, k 2, µ). Then, in spectral coordinates, the residues of the matrix Lz) are given by [ )] L = µ γ z µ ρ k γ z2) γ z ϕζ 2, z ) [ 4.4) µ ρ2 k 2 + ρ2 π z 2 z ϕz 2, ζ ) )], L 2 = µ γ z [ 2 µ ρ ] k ϕζ 2, z )) [ )] 4.5) z z 2 µ ρ 2 k 2 + ρ2γ z) πγ z ϕz 2) 2, ζ ). [ ] ai [ ] Proof. Let L i = α i bi. Then 4.6) Lz) 2 = α + α 2 µz γ) = z z z z 2 z z )z z 2 ), and so α = µγ z )/z 2 z ), α 2 = µγ z 2 )z z 2 ), and α + α 2 = µ. 6

7 The normalization at infinity and the definition of π, 4.7) 4.8) immediately give 4.9) 4.0) ρ k = α a + α 2 a 2 = µa + α 2 a 2 a ) = α a a 2 ) + µa 2, Lγ) = π γ ζ 2) γ z ) = ρ + α a a 2 ) γ z ) a = µ ρ k α 2 a 2 a )) = µ a 2 = µ ρ k ϕζ 2, z )). = ρ + α 2a 2 a ), γ z 2 ) ρ k γ z ) 2) γ z ) ϕζ 2, z ), γ ζ ) Using the equation Lγ) Lγ) 22 = det Lγ) we get Lγ) 22 = ρρ2 π γ z. This, and the condition ρ 2) 2k 2 = α b + α 2 b 2, allows us to find the expressions for b, b 2 in exactly the same way. Corollary 4.4. In the same gauge, the residues M i of the inverse matrix 4.) Mz) = Lz) = diag{/ρ, /ρ 2 } M z ζ M 2 z ζ 2 are given by 4.2) 4.3) M = µ [ γ ζ µ ρ ρ 2 ζ 2 ζ M 2 = µ ρ ρ 2 γ ζ 2 ζ ζ 2 ρ2 k ρ2 π ϕζ 2, z ) ) [ µ ρ 2 k ρ2γ ζ) ] [ πγ ζ 2) ϕζ 2, z ) )] µ ρ k 2 + γ ζ2 γ ζ ϕz 2, ζ ), )] [ µ ρ k 2 + ϕz 2, ζ )) ]. Proof. Since Mz) has the same form as Lz), we only have to determine the type and spectral coordinates of Mz) in terms of those of Lz), Lz): z z 2 ζ ζ 2 ρ ρ 2 k k 2 µ γ π Mz): ζ ζ 2 z z 2 ρ ρ 2 k k 2 µ ρ ρ 2 γ γ z )γ ζ ) πγ z 2)γ ζ 2) and take the negative sign in the definitions of M i into account. Note that in computing the type of Mz) we used the equation M = L 0 L L 0, which follows from the condition res Lz)Mz) = Difference Painlevé V. Consider now the isomonodromy transformation given by Rz) = B r z): 4.4) Lz) = B l 2z)L 0 B r z) Lz) = B r z + )B l 2z)L 0 = B l z)l 0 Br 2 z). Theorem 4.5. The type and spectral coordinates of Lz) in terms of those of Lz) are given by 4.5) Lz): z z 2 ζ ζ 2 ρ ρ 2 k k 2 µ γ π Lz): z = z z 2 = z 2 ζ = ζ ζ2 = ζ 2 ρ ρ 2 k k 2 µ γ π, where 4.6) 4.7) 4.8) µ = µ ρ π ρ 2 ) ρ 2 π ρ ) γ + γ = z 2 + ζ 2 + ρ k z + ζ 2 ) π ρ + ρ 2k 2 z + ζ 2 + ) π ρ 2 ππ = ρ ρ 2 γ z )γ ζ ) γ z 2 )γ ζ 2 ) Equations ) are the difference Painlevé V equations of Sakai s hierarchy [0], first obtained in this setting in [9] Theorem B). 7

8 Proof. First note that 4.9) L = G l 2L 0 + L 0 G r, L = G r + G l 2)L 0. Thus, using 3.), 4.4), and 4.2), we get µ = L ) 2 = µ + [G r, L 0 ] 2 = µ + ρ ρ 2 )G r ) 2 = µ + ρ ρ 2 ) z ζ b c = µ + πρ ) ρ 2 ) = µ ρ π ρ 2 ) ρ 2 π ρ ) ρ 2 π ρ ). From the uniqueness of the left and right divisors we see that B r z + ) = B l z) and B l 2z) = L0 B r 2z). Using the first equation, we see that G r = G l and so c b )/b c ) = ã d )/d ã). In particular, µ 4.20) µ b c = z ζ ) ρ π ρ 2 ) = π d ã = z ζ ) ρ πγ z 2)γ ζ ) 2 ) γ z )γ ζ, ) which gives 4.8). Similarly, comparing the first components of the normalized vectors c and ã gives 4.7) and completes the proof. I am very grateful to I. Krichever, A. Borodin, M. Gekhtman, H. Sakai, F. Soloviev and T. Takenawa for interesting and helpful discussions. I also want to thank the organizers of the SIDE 8 conference for the invitation to participate at the conference and the opportunity to present a part of this work. Finally, I thank the referees for many useful suggestions. References References [] Veselov A 2003 Yang-Baxter maps and integrable dynamics Phys. Lett. A 343) [2] Veselov A 2007 Yang-Baxter maps: dynamical point of view. Combinatorial aspect of integrable systems MSJ Mem.), vol 7 Math. Soc. Japan, Tokyo) [3] Krichever I 2002 Isomonodromy equations on algebraic curves, canonical transformations and Whitham equations Mosc. Math. J. 24) [4] Krichever I 2002 Vector bundles and Lax equations on algebraic curves. Comm. Math. Phys. 2292) [5] Moser J and Veselov A 99 Discrete versions of some classical integrable systems and factorization of matrix polynomials. Comm. Math. Phys. 392) [6] Suris Y 2004 Discrete Lagrangian models Discrete integrable systems Lecture Notes in Phys. vol. 644) 84 [7] Borodin A 2004 Isomonodromy transformations of linear systems of difference equations. Ann. of Math. 2) 603) 4 82 [8] Krichever I 2004 Analytic theory of difference equations with rational and elliptic coefficients and the Riemann-Hilbert problem. Uspekhi Mat. Nauk 596) 7 54 [9] Arinkin D and Borodin A 2006 Moduli spaces of d-connections and difference Painlevé equations Duke Math. J. 343) [0] Sakai H 200 Rational surfaces associated with affine root systems and geometry of the Painlevé equations. Comm. Math. Phys. 220) [] Arinkin D and Borodin A 2007 Tau-function of discrete isomonodromy transformations and probability Preprint math/ [2] Krichever I and Phong D 998 Symplectic forms in the theory of solitons. Surveys in differential geometry: integral systems [integrable systems] Surv. Differ. Geom.) vol. IV Int. Press Boston, MA) [3] Dzhamay A 2008 On the Lagrangian structure of the discrete isospectral and isomonodromic transformations. Int. Math. Res. Not. IMRN rnn02 [4] Soloviev F 2008 Quadratic algebras and integrable chains Preprint math-ph/ [5] Veselov A 99 Integrable Lagrangian relations and factorization of matrix polynomials Funktsional. Anal. i Prilozhen. 252) [6] Jimbo M and Sakai S 996 A q-analog of the sixth Painlevé equation Lett. Math. Phys. 382)