arxiv:nlin/ v1 [nlin.si] 16 Nov 2003
|
|
- Wilfred Sullivan
- 5 years ago
- Views:
Transcription
1 Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/ v1 [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 ( ) 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 c 2003 Willy Hereman
2 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], = ( )+ ( ), = ( ) ( ), 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
3 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) = α( )+ ( ), = α( ) ( ), 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 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 ( ++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 + ), J (3) n = 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, 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, }
5 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 ( ) 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 ) 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 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, )
7 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 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 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 ) = 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 c 2 v 2 n 1 c 3 +1 c 3 v 2 n
9 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 + ) 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 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) Table 1 Computation of j 1,n ( )( Ẽn) i L (i+1) ( Ẽn) (D I) i ( L v (i+1) n ( Ẽn)) v 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 = = (59) = +1 +vn 2 0 (j 1,n ( )( Ẽn)[λ ]+j 2,n ( )( Ẽn)[λ ]) dλ λ (3λ λv 2 n ) dλ
10 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 c 3 1 +c 4 +c 5 +1, G 2 = c 6 u 4 n +c 7 u 2 n 1 1 +c c 9 u 2 n 1 +c c 11 v 2 n 1 +c 12 u 2 n +c c 14 u 2 n+1 +c 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 )
11 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 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, )
13 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 c 7 u 2 n+1 +c 8 1 +c 9 )I +(c 10 u 2 n +c c 12 u 2 n+1 +c c 14 )D, (R 0 ) 22 = (c 15 +c )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 c 7 u 2 n+1 +c 8 1 +c 9 )I +(c 10 u 2 n +c c 12 u 2 n+1 +c c 14 )D, R 22 = (c 15 +c )I+c 17 ( +1 )(D I) 1 1 I )
14 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 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 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 I++1 P n 1 P n I, R 22 = ( )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), [ASY00] V E Adler, A B Shabat and R I Yamilov, Symmetry approach to the integrability problem, Theor Math Phys 125 (2000), I
15 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), [GH97b], Mathematica package InvariantsSymmetriesm is available since 1997 from [GH98], Computation of conserved densities for nonlinear lattices, Physica D 123 (1998), [GH99], Algorithmic computation of higher-order symmetries for nonlinear evolution and lattice equations, Adv Comp Math 11 (1999), [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 [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), [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), [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 [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), [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), [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 address: whereman@minesedu Department of Mathematics, Faculty of Sciences, Free University, De Boelelaan 1081a, 1081 HV Amsterdam, The Netherlands address: jansa@csvunl Physics Department , California Institute of Technology, Pasadena, California address: jack@itscaltechedu Department of Mathematics, Brock University, St Catharines, Ontario, Canada L2S 3A1 address: jwang@spartanacbrockuca
Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations
Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations Ünal Göktaş and Willy Hereman Department of Computer Engineering Turgut Özal University, Keçiören, Ankara, Turkey
More informationSymbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations
Symbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations Ünal Götaş and Willy Hereman Abstract Algorithms for the symbolic
More informationContinuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman
Continuous and Discrete Homotopy Operators with Applications in Integrability Testing Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs
More informationThe first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3
Invariants and Symmetries for Partial Differential Equations and Lattices Λ Ünal Göktaοs y Willy Hereman y Abstract Methods for the computation of invariants and symmetries of nonlinear evolution, wave,
More informationSymbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN
. Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN Joint work with ÜNAL GÖKTAŞ and Grant Erdmann Graduate Seminar Department
More informationSymbolic Computation of. Conserved Densities and Fluxes for. Nonlinear Systems of Differential-Difference. Equations
Symbolic Computation of Conserved Densities and Fluxes for Nonlinear Systems of Differential-Difference Equations by Holly Eklund Copyright c 003 Holly Eklund, All rights reserved. ii A thesis submitted
More informationarxiv:nlin/ v1 [nlin.si] 16 Nov 2003
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/0311028v1 [nlin.si] 16 Nov 2003 Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for
More informationON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS
Proceedings of the International Conference on Difference Equations, Special Functions and Orthogonal Polynomials, World Scientific (2007 ON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS ANASTASIOS
More informationRESEARCH ARTICLE. A symbolic algorithm for computing recursion operators of nonlinear PDEs
International Journal of Computer Mathematics Vol. 00, No. 00, Month 200, 1 27 RESEARCH ARTICLE A symbolic algorithm for computing recursion operators of nonlinear PDEs D.E. Baldwin and W. Hereman Department
More informationApplication of Euler and Homotopy Operators to Integrability Testing
Application of Euler and Homotopy Operators to Integrability Testing Willy Hereman Department of Applied Mathematics and Statistics Colorado School of Mines Golden, Colorado, U.S.A. http://www.mines.edu/fs
More informationAlgorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of. Nonlinear Evolution and Differential-Difference.
Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of Nonlinear Evolution and Differential-Difference Equations by Ünal GÖKTAŞ A thesis submitted to the Faculty and the
More informationSymbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs.
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs. Douglas Baldwin, Willy Hereman,
More informationAbout Integrable Non-Abelian Laurent ODEs
About Integrable Non-Abelian Laurent ODEs T. Wolf, Brock University September 12, 2013 Outline Non-commutative ODEs First Integrals and Lax Pairs Symmetries Pre-Hamiltonian Operators Recursion Operators
More informationSymmetries and Group Invariant Reductions of Integrable Partial Difference Equations
Proceedings of 0th International Conference in MOdern GRoup ANalysis 2005, 222 230 Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations A. TONGAS, D. TSOUBELIS and V. PAPAGEORGIOU
More informationINTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION
v..1r0180507 *018.11.15#5f9cb4 INTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION CORINA N. BABALIC University of Craiova 13 A.I. Cuza, 00585, Craiova, Romania E-mail:
More informationSolving Nonlinear Wave Equations and Lattices with Mathematica. Willy Hereman
Solving Nonlinear Wave Equations and Lattices with Mathematica Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs home/whereman/
More information2. Examples of Integrable Equations
Integrable equations A.V.Mikhailov and V.V.Sokolov 1. Introduction 2. Examples of Integrable Equations 3. Examples of Lax pairs 4. Structure of Lax pairs 5. Local Symmetries, conservation laws and the
More informationFibonacci tan-sec method for construction solitary wave solution to differential-difference equations
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 7 (2011) No. 1, pp. 52-57 Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations
More informationThe Complete Set of Generalized Symmetries for the Calogero Degasperis Ibragimov Shabat Equation
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 209 214 The Complete Set of Generalized Symmetries for the Calogero Degasperis Ibragimov Shabat Equation Artur SERGYEYEV
More informationA MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE
Pacific Journal of Applied Mathematics Volume 1, Number 2, pp. 69 75 ISSN PJAM c 2008 Nova Science Publishers, Inc. A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE Wen-Xiu Ma Department
More informationarxiv: v2 [math-ph] 24 Feb 2016
ON THE CLASSIFICATION OF MULTIDIMENSIONALLY CONSISTENT 3D MAPS MATTEO PETRERA AND YURI B. SURIS Institut für Mathemat MA 7-2 Technische Universität Berlin Str. des 17. Juni 136 10623 Berlin Germany arxiv:1509.03129v2
More informationarxiv: v1 [nlin.si] 25 Mar 2009
Linear quadrilateral lattice equations and multidimensional consistency arxiv:0903.4428v1 [nlin.si] 25 Mar 2009 1. Introduction James Atkinson Department of Mathematics and Statistics, La Trobe University,
More informationB.7 Lie Groups and Differential Equations
96 B.7. LIE GROUPS AND DIFFERENTIAL EQUATIONS B.7 Lie Groups and Differential Equations Peter J. Olver in Minneapolis, MN (U.S.A.) mailto:olver@ima.umn.edu The applications of Lie groups to solve differential
More informationOn Local Time-Dependent Symmetries of Integrable Evolution Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2000, Vol. 30, Part 1, 196 203. On Local Time-Dependent Symmetries of Integrable Evolution Equations A. SERGYEYEV Institute of Mathematics of the
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions
. Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden,
More informationBreaking soliton equations and negative-order breaking soliton equations of typical and higher orders
Pramana J. Phys. (2016) 87: 68 DOI 10.1007/s12043-016-1273-z c Indian Academy of Sciences Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders ABDUL-MAJID
More informationSymmetry Methods for Differential and Difference Equations. Peter Hydon
Lecture 2: How to find Lie symmetries Symmetry Methods for Differential and Difference Equations Peter Hydon University of Kent Outline 1 Reduction of order for ODEs and O Es 2 The infinitesimal generator
More informationCanonical Forms for BiHamiltonian Systems
Canonical Forms for BiHamiltonian Systems Peter J. Olver Dedicated to the Memory of Jean-Louis Verdier BiHamiltonian systems were first defined in the fundamental paper of Magri, [5], which deduced the
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 6: Scattered Data Interpolation with Polynomial Precision Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationThe Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations
Symmetry in Nonlinear Mathematical Physics 1997, V.1, 185 192. The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations Norbert EULER, Ove LINDBLOM, Marianna EULER
More informationThe Modified (G /G)-Expansion Method for Nonlinear Evolution Equations
The Modified ( /-Expansion Method for Nonlinear Evolution Equations Sheng Zhang, Ying-Na Sun, Jin-Mei Ba, and Ling Dong Department of Mathematics, Bohai University, Jinzhou 11000, P. R. China Reprint requests
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions
Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions Willy Hereman and Douglas Poole Department of Mathematical and Computer Sciences Colorado
More informationEuler characteristic of the truncated order complex of generalized noncrossing partitions
Euler characteristic of the truncated order complex of generalized noncrossing partitions D. Armstrong and C. Krattenthaler Department of Mathematics, University of Miami, Coral Gables, Florida 33146,
More informationarxiv:nlin/ v2 [nlin.si] 9 Oct 2002
Journal of Nonlinear Mathematical Physics Volume 9, Number 1 2002), 21 25 Letter On Integrability of Differential Constraints Arising from the Singularity Analysis S Yu SAKOVICH Institute of Physics, National
More informationarxiv:nlin/ v2 [nlin.si] 15 Sep 2004
Integrable Mappings Related to the Extended Discrete KP Hierarchy ANDREI K. SVININ Institute of System Dynamics and Control Theory, Siberian Branch of Russian Academy of Sciences, P.O. Box 1233, 664033
More informationSymmetry Classification of KdV-Type Nonlinear Evolution Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 1, 125 130 Symmetry Classification of KdV-Type Nonlinear Evolution Equations Faruk GÜNGÖR, Victor LAHNO and Renat ZHDANOV Department
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationCompatible Hamiltonian Operators for the Krichever-Novikov Equation
arxiv:705.04834v [math.ap] 3 May 207 Compatible Hamiltonian Operators for the Krichever-Novikov Equation Sylvain Carpentier* Abstract It has been proved by Sokolov that Krichever-Novikov equation s hierarchy
More informationCLASSIFICATION OF NON-ABELIAN CHERN-SIMONS VORTICES
CLASSIFICATION OF NON-ABELIAN CHERN-SIMONS VORTICES arxiv:hep-th/9310182v1 27 Oct 1993 Gerald V. Dunne Department of Physics University of Connecticut 2152 Hillside Road Storrs, CT 06269 USA dunne@hep.phys.uconn.edu
More informationComplete integrability of a modified vector derivative nonlinear Schrödinger equation
Complete integrability of a modified vector derivative nonlinear Schrödinger equation arxiv:solv-int/9407002v1 14 Jul 1994 Ralph Willox 1, Willy Hereman 2 and Frank Verheest 3 1 Dienst Theoretische Natuurkunde,
More informationEvolutionary Hirota Type (2+1)-Dimensional Equations: Lax Pairs, Recursion Operators and Bi-Hamiltonian Structures
Symmetry, Integrability and Geometry: Methods and Applications Evolutionary Hirota Type +-Dimensional Equations: Lax Pairs, Recursion Operators and Bi-Hamiltonian Structures Mikhail B. SHEFTEL and Devrim
More informationLax Representations and Zero Curvature Representations by Kronecker Product
Lax Representations and Zero Curvature Representations by Kronecker Product arxiv:solv-int/9610008v1 18 Oct 1996 Wen-Xiu Ma and Fu-Kui Guo Abstract It is showed that Kronecker product can be applied to
More informationTHE NUMBER OF ORTHOGONAL CONJUGATIONS
THE NUMBER OF ORTHOGONAL CONJUGATIONS Armin Uhlmann University of Leipzig, Institute for Theoretical Physics After a short introduction to anti-linearity, bounds for the number of orthogonal (skew) conjugations
More informationarxiv:nlin/ v1 [nlin.si] 13 Apr 2005
Diophantine Integrability R.G. Halburd Department of Mathematical Sciences, Loughborough University Loughborough, Leicestershire, LE11 3TU, UK (Dated: received January 1, revised March 7, 005) The heights
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationOn bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy
arxiv:nlin/3139v2 [nlin.si] 14 Jan 24 On bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy S. Yu. Sakovich Mathematical Institute, Silesian University, 7461 Opava, Czech
More informationExact solutions of nonlinear partial differential equations by the method of group foliation reduction
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany STEPHEN S. ANCO 1, S. ALI, THOMAS WOLF 3 Exact solutions of nonlinear partial differential equations by the
More informationA new integrable system: The interacting soliton of the BO
Phys. Lett., A 204, p.336-342, 1995 A new integrable system: The interacting soliton of the BO Benno Fuchssteiner and Thorsten Schulze Automath Institute University of Paderborn Paderborn & Germany Abstract
More informationMATH 320, WEEK 7: Matrices, Matrix Operations
MATH 320, WEEK 7: Matrices, Matrix Operations 1 Matrices We have introduced ourselves to the notion of the grid-like coefficient matrix as a short-hand coefficient place-keeper for performing Gaussian
More informationOn a WKB-theoretic approach to. Ablowitz-Segur's connection problem for. the second Painleve equation. Yoshitsugu TAKEI
On a WKB-theoretic approach to Ablowitz-Segur's connection problem for the second Painleve equation Yoshitsugu TAKEI Research Institute for Mathematical Sciences Kyoto University Kyoto, 606-8502, JAPAN
More informationarxiv: v1 [math.co] 13 May 2016
GENERALISED RAMSEY NUMBERS FOR TWO SETS OF CYCLES MIKAEL HANSSON arxiv:1605.04301v1 [math.co] 13 May 2016 Abstract. We determine several generalised Ramsey numbers for two sets Γ 1 and Γ 2 of cycles, in
More informationDiscriminative Direction for Kernel Classifiers
Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering
More informationOn recursion operators for elliptic models
On recursion operators for elliptic models D K Demskoi and V V Sokolov 2 School of Mathematics and Statistics, University of New South Wales, Sydney, NSW 2052, Australia E-mail: demskoi@maths.unsw.edu.au
More informationA Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations
A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations Kurt Schlacher, Andreas Kugi and Kurt Zehetleitner kurt.schlacher@jku.at kurt.zehetleitner@jku.at,
More informationarxiv: v1 [math.ap] 11 Jan 2014
THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for
More informationarxiv:math-ph/ v1 25 Feb 2002
FROM THE TODA LATTICE TO THE VOLTERRA LATTICE AND BACK arxiv:math-ph/0202037v1 25 Feb 2002 (1) PANTELIS A DAMIANOU AND RUI LOJA FERNANDES Abstract We discuss the relationship between the multiple Hamiltonian
More informationSymmetry Reductions of Integrable Lattice Equations
Isaac Newton Institute for Mathematical Sciences Discrete Integrable Systems Symmetry Reductions of Integrable Lattice Equations Pavlos Xenitidis University of Patras Greece March 11, 2009 Pavlos Xenitidis
More informationBäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations
arxiv:nlin/0102001v1 [nlin.si] 1 Feb 2001 Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations Ayşe Karasu (Kalkanli) and S Yu Sakovich Department
More informationON POINCARE POLYNOMIALS OF HYPERBOLIC LIE ALGEBRAS
ON POINCARE POLYNOMIALS OF HYPERBOLIC LIE ALGEBRAS Meltem Gungormez Dept. Physics, Fac. Science, Istanbul Tech. Univ. 34469, Maslak, Istanbul, Turkey e-mail: gungorm@itu.edu.tr arxiv:0706.2563v2 [math-ph]
More informationData representation and classification
Representation and classification of data January 25, 2016 Outline Lecture 1: Data Representation 1 Lecture 1: Data Representation Data representation The phrase data representation usually refers to the
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationSymmetry reductions and travelling wave solutions for a new integrable equation
Symmetry reductions and travelling wave solutions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX 0, 50 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationMoving Boundary Problems for the Harry Dym Equation & Reciprocal Associates
Moving Boundary Problems for the Harry Dym Equation & Reciprocal Associates Colin Rogers Australian Research Council Centre of Excellence for Mathematics & Statistics of Complex Systems & The University
More informationAutomata on linear orderings
Automata on linear orderings Véronique Bruyère Institut d Informatique Université de Mons-Hainaut Olivier Carton LIAFA Université Paris 7 September 25, 2006 Abstract We consider words indexed by linear
More informationWRONSKIANS AND LINEAR INDEPENDENCE... f (n 1)
WRONSKIANS AND LINEAR INDEPENDENCE ALIN BOSTAN AND PHILIPPE DUMAS Abstract We give a new and simple proof of the fact that a finite family of analytic functions has a zero Wronskian only if it is linearly
More informationRational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation
Commun. Theor. Phys. (Beijing, China) 43 (005) pp. 975 98 c International Academic Publishers Vol. 43, No. 6, June 15, 005 Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional
More informationExact Solutions of Nonlinear Partial Dif ferential Equations by the Method of Group Foliation Reduction
Symmetry, Integrability and Geometry: Methods and Applications Exact Solutions of Nonlinear Partial Dif ferential Equations by the Method of Group Foliation Reduction SIGMA 7 (2011), 066, 10 pages Stephen
More informationSYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman
. SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado
More informationGauge Fixing and Constrained Dynamics in Numerical Relativity
Gauge Fixing and Constrained Dynamics in Numerical Relativity Jon Allen The Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation is reviewed. Gauge freedom is discussed and
More informationSymmetries and solutions of field equations of axion electrodynamics
Symmetries and solutions of field equations of axion electrodynamics Oksana Kuriksha Petro Mohyla Black Sea State University, 10, 68 Desantnukiv Street, 54003 Mukolaiv, UKRAINE Abstract The group classification
More informationLINEAR EQUATIONS WITH UNKNOWNS FROM A MULTIPLICATIVE GROUP IN A FUNCTION FIELD. To Professor Wolfgang Schmidt on his 75th birthday
LINEAR EQUATIONS WITH UNKNOWNS FROM A MULTIPLICATIVE GROUP IN A FUNCTION FIELD JAN-HENDRIK EVERTSE AND UMBERTO ZANNIER To Professor Wolfgang Schmidt on his 75th birthday 1. Introduction Let K be a field
More informationRepresentations of Sp(6,R) and SU(3) carried by homogeneous polynomials
Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials Govindan Rangarajan a) Department of Mathematics and Centre for Theoretical Studies, Indian Institute of Science, Bangalore 560 012,
More informationConditional symmetries of the equations of mathematical physics
W.I. Fushchych, Scientific Works 2003, Vol. 5, 9 16. Conditional symmetries of the equations of mathematical physics W.I. FUSHCHYCH We briefly present the results of research in conditional symmetries
More information(Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1. Lisa Carbone Rutgers University
(Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1 Lisa Carbone Rutgers University Slides will be posted at: http://sites.math.rutgers.edu/ carbonel/ Video will be
More informationClassification of cubics up to affine transformations
Classification of cubics up to affine transformations Mehdi Nadjafikah and Ahmad-Reza Forough Abstract. Classification of cubics (that is, third order planar curves in the R 2 ) up to certain transformations
More information-state problems and an application to the free particle
-state problems and an application to the free particle Sourendu Gupta TIFR, Mumbai, India Quantum Mechanics 1 2013 3 September, 2013 Outline 1 Outline 2 The Hilbert space 3 A free particle 4 Keywords
More informationJOINING OF CIRCUITS IN P SL(2, Z)-SPACE
Bull. Korean Math. Soc. 5 (015), No. 6, pp. 047 069 http://dx.doi.org/10.4134/bkms.015.5.6.047 JOINING OF CIRCUITS IN P SL(, Z)-SPACE Qaiser Mushtaq and Abdul Razaq Abstract. The coset diagrams are composed
More informationToward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach
Commun. Theor. Phys. 57 (2012) 5 9 Vol. 57, No. 1, January 15, 2012 Toward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach G. Darmani, 1, S. Setayeshi,
More informationExact solutions through symmetry reductions for a new integrable equation
Exact solutions through symmetry reductions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX, 1151 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es
More informationarxiv:hep-th/ v1 16 Jul 1992
IC-92-145 hep-th/9207058 Remarks on the Additional Symmetries and W-constraints in the Generalized KdV Hierarchy arxiv:hep-th/9207058v1 16 Jul 1992 Sudhakar Panda and Shibaji Roy International Centre for
More informationEXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL
Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 05, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE OF NONTRIVIAL
More informationPeriodic and Solitary Wave Solutions of the Davey-Stewartson Equation
Applied Mathematics & Information Sciences 4(2) (2010), 253 260 An International Journal c 2010 Dixie W Publishing Corporation, U. S. A. Periodic and Solitary Wave Solutions of the Davey-Stewartson Equation
More informationBLOCKS IN THE CATEGORY OF FINITE-DIMENSIONAL REPRESENTATIONS OF gl(m n)
BLOCKS IN THE CATEGORY OF FINITE-DIMENSIONAL REPRESENTATIONS OF gl(m n) VERA SERGANOVA Abstract. We decompose the category of finite-dimensional gl (m n)- modules into the direct sum of blocks, show that
More informationGravitational Waves: Just Plane Symmetry
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 2006 Gravitational Waves: Just Plane Symmetry Charles G. Torre Utah State University Follow this and additional works at:
More informationNew approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations
Volume 28, N. 1, pp. 1 14, 2009 Copyright 2009 SBMAC ISSN 0101-8205 www.scielo.br/cam New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations HASSAN A. ZEDAN Mathematics
More informationSymmetry Reductions of (2+1) dimensional Equal Width. Wave Equation
Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,
More informationarxiv:math/ v1 [math.co] 6 Dec 2005
arxiv:math/05111v1 [math.co] Dec 005 Unimodality and convexity of f-vectors of polytopes Axel Werner December, 005 Abstract We consider unimodality and related properties of f-vectors of polytopes in various
More informationA multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system
Chaos, Solitons and Fractals 30 (006) 197 03 www.elsevier.com/locate/chaos A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Qi Wang a,c, *,
More informationA determinant characterization of moment sequences with finitely many mass-points
To appear in Linear and Multilinear Algebra Vol 00, No 00, Month 20XX, 1 9 A determinant characterization of moment sequences with finitely many mass-points Christian Berg a and Ryszard Szwarc b a Department
More informationarxiv: v1 [math-ph] 4 May 2016
arxiv:1605.01173v1 [math-ph] 4 May 2016 On the Classifications of Scalar Evolution Equations with Non-constant Separant Ayşe Hümeyra BİLGE Faculty of Engineering and Natural Sciences, Kadir Has University,
More informationNew explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation
Physics Letters A 07 (00) 107 11 www.elsevier.com/locate/pla New explicit solitary wave solutions for ( + 1)-dimensional Boussinesq equation and ( + 1)-dimensional KP equation Yong Chen, Zhenya Yan, Honging
More informationConservation laws for difference equations. Peter Hydon (Surrey) Elizabeth Mansfield (Kent) Olexandr Rasin (Bar-Ilan) Tim Grant (Surrey)
Conservation laws for difference equations Peter Hydon (Surrey) Elizabeth Mansfield (Kent) Olexandr Rasin (Bar-Ilan) Tim Grant (Surrey) Introduction: What is a conservation law? What is a conservation
More informationA supersymmetric Sawada-Kotera equation
A supersymmetric Sawada-Kotera equation arxiv:0802.4011v2 [nlin.si] 7 Dec 2008 Kai Tian and Q. P. Liu Department of Mathematics, China University of Mining and Technology, Beijing 100083, P.R. China Abstract
More informationTHE SEMISIMPLE SUBALGEBRAS OF EXCEPTIONAL LIE ALGEBRAS
Trudy Moskov. Matem. Obw. Trans. Moscow Math. Soc. Tom 67 (2006) 2006, Pages 225 259 S 0077-1554(06)00156-7 Article electronically published on December 27, 2006 THE SEMISIMPLE SUBALGEBRAS OF EXCEPTIONAL
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationExact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation
Exact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation Francesco Demontis (based on a joint work with C. van der Mee) Università degli Studi di Cagliari Dipartimento di Matematica e Informatica
More informationNILPOTENCY INDEX OF NIL-ALGEBRA OF NIL-INDEX 3
J. Appl. Math. & Computing Vol. 20(2006), No. 1-2, pp. 569-573 NILPOTENCY INDEX OF NIL-ALGEBRA OF NIL-INDEX 3 WOO LEE Abstract. Nagata and Higman proved that any nil-algebra of finite nilindex is nilpotent
More informationOn the Linearization of Second-Order Dif ferential and Dif ference Equations
Symmetry, Integrability and Geometry: Methods and Applications Vol. (006), Paper 065, 15 pages On the Linearization of Second-Order Dif ferential and Dif ference Equations Vladimir DORODNITSYN Keldysh
More informationDegenerate Perturbation Theory. 1 General framework and strategy
Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying
More information