Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/0311029v1 [nlinsi] 16 Nov 2003 Symbolic Computation of Conserved Densities, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations Willy Hereman, Jan A Sanders, Jack Sayers, and Jing Ping Wang This paper is dedicated to Ryan Sayers (1982-2003) Abstract Algorithms for the symbolic computation of conserved densities, fluxes, generalized symmetries, and recursion operators for systems of nonlinear differential-difference equations are presented In the algorithms we use discrete versions of the Fréchet and variational derivatives, as well as discrete Euler and homotopy operators The algorithms are illustrated for prototypical nonlinear lattices, including the Kac-van Moerbeke (Volterra) and Toda lattices Results are shown for the modified Volterra and Ablowitz-Ladik lattices 1 Introduction The study of complete integrability of nonlinear differential-difference equations (DDEs) largely parallels that of PDEs [AC91, ASY00, W98] Indeed, as in the continuous case, the existence of sufficiently many conserved densities and generalized symmetries is a predictor for complete integrability Based on the first few densities and symmetries a recursion operator (which maps symmetries to symmetries) can be constructed The existence of a recursion operator, which allows one to generate an infinite sequence of symmetries, confirms complete integrability There is a vast body of work on the complete integrability of DDEs Consult eg [ASY00, HH03] for various approaches and references In this article we describe algorithms to symbolically compute polynomial conservation laws, fluxes, generalized symmetries, and recursion operators for DDEs The design of these algorithms heavily relies on related work for PDEs [GH97a, HG99, HGCM98] and work by Oevel et al [OZF93] There exists a close analogy between the continuous and discrete (in space) cases As shown in [HM02, MH02], this analogy can be completely formalized 1991 Mathematics Subject Classification Primary 37J35, 37K10; Secondary 35Q58, 37K05 Key words and phrases Homotopy operators, conserved densities, fluxes, recursion operators, integrability, differential-difference, DDEs, semi-discrete lattices This material is based upon work supported by the National Science Foundation under Grant No CCR-9901929 1 c 2003 Willy Hereman
2 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG andboththeoriescanbeformulatedintermsofcomplexes Thesameappliesforthe formulation in terms of Lie algebra complexes in[w98] This allows one to translate by analogy the existing algorithms immediately (although complications arise when there is explicit dependence on the space variable in the discrete case) One of the more useful tools following from the abstract theory is the homotopy operator, which is based on scaling vectorfields, and goes back to Poincaré in the continuous case This operator allows one to directly integrate differential forms and can be straightforwardly implemented in computer algebra packages, since it reduces to integration over one scaling parameter The discrete analogue as described in [HM02, MH02] does the corresponding job in the discrete case and we use it to compute fluxes In this paper we do not explicitly use the abstract framework, yet, it has been a motivating force for the development of our algorithms The algorithms in this paper can be implemented in any computer algebra system Our Mathematica package InvariantsSymmetriesm [GH97b] computes densities and symmetries, and therefore aids in automated testing of complete integrability of semi-discrete lattices Mathematica code to automatically compute recursion operators is still under development The paper is organized as follows In Section 2, we give key definitions and introduce the Kac-van Moerbeke(KvM)[KvM75], Toda[T81] and Ablowitz-Ladik (AL) lattices [AL75], which will be used as prototypical examples throughout the paper The discrete higher Euler operators(variational derivatives) and the discrete homotopy operator are introduced in Section 3 These operators are applied in the construction of densities and fluxes in Section 5 The algorithm for higher-order symmetries is outlined in Section 6 The algorithms for scalar and matrix recursion operators are given in Section 7 and Section 8 The paper concludes with two more examples in Section 9 2 Key Definitions and Prototypical Examples Definition 21 A nonlinear (autonomous) DDE is an equation of the form (21) = F(, 1,,+1,), where and F are vector-valued functions with m components The integer n corresponds to discretization in space; the dot denotes one derivative with respect to the continuous time variable (t) For simplicity, we will denote the components of by (,,w n, ) For brevity, we write F( ), although F typically depends on and a finite number of its forward and backward shifts We assume that F is polynomial with constant coefficients If present, parameters are denoted by lower-case Greek letters No restrictions are imposed on the level of the forward or backward shifts or the degree of nonlinearity in F (22) Example 22 The AL lattice [AL75], = (+1 2 + 1 )+ (+1 + 1 ), = (+1 2 + 1 ) (+1 + 1 ), is a completely integrable discretization of the nonlinear Schrödinger equation Definition 23 A DDE is said to be dilation invariant if it is invariant under a scaling (dilation) symmetry
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 3 Example 24 The KvM lattice [KvM75], (23) = (+1 1 ), is invariant under (t, ) (λ 1 t,λ ), where λ is an arbitrary scaling parameter Example 25 The Toda lattice [T81] in polynomial form [GH98], (24) = 1, = ( +1 ), is invariant under the scaling symmetry (25) (t,, ) (λ 1 t,λ,λ 2 ) Thus, and correspond to one, respectively two derivatives with respect to t, (26) d dt, d2 dt 2 Definition 26 We define the weight, w, of a variable as the number of t- derivatives the variable corresponds to Since t is replaced by t λ, we set w( d dt ) = 1 In view of (26), we have w() = 1, andw( ) = 2fortheTodalattice Weightsofdependentvariablesarenonnegative, rational, and independent of n Definition 27 The rank of a monomial is defined as the total weight of the monomial An expression is uniform in rank if all its terms have the same rank Example 28 In the first equation of (24), all the monomials have rank 2; in the second equation all the monomials have rank 3 Conversely, requiring uniformity inrankforeachequationin(24)allowsonetocomputetheweightsofthedependent variables with simple linear algebra Indeed, (27) w( )+1=w( ), w( )+1=w( )+w( ), yields (28) w( ) = 1, w( ) = 2, which is consistent with (25) Dilation symmetries, which are special Lie-point symmetries, are common to many lattice equations Lattices that do not admit a dilation symmetry can be made scaling invariant by extending the set of dependent variables using auxiliary parameters with scaling Example 29 The AL lattice is not dilation invariant Introducing an auxiliary parameter α, hence replacing (22) by (29) = α(+1 2 + 1 )+ (+1 + 1 ), = α(+1 2 + 1 ) (+1 + 1 ), and requiring uniformity in rank, gives (210) w( )+1 = w(α)+w( ) = 2w( )+w( ), w( )+1 = w(α)+w( ) = 2w( )+w( ) Obviously, (211) w( )+w( ) = w(α) = 1
4 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Several scales are possible The choice w( ) = w( ) = 1 2,w(α) = 1, corresponds to the scaling symmetry (212) (t,,,α) (λ 1 t,λ 1 2 un,λ 1 2 vn,λα) Definition 210 A scalar function ρ n ( ) is a conserved density of (21) if there exists a scalar function J n ( ) called the associated flux, such that (213) D t ρ n + J n = 0 is satisfied on the solutions of (21) [O93] In (213), we used the (forward) difference operator, = D I, defined by (214) J n = (D I)J n = J n+1 J n, where D denotes the up-shift (forward or right-shift) operator, DJ n = J n+1, and I is the identity operator Operator takes the role of a spatial derivative on the shifted variables as many DDEs arise from discretization of a PDE in (1 +1) variables Most, but not all, densities are polynomial in (215) (216) Example 211 The first three conservation laws for (23) are D t (ln( ))+(+1 + ) ( + 1 ) = 0, D t ( )+(+1 ) ( 1 ) = 0, (217) D t ( 1 2 u2 n ++1 )+ +1 (+1 ++2 ) 1 ( ++1 ) = 0 (218) (219) (220) (221) Example 212 For the Toda lattice (24) the first four density-flux pairs are ρ (0) n = ln( ), J n (0) n, ρ (1) n =, J n (1) n 1, ρ (2) n = 1 2 u2 n +, J (2) n = 1, ρ (3) n = 1 3 3 + ( 1 + ), J (3) n = 1 1 + 1 2 The above densities are uniform of ranks 0 through 3 The corresponding fluxes are also uniform in rank with ranks 1 through 4 In general, if in (213) rankρ n = R then rankj n = R+1, since w(d t ) = 1 This comes as no surprise since the conservation law (213) holds on solutions of (21), hence it inherits the dilation symmetry of (21) In Section 5 we will give an algorithm to compute polynomial conserved densities and fluxes and use(24) to illustrate the steps Non-polynomial densities(which are easy to find by hand) can be computed with the method given in [HH03] Definition 213 Compositions of D and D 1 define an equivalence relation ( ) on monomials in the components of All shifted monomials are equivalent Example 214 For example, 1 +1 +2 +4 3 1 Definition 215 The main representative of an equivalence class is the monomial of the class with n as lowest label on any component of u (or v if u is missing) Example 216 For example, +2 is the main representative of the class {, 2, 1 +1, +2,+1 +3, } Use lexicographical ordering to resolve conflicts For example, +2 (not 2 ) is the main representative of the class {, 3 1,,+2 +4, }
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 5 Definition 217 A vector function G( ) is called a generalized symmetry of (21) if the infinitesimal transformation +ǫg leaves (21) invariant up to order ǫ Consequently, G must satisfy [O93] (222) D t G = F ( )[G], on solutions of (21) F ( )[G] is the Fréchet derivative of F in the direction ofg For the scalar case (N = 1), the Fréchet derivative in the direction of G is (223) F ( )[G] = ǫ F( +ǫg) ǫ=0 = k which defines the Fréchet derivative operator (224) F ( ) = k F +k D k F +k D k G, In the vector case with two components, and, the Fréchet derivative operator is (225) F ( ) = Applied to G = (G 1 k F 1 +k D k k k F2 +k D k F 1 +k D k k F2 +k D k G 2 ) T, where T is transpose, one obtains (226) F i ( )[G] = k F i D k G 1 + +k k F i +k D k G 2, i = 1,2 In (223) and (226) summation is over all positive and negative shifts (including k = 0) For k > 0, D k = D D D (k times) Similarly, for k < 0 the down-shift operator D 1 is applied repeatedly The generalization of (225) to a system with N components should be obvious (227) (228) Example 218 The first two symmetries of (23) are G (1) = (+1 1 ), G (2) = +1 ( ++1 ++2 ) 1 ( 2 + 1 + ) These symmetries are uniform in rank (rank 2 and 3, respectively) The symmetries of ranks 0 and 1 are both zero Example 219 The first two non-trivial symmetries of (24) are ( ) (229) G (1) v = n 1, ( +1 ) and (230) G (2) = ( vn ( ++1 ) 1 ( 1 + ) (u 2 n+1 u 2 n ++1 1 ) The above symmetries are uniform in rank For example, rankg (2) 1 = 3 and rankg (2) 2 = 4 Symmetries of lower ranks are trivial An algorithm to compute generalized symmetries will be outlined in Section 6 and applied to (24) )
6 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Definition 220 A recursion operator R connects symmetries (231) G (j+s) = RG (j), where j = 1,2,, and s is the seed In most cases the symmetries are consecutively linked (s = 1) For N-component systems, R is an N N matrix operator The defining equation for R [O93, W98] is (232) D t R+[R,F ( )] = R t +R [F]+R F ( ) F ( ) R = 0, where [, ] denotes the commutator and the composition of operators The operator F ( ) was defined in (225) R [F] is the Fréchet derivative of R in the direction of F For the scalar case, the operator R is often of the form (233) R = U( )O((D I) 1,D 1,I,D)V( ), and then (234) R [F] = k (D k U F) OV + +k k UO(D k F) V +k For the vector case, the elements of the N N operator matrix R are often of the form R ij = U ij ( )O ij ((D I) 1,D 1,I,D)V ij ( ) For the two-component case R [F] ij = k (D k F 1 ) U ij +k O ij V ij + k (D k F 2 ) U ij +k O ij V ij (235) + k U ij O ij (D k F 1 ) V ij +k + k U ij O ij (D k F 2 ) V ij +k (236) Example 221 The KvM lattice (23) has recursion operator R = (I+D)( D D 1 )(D I) 1 1 I = D 1 +( ++1 )I+ D+ (+1 1 )(D I) 1 1 I Example 222 The Toda lattice (24) has recursion operator ( I D 1 I+( 1 )(D I) 1 1 (237) R = I D +1 I+ ( +1 )(D I) 1 1 I I In Section 7 we will give an algorithm for the computation of scalar recursion operators like (236) In Section 8 we cover the matrix case and show how (237) is computed The algorithms complement those for recursion operators of PDEs presented in [HG99] and elsewhere in these proceedings [BHS04] We now introduce two powerful tools which will be used in the computation of densities and fluxes 3 The Discrete Variational Derivative (Euler Operator) Definition 31 A function E n ( ) is a total difference if there exists another function J n ( ), such that E n = J n = (D I)J n Theorem 32 A necessary and sufficient condition for a function E n with positive shifts up to level p 0, to be a total difference is that (31) L (0) (E n ) = 0, )
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 7 where L (0) is the discrete variational derivative (Euler operator) [ASY00] defined by (32) L (0) = p 0 ( D k ) = (I+D 1 +D 2 + +D p0 ) k=0 A proof is given in eg [HH03] Remark 33 ToverifythatanexpressionE( q,,,,+p )involving negative shifts is a total difference, one must first remove the negative shifts by replacing E n by Ẽn = D q E n Applied to Ẽn, (32) terminates at p 0 = p+q We now introduce a tool to invert the total difference operator = D I 4 The Discrete Homotopy Operator Given is an expression E n (free of negative shifts) We assume that one has verified that E n KerL (0), ie L (0) (E n ) = 0 Consequently, E n Im So, J n such that E n = J n To compute J n = 1 (E n ) one must invert the operator = D I Working with the formal inverse, (41) 1 = D 1 +D 2 +D 3 +, is impractical, perhaps impossible We therefore present the (discrete) homotopy operator which circumvents the above infinite series In analogy to the continuous case [O93], we first introduce the discrete higher Euler operators Definition 41 The discrete higher Euler operators are defined by (42) L u (i) n = p i ( ) k ( D k ) i (43) (44) (45) The higher Euler operators all terminate at some maximal shifts p i k=i Example 42 For scalar component, the first higher Euler operators are: L (0) L (1) L (2) = (I+D 1 +D 2 +D 3 + +D p0 ), = (D 1 +2D 2 +3D 3 +4D 4 + +p 1 D p1 ), = (D 2 +3D 3 +6D 4 +10D 5 + + 1 2 p 2(p 2 1)D p2 ) Note that (43) coincides with (32) Similar formulae hold for L (i) Next, we introduce the homotopy operators For notational simplicity we show the formulae for the two component case = (, ) Definition 43 The total homotopy operator is defined as (46) H = 1 where the homotopy operators are given by 0 (j 1,n ( )[λ ]+j 2,n ( )[λ ]) dλ λ, p 1 1 p 2 1 (47) j 1,n ( ) = (D I) i ( L u (i+1) n ), j 2,n ( ) = (D I) i ( L u (i+1) n ) i=0 i=0
8 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Theorem 44 J n = 1 (E n ) can be computed as J n = H(E n ) A similar theorem and proof for continuous homotopy operators is given in [O93] By constructing a similar variational bicomplex the theorem still holds in the discrete case See [HM02, MH02] and elsewhere in these proceedings [M04] Remark 45 Since the theory for the recursion operator can be formulated in the language of Lie algebra complexes, the results in [SW01a] are immediately applicable This gives one explicit conditions from which the existence of infinite hierarchies of local symmetries and cosymmetries can be concluded Also, the problems with the definition of a recursion operator as treated in [SW01b] also play a role in the discrete case Remark 46 Note that p 1 and p 2 in (47) are the highest shifts for and in E n Furthermore, j r,n ( )[λ ] means that in j r,n ( ) one replaces λ, and, of course, +1 λ+1, +2 λ+2, etc Remark 47 In practice, one does not compute the definite integral(46) since the evaluation at boundary λ = 0 may cause problems Instead, one computes the indefinite integral and evaluates the result at λ = 1 5 Algorithm to Compute Densities and Fluxes As an example, we compute the density-flux pair (221) for (24) Assuming that the weights (28) are computed and the rank of the density is selected (R = 3 here), the algorithm has three steps Step 1: Construct the form of the density List all monomials in and of rank 3 or less: G={ 3, 2,,, } Next, for each monomial in G, introduce the correct number of t-derivatives so that each term has weight 3 Thus, using (24), d 0 dt 0( 3 ) = u 3 d 0 n, dt 0( ) =, d dt ( 2 d ) = 2 = 2 1 2, dt () = = +1, d 2 dt 2() = d dt ( ) = d dt ( 1 ) = 1 1 1 ++1 Gather all terms in H = { 3, 1,, 1 1,+1 } Identify members belonging to the same equivalence classes and replace them by their main representatives For example, 1 +1, so both are replaced by 1 So, H is replaced by I = { 3, 1, }, which has the building blocks of the density Linear combination of the monomials in I with constant coefficients c i gives (51) ρ n = c 1 3 +c 2 1 +c 3 Step 2: Determine the coefficients Require that (213) holds Compute D t ρ n Use (24) to remove and and their shifts Thus, (52) E n = D t ρ n = (3c 1 c 2 )u 2 n 1 +(c 3 3c 1 )u 2 n +(c 3 c 2 ) 1 +c 2 1 1 +c 2 v 2 n 1 c 3 +1 c 3 v 2 n
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 9 Compute Ẽn = DE n First, apply (32) for component to Ẽn : (53) L (0) (Ẽn) = (I+D 1 +D 2 )(Ẽn) Second, apply (32) for component to Ẽn : (54) = 2(3c 1 c 2 ) 1 +2(c 3 3c 1 ) +(c 2 c 3 ) 1 1 +(c 2 c 3 )+1 L (0) (Ẽn) = (I+D 1 )(Ẽn) = (3c 1 c 2 )u 2 n+1 +(c 3 c 2 )+1 +(c 2 c 3 ) +1 +2(c 2 c 3 ) +(c 3 3c 1 )u 2 n +(c 3 c 2 ) 1 Both (53) and (54) must vanish identically Solve the linear system (55) S = {3c 1 c 2 = 0,c 3 3c 1 = 0,c 2 c 3 = 0} The solution is 3c 1 = c 2 = c 3 Substituting c 1 = 1 3,c 2 = c 3 = 1, into (51) (56) ρ n = 1 3 3 + ( 1 + ) Step 3: Compute the flux In view of (213), we will compute J n = 1 ( E n ) with the homotopy operator introduced in Section 4 Alternative methods are described in [GH98, HH03] Insert c 1 = 1 3,c 2 = c 3 = 1 into (52) and compute (57) Ẽ n = DE n = +1 +v 2 n +1+2 +1 v 2 n+1 We apply formulae (47) to Ẽn The pieces are listed in Tables 1 and 2 i ( Ẽn) (D I) i ( L (i+1) ( Ẽn)) L (i+1) 0 1 1 ++1 1 1 + +1 1 1 1 +1 1 1 Table 1 Computation of j 1,n ( )( Ẽn) i L (i+1) ( Ẽn) (D I) i ( L v (i+1) n ( Ẽn)) 0 +1 +2 +1 +2v 2 n Table 2 Computation of j 2,n ( )( Ẽn) Adding the terms in the right columns in Table 1 and Table 2, (58) j 1,n ( )( Ẽn) = 2 +1, j 2,n ( )( Ẽn) = +1 +2vn 2 Thus, the homotopy operator (46) gives J n = = 1 0 1 (59) = +1 +vn 2 0 (j 1,n ( )( Ẽn)[λ ]+j 2,n ( )( Ẽn)[λ ]) dλ λ (3λ 2 +1 +2λv 2 n ) dλ
10 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG After a backward shift, J n = D 1 Jn, we obtain the final result: (510) ρ n = 1 3 u3 n + ( 1 + ), J n = 1 1 +v 2 n 1 6 Algorithm to Compute Generalized Symmetries As an example, we compute the symmetry (230) of rank (3,4) for (24) The two steps of the algorithm [GH99] are similar to those in Section 5 Step 1: Construct the form of the symmetry Start by listing all monomials in and of ranks 3 and 4, or less: K 1 = {u 3 n,u2 n,,, }, K 2 = {u 4 n,u3 n,u2 n,u 2 n,,,v 2 n,} As in Step 1 in Section 5, for each monomial in K 1 and K 2, introduce the necessary t-derivatives so that each term has rank 3 and 4, respectively At the same time, use (24) to remove all t derivatives Doing so, based on K 1 we obtain (61) L 1 = {u 3 n, 1 1, 1,,+1 } Similarly, based on K 2, we get (62) L 2 = {u 4 n,u2 n 1 1, 1 1,u 2 n 1, 2 1,v 2 n 1,u2 n, +1,u 2 n+1, 1,v 2 n, +1 } In contrast to the strategy for densities, we do not introduce the main representatives, but linearly combine the monomials in L 1 and L 2 to get the form of the symmetry: (63) G 1 = c 1 u 3 n +c 2 1 1 +c 3 1 +c 4 +c 5 +1, G 2 = c 6 u 4 n +c 7 u 2 n 1 1 +c 8 1 1 +c 9 u 2 n 1 +c 10 2 1 +c 11 v 2 n 1 +c 12 u 2 n +c 13 +1 +c 14 u 2 n+1 +c 15 1 +c 16 v 2 n +c 17 +1, with constant coefficients c i Step 2: Determine the coefficients To compute the coefficients c i we require that (222) holds on solutions of (21) Compute D t G 1 and D t G 2 and use (24) to remove,, and their shifts This gives the left hand sides of (222) Use (226) to get the right hand sides of (222): (64) F 1 ( )[G] = D 1 G 2 IG 2, F 2 ( )[G] = IG 1 DG 1 +( +1 )IG 2 Substitute (63) into (64) and equate the corresponding left and right hand sides Since all monomials in, and their shifts are independent, one obtains the linear system that determines the coefficients c i The solution is (65) c 1 = c 6 = c 7 = c 8 = c 9 = c 10 = c 11 = c 13 = c 16 = 0, c 2 = c 3 = c 4 = c 5 = c 12 = c 14 = c 15 = c 17 With the choice c 17 = 1, the symmetry (63) finally becomes (66) G 1 = ( ++1 ) 1 ( 1 ), G 2 = (u 2 n+1 u2 n ++1 1 )
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 11 7 Algorithm to Compute Scalar Recursion Operators In this section we show how to compute the recursion operator (236) of (24) Again, we will use the concept of rank invariance to construct a candidate recursion operator The defining equation (232) is then used to determine the coefficients We observe that (236) in expanded form naturally splits into two pieces: (71) R = R 0 +R 1, where R 0 contains only terms with shift operators D 1,I, and D, and R 1 has terms involving (D I) 1 Step 1: Determine the rank of the recursion operator Inviewof(231)andassumingthat thesymmetriesarelinkedconsecutively(s = 1), the recursion operator R has rank (72) R = rank R = rank G (2) rank G (1) = 3 2 = 1, where we used (227) and (228) If the assumption turns out to be correct, then the recursion operator has rank 1 If the assumption were falls because symmetries are not linked consecutively, then R must be adjusted and the subsequent steps must be repeated See [BHS04] for examples of PDEs for which that happens Step 2: Determine the form of the recursion operator We split this into two sub-steps (i) Determine the form of R 0 The candidate R 0 is a sum of terms involving D 1,I, and D The coefficients of these terms are linear combinations of 1,, and +1, so that all terms have the correct rank (73) R 0 = (c 1 1 +c 2 +c 3 +1 )D 1 +(c 4 1 +c 5 +c 6 +1 )I +(c 7 1 +c 8 +c 9 +1 )D, where the c i are constant coefficients A few remarks are in place First, in R 0 we moved the operators D 1,I, and D all the way to the right Second, the maximum up-shift and down-shift operator that should be included can be determined by comparing two consecutive symmetries Indeed, if the maximum up-shift in the first symmetry is +p, and the maximum up-shift in the next symmetry is +p+r, then R 0 must have D,D 2,,D r The same line of reasoning determines the minimum down-shift operator to be included In our example, there is no need to include terms in D 2,D 2, etc Third, the coefficients of the operators can be restricted to linear combinations of the terms appearing in F Hence, no terms in ±2,±3 and so on occur in (73) (ii) Determine the form of R 1 As in the continuous case [HG99], R 1 is a linear combination (with constant coefficients c jk ) of sums of all suitable products of symmetries and covariants (Fréchet derivatives of densities) sandwiching (D I) 1, ie (74) R 1 = j k c jk G (j) (D I) 1 ρ (k) n For (23), G (1) in (227) and ρ (0) n = ln( ) in (215) are the only suitable pair Indeed, using (224) we have ρ (0) n = ln( ) = 1 I, which has rank -1 Combined
12 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG with G (1) of rank 2, we have a term of rank 1 Other combinations of symmetries and covariants would exceed rank 1 Therefore, (75) R 1 = c 10 (+1 1 )(D I) 1 ( 1 )I where c 10 a constant coefficient Using (71) and renaming c 10 to c 10, R = (c 1 1 +c 2 +c 3 +1 )D 1 +(c 4 1 +c 5 +c 6 +1 )I (76) +(c 7 1 +c 8 +c 9 +1 )D+c 10 (+1 1 )(D I) 1 ( 1 )I is the candidate recursion operator for (23) Step 3: Determine the coefficients Starting from (76), we use (234) with F = (+1 1 ) to compute R [F] The partial results will not be shown due to length Using (223), we compute (77) F = D 1 +(+1 1 )I+ D Then we compute R F and F R After substituting the pieces into (232) we simplify the resulting expression using rules such as (78) (D I) 1 D = D(D I) 1 = I+(D I) 1, (D I) 1 D 1 = D 1 (D I) 1 = D 1 +(D I) 1 We further simplify by recursively using formulas like (79) DU( )(D I) 1 V( )I = U(+1 )V( )I+U(+1 )(D I) 1 V( )I, (D I) 1 U( )V( )D = U( 1 )V( 1 )I+(D I) 1 U( 1 )V( 1 )I Finally, we equate like terms to obtain a linear system for the c i Substituting the solution (710) c 1 = c 3 = c 4 = c 7 = c 9 = 0, c 2 = c 5 = c 6 = c 8 = c 10 = 1 into (76) we obtain the final result (711) R = D 1 +( ++1 )I+ D+ (+1 1 )(D I) 1 1 I A straightforward computation confirms that RG (1) = G (2) with G (1) in (227) and G (2) in (228) 8 Algorithm to Compute Matrix Recursion Operators We construct the recursion operator(237) for(24) Now all the terms in (232) are 2 2 matrix operators Step 1: Determine the rank of the recursion operator The difference in rank of symmetries is again used to compute the rank of the elements of the recursion operator Using (28), (229) and (230), ( (81) rankg (1) 2 = 3 ), rankg (2) = Assuming that RG (1) = G (2), we use the formula (82) rankr ij = rankg (k+1) i ( 3 4 rankg (k) j, )
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 13 to compute a rank matrix associated to the operator ( ) 1 0 (83) rank R = 2 1 Step 2: Determine the form of the recursion operator As in the scalar case, we build a candidate R 0 : ( ) (R0 ) R 0 = 11 (R 0 ) (84) 12, (R 0 ) 21 (R 0 ) 22 with (R 0 ) 11 = (c 1 +c 2 +1 )I, (R 0 ) 12 = c 3 D 1 +c 4 I, (R 0 ) 21 = (c 5 u 2 n +c 6 +1 +c 7 u 2 n+1 +c 8 1 +c 9 )I +(c 10 u 2 n +c 11 +1 +c 12 u 2 n+1 +c 13 1 +c 14 )D, (R 0 ) 22 = (c 15 +c 16 +1 )I Analogous to the scalar case, the elements of matrix R 1 are linear combinations with constant coefficients of all suitable products of symmetries and covariants sandwiching (D I) 1 Hence, (85) j k c jk G (j) (D I) 1 ρ (k) n, where denotes the matrix outer product, defined as ( ) ( G (j) ) 1 G (j) (D I) 1 (ρ (k) n,1 ρ (k) n,2 = 2 G (j) 1 (D I) 1 ρ (k) n,1 G (j) 1 (D I) 1 ρ (k) n,2 G (j) 2 (D I) 1 ρ (k) n,1 G (j) 2 (D I) 1 ρ (k) n,2 Only the pair (G (1),ρ (0) n ) can be used, otherwise the ranks in (83) would be exceeded Using (226) and (218) we compute (86) ρ (0) n = ( 1 0 I ), Therefore, using (85) and renaming c 10 to c 17, ( 0 c17 ( 1 )(D I) 1 1 (87) R 1 = I 0 c 17 ( +1 )(D I) 1 1 I Adding (84) and (87) we obtain ( ) R11 R (88) R = 12, R 21 R 22 with R 11 = (c 1 +c 2 +1 )I ) R 12 = c 3 D 1 +c 4 I+c 17 ( 1 )(D I) 1 1 I, R 21 = (c 5 u 2 n +c 6 +1 +c 7 u 2 n+1 +c 8 1 +c 9 )I +(c 10 u 2 n +c 11 +1 +c 12 u 2 n+1 +c 13 1 +c 14 )D, R 22 = (c 15 +c 16 +1 )I+c 17 ( +1 )(D I) 1 1 I )
14 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Step 3: Determine the coefficients All the terms in (232) need to be computed The strategy is similar to the scalar case, yet the computations are much more cumbersome Omitting the details, the result is: c 2 = c 5 = c 6 = c 7 = c 8 = c 10 = c 11 = c 12 = c 13 = c 15 = 0, c 1 = c 3 = c 4 = c 9 = c 14 = c 16 = 1, and c 17 = 1 Substitution the constants into (88) gives ( I D 1 I+( 1 )(D I) 1 1 ) (89) R = I I D +1 I+ ( +1 )(D I) 1 1 It is straightforward to verify that RG (1) = G (2) with G (1) in (229) and G (2) in (230) 9 More Examples Example 91 The modified Volterra lattice [ASY00, HH03], (91) = u 2 n (+1 1 ), has two non-polynomial densities ρ (0) n = 1 and ρ (1) n = ln( ), and infinitely many polynomial densities The first two symmetries, (92) (93) G (1) = u 2 n(+1 1 ), G (2) = u 2 nu 2 n+1( ++2 ) u 2 n 1u 2 n( 2 + ), are linked by the recursion operator (94) R = u 2 n D 1 +2 +1 I+u 2 n D+2u2 n (+1 1 )(D I) 1 1 I Example 92 The AL lattice (22) has infinitely many densities [GH98] and symmetries [GH99] The recursion operator is of the form (88) with R 11 = P n D 1 1 +1 I 1 P n 1 P n I, R 12 = 1 I 1 1 I 1 P n 1 P n I, R 21 = +1 I+ 1 +1 I++1 P n 1 P n I, R 22 = ( +1 + 1 )I+P n D+ 1 1 I++1 P n 1 (95) I, P n where P n = 1+ and = D I This recursion operator has an inverse, which is quite exceptional The investigation of the recursion operator structure of the AL lattice is work in progress Acknowledgements M Hickman and B Deconinck are gratefully acknowledged for valuable discussions D Baldwin is thanked for proof reading the manuscript References [AC91] M J Ablowitz and P A Clarkson, Solitons, nonlinear evolution equations and inverse scattering, London Math Soc Lect Note Ser 149, Cambridge University Press, Cambridge, U K, 1991 [AL75] M J Ablowitz and J F Ladik, Nonlinear differential-difference equations, J Math Phys 16 (1975), 598 603 [ASY00] V E Adler, A B Shabat and R I Yamilov, Symmetry approach to the integrability problem, Theor Math Phys 125 (2000), 1603 1661 I
DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 15 [BHS04] D Baldwin, W Hereman and J Sayers, Symbolic algorithms for the Painlevé test, special solutions and recursion operators for nonlinear PDEs, present proceedings [GH97a] Ü Göktaş and W Hereman, Symbolic computation of conserved densities for systems of nonlinear evolution equations, J Symb Comp 24 (1997), 591 621 [GH97b], Mathematica package InvariantsSymmetriesm is available since 1997 from http://wwwmathsourcecom/cgi-bin/msitem?0208-932 [GH98], Computation of conserved densities for nonlinear lattices, Physica D 123 (1998), 425 436 [GH99], Algorithmic computation of higher-order symmetries for nonlinear evolution and lattice equations, Adv Comp Math 11 (1999), 55 80 [HG99] W Hereman and Ü Göktaş, Integrability Tests for Nonlinear Evolution Equations Computer Algebra Systems: A Practical Guide, (M Wester, ed), Wiley, New York, 1999, pp 211 232 [HGCM98] W Hereman, Ü Göktaş, M D Colagrosso and A J Miller, Algorithmic integrability tests for nonlinear differential and lattice equations, Comp Phys Comm 115 (1998), 428 446 [HH03] M S Hickman and W A Hereman, Computation of densities and fluxes of nonlinear differential-difference equations, Proc Roy Soc Lond A (2003), in press [HM02] P E Hydon and E L Mansfield, A variational complex for difference equations, Found Comp Math (2003), in press [KvM75] M Kac and P van Moerbeke, On an explicitly soluble system of nonlinear differential equations related to certain Toda lattices, Adv Math 16 (1975), 160 169 [MH02] E L Mansfield and P E Hydon, On a variational complex for difference equations The Geometrical Study of Differential Equations Proc NSF-CBMS Conf (J A Leslie and T P Robart, eds), American Mathematical Society, Providence, Rhode Island, 2002, pp 121 129 [M04] E L Mansfield and R Quispel, Towards a variational complex for the finite element method, present proceedings [OZF93] W Oevel, H Zhang and B Fuchssteiner, Mastersymmetries and multi-hamiltonian formulations for some integrable lattice systems, Prog Theor Phys 81 (1989), 294 308 [O93] P J Olver, Applications of lie groups to differential equations, 2nd ed, Graduate Texts in Mathematics 107, Springer, New York, 1993 [SW01a] J A Sanders and J P Wang, Integrable systems and their recursion operators, Nonlinear Anal 47 (2001), 5213 5240 [SW01b], On recursion operators, Physica D 149 (2001), 1 10 [T81] M Toda, Theory of nonlinear lattices, Springer, Berlin, 1981 [W98] J P Wang, Symmetries and conservation laws of evolution equations, PhD Thesis, Thomas Stieltjes Institute for Mathematics, Amsterdam, 1998 Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, Colorado 80401-1887 E-mail address: whereman@minesedu Department of Mathematics, Faculty of Sciences, Free University, De Boelelaan 1081a, 1081 HV Amsterdam, The Netherlands E-mail address: jansa@csvunl Physics Department 103-33, California Institute of Technology, Pasadena, California 91125 E-mail address: jack@itscaltechedu Department of Mathematics, Brock University, St Catharines, Ontario, Canada L2S 3A1 E-mail address: jwang@spartanacbrockuca