arxiv: v2 [nlin.si] 28 Jan 2013

Transcription

1 arxiv:29.472v2 [nlin.si] 28 Jan 23 Discrete Painlevé equations and their Lax pairs as reductions of integrable lattice equations C. M. Ormerod Peter H. van der Kamp and G.R.W. Quispel Abstract. We describe a method to obtain Lax pairs for periodic reductions of a rather general class of integrable non-autonomous lattice equations. The method is applied to obtain reductions of the non-autonomous discrete Korteweg-de Vries equation and non-autonomous discrete Schwarzian Korteweg-de Vries equation which yield a discrete analogue of the fourth Painlevé equation a q-analogue of the sixth Painlevé equation and the q- Painlevé equation with a symmetry group of affine Weyl type E ) 6. Integrable partial difference equations are discrete time and discrete space analogues of integrable partial differential equations which often admit classical integrable partial differential equations as continuum limits [ ]. Integrable ordinary difference equations are discrete analogues of integrable ordinary differential equations. Integrable ordinary difference equations admit integrable ordinary differential equations as continuum limits [7]. Integrable ordinary and partial difference equations possess discrete analogues of many of the properties associated to the integrability of their continuous counterparts [ ]. We consider partial difference equations whose evolution on a lattice of points w lm is determined by the equation.) Qw lm w l+m w lm+ w l+m+ ;αβ) = where α and β are parameters associated with the horizontal and vertical edges respectively. The equation is imposed on each square on the space of independent variables lm) Z 2. From a suitable staircase of initial conditions [3] one may determine w lm for all lm) Z 2. Imposing the periodic constraint that.2) w l+sm+s 2 = w lm defines a periodic reduction [ ]. We will assume for simplicity that s and s 2 are both positive. In an analogous way to how similarity reductions of partial differential equations yield ordinary differential equations [6] periodic reductions given by.2) yield ordinary difference equations [ ]. Given a partial differential equation with some similarity reduction there is a procedure that allows one to obtain a Lax representation of the resulting ordinary differential equation from the Lax representation of the partial differential 2 Mathematics Subject Classification. 39A4; 37K5; 35Q5.

2 2 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL equation. This holds for autonomous and non-autonomous reductions [9]. The discrete analogue of this procedure is fairly straightforward for autonomous reductions [ ] however there has been no direct method for determining the Lax representation for non-autonomous reductions [3 4]. Given a reduction another task is to determine whether the reduction is a known system of difference equations. For autonomous reductions one may be able to find a certain parameterisation which identifies the system as a known QRT mapping [34 35] which may be classified in terms of elliptic surfaces [8]. For nonautonomous reductions one may be able find a parameterisation of the equation that identifies the system as one of the Painlevé equations which are classified by the group of symmetries of their surface of initial conditions [4]. The aim of this note is to demonstrate a method which we outline in by which we may directly obtain a Lax representation of both autonomous and nonautonomous reductions from a Lax representation of partial difference equations in an algorithmic manner. The method gives Lax representations in a manner that is general and concise enough to directly provide the Lax integrability of entire hierarchies of reductions. As an application of this method we present a reduction of the non-autonomous discrete Schwarzian Korteweg-de Vries equation which is a non-autonomous version of Q δ= in the classification of Adler Bobenko and Suris [2 3]) to the q-painlevé equation with E ) 6 symmetry which is associated with a surface with A ) 2 symmetry or q-pa ) 2 )):.3a).3b) y z )y z ) = a y )a 2 y )a 3 y )a 4 y ) b q 4 ty )b 2 q 4 ty ) yz )y z ) = θ z a )z a 2 )z a 3 )z a 4 ) b b 2 tz +θ )a a 2 a 3 a 4 +θ q 4 tz) where t = q 4 t the a i b i and θ are fixed parametersand q is some complex number whose modulus is not. This is the q-painlevé equation whose group of Bäcklund transformations is an affine Weyl group of type E ) 6 [4]. Our method stands in contrast to two methods of performing reductions of partial difference equations in the literature namely the method of Hydon et al. [38] which is based on the existence of certain Lie point symmetries and the method of Grammaticos and Ramani who perform autonomous reductions then deautonomize the equation via singularity confinement [36]. While the first method seems to rely on a similar approach to ours neither method gives rise to the associated linear problem for the reduced equation. The approach most similar to our method has been discussed by Hay et al. [4] in which the form of the monodromy matrix for autonomous reductions and its properties are used as an ansatz for an associated linear problem of the non-autonomous reductions of the lattice modified Kortewegde Vries equation. A further extension to this work successfully determined the associated linear problem for a hierarchy of systems [3]. To demonstrate our method we first provide some simple examples in 2. We present an autonomous reduction of the discrete potential Korteweg-de Vries equation dkdv) [24] then present the non-autonomous generalization of this example. In 3 we first present the q-painlevé equation associated with the A ) 3 surface otherwise known as q-p VI [7]) as a reduction of 3.) before going to the higher case where we present the above-mentioned reduction of 3.) to.3) which we believe to be the first known reduction to this equation.

3 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS3. The method We start by imposing.2) as a constraint on our initial conditions then the periodicity gives us that there are s +s 2 independent initial conditions to define. We solve this periodicity constraint by a specific labelling following [3 43]; let s = ag and s 2 = bg where g = gcds s 2 ) then the direction of the generating shift cd) associated with the increment n n+ is chosen so that ) a b det =. c d We specify an n Z and a p Z g by letting ) ) a b l m n = det p det l m c d where the labelling of variables is specified by.) w lm w p n. mod g In the case in which g = the superscript will be omitted. The reduction in the autonomous case is a system of g equations given by.2) Qw p n wp+d n b wp c n+a wp c+d n+a b ;αβ) = p =...g where α and β are constants. In the nonautonomous setting we have.3) Qw p n wp+d n b wp c n+a wp c+d n+a b ;α lβ m ) = p =...g whereα l and β m will be aposteriori constrainedfunctions ofland m. We will now outline how to obtain Lax representations for the autonomous and nonautonomous reductions respectively... Autonomous reductions. It is known that multilinear partial difference equations that are consistent around a cube are in a sense their own Lax pair [2 4 5]. For a generic multilinear equation.) that is consistent around a cube a Lax pair may be written as.4a).4b) where.5a).5b) φ l+m = L lm φ lm φ lm+ = M lm φ lm Qxuv;αγ) L lm = λ lm v 2 Qxuvy;αγ) v y Qxuv;βγ) M lm = µ lm u 2 Qxuvy;βγ) u y Qxu;αγ) Qxuy;αγ) y x = w lm u = w l+m Qxv;βγ) Qxvy;βγ) y x = w lm v = w lm+ where γ is a spectral parameter. The compatibility condition is.6) M l+m L lm = L lm+ M lm

4 4 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL forcing the prefactors λ lm and µ lm to be chosen in a manner that satisfies the equation detl lm+ = detm l+m. detl lm detm lm When the prefactors are appropriately chosen imposing.6) is equivalent to.). In practice it is often computationally convenient to deal with some transformation of this Lax pair. To obtain a Lax representation for the system of ordinary difference equations.2) we define two operators A n and B n associated with the shifts lm) l+s m+s 2 ) and the generating shift lm) l+cm+d) respectively. These operators have the effect.7a).7b) φ n = A n φ n φ n+ = B n φ n where one representation that is simple to write is as follows: A n B n s 2 j= d j= M l+sm+j M l+cm+j c s i= L l+im L l+im i= where the dependence on n and p is specified by The compatibility condition L lm w lm w l+m ;γ) L p n γ) = Lp n wp n wp+d n b ;γ) M lm w lm w l+m ;γ) M p n γ) = Mp n wp n wp c n+a ;γ)..8) A n+ B n B n A n = is equivalent to imposing.2). We call A n the monodromy matrix for the following reason: by identifying all points in Z 2 that are multiples of s s 2 ) apart we may consider the space in which the new system exists as being cylindrical. We wrap around in a manner that connects the points that are identified by the periodic reduction. The monodromy matrix rather than presenting a trivial action as.7a) suggests expresses the action of wrapping around the cylinder as in figure. The monodromy matrix can be expressed as a function of the s + s 2 initial conditions wnw n+...w g a+b ) by following the standard staircase. Geometrically the standard staircase is the path between two lines which squeeze a set of squares with the same values i.e. a set of squares shifted by s s 2 ) [43]. One advantage of the generating shift is that every other shift in n may be expressed as some power of the generating shift by construction [44]. Furthermore this generating shift allows us to constrain the non-linear component where we need to use.2) to just g places. We have illustrated the standard staircase and generating shift in figure 2. In practice the product follows the path of a standard staircase [44].

5 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS5 w 2 n w n w n Figure. A pictorial representation of the way in which the monodromy matrix wraps around to similar points for a 96)- reduction. The points w p given for reference to figure 2. w 2 n+ w n w n w n+ w 2 n w 2 n w n+ w n w n w n w n Figure 2. The full labelling of variables in the 9 6)-reduction of figure. Inthisexampletheshiftpn) p+n)correspondsto the shift ab) = 32) and the shift pn) pn+) corresponds to the shift cd) = ). In the example defined by figure 2 if we allow our monodromy matrix to follow the standard staircase the monodromy matrix is A n L l+8m+6 M l+8m+5 L l+7m+5 L l+6m+5 M l+6m+4 L l+5m+4 M l+5m+3 L l+4m+3 L l+3m+3 M l+3m+2 L l+2m+2 M l+2m+ L l+m+ L lm+ M lm

6 6 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL and the other half of the Lax pair is B n M l+m L lm. The resulting compatibility condition.8) gives the evolution equations for the w i n+ i = 2: Qw n 2w 2 n 4w n+w n ;αβ) = Qw 2 n 2w n 4w n+w 2 n ;αβ) = Qwn 2 w n 4 w2 n+ w n ;αβ) =. In general this procedure gives us an s +s 2 dimensional mapping φ : C s+s2 C s+s2 which applied to wnw n+...wg n+a+b ) gives w n+ w n+2...wg n+a+b ). This new set of values forms a new standard staircase. As a matter of fact this new standard staircase is the old one translated by the generating shift..2. Nonautonomous reductions. To deautonomize this theory we consider the α and β to be functions of l and m. As L lm and M lm are shifted in only m and l respectively in the compatibility condition.6) replacing α and β with α l and β m which are arbitrary functions of l and m respectively preserves the Lax integrability. Hence our basic non-autonomous lattice equations may be considered to be of the form.9) Qw lm w l+m w lm+ w l+m+ ;α l β m ) = where the Lax representation is specified by.4) where.a) Qxuv;α lγ) L lm = λ lm v 2 Qxuvy;α l γ) v y.b) Qxuv;βmγ) M lm = µ lm u 2 Qxuvy;β m γ) u y Qxu;α l γ) Qxuy;α l γ) y x = w lm u = w l+m Qxv;β m γ) Qxvy;β m γ) y x = w lm v = w lm+ where γ is a spectral parameter and the prefactors λ lm and µ lm are chosen to satisfy the compatibility conditions in an analogous manner to the autonomous case. Ifoneassumesthatthe αandβ arefunctionsofboth l andm i.e. α = α lm and β = β lm then demanding that α lm is independent of m and β lm is independent of l has also been shown to be a necessary condition for singularity confinement for equations in the ABS list []. The above constitutes a Lax pair interpretation of this constraint. Let us now specialise our choice of systems to those that admit representations of the additive form.) Qw lm w l+m w lm+ w l+m+ ;α l β m ) =

7 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS7 or the multiplicative form.2) Q w lm w l+m w lm+ w l+m+ ; α ) l = β m with a possible additional dependence on α l+ α l and β m+ β m in the additive case or α l+ /α l and β m+ /β m in the multiplicative case. A list of transformed equations appears in table where subscripts m and a denote those functions.) dependent on a multiplicative or additive combination of α l and β m respectively. This list is restricted to equations that we could find transformations to forms admitting additive or multiplicative reductions. This does not include the equations Q3 or Q4. ABS Qxuvy;α l β m ) H a w lm w l+m+ )w l+m w lm+ )+β m α l H m w lm β ) m+ w l+m+ w l+m β ) m+ w lm+ + α2 l β m β m βm 2 H2 m w lm β ) m+ w l+m+ w l+m β ) m+ w lm+ α2 l β m β m βm 2 + α ) l w lm +w l+m + β ) m+ w lm+ +w l+m+ ) + β m β m H3 δ= m H3 δ m Q δ= m Q δ m Q2 m α l w lm w l+m +w lm+ w l+m+ ) w lm w lm+ +w l+m w l+m+ ) β m α 2 ) l βm 2 w lm w l+m + β2 m+ βm 2 w lm+ w l+m+ β ) m+ α 4 w lm w lm+ +w l+m w l+m+ )+δ l β m βm 4 α l w lm w lm+ )w l+m w l+m+ ) β m w lm w l+m )w lm+ w l+m+ ) α 2 l βm 2 w lm β ) m+ w lm+ w l+m β ) m+ w l+m+ β m β m β m+ w lm w l+m )w lm+ w l+m+ )+ δα2 l α 2 l β m βm 2 ) ) α l w l+m β2 m+ β m βm 2 w l+m+ w lm β2 m+ βm 2 w lm+ β2 m+ βm 2 αl β 2 m ) w lm w l+m )w lm+ w l+m+ ) α ) l w lm +w l+m + β2 m+ β m β m βm 2 α ) l αl α 2 l β m β m βm 2 α ) l + β m Table. A list of various lattice equations taken from [2 3]) in a suitable form for non-autonomous reductions. ) w lm+ + β2 m+ βm 2 w l+m+

8 8 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL With each lattice equation written in terms of α l β m or α l /β m the necessary requirement for.2) to be consistent is the requirement that in the additive case and α l β m = α l+s β m+s2 α l β m = α l+s β m+s2 in the multiplicative case. By a separation of variables argument we define h and q by letting.3a) α l+s α l = β m+s2 β m := habg.3b) α l+s α l = β m+s 2 β m := q abg in the additive and multiplicative cases respectively. Although it is not a technical requirement we will assume that h is not and that q is not a root of unity. We solve the additive and multiplicative case by letting α l = hlb+a l β m = hma+b m α l = a l q bl β m = b m q am where a l and b m are sequences that are periodic of order s and s 2 respectively not related to the constants a and b). This choice of α l and β m ensures the consistency of the reduction with as many degrees of freedom as the sum of the orders of the difference equations satisfied by α l and β m.3a) and.3b) i.e. s +s 2. To provide a non-autonomous Lax pair for the non-autonomous reduction we need to choose a spectral variable x in a manner that couples a linearly independent direction with the spectral variable γ. While any linearly independent direction may be considered a valid choice we present a simple choice. Our choice of spectral parameter is specified by introducing the variable k = l and x = hbk γ in the additive case and x = q bk /γ in the multiplicative case. In the additive case L lm = L lm α l γ) L n a l +x) M lm = M lm β m α l )+α l γ)) M n x+hn+b m ) and in the multiplicative case L lm = L lm α l /γ) = L lm a l x) M lm = M lm β m /α l )α l /γ)) = M lm b m xq n ). This gives us a non-standard Lax pair which in the additive case reads and in the multiplicative case reads Y n x+abgh) = A n x)y n x) Y n+ x+cbh) = B n x)y n x) Y n q abg x) = A n x)y n x) Y n+ q cb x) = B n x)y n x)

9 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS9 where.4a).4b) The compatibility conditions s 2 A n x) j= d B n x) j= M l+sm+j M l+cm+j c s i= L l+im L l+im. i= A n+ x+cbh)b n x) = B n x+abgh)a n x) A n+ q cb x)b n x) = B n q abg x)a n x) in the additive and multiplicative cases respectively gives us.3). This choice of spectral variable has the advantage that the spectral matrix and deformation matrix A n x) and B n x) have a simple dependence on the independent variable n. 2. Some simple examples In this section wepresentsomeexamples ofthe theoryabove. An examplethat hasappearedrecentlyistheexampleofq-p VI asareductionofthediscretemodified Korteweg-de Vries equation [28] which also gave rise via ultradiscretization to the first known Lax representation of u-p VI. 2.. Autonomous example. We consider some additive examples in particular we will consider reductions of the discrete potential Korteweg-de Vries equation 2.) w lm w l+m+ )w l+m w lm+ ) = α β labelled as H a in table which possesses a Lax representation of the form.4) where L lm and M lm are specified by ) wlm α γ w L lm = lm w 2.2a) l+m w l+m ) wlm β γ w M lm = lm w 2.2b) lm+. w lm+ Let us consider a reduction.2) where s = 2 and s 2 = with a labelling indicated in figure 3. This gives us g = a = 2 and b = hence n = 2m l and the direction that characterises the generating shift c d) is chosen to be ). w 4 w 3 w 4 w 3 w 2 w w 3 w 2 w w w w Figure 3. The labelling of initial conditions with 2) periodicity and an evolution in the )-direction.

10 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL The product formula for the monodromy matrix A n and the matrix that is related to the generating shift B n are A n =M l+2m L l+m L lm ) ) wn β γ w = n w n+2 wn+ α γ w n w n+ w n+2 w n ) wn+2 α γ w n+ w n+2 w n+ B n =M l+m L lm ) ) wn+ β γ w = n+ w n+3 wn+2 α γ w n+ w n+2. w n+3 w n+ The compatibility condition given by.8) reads This simplifies to M l+3m+ L l+2m+ L l+m+ M l+m L lm = M l+3m+ L l+2m+ M l+2m L l+m L lm. L l+m+ M l+m = M l+2m L l+m ) ) wn+3 α γ w n+2 w n+3 wn+ β γ w n+ w n+3 w n+2 w n+3 ) ) wn β γ w = n w n+2 wn+ α γ w n w n+ w n+2 w n which defines the evolution of this autonomous reduction to be given by the equation 2.3) w n w n+3 )w n+ w n+2 ) = α β. If we let y n = w n w n+ this equation is equivalent to 2.4) y n +y n +y n+ = α β y n which is a well known example of a second order difference equation of QRT type Nonautonomous example. The autonomous equation and Lax representation generalise naturally to the non-autonomous case by replacing α and β by α l and β m respectively. Furthermore we may satisfy the periodicity constraint We solve this constraint by letting α l+2 α l = β m+ β m := 2h. α l = hl+a l β m = 2hm+b m where a l is periodic of order two and b m is constant and hence may be taken to be without loss of generality. The evolution equation for this system may be represented as an application of the nonautonomous version of 2.) translated by the vector ); 2.5) w n w n+3 )w n+ w n+2 ) = α l+ β m+. recalling that n = 2m l. The increment in l and m by directly corresponds to the increment in n by. In the simplest case where a l = a is constant rather than periodic) letting y n = w n w n+

11 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS results in the evolution equation or alternatively y n +y n+ +y n+2 = α l+ β m+ y n+ = hn h+a y n+ 2.6) y n +y n +y n+ = hn+a y n. To form the Lax pair for this reduction we choose a spectral variable to be x = hl γ. Since α l and β m appear with γ in L lm and M lm this gives us α l γ = hl γ +a = a +x β m γ = 2hm γ = hn+x. which means we may interpret our l and m variables in terms of x and n variables. Thelm) l+2m+)shiftandlm) l+m+)givesusnx) nx+2h) and nx) n+x+h) hence we have a linear system of the form where Y n x+2h) = A n x)y n x) Y n+ x+h) = B n x)y n x) A n x) M l+2m L l+m L lm wn hn+x w = n w n+2 ) ) wn+ a 2 +x w n w n+ w n+2 w n ) wn+2 a +x w n+ w n+2 w n+ B n x) M l+m L lm ) ) wn+ hn+x w = n+ w n+3 wn+2 a +x w n+ w n+2. w n+3 w n+ The compatibility condition is 2.7) A n+ x+h)b n x) = B n x+2h)a n x) which gives 2.5). To express everything in terms of y n = w n w n+ we use the matrices obtained from applying a gauge transformation L lm S l+m L lms lm and M lm S lm+ M lms lm where S lm = ) wlm

12 2 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL in which case A n x) and B n x) are given explicitly in terms of products of matrices with entries expressible in terms of the y n variables: yn y A n x) = n+ y n +y n+ ) 2 ) +hn+x yn a +h+yn 2 +x ) y n y n+ y n yn+ a +yn+ 2 +x ) y n+ yn+ y B n x) = n+2 y n+ +y n+2 ) 2 ) +hn+x y n+ y n+2 ) yn+ a +yn+ 2 +x. y n+ for which 2.7) now gives 2.6) as required. We recover Lax matrices for 2.4) by letting A n := lim A n x) h B n := lim B n x) h whose compatibility.8) gives 2.4) where α β = a. This is a way in which the autonomous and nonautonomous reductions are related. While the simple case above demonstrates the basic mechanisms in the method the periodicity constraint allows us to build in an extra variable. If we allow the full generality of a periodic value of a l by letting { hl+a where l is odd α l = hl+a 2 where l is even we obtain A n x) M l+2m L l+m L lm wn hn+x w A n x) = n w n+2 ) ) wn+ a 2 +x w n w n+ w n+2 w n ) wn+2 a +x+w n+ w n+2 w n+ B n x) M l+m L lm ) ) wn+ hn+x w B n x) = n+ w n+3 wn+2 a +x w n+ w n+2. w n+3 w n+ The compatibility condition needs to take into account that the shift n x) n + x + h) also shifts the position on the lattice hence we need to couple the non-linear component with a swapping of the roles of a and a 2 hence the compatibility condition given by 2.7) results in the evolution w n w n+3 )w n+ w n+2 ) = hn+a 2 a a 2 +h a 2 a h. If we restrict our attention to even powers of this map the a i are constant. We proceed to specify a change of variables which is motivated by some historical context. In the 98 s Novikov and Veselov formalized the derivation of the Hamiltonian structure of hierarchies of soliton equations from their Lax representations [45 9]. The Darboux coordinates are the poles of the Baker-Akhiezer function

13 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS 3 and the eigenvalues of the spectral matrix at these poles which for 2 2 matrices coincides with the roots in the spectral variable of the off-diagonal elements and the diagonal elements evaluated at those roots. The Hamiltonian description of the isomonodromic deformations of [6] are expressed in these coordinates. While the link between the symplectic structure for discrete Painlevé equations and the discrete isomonodromic deformations is not as well developed it is interesting to note that the parameterizations of many known discrete Lax pairs are provided in terms of these coordinates [2 4 7]. Perhaps a discrete analogue of the Baker- Akhiezer function holds the key to linking the geometric theory and the theory of discrete isomonodromic deformations. It may also provide an algorithmic manner of describing the symplectic structure for hierarchies of equations. Our new coordinates are y n z n ) where y n is the root in the spectral variable of the 2)-element and z n is a variable which parameterizes the diagonal elements at this root. In the lattice variables w n these are y n = w n+2 w n )w n+2 w n+ ) a 2 z n = w n+2 w n )y n +a ) y n +nh from which we extricate the second-order system This system is equivalent to y n +y n+2 = z 2 n a +a 2 ) z n z n+2 = y n+2 +a )y n+2 +a 2 ). y n+2 +n+2)h) y n +y n 2 +a +a 2 )y n +y n+2 +a +a 2 ) = y n +a ) 2 y n +a 2 ) 2 y n +nh) 2 which is a special version of d-p IV [37]. The second power of the generating shift is equivalent to the shift lm) lm+) hence a simplified Lax representation for this system is Y n x+2h) = A n x)y n x) Y n+2 x) = B n x)y n x) where B n x) M lm. We may simplify the spectral matrix via a gauge transformation to be A n x) y n +a )y n +a 2 ) x 2 +δx+ǫ z n x y n y n +nh)z n where δ =y n +a +a 2 +hn ǫ =y n +hn)y n +a )+hn+y n +a )a 2 and a deformation matrix under the same transformation becomes ) zn x+hn+zn B n x) = 2. z n Using this method one is able to provide product formulas for the Lax representations of q-p II q-p V of Hay et al. [4] and the Lax pair for a version of q-p III

14 4 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL of Joshi et al. [8]. Furthermore this shows and explains that the Lax pairs for the reductions of [4] factorize in a nice way. The resulting factorizations provide a simple way to compute the compatibility. 3. Reductions of dskdv We consider periodic reductions of the nonautonomous discrete Schwarzian Korteweg-de Vries equation; ) 3.) α l + w lm+ w l+m+ w l+m w lm ) = β m + w l+m w l+m+ w lm+ w lm which possesses a Lax representation of the form.) where ) w lm w l+m 3.2a) L lm = α l γw lm w l+m ) w lm w lm+ 3.2b) M lm = β m. γw lm w lm+ ) 3.. q-pa ) 3 ). Recently one of the authors derived a q-analogue of the sixth Painlevé equation or q-pa ) 3 )) as a reduction of the non-autonomous modified Korteweg-de Vries equation [28]. This work demonstrated the specified method where the transformation between the lattice variables the w i and the Painlevé variables was relatively simple. We provide a similar but more complicated relation between the lattice variables for the discrete Schwarzian Korteweg-de Vries equation and the Painlevé variables of q-pa ) 3 ). We now provide a new reduction of the lattice equation to a version of q-p VI given by b b 2 z ) q 2 z ) 3.3a) 3.3b) yy = q 2 b b 2 t θ z)b b 2 t θ 2 z) b q 2 ty ) b 2 q 2 ty ) zz = q 2 a y )a 2 y ) t = q 2 t where θ θ 2 = a a 2 b b 2. This extends our previous result in [28] by the addition of an extra parameter which represents an integral that is used to reduce the order of the map given by θ or θ 2 ). Obtaining a q-analogue of the sixth Painlevé equation from a discrete analogue of the Schwarzian KdV equation has some historic significance as it could be considered a discrete analogue of the reduction of the Schwarzian KdV equation to the sixth Painlevé equation [25 26]. We consider a reduction of 3.) given by.2) where s = s 2 = 2. The labelling given by.) is depicted in figure 4. The constraint.3b) becomes 3.4) which introduces the parameter q. α l+2 α l = β m+2 β m := q 2

15 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS 5 w 3 w 2 w 3 w 2 w w 2 w w w w w Figure 4. The reduction and the labelling of variables. We satisfy the periodicity constraints of 3.4) by explicitly setting α l = { a q l if l is even a 2 q l if l is odd β m = { b q m if m is even b 2 q m if m is odd. Naturally the dependent variable n is invariant along the direction of the reduction. The correspondence between the associated linear problem here and that of Sakai [4] is made more natural by specifying an independent variable t and our spectral variable x following previous sections: t = q m l = q n x = ql γ = qk γ. For the sake of clarity the reduced lattice variables may be regarded as functions of t or n under the identification w i n+k = w i q k t) so that the shift n n+ is equivalent to the shift t qt. We know that the shift lm) l +2m+2) is equivalent to the shift xt) q 2 xt). Since a = and b = the direction of our generating shift is lm) lm+) which is equivalent to the shift xt) xqt). This gives us a Lax pair of the form 3.5a) 3.5b) Yq 2 xt) = Axt)Yxt) Yxqt) = Bxt)Yxt) where the matrix Axt) is given by the product 3.6) 3.7) Axt) L l+m+2 M l+m+ L lm+ M lm wn+ Axt) = w w n n a 2 x w n+ txb wn+ 2 w n wn wn+ wn+ wn wn wn+ a x txb wn+. w n wn wn+

16 6 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL The deformation matrix Bxt) corresponds to the shift m m + hence is given by 3.8) 3.9) The compatibility condition written as Bxt) M lm wn Bxt) = w n+ txb. wn wn+ Axqt)Bxt) = Bq 2 xt)axt) yeilds the evolution equation ) ) wn+ =a wn w n +tb w n wn w 3.a) ) n a w n wn) +tb w n wn b b 2 q ) wn+ =a 2 w n wn w n +tb 2 wn w n wn ) 3.b) ) a 2 w n wn) +tb2 w n wn b 2 qb. Notice that the generating shift swaps the roles of b and b 2 coupled with multiplicative factors introduced to compensate for the dependence of β m on q m. We claim that the second iterate of this mapping is q-pa ) 3 ). To make the full correspondence with 3.3) let us expand out 3.7) to give a matrix of the form Axt) = A +A x+a 2 x 2 where A = I and the eigenvalues of A 2 are θ t and θ 2 t where ) w n wn+) where θ i = θ i n) satisfies θ t = tb b 2 w n+ wn ) w n wn+ w n wn+) θ 2 t = a ) a 2 w n wn+ w n wn+) t ) wn+ n) w w n wn+ θ n+)θ n) = a a 2 b b 2 q θ 2 n+)θ 2 n) = qa a 2 b b 2 hence θ i n + 2) = θ i n). That is to say that θ and θ 2 are 2-integrals of the generating shift [2] with the additional constraint θ θ 2 = a a 2 b b 2. Secondly from the product form we have that detaxt) = a x )a 2 x )b tx )b 2 tx ). In accordance with the motivation given in the previous section we choose to parameterize the spectral matrix Ax t) in terms of the x-root of 2 )-element of Axt) and the diagonal entries at that root. This provides us with a spectral matrix in the curious form 3.) Ax t) = txx y n +ζ n )θ + txδ n ω n θ 2 txx y n )θ txx y n +η n )θ 2 + ω n

17 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS 7 where ζ n = y na )y n a 2 ) η n = b ty n )b 2 ty n )z n. z n Fixing the determinant requires that δ n = + + a a 2 b t + b 2 t 2y+ζ n +η n. Under this identification the y n and z n are specified by the lattice variables 3.2a) [ y n = wn wn+)w n+ wn)a 2 wn+ wn)+b twn wn+)) ] +a wn w n+ )w n w n+ ))+b 2twn w n+ )w n+ w n )w n w n+ ) [ b twn+ w n )a 2wn+ w n+ )w n+ w n )+b 2twn+ w n ) w n w n+)) a a 2 w n w n+)w n w n+) 2] 3.2b) [ z n = wn wn+)w n wn+)a wn+ wn)w n wn+) ] tb 2 wn wn+)w n+ wn) b twn+ wn+)w n+ wn)) [ wn w n+ )w n+ w n )w n w n+ )a wn a 2wn+ ] +a 2 a )wn+)+b 2 twn+ wn)w n wn+)). We may now make the correspondence with 3.3) via the identification that y = y n y = y n+2 z = z n and z = z n+2. One may verify 3.3) remarkably using the evolution 3.) alone. We also have the gauge factor given by [ ω n = b b 2 t 2 wn+ 2 n) w w n wn+ ) 2 ] [ b t wn+ ) w n a2 w n+ wn+ ) w n+ wn ) +b2 t wn+ ) w n w n wn+)) a a 2 w n wn+) w n wn+ ) ] 2 which satisfies the equation ω n+2 z n )θ z n a 2 y n )+) a 2 tb b 2 y n ) = ω n y n θ z n tb b 2 )a 2 2 y nz n +θ ty n a 2 tb +b 2 )y n +z n )). We may now parameterize the deformation matrix for the double shift in terms of y n z n and w n as Bxt) = I+ xθ tb b 2) δ nω nθ 2 b +b 2+ζ n y n)θ b +b 2+ζ n y n)θ xtb y n ζ n) )θ tb 2 )tb 2y ζ n) )θ tb 2 2) xb b 2b +b 2+ζ n y n)θ )t 2 δ nω nθ tb b 2)θ 2 θ tb b 2 Curiously the coefficient of x is lower triangular and the constant coefficient is upper-triangular. This rather simplified Lax pair comes at the expense of the requirement that the variables y n and z n explicitly lie on the biquadratic Vyzt) =θ z ) 2 θ yz )a +a 2 )z tb +b 2 )) +y 2 a a 2 z θ t)θ z tb b 2 ) =.

18 8 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL which makes the computation of the compatibility condition Axq 2 t)bxt) = Bq 2 xt)axt) slightly more difficult. One can verify directly that Vy z t) = implies that Vy q 2 zq 2 t) = and Vy z q 2 t) = under the evolution defined by 3.3). This forms an explicit parameterization of the member of the pencil of biquadratics for each t. The existence of such a parameterization is not without precedent and has appeared in the work of Yamada [47] and Noumi et al. [27]. We do remark that it is interesting that this explicit dependence on the bi-quadratic curve essentially came from a condition on the Lax matrices q-pa ) 2 ) as a reduction of dskdv. While a reduction to q-pa) 3 ) has been provided by one of the authors previously [28] we know of no reduction from anintegrablelattice equationto anymember ofthe hierarchyaboveq-pa ) 3 ). We wish to extend this further and consider reductions of a Painlevé equation that is higher up in the classification scheme; namely the q-painlevé equation with E ) 6 symmetry which is associated with a surface with A ) 2 symmetry or q-pa ) 2 )) given by.3). We will find that this equation appears as a reduction of 3.) or Q δ= m [23]. or the non-autonomous discrete Schwarzian Korteweg-de Vries equation w 6 w 5 w 6 w 5 w 4 w 3 w6 w5 w4 w3 w 2 w w 5 w 4 w 3 w 2 w w w 3 w 2 w w w w Figure 5. The reduction and the labelling of variables. We impose.2) on 3.) with s = 4 and s 2 = 2 hence our constaint.3b) becomes α l+4 = β m+2 := q 4 α l β m which we solve in a similar manner as before; a q l if l = mod 4 a 3.3) α l = 2 q l if l = mod 4 a 3 q l β if l = 2 mod 4 m = a 4 q l if l = 3 mod 4 { b q 2m if m = mod 2 b 2 q 2m if m = mod 2 We use the same Lax pair for 3.) as the previous section namely.4) where L lm and M lm are specified by 3.2). The t-direction is also chosen to be constant in the direction of the reduction. We choose our spectral parameter x and independent variable t in a manner in which correspondence with q-pa ) 2 ) comes more.

19 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS 9 naturally that is t = q 2m l = q n x = qk γ = ql γ. We note that the generating shift is no longer just a shift in m alone but a simultaneous shift in l and m i.e. lm) l +m+) hence we present a Lax pair that represents the shift lm) l + 4m + 2) which is now equivalent to xt) q 4 xt) and the shift lm) l + m + ) which is equivalent to the shift xt) qxqt). Hence our Lax pair is a linear system satisfying 3.4a) 3.4b) Yq 4 xt) = Axt)Yxt) Yqxqt) = Bxt)Yxt) where we compose A and B in terms of L and M in the following way 3.5a) Axt) L l+3m+2 L l+2m+2 M l+2m+ L l+m+ L lm+ M lm wn+ Axt) = w n wn+2 xa 4 w n+ 3.5b) xa 3 wn+ w n wn+2 w n+ wn wn+2 wn+ wn txb 2 xa 2 wn wn+2 wn+ w n wn+2 wn+ wn wn+2 xa txb wn+2. w n+ wn w n+2 The deformation matrix may be written as 3.5c) Bxt) L lm+ M lm 3.5d) wn+2 Bxt) = w w n+ n wn+2 xa txb wn+2 w n+ wn w n+2 where we have written these as functions of x and t. The compatibility of the system defined by 3.4) is given by 3.6) Aqxqt)Bxt) = Bq 4 xt)axt). One may use 3.6) to derive the required evolution equations; the equations defining the evolution may be written as wn+ = a ) 4 w n 2 wn w n +b 2 twn w n wn 2) 3.7a) a 4 w n 2 wn) +b2 t ) wn w n 2 3.7b) 3.7c) 3.7d) wn+ = a ) 2wn w n 2 wn +b t ) wn wn 2 w n wn 2 a 2 b t) a 2 wn +b twn a a 2 q a 2 a 3 q a 3 a 4 q a 4 a q 3 b b 2 q 2 b 2 q 2 b.

20 2 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL Thedifficult step is toextract asecondordersystemfrom this seeminglysixth order system via some special parameterization. To do this we observe the properties of 3.4) as an associated linear problem. Firstly expanding out 3.5b) in the spectral parameter we find Axt) = A +A x+a 2 x 2 +A 3 x 3 where A is the identity matrix and A 3 is a lower triangular matrix with diagonal entries θ t and θ 2 t where θ = a ) ) 2a 3 b w n wn+ w n+2 wn+ w n wn+2) ) ) w n wn+ w n+2 wn w n+ wn+2) θ 2 = a ) ) a 4 b 2 w n wn+ w n wn+2 w n+ wn+2) ) ) w n+ wn w n+2 wn+ w n wn+2) where θ and θ 2 are invariants in accordance with the evolution equations given by 3.7). The simplicity of the individual factors of Axt) from 3.5b) give us that the determinant is 3.8) detaxt) = xa )xa 2 )xa 3 )xa 4 )txb )txb 2 ). These properties should remind us of the properties of a special case of the spectral matrix in the Lax pair of Sakai [4]. We may simply use a transformation of the form Yxt) SYxt) where S is a lower triangular matrix that is constant in x such that the transformation A SAS diagonalises A 3. While this matrix S isn t particularly nice to write down the resulting matrix is in the general form Axt) = x θ tx y n )x ǫ n )+ζ n ) θ 2 tω n x y n ) θ txγ n +δ n ) +I. θ 2 tx y n )x χ n )+η n ) ω n The terms ǫ n χ n γ n and δ n may be determined from 3.8) in terms of y n and ζ n η n the a i s and the b i s. The relation between this form of Axt) and the known Lax pairs of Sakai [4] and Yamada [47] has recently been found by one of the authors [46]. Following the motivation from the previous section this parameterization defines y n to be the root of the 2)-entry leaving a choice of ζ n and η n such that θ tζ n y n +)θ 2 tη n y n +) = detay n t) where the right hand side is defined by 3.8). With this in mind the defining equations for ζ n and η n are θ tζ n y n + = y n z n )b ty n )b 2 ty n ) θ 2 tη n y n + = a y n )a 2 y n )a 3 y n )a 4 y n ). y n z n )

21 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS 2 This specifies y n and z n in terms of the lattice variables: 3.9) [ y n = wn+2 wn+)w n wn+)w n wn+2)b 2 twn wn)w n+ wn+2) a 3 w n+ w n)w n w n+2))+a 2 w n w n)w n w n+)w n w n+2) wn+2 wn+))+a wn wn+)w n wn+)w n+2 wn)w n wn+2) ] [ wn+ wn+2) a wn wn+)w n+2 wn)a 3 wn+ wn) w n w n+2 )w n+2 w n+ )+b 2tw n w n+ )w n w n+2 ) w n+ w n+2 )) a 2a 3 w n+ w n )w n w n+ )w n+2 w n+ ) w n w n+2 )2] 3.2) z n =[a wn w n+ )w n+2 w n+2 )a 3wn+ w n )w n w n+2 )w n+2 w n+ ) +b 2 tw n w n+ )w n w n+2 )w n+ w n+2 )) a 2w n+2 w n+ ) w n w n+2 )w n+2 w n )w n+2 w n+ ))+a 2a 3 w n+ w n )w n+2 w n+ ) w n+ w n+2 )w n w n+2 )2 ] [w n w n+2 )w n+2 w n+ )w n w n+ ) a w n+2 w n )w n+ w n+2 ) b 2tw n+2 w n+ )w n w n+2 )) +a 2 w n w n+ )w n+2 w n+ )w n w n+2 ))] which satisfy.3a) and.3b) under the identification of y = y n z = z n y = y n+4 and z = z n+4. While these expressions may fail to be succinct they do succeed in being very explicit. Furthermore using this identification it is possible to directly verify.3a) and.3b) from 3.7) alone. 4. Conclusion Given a nonautonomous partial difference equation which admits an additive or multiplicative form.) or.2) we have outlined a direct method for finding a Lax representation for any periodic travelling wave) reduction of the form.2). The method outlines how the Lax matrices may be expressed in terms of products of the Lax matrices of the partial difference equation from which it was derived. We have shown that the method applies to deautonomized versions of the equations in the ABS list with the exception of Q3 and Q4. We have concentrated on two cases the discrete Korteweg-de Vries equation and the discrete Schwarzian Korteweg-de Vries equation. We have also found the relation between the derived reductions and discrete Painlevéequations d-p I and d-p IV respectively reductions from the discrete Schwarzian Korteweg-de Vries equation and q-p VI and q-pa ) 2 ). The latter is to ourbest knowledge the highest full parametermember ofthe Sakai classification derived as a reduction so far. Acknowledgments This research is supported by Australian Research Council Discovery Grant #DP77.

22 22 C. M. ORMEROD PETER H. VAN DER KAMP AND G.R.W. QUISPEL References [] M. J. Ablowitz and J. F. Ladik. On the solution of a class of nonlinear partial difference equations. Studies in Appl. Math. 57): 2 976/77. [2] V. E. Adler A. I. Bobenko and Yu. B. Suris. Classification of integrable equations on quadgraphs. The consistency approach. Comm. Math. Phys. 2333): [3] V. E. Adler A. I. Bobenko and Yu. B. Suris. Discrete nonlinear hyperbolic equations: classification of integrable cases. Funct. Anal. Appl. 43): [4] Alexander I. Bobenko and Yuri B. Suris. Integrable systems on quad-graphs. Int. Math. Res. Not. ): [5] T. Bridgman W. Hereman G.R.W. Quispel and P.H. Kamp. Symbolic Computation of Lax pairs of Partial Difference Equations using Consistency Around the Cube. Foundations of Computational Mathematics pages [6] Peter A. Clarkson and Martin D. Kruskal. New similarity reductions of the Boussinesq equation. J. Math. Phys. 3): [7] Robert Conte editor. The Painlevé property. CRM Series in Mathematical Physics. Springer- Verlag New York 999. [8] Johannes J. Duistermaat. Discrete integrable systems. Springer Monographs in Mathematics. Springer New York 2. [9] Hermann Flaschka and Alan C. Newell. Monodromy- and spectrum-preserving deformations. I. Comm. Math. Phys. 76): [] B. Grammaticos A. Ramani J. Satsuma R. Willox and A. S. Carstea. Reductions of integrable lattices. J. Nonlinear Math. Phys. 2suppl. ): [] Basil Grammaticos and Alfred Ramani. Singularity confinement property for the nonautonomous) Adler-Bobenko-Suris integrable lattice equations. Lett. Math. Phys. 92): [2] F. A. Haggar G. B. Byrnes G. R. W. Quispel and H. W. Capel. k-integrals and k-lie symmetries in discrete dynamical systems. Physica A: Statistical Mechanics and its Applications 2332): [3] Mike Hay. Hierarchies of nonlinear integrable q-difference equations from series of Lax pairs. J. Phys. A 434): [4] Mike Hay Jarmo Hietarinta Nalini Joshi and Frank W. Nijhoff. A Lax pair for a lattice modified KdV equation reductions to q-painlevé equations and associated Lax pairs. J. Phys. A 42):F6 F [5] Ryogo Hirota. Nonlinear partial difference equations. I. A difference analogue of the Kortewegde Vries equation. J. Phys. Soc. Japan 434): [6] Michio Jimbo and Tetsuji Miwa. Monodromy preserving deformation of linear ordinary differential equations with rational coefficients. II. Physica D 23): [7] Michio Jimbo and Hidetaka Sakai. A q-analog of the sixth Painlevé equation. Lett. Math. Phys. 382): [8] N. Joshi F. W. Nijhoff and C. Ormerod. Lax pairs for ultra-discrete Painlevé cellular automata. J. Phys. A 3744):L559 L [9] I. M. Krichever and D. H. Phong. Symplectic forms in the theory of solitons. In Surveys in differential geometry: integral systems Surv. Differ. Geom. IV pages Int. Press Boston MA 998. [2] Mikio Murata. Lax forms of the q-painlevé equations. Journal of Physics A: Mathematical and Theoretical 42):52 7pp) 29. [2] F. W. Nijhoff. Lax pair for the Adler lattice Krichever-Novikov) system. Phys. Lett. A 297-2): [22] F. W. Nijhoff and V. G. Papageorgiou. Similarity reductions of integrable lattices and discrete analogues of the Painlevé II equation. Phys. Lett. A 536-7): [23] F. W. Nijhoff G. R. W. Quispel and H. W. Capel. Direct linearization of nonlinear differencedifference equations. Phys. Lett. A 974): [24] Frank W. Nijhoff and Hans Capel. The discrete Korteweg-de Vries equation. Acta Appl. Math. 39-3): KdV 95 Amsterdam 995). [25] Frank W. Nijhoff Andrew Hone and Nalini Joshi. On a Schwarzian PDE associated with the KdV hierarchy. Phys. Lett. A ):

23 DISCRETE PAINLEVÉ EQUATIONS AND THEIR LAX PAIRS AS REDUCTIONS OF INTEGRABLE LATTICE EQUATIONS 23 [26] Frank W. Nijhoff Nalini Joshi and Andrew Hone. On the discrete and continuous Miura chain associated with the sixth Painlevé equation. Phys. Lett. A 2645): [27] M. Noumi S. Tsujimoto and Y. Yamada. Padé interpolation for elliptic Painlevé equation. ArXiv e-prints April 22. arxiv: [28] Christopher M. Ormerod. Reductions of lattice mkdv to q-p V I. Phys. Lett. A 22. In Press. [29] V. G. Papageorgiou F. W. Nijhoff and H. W. Capel. Integrable mappings and nonlinear integrable lattice equations. Phys. Lett. A 472-3): [3] V. G. Papageorgiou F. W. Nijhoff B. Grammaticos and A. Ramani. Isomonodromic deformation problems for discrete analogues of Painlevé equations. Phys. Lett. A 64): [3] G. R. W. Quispel H. W. Capel V. G. Papageorgiou and F. W. Nijhoff. Integrable mappings derived from soliton equations. Physica A 73-2): [32] G. R. W. Quispel H. W. Capel and J. A. G. Roberts. Duality for discrete integrable systems. J. Phys. A 388): [33] G. R. W. Quispel F. W. Nijhoff H. W. Capel and J. van der Linden. Linear integral equations and nonlinear difference-difference equations. Phys. A 252-3): [34] G. R. W. Quispel J. A. G. Roberts and C. J. Thompson. Integrable mappings and soliton equations. Phys. Lett. A 267): [35] G. R. W. Quispel J. A. G. Roberts and C. J. Thompson. Integrable mappings and soliton equations. II. Physica D 34-2): [36] A. Ramani A. S. Carstea and B. Grammaticos. On the non-autonomous form of the Q 4 mapping and its relation to elliptic Painlevé equations. J. Phys. A 4232): [37] A. Ramani and B. Grammaticos. Discrete Painlevé equations: coalescences limits and degeneracies. Phys. A 228-4): [38] Olexandr G. Rasin and Peter E. Hydon. Symmetries of integrable difference equations on the quad-graph. Stud. Appl. Math. 93): [39] O. Rojas P.H. van der Kamp and G. R. W. Quispel. Lax representation for integrable O Es. Proceedings Symmetry And Perturbation Theory 27. pages Editted by R. Vitolo Published by World Scientific. [4] Hidetaka Sakai. Rational surfaces associated with affine root systems and geometry of the Painlevé equations. Comm. Math. Phys. 22): [4] Hidetaka Sakai. Lax form of the q-painlevé equation associated with the A ) 2 surface. J. Phys. A 3939): [42] Dinh T Tran Peter H van der Kamp and G R W Quispel. Involutivity of integrals of sinegordon modified KdV and potential KdV maps. Journal of Physics A: Mathematical and Theoretical 4429): [43] Peter H. van der Kamp. Initial value problems for lattice equations. J. Phys. A 424): [44] Peter H. van der Kamp and G. R. W. Quispel. The staircase method: integrals for periodic reductions of integrable lattice equations. J. Phys. A 4346): [45] A. P. Veselov and S. P. Novikov. Poisson brackets that are compatible with the algebraic geometry and the dynamics of the Korteweg-de Vries equation on the set of finite-gap potentials. Dokl. Akad. Nauk SSSR 2663): [46] Nicolas S. Witte and Christopher M. Ormerod. Construction of a Lax pair for the E ) 6 q-painlevé system. arxiv:27.4. [47] Yasuhiko Yamada. Lax formalism for q-painlevé equations with affine Weyl group symmetry of type E n ). Int. Math. Res. Not. IMRN 7): C. M. Ormerod La Trobe University Bundoora VIC 386 Australia address: C.Ormerod@latrobe.edu.au Peter H. van der Kamp La Trobe University Bundoora VIC 386 Australia address: P.vanderKamp@latrobe.edu.au G.R.W. Quispel La Trobe University Bundoora VIC 386 Australia address: R.Quispel@latrobe.edu.au