Characteristic Lie Algebras of Integrable Differential-Difference Equations in 3D
aa r X i v : . [ n li n . S I] F e b Characteristic Lie Algebrasof Integrable Differential-Difference Equations in 3D
I. Habibullin and A. Khakimova
Abstract
The purpose of this article is to develop an algebraic approach to the problem of integrable classification ofdifferential-difference equations with one continuous and two discrete variables. As a classification criterion,we put forward the following hypothesis. Any integrable equation of the type under consideration admitsan infinite sequence of finite-field Darboux-integrable reductions. The property of Darboux integrability of afinite-field system is formalized as finite-dimensionality condition of its characteristic Lie-Rinehart algebras.That allows one to derive effective integrability conditions in the form of differential equations on the righthand side of the equation under study. To test the hypothesis, we use known integrable equations from theclass under consideration. In this article, we show that all known examples do have this property.
The theory of integrability is an important component of modern mathematical physics. At present, theproblem of classifying integrable equations with three and more independent variables is being activelystudied. For these classes of equations, in contrast to the case of two variables, there are several approachesto the problem of integrable classification (see [1–5, 7–9, 11–17, 19, 20, 31, 32, 35–38, 44, 45]) mostly based ongeometric ideas and constructions. Note that the well-known generalized symmetry method is not efficientin higher dimensions due to the non-locality problem.The hydrodynamic reduction method is one of the most widespread among specialists in mathematicalphysics. Let us briefly explain the essence of this method for the classification of three-dimensional integrableequations, based on the following two-stage procedure. The first consists in classification of dispersionless in-tegrable systems by using the method of hydrodynamic reductions (initiated by Gibbons, Kodama and Tsarev(see [19, 20]) and developed into an efficient integrability test by the group of Ferapontov, Khusnutdinova,Novikov, Pavlov, Odesskii and Sokolov ( [12, 35–37]). Second is reconstruction of dispersive deformations ofdispersionless integrable systems based on the method of dispersive deformations of hydrodynamic reductionsdeveloped by Ferapontov, Moro and Novikov ( [14, 15]).The article discusses a problem of developing an algebraic integrability criterion based on the Darbouxintegrable reductions for the differential-difference equations of the form: u jn +1 ,x = F ( u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 ) , −∞ < n, j < ∞ , (1)where the sought function u depends on three variables, on a real x and two integers j and n . We assumethat ∂F∂u jn,x = 0 , ∂F∂u j +1 n = 0 , ∂F∂u j − n +1 = 0 . The lattice is defined on a quadrilateral graph. To preserve the parity of the forward and backward directionsin n we assume that equation (1) can be uniquely rewritten as u jn − ,x = G ( u jn,x , u j +1 n − , u jn , u jn − , u j − n ) . Integrable equations of the form (1) have been studied by many authors (see [17] and the referencestherein). Note that, in some particular cases, equations of the form (1) can be obtained as the Backlundtransformation for the two-dimensional Toda-type chains u n,xy = f ( u n +1 , u n , u n − , u n,x , u n,y ) , −∞ < n < ∞ . (2)Earlier the lattices (1) have been investigated by means of the over mentioned method of hydrodynamicreductions [17]. Within the framework of this approach, a certain class of lattices of the form (1) is studiedand integrable cases are selected. n the studies [23, 24, 26, 39] a classification algorithm was developed for integrable equations with threeindependent variables based on the idea of the Darboux integrable reductions and characteristic Lie algebras.Efficiency of the method was illustrated by applying to the equations from the class (2) (see [23]).It was observed in the article [10] that it is very convenient to study the classification problem for (2) bycombining different approaches. At the first step, equations are selected that satisfy the requirement thatthe dispersionless limit of the equation is integrable, that is its characteristic variety defines a conformalstructure, which is Einstein-Weyl, on every solution. As a result one obtaines a class of equations havingthe freedom in only one function of one variable and a few constant parameters. And then, to the obtainedequations from this class the test of Darboux integrable reductions is applied.In this paper, we show at the level of examples that integrable differential-difference lattices of the form(1) admit Darboux integrable reductions. More precisely, the purpose of the article is to approve the following Conjecture 1.1
A lattice of the form (1) is integrable if and only if there exists a pair of functions H (1) and H (2) of four variables such that for any choice of the integer N , a system of the hyperbolic type differential-difference equations u n +1 ,x = H (1) ( u n,x , u n , u n +1 , u n ) ,u jn +1 ,x = F j ( u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 ) , < j < N, (3) u Nn +1 ,x = H (2) ( u Nn,x , u Nn +1 , u Nn , u N − n +1 ) , −∞ < n < ∞ obtained from (1) is integrable in the sense of Darboux. Let us comment on the essence of the Conjecture 1.1. Finite-field reductions of integrable lattices thatinherit integrability are quite common, for example, they are used to construct particular solutions of chains.Such reductions are provided by imposing boundary conditions (cutoff conditions) compatible with the in-tegrability property. The most popular of this kind conditions is the periodicity closure constraint. Usually,the finite-field reductions obtained as a result of integrable cutoffs are soliton systems, and only in the excep-tional case, which is obtained by imposing degenerate cutoff conditions, do we arrive at systems integrable inthe sense of Darboux. Our hypothesis is that, on the one hand, any integrable chain of the form (1) admitsdegenerate boundary conditions, which reduce the chain to a Darboux integrable system, and on the otherhand, a chain admitting degenerate boundary conditions with this property, certainly, is integrable. Exis-tence of Darboux integrable reductions allows one to apply the technique of the characteristic Lie-Rinehartalgebras for the classifying the lattices (1).The systematic study of partial differential equations of hyperbolic type of the form u x,y = f ( x, y, u, u x , u y ) , the integration problem of which is reduced to solving ordinary differential equations began apparently inthe second half of the nineteenth century. The problem of a complete description of this class of equations,called Darboux integrable equations or equations of Liouville type was studied in the well-known works ofE. Goursat. The task of complete classification turned out to be very difficult, and later many researchersworked in this field. The reader can find an overview of the results and a description of the current state ofthe problem in a detailed article [47].Systems of nonlinear hyperbolic equations remain less studied from the point of view of Darboux inte-grability. The exception is a class of the so-called exponential type systems [18, 21, 33, 34, 41] u ix,y = exp X j a i,j u j . It is known that this system is Darboux integrable if and only if A = { a i,j } is the Cartan matrix of a semi-simple Lie algebra. It is remarkable that general solution for Darboux integrable systems can be constructedin a closed form (see [6, 18, 33, 43, 47]).Characteristic algebras for the Darboux integrable scalar differential-difference equations of the form u n +1 ,x = u n,x + f ( u n , u n +1 ) are introduced and investigated in [25, 27, 28]. In [28] the complete list of such equations is given. Darbouxintegrable systems of differential-difference equations of the exponential type are studied in [30] and [42].We briefly discuss the content of the article. In §2 we consider finite systems of the differential-differenceequations of the form (4). We recall the definition of the x - and n -integrals and define the integrability in thesense of Darboux. We introduce the characteristic vector fields and the characteristic Lie-Rinehart algebras n the directions of n and x . In §3, we study the structure of characteristic algebras and describe the action oftheir automorphisms, which will be needed when solving the classification problem. We then in §4 illustratehow to use characteristic algebra to study the problem of integrable classification of the systems of the form(3). As an example, we derive one of the necessary conditions of integrability for the system (3). First,we derive a system of nonlinear functional equations from the requirement that the characteristic algebra isfinite-dimensional, and then we reduce these functional equations to a system of differential equations. InSection 5, we study all seven known differential-difference integrable lattices with one continuous and twodiscrete independent variables. In each of these equations, we make a change of variables in order to obtainan equation that admits a degenerate cutoff. By imposing the cutoff conditions at two points, we reducethe equation to a finite-field system of differential-difference equations. Further, we show that the resultingfinite-field systems are Darboux integrable by presenting for these systems complete sets of integrals in bothcharacteristic directions (among them there are scalar equations, they are presented separately in section5.8). This circumstance confirms the correctness of the approach to the classification problem formulated inthe Conjecture 1.1 above. We use two different methods to construct integrals. In some cases, when the Laxpair satisfies the conditions of Lemma 5.1 or Lemma 5.2, it is convenient to use the Lax pair. In other cases,we use the method of characteristic Lie-Reinhart algebras. This method is always applicable for constructingintegrals, although it is more labor-consuming. In this section we concentrate on a system of differential-difference equations of general form u jn +1 ,x = F j ( x, n, u jn,x , u n , u n +1 ) , j = 1 , , . . . , N, (4)where u n = ( u n , u n , . . . , u Nn ) . We request that functions F j are analytic in a domain in C + , where C isthe complex plane.We assume that system (4) can be rewritten in the converse way u jn − ,x = G j ( x, n, u jn,x , u n , u n − ) , j = 1 , , . . . , N. Variables u jn , u jn ± , u jn ± , . . . and u jn,x , u jn,xx , u jn,xxx , . . . are called dynamical variables, we treat them as in-dependent ones.In what follows we will be interested in the cases when system (4) admits nontrivial x - and n -integrals.Let us give the necessary definitions. We call a function I = I ( x, n, u, u n,x , u n,xx , . . . ) n -integral of (4) if itsatisfies the condition D n I = I by means of the system, where D n is the shift operator acting due to therule D n y ( n ) = y ( n + 1) . Similarly function J = J ( x, n, u n , u n ± , u n ± , . . . ) is an x -integral if the equation D x J = 0 holds, where D x stands for the operator of the total derivative with respect to x . Integrals I = I ( x ) and J = J ( n ) depending only on x and correspondingly only on n are called trivial integrals. System(4) is called integrable in sense of Darboux if it possesses the complete set of nontrivial integrals in bothcharacteristic directions x and n . Completeness means that the number of functionally independent integralsin each direction coincides with the order N of the system.In order to find the integrals one can use the so-called characteristic operators. Let us discuss howthese operators are derived. Assume that J is an x -integral, then according to the definition we have D x J ( x, n, u n , u n ± , u n ± , . . . ) = 0 , we use the chain rule and get K J = ∂∂x + p X j =1 u jn,x ∂∂u jn + u jn +1 ,x ∂∂u jn +1 + u jn − ,x ∂∂u jn − + . . . J = 0 . (5)We specify the operator K by means of the system (4) and obtain K = ∂∂x + p X j =1 u jn,x ∂∂u jn + F jn ∂∂u jn +1 + G jn ∂∂u jn − + F jn +1 ∂∂u jn +2 + G jn − ∂∂u jn − + . . . , (6)where F jn + i = F j ( x, n + i, u jn + i,x , u n + i , u n + i +1 ) and G jn − i = G j ( x, n − i, u jn − i,x , u n − i , u n − i − ) . We call K the characteristic operator in the x direction.An important feature of the equation (5) is that the solution J does not depend on the variables u n,x , u n,x , . . . , u Nn,x despite the fact that the coefficients of the equation depend on them. Therefore inaddition to (5) any x -integral solves also a set of the following equations X k J = 0 , k = 1 , , . . . , N, (7) here X k = ∂∂u kn,x . Evidently x -integral is nullified by the multiple commutators of the operators X , X , . . . , X N , K (8)as well as linear combinations of the obtained operators with variable coefficients, depending on the dynamicalvariables. In other words, any x -integral belongs to the kernel of the operators in the Lie-Rinehart algebra L x , generated by the operators (8) over the ring A of locally analytic functions depending on the dynamicalvariables. We call it characteristic algebra in the direction of x . Actually algebra L x is a finitely generatedmodule over the ring A . Here the elements in L x can be multiplied by functions from A . The consistencyconditions:1) [ W , aW ] = W ( a ) W + a [ W , W ] ,2) ( aW ) b = aW ( b ) are assumed to be valid for any W , W ∈ L x and a, b ∈ A . In other words we request that, if W ∈ L x and a ∈ A then aW ∈ L x (see [40]).The algebra L x is of a finite dimension if it admits a basis consisting of a finite number of the elements W , W , . . . , W k ∈ L x such that an arbitrary operator W ∈ L x is represented as a linear combination of theform W = a W + a W + · · · + a k W k , (9)where the coefficients are functions a , a , . . . , a k ∈ A . If in (9) W = 0 then we have a = 0 , a = 0 , . . . , a k = 0 .Let us now discuss the n -integrals. Assume that a function I = I ( x, n, u n , u n,x , u n,xx , . . . ) is an n -integralfor the system (4), i.e. the following relation holds D n I = I , or the same I ( x, n + 1 , u n +1 , u n +1 ,x , u n +1 ,xx , ... ) = I ( x, n, u n , u n,x , u n,xx , . . . ) . Due to the equation (4), that can be written in a shortened way as u n +1 ,x = F ( u n,x , u n , u n +1 ) where F = ( F , F , . . . , F N ) we rewrite the latter as I ( x, n + 1 , u n +1 , F, D x F, D x F, . . . ) = I ( x, n, u n , u n,x , u n,xx , . . . ) . (10)The problem of finding n -integrals is much more difficult than finding x -integrals, indeed, in this case, oneneeds to solve not a differential but a functional equation (10). Below we introduce characteristic operatorsthat allow to study completely this equation. It is easily observed that the right hand side of (10) does notcontain dependence on the variable u n +1 = ( u n +1 , u n +1 , . . . , u Nn +1 ) therefore the left hand side satisfies therelations ∂∂u jn +1 D n I = 0 , which evidently implies Y j, I := D − n ∂∂u jn +1 D n I = 0 for j = 1 , , . . . , N. (11)It can be proved by a direct computation that characteristic operators defined in (11) act on the dynamicalvariables u n , u n,x , u n,xx , . . . as vector fields of the form Y j, = ∂∂u jn + N X i =1 D − n ∂F in ∂u jn +1 ! ∂∂u in,x + D − n ∂F in,x ∂u jn +1 ! ∂∂u in,xx + . . . , (12)where F in,x := D x F in . We set Y j, = ∂∂u jn +1 and then rewrite representation as follows Y j, = ∂∂u jn + N X i =1 D − n (cid:0) Y j, (cid:0) F in (cid:1)(cid:1) ∂∂u in,x + D − n (cid:0) Y j, (cid:0) F in,x (cid:1)(cid:1) ∂∂u in,xx + . . . . In contrast to the case of the previously studied x -integrals, in this case there are additional characteristicoperators, which also annul I . For any natural number k , the following relation holds: Y j,k I := D − kn ∂∂u jn +1 D kn I = 0 for j = 1 , , . . . , N. (13)We give uniform representations for all characteristic operators (13) for k ≥ Y j,k = N X i =1 D − n (cid:0) Y j,k − (cid:0) F in (cid:1)(cid:1) ∂∂u in,x + D − n (cid:0) Y j,k − (cid:0) F in,x (cid:1)(cid:1) ∂∂u in,xx + . . . . (14) ince the coefficients of the equations (13) depend on the variables u jn − s for ≤ s ≤ k and ≤ j ≤ N , whilethe solution to the equation does not depend on them, then such equations ¯ X i,s I = 0 , where ¯ X i,s = ∂∂u in − s should hold for these values of i and s .It is important that the operators Y j,i commute with each other. Lemma 2.1
For any i, i ′ ≥ and j, j ′ such that ≤ j, j ′ ≤ N the relation [ Y j,i , Y j ′ ,i ′ ] = 0 holds. Proof 2.1
Let us first assume that i ′ = 0 , then obviously [ Y j,i , Y j ′ , ] = " Y j,i , ∂∂u j ′ n +1 = 0 since the coefficients of the vector-field Y j,i do not depend on the variables u jn +1 for any j ∈ [1 , N ] due toexplicit formula (14). Now consider general case assuming that i = i ′ + k , k ≥ : [ Y j,i ′ + k , Y j ′ ,i ′ ] = h D − i ′ n Y j,k D i ′ n , D − i ′ n Y j ′ , D i ′ n i = D − i ′ n [ Y j,k , Y j ′ , ] D i ′ n . Latter vanishes due to the previous step.
Explicit expressions for the commutators of the basic operators Y j,k and X j,k with the operator D x ofthe total derivative with respect to x provide an effective tool for studying the system (4). It can easily beapproved that the following operator identities hold [ D x , Y j, ] = − N X i =1 ∞ X k =1 Y j, (cid:0) D k − n F in (cid:1) ∂∂u in + k , (15) [ D x , Y j, ] = − N X i =1 D − n (cid:0) Y j, (cid:0) F in (cid:1)(cid:1) Y i, + ∞ X k =1 Y j, (cid:0) D k − n F in (cid:1) ∂∂u in + k + ∞ X k =1 Y j, (cid:0) D − kn G in (cid:1) ∂∂u in − k ! , (16) [ D x , Y j,m ] = − N X i =1 m X s =1 D − sn (cid:0) Y j,m − s (cid:0) F in (cid:1)(cid:1) Y i,s + ∞ X k =1 Y j,m (cid:0) D k − n F in (cid:1) ∂∂u in + k + ∞ X k =1 Y j,m (cid:0) D − kn G in (cid:1) ∂∂u in − k ! , (17) (cid:2) D x , ¯ X j, (cid:3) = − N X i =1 ∞ X k =1 ¯ X j, (cid:0) D − kn G in (cid:1) ∂∂u in − k . (18)In a particular scalar case of the system (4) when N = 1 formulas (15)–(18) have been derived in [28].Let A k be a ring of locally analytic functions of the variables u n − k , u n − k +1 , . . . , u n ; u nx , u n,xx , u n,xxx , . . . .Here we formulate two theorems relating integrability in the sense of Darboux and properties of the charac-teristic Lie-Rinehart algebras. Theorem 2.1
System (4) admits a complete set of n -integrals if and only if the following two conditionshold:1) The linear space V spanned by the operators { Y i,s } has finite dimension, which we denote by N .Assume that Z , Z , . . . , Z N constitute a basis in V , such that for any Z ∈ V we have an expansion Z = λ Z + λ Z + . . . + λ N Z N , We emphasize that the coefficients in this expansion, as a rule, are not constant, they are analytic functions ofdynamic variables. Then due to construction of the space V there exists a number N , such that [ Z j , ¯ X i,s ] = 0 for all j = 1 , , . . . , N , i = 1 , , . . . , N and s > N .2) The Lie-Rinehart algebra L n generated by the operators { Z j } N j =1 and { ¯ X i,s } N ,Ns =1 ,i =1 over the ring A N has a finite dimension. Theorem 2.2
System (4) admits a complete set of x -integrals if and only if the characteristic Lie-Rinehartalgebra L x has a finite dimension. The proof of these two theorems is beyond the scope of this paper. In the particular case when N = 1 Theorems 2.1, 2.2 are proved in [29]. For systems of differential equations of hyperbolic type, similar statementis proved in [46]. From these theorems we derive
Corollary 2.1
System (4) is integrable in the sense of Darboux if and only if both characteristic algebras L x and L n have finite dimension. Some properties of the characteristic algebras and their potentialapplications
In this section we will discuss some important properties of the characteristic algebras. Here we outline anapproach to the problem of integrable classification of the lattices from the class (1) based on algebraic ideas.Let us first concentrate on the algebra L x for the system (4), generated as it was discussed in §2 above bythe operators X , X , . . . , X N , K . It is easy to see that the operator D x of the total derivative with respectto x acts on the functions of the dynamical variables u n,x , u n , u n ± , u n ± , . . . as follows D x H ( u n,x , u n , u n ± , u n ± , . . . ) = N X j =1 u jn,xx X j + K H ( u n,x , u n , u n ± , u n ± , . . . ) , where the operators K and X j have already been defined above (see (6), (7)). Since the operators D n and D x commute with each other we have a relation D n D x D − n = D x or, the same D n N X j =1 u jn,xx X j + K D − n = N X j =1 u jn,xx X j + K . (19)We simplify the left hand side of (19) due to the relations D n u jn,xx = u jn +1 ,xx D n and due to u jn +1 ,xx = D x F jn = ∂F jn ∂u jn,x · u jn,xx + K ( F jn ) . As a result we obtain the relation N X j =1 ∂F jn ∂u jn,x · u jn,xx + K (cid:0) F jn (cid:1)! D n X j D − n + D n K D − n = N X j =1 u jn,xx X j + K . (20)Let us define an automorphism of the algebra L x , acting according to the formula Z → D n ZD − n . (21)It is remarkable that equation (20) allows one to describe the action of the automorphism on the basicoperators. Indeed, since variables u n,xx , u n,xx , . . . , u Nn,xx are regarded as independent ones, we can comparethe coefficients in front of these variables in (20) and obtain the formulas D n X j D − n = 1 ∂F jn ∂u jn,x X j , (22) D n K D − n = K − N X j =1 K ( F jn ) ∂F jn ∂u jn,x X j . (23)Elements of the algebra L x are vector fields with infinite number of components. Therefore, to proverelations of the form Y = 0 for Y ∈ L x , it is necessary to check an infinite set of conditions. The followinglemma provides a convenient tool for exploring such questions. Lemma 3.1
Assume that the vector field K = N X j =1 ∞ X k =1 α j ( k ) ∂∂u jn + k + α j ( − k ) ∂∂u jn − k ! , solves the equation D n KD − n = hK, (24) where the factor h is a function of the dynamical variables, then K = 0 . Proof 3.1
First, we find an explicit representation of the operator D n KD − n and then by substituting theresult into equation (24) we obtain D n ( α j ( − N X s =1 D n ∂∂u jn − G sn ! X s + D n ( α j ( − ∂∂u jn + D n ( α j ( − ∂∂u jn − ∞ X k =2 D n ( α j ( k − ∂∂u jn + k + D n ( α j ( − k − ∂∂u jn − k (25) = h ∞ X k =1 α j ( k ) ∂∂u jn + k + α j ( − k ) ∂∂u jn − k ! for any j from the segment ≤ j ≤ N .By comparing the coefficients before the operators ∂∂u jn + k and ∂∂u jn − k in the equation (25) one can easilyprove that α j ( k ) = 0 for all values of k . For the case of the algebra L n we have a similar statement. Lemma 3.2
Suppose that the vector field Y = N X j =1 α j (1) ∂∂u jn,x + α j (2) ∂∂u jn,xx + α j (3) ∂∂u jn,xxx + · · · ! (26) solves the equation [ D x , Y ] = hY, (27) where h is a function of the dynamical variables, then Y = 0 . Proof 3.2
From (26) and (27) we easily obtain the equation N X j =1 − α j (1) ∂F jn ∂u jn,x ∂∂u jn +1 − α j (1) ∂∂u jn + ( D x α j (1) − α j (2)) ∂∂u jn,x + ( D x α j (2) − α j (3)) ∂∂u jn,xx + · · · ! = h N X j =1 α j (1) ∂∂u jn,x + α j (2) ∂∂u jn,xx + α j (3) ∂∂u jn,xxx + · · · ! . Now we compare the coefficients before the operators ∂∂u jn +1 , ∂∂u jn , ∂∂u jn,x , ∂∂u jn,xx , . . . to make sure that α j ( k ) = 0 for all k . That completes the proof of Lemma 3.2. In this section we concentrate on the differential-difference system (3) being a reduction of the differential-difference equation (1). System (3) is a particular case of the general system (4), where assumed that F = H (1) ( u n,x , u n , u n +1 , u n ) ,F j = F ( u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 ) for ≤ j ≤ N − ,F N = H (2) ( u Nn,x , u Nn +1 , u Nn , u N − n +1 ) . Suppose that system (3) is integrable, then its characteristic algebra L x has a finite dimension. The problemis figuring out what this means in terms of the right hand side of the differential-difference equation (1),i.e. in terms of the function F ( u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 ) . Below we make the first step towards solving theproblem. We define a sequence of the operators in L x according to the rule K = [ X j , K ] , K = [ X j , K ] , . . . , K m = [ X j , K m − ] , . . . . (28)Since the algebra L x is finite dimensional then there exists an integer M such that operator K M is representedas a linear combination of the preceding members of the sequence: K M = λK M − + λ K M − + · · · + λ M − K , (29)where the operators K M − , K M − , . . . , K are linearly independent.Explicit formula (6) shows that operator K = [ X j , K ] has the following coordinate representation K = ∂∂u jn + X j ( F jn ) ∂∂u jn +1 + X j ( G jn ) ∂∂u jn − + X j ( F jn +1 ) ∂∂u jn +2 + X j ( G jn − ) ∂∂u jn − + . . . , hile for m ≥ we have K m = X mj ( F jn ) ∂∂u jn +1 + X mj ( G jn ) ∂∂u jn − + X mj ( F jn +1 ) ∂∂u jn +2 + X mj ( G jn − ) ∂∂u jn − + . . . . (30)Since operator K contains the term ∂∂u jn and the other operators in (29) do not, then the coefficient λ M − vanishes.Automorphism (21) provides an effective tool for studying the problem of decomposition (29). Above in(22), (23) we have already described the action of the automorphism on the basic operators. We can derivesimilar formulas for the members of the sequence (28) D n K D − n = q j K − q j K ( F jn ) X j − q j − F j − n,u jn q j X j − − q j +1 F j +1 n,u jn X j +1 , (31) D n K D − n = q j K + q j X j ( q j ) K + r X j − q j − F j − u jn q j X j ( q j ) X j − , (32) D n K D − n = q j K + 3 q j X ( q j ) K + ( q j X j ) ( q j ) K + r X j − q j − F j − n,u jn (cid:8) ( q j X j ) ( q j ) (cid:9) X j − . (33)Here q j = (cid:16) ∂F jn ∂u jn,x (cid:17) − , and r , r are functions on q j and its derivatives. For the general case we can proveby induction that D n K m D − n = q mj K m + p m K m − + s m K m − + . . . + ( q j X j ) m − ( q j ) K + r m X j + q j − F j − n,u jn (cid:8) ( q j X j ) m − ( q j ) (cid:9) X j − , (34)where p m = 12 m ( m − q m − j X j ( q j ) ,s m = 16 m ( m − m − q m − j X j ( q j ) + 124 m ( m − m − m − q m − ( X j ( q j )) . Explicit expression for r m is rather complicated and we do not specify it.Let us apply automorphism (21) to both sides of the equation (30) and get q Mj K M + p M K M − + s M K M − + . . . + ( q j X j ) M − ( q j ) K + r M X j + q j − F j − n,u jn (cid:8) ( q j X j ) M − ( q j ) (cid:9) X j − = D n ( λ ) (cid:0) q M − j K M − + . . . (cid:1) + D n ( λ ) (cid:0) q M − j + . . . (cid:1) + . . . . (35)We first replace K M due to (30) and then collect in (35) the coefficients in front of the linearly independentoperators K M − , K M − , . . . , K : K M − : q Mj λ + p M = D n ( λ ) q M − j , (36) K M − : q Mj λ + s M = D n ( λ ) p M − + D n ( λ ) q M − j , (37) . . . . . . . . . . . . . . .K : ( q j X j ) M − q j = D n ( λ )( q j X j ) M − q j + . . . + D n ( λ ) q j X j ( q j ) . (38)Due to the explicit expression (30) the coefficients of the operators K , K , . . . , K M depend on u jn,x andon the variables u j +1 n , u jn , u j − n and their shifts on n . The factors λ, λ , . . . , λ M − might depend only onthese variables. However we have relations (36)–(38) which provide an additional restriction for the factors.Indeed the coefficients q j , p m , s m in (36)–(38) depend only on the variables u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 , hence λ, λ , . . . , λ M − might depend only on the variables u jn,x , u jn .Relations (36)–(38) determine a system of functional equations with the set of unknown functions λ , λ ,. . . , λ M − depending on two variables u jn,x , u jn . Indeed, each of these equations contains the unknowns takenat different points ( u jn,x , u jn ) and (cid:0) D n u jn,x , D n u jn (cid:1) = (cid:16) F (cid:16) u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 (cid:17) , u jn +1 (cid:17) . An importantpeculiarity of the system (36)–(38) is that it is highly overdetermined since the coefficients of the equationsdepend on five independent variables while the solutions depend only on two of them.Let us finalize the reasoning above as a statement. Theorem 4.1 (Necessary condition of integrability of system (3) and hence due to Conjec-ture 1.1 of the equation (1).)
If system (3) admits the complete set of x -integrals then the overdeterminedsystem of equations (36)–(38) admits a solution ( λ, λ , . . . , λ M − ) depending only on u jn,x , u jn for any j fromthe segment ≤ j ≤ N . elow we show how to reduce the system (36)–(38) to a system of differential equations by using thecharacteristic operators in the direction of n . We start with the equation (36) that is specified to the form λ + ε q j X j ( q j ) = D n ( λ ) 1 q j , where ε = M ( M − . Since q j = (cid:16) ∂F jn ∂u jn,x (cid:17) − equation reduces to λ ∂F jn ∂u jn,x − ε ∂ F jn ∂ ( u jn,x ) = D n ( λ ) ∂F jn ∂u jn,x ! . (39)We can rewrite the latter equation as λ − ε ∂∂u jn,x log ∂F jn ∂u jn,x = D n ( λ ) ∂F jn ∂u jn,x . (40)Let us apply the operator ∂∂u jn +1 to both sides of the equation (40). By taking into account that λ = λ (cid:0) u jn,x , u jn (cid:1) we obtain − ε ∂∂u jn +1 ∂∂u jn,x log ∂F jn ∂u jn,x ! = ∂∂u jn +1 D n ( λ ) ∂F jn ∂u jn,x ! . (41)As it was discussed above we have the relation D − n ∂∂u jn +1 D n = Y j, (see formula (12)) when the operatorsact on the variables u jn , u jn − , u jn,x . Since ∂∂u jn,x = X j we can represent (41) in the form − εY j, (cid:8) D − n X j log X j ( F jn ) (cid:9) = Y j, (cid:8) λD − n ( X j ( F jn )) (cid:9) . (42)In what follows we need some rather simple formulas which are produced by the obvious identity u jn,x = G j (cid:16) F j (cid:16) u jn,x , u j +1 n , u jn +1 , u jn , u j − n +1 (cid:17) , u j +1 n , u jn +1 , u jn , u j − n +1 (cid:17) . Differentiation of the identity with respect to the variable u jn,x gives rise to equation ∂D n G jn ∂F jn × ∂F jn ∂u jn,x which implies D n ∂G jn ∂u jn,x ! ∂F jn ∂u jn,x ! or the same D − n X j F jn = 1 X j G jn . (43)By combining (22) and (43) we find D − n X j D n = 1 X j G jn X j . (44)By differentiating identity with respect to u jn +1 we obtain an equation D n (cid:18) ∂G jn ∂u jn (cid:19) + D n ∂G jn ∂u jn,x ! ∂F jn ∂u jn +1 , which implies D − n ∂F jn ∂u jn +1 = − ∂G jn /∂u jn ∂G jn /∂u jn,x . Now equation (42) is easily converted due to (43), (44) to the form εY j, X j G jn ( X j G jn ) ! = Y j, λX j G jn ! . (45) et as replace the characteristic operator Y j, by its explicit representation (12). Since in (45) the operatoris applied to functions depending on the variables u jn,x , u j +1 n − , u jn , u jn − , u j − n only we omit the terms whichwe do not use. More precisely we use the relation that holds for the function h depending on these variables Y j, h = ∂∂u jn − D − n ∂F jn ∂u jn +1 ! ∂∂u jn,x ! h. Therefore equation (45) gets the form ∂∂u jn − ∂G jn /∂u jn ∂G jn /∂u jn,x ∂∂u jn,x ! λ∂G jn /∂u jn,x − ε ∂ G jn /∂ ( u jn,x ) (cid:16) ∂G jn /∂u jn,x (cid:17) = 0 . (46)Thus the functional equation (36) produces a differential consequence (46). In a similar way we can derivedifferential consequences for the other functional equations of the system (36)–(38). For example (37) implies ∂∂u jn − ∂G jn /∂u jn ∂G jn /∂u jn,x ∂∂u jn,x ! δ X j ( G jn )( X j ( G jn )) + γ ( X j ( G jn )) ( X j ( G jn )) − ε λX j ( G jn )( X j ( G jn )) + λ ( X j ( G jn )) ! = 0 , where ε = ( M − M − , δ = M ( M − M − , γ = M ( M − M − M − .By using the identity ( q j X j ) m q j = X m +1 j ( G jn ) , that is easily proved by induction one can rewrite the lastequation of the system as follows (cid:0) X Mj − λX M − j − λ X M − j − . . . − λ M − X j (cid:1) G jn = 0 . Finally we can state that the overdetermined system of functional equations (36)–(38) is reduced to anoverdetermined system of differential equations, that is certainly more simple. An equation (46) is of a specialinterest, it should admit a solution λ = λ (cid:0) u jn , u jn,x (cid:1) depending only on two variables u jn and u jn,x . It is asevere requirement since the coefficients of the equation (46) depend generally on five variables u jn,x , u j +1 n − , u jn , u jn − , u j − n . Example 4.1
We consider an illustrative example by taking a scalar equation of the form u n +1 ,x = u n,x u n +1 . (47) Here F n = u n,x u n +1 , ∂F n ∂u n,x = u n +1 and ∂ F n ∂ ( u n,x ) = 0 . Equation (39) turns into D n ( λ ) λ = 1 u n +1 (48) Due to equation (47) we can rewrite (48) as D n ( λ ) λ = u n,x u n +1 ,x . Evidently it has a solution of the necessary form λ = cu n,x where c is a constant. The aim of this section is discussing the integrable lattices of the form (1) given in [17]. Since these equationsare written as v sk +1 ,x = ¯ F ( v s +1 k,x , v sk , v sk +1 , v s +1 k , v s +1 k +1 ) , −∞ < s, k < ∞ , at first we convert them into the form (1) by a linear transformation of the independent variables: v sk = u j +1 n , n = k , j = − k − s . Below we give the list of integrable differential-difference equations from [17] rewrittenin the new variables u jn +1 ,x = u jn,x + e u jn +1 − u j +1 n − e u j − n +1 − u jn ; (49) u jn +1 ,x = u jn,x + e u j − n +1 e u jn +1 − e u j − n +1 e u jn − e u jn e u j +1 n + e u jn +1 e u j +1 n ; (50) ) u jn +1 ,x = u jn,x (cid:16) u jn +1 (cid:17) u j − n +1 u j +1 n ; (51) u jn +1 ,x = u jn,x (cid:16) u jn +1 − u j +1 n (cid:17)(cid:16) u j − n +1 − u jn (cid:17) ; (52) u jn +1 ,x = u jn,x u jn +1 (cid:16) u jn +1 − u j +1 n (cid:17) u j +1 n (cid:16) u j − n +1 − u jn (cid:17) ; (53) u jn +1 ,x = u jn,x (cid:16) u j − n +1 − u jn +1 (cid:17) (cid:16) u jn +1 − u j +1 n (cid:17)(cid:16) u j − n +1 − u jn (cid:17) (cid:16) u jn − u j +1 n (cid:17) ; (54) u jn +1 ,x = u jn,x sinh (cid:16) u j − n +1 − u jn +1 (cid:17) sinh (cid:16) u jn +1 − u j +1 n (cid:17) sinh (cid:16) u j − n +1 − u jn (cid:17) sinh (cid:16) u jn − u j +1 n (cid:17) (55)their Lax pairs, also found in [17] are respectively, of the form ( ϕ jn +1 = − ϕ j +1 n + e u jn − u jn +1 ϕ jn ,ϕ jn,x = − e u j − n − u jn ϕ j − n ; (56) ( ϕ jn +1 = − e − u jn +1 + u jn (cid:0) ϕ j +1 n − ϕ jn (cid:1) + ϕ j +1 n ,ϕ jn,x = e u j − n − u jn (cid:0) ϕ j − n − ϕ jn (cid:1) ; (57) ϕ jn +1 = − u jn +1 u j +1 n (cid:0) ϕ jn − ϕ j +1 n (cid:1) ,ϕ jn,x = − u jn,x u jn (cid:0) ϕ j − n − ϕ jn (cid:1) ; (58) ( ϕ jn +1 = (cid:16) u j +1 n − u jn +1 (cid:17) ϕ jn + ϕ j +1 n ,ϕ jn,x = − u jn,x ϕ j − n ; (59) ϕ jn +1 = (cid:18) − u jn +1 u j +1 n (cid:19) ϕ jn − u jn +1 u j +1 n ϕ j +1 n ,ϕ jn,x = u jn,x u jn (cid:0) ϕ j − n + ϕ jn (cid:1) ; (60) ϕ jn +1 = u jn +1 − u j +1 n u jn − u j +1 n ϕ jn + (cid:18) − u jn +1 − u j +1 n u jn − u j +1 n (cid:19) ϕ j +1 n ,ϕ jn,x = u jn,x u j − n − u jn (cid:0) ϕ j − n + ϕ jn (cid:1) ; (61) ϕ jn +1 = e ( ujn +1 − uj +1 n ) − e ( ujn − uj +1 n ) − ϕ jn + (cid:18) − e ( ujn +1 − uj +1 n ) − e ( ujn − uj +1 n ) − (cid:19) ϕ j +1 n ,ϕ jn,x = u jn,x e ( uj − n − ujn ) − (cid:0) ϕ j − n − ϕ jn (cid:1) . (62)It can be shown that equations (49)-(55) admit infinite sequences of suitable finite field reductions beingintegrable in the sense of Darboux systems of differential-difference hyperbolic type equations. Simultane-ously, the Lax pairs (56)-(62) of equations also pass into the Lax pairs of the corresponding reductions. Thecrucial point of our algorithm is finding a special constraint consistent with the lattice that divides the latticeinto two independent parts. We call such kind of constraint a degenerate boundary condition. Below we lookfor the degenerate boundary conditions for the equations (49)-(55) and appropriate boundary conditions fortheir Lax pairs.The Lax pair is an important attribute of the integrable soliton equations. They also have useful appli-cations for Darboux integrable systems. For example, in some cases, Lax pairs provide an effective tool forconstructing characteristic integrals. We will now outline a simple algorithm suitable for this purpose.Assume that system of equations (4) admits a Lax pair, i.e. it is a compatibility condition of a pair ofsystems of linear equations Φ n +1 = U n Φ n , (63) Φ n,x = V n Φ n , (64) here the potentials U n and V n of the linear systems are matrices with the following triangular structure U n = a ,n a ,n . . . a N,n a ,n . . . a N,n . . . . . . . . . . . . . . . a
NN,n , (65) V n = b ,n . . . b ,n b ,n . . . . . . . . . . . . . . .b N ,n b N ,n . . . b NN,n . (66)Let P ( k ) n be a product of the matrices P ( k ) n = U n + k U n + k − · · · U n , then obviously we have D x P ( k ) n = V n + k +1 P ( k ) n − P ( k ) n V n . (67)Evidently P ( k ) n is an upper triangular matrix P ( k ) n = p ( k )11 ,n p ( k )12 ,n . . . p ( k )1 N,n p ( k )22 ,n . . . p ( k )2 N,n . . . . . . . . . . . . . . . p ( k ) NN,n . In the equality (67), we select the matrix elements located at the intersection of the first row and the lastcolumn. As a result, we obtain a scalar equality of the form D x p ( k )1 N,n = p ( k )1 N,n ( D k +1 n b − b NN ) . (68)Now we can conclude: Lemma 5.1
If the entries of the matrix (66) satisfy the relations D n b = b and b = b NN then function J = p ( k )1 N,n is an x -integral. Proof 5.1
In such a case (68) implies D x p ( k )1 N,n = 0 . For constructing n -integrals we use the higher order derivatives of the equation system (64) with respectto x . Evidently they solve linear systems as Φ n,xx = ( V n,x + V n )Φ n , Φ n,xxx = ( V n,xx + 2 V n V n,x + V n V n,x + V n ))Φ n . For arbitrary k we find D kx (Φ n ) = R ( k ) n (cid:16) V n , V n,x , . . . , D ( k − x ( V n ) (cid:17) Φ n , where R ( k ) n is a polynomial with constant coefficients on all of its arguments R ( k ) n = r ( k )11 ,n . . . r ( k )21 ,n r ( k )22 ,n . . . . . . . . . . . . . . .r ( k ) N ,n r ( k ) N ,n . . . r ( k ) NN,n . By applying the shift operator D n to the latter equation we obtain D kx Φ n +1 = D n ( R ( k ) n ) U n Φ n . (69)On the other hand by differentiating (63) with respect to x k times we find D kx Φ n +1 = k X j =0 c jk D jx ( U n )( R ( k − j ) n )Φ n . (70) quations (69) and (70) evidently imply D n ( R ( k ) n ) U n = k X j =0 c jk D jx ( U n )( R ( k − j ) n ) . (71)Let us pass in the equality (71), to the matrix elements located at the left lower corner. As a result we arriveat the equation D n ( r ( k ) N ,n ) a = a NN r ( k ) N ,n + a NN,x r ( k − N ,n + a NN,xx r ( k − N ,n + · · · + D ( k ) x ( a NN ) . (72) Lemma 5.2
If the entries of the matrix (65) satisfy the conditions D x a = 0 and a = a NN then function I = r ( k ) N ,n is an n -integral. Proof 5.2
Proof follows right away from equation (72).
To find the desired degenerate boundary condition for (49) we make the following change of the variables: u jn = v jn − log ε for j > , ε > ,u n = v n ,u jn = v jn + log ε for j < . Then (49) converts into v n +1 ,x = v n,x + εe v n +1 − v v − εe v − n +1 − v n ,v n +1 ,x = v n,x + e v n +1 − v n − εe v n +1 − v n ,v − n +1 ,x = v − n,x + εe v − n +1 − v n − e v − n +1 − v − n ,v jn +1 ,x = v jn,x + e v jn +1 − v j +1 n − e v j − n +1 − v jn , for | j | > . Now we take the limit for ε → and obtain an equation v n +1 ,x = v n,x and two semi-infinite lattices whichare not related to each other: v n +1 ,x = v n,x + e v n +1 − v n ,v jn +1 ,x = v jn,x + e v jn +1 − v j +1 n − e v j − n +1 − v jn , for j ≥ , and v − n +1 ,x = v − n,x − e v − n +1 − v − n ,v jn +1 ,x = v jn,x + e v jn +1 − v j +1 n − e v j − n +1 − v jn , for j ≤ − . By applying this manipulation to the lattice at two fixed points j = − and j = N + 1 we obtain a finitefield system of the form (3) which is a desired reduction of the equation (49): u n +1 ,x = u n,x + e u n +1 − u n ,u jn +1 ,x = u jn,x + e u jn +1 − u j +1 n − e u j − n +1 − u jn , < j < N − ,u Nn +1 ,x = u Nn,x − e u N − n +1 − u Nn . (73)In [42] it was proved that system (73) is integrable in the sense of Darboux for arbitrary natural N . Werewrite the Lax pair (56) by means of the obtained boundary conditions ϕ n +1 = − ϕ n + e u n − u n +1 ϕ n ,ϕ jn +1 = − ϕ j +1 n + e u jn − u jn +1 ϕ jn , < j < N − ,ϕ Nn +1 = e u Nn − u Nn +1 ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = − e u j − n − u jn ϕ j − n , < j < N − ,ϕ Nn,x = − e u N − n − u Nn ϕ N − n . When deriving it from (56) we set ϕ − n = 0 , ϕ N +1 n = 0 . Let us study in more details the case N = 2 thatcorresponds to the system u n +1 ,x = u n,x + e u n +1 − u n ,u n +1 ,x = u n,x + e u n +1 − u n − e u n +1 − u n ,u n +1 ,x = u n,x − e u n +1 − u n (74) dmitting the Lax pair ϕ n +1 = − ϕ n + e u n − u n +1 ϕ n ,ϕ n +1 = − ϕ n + e u n − u n +1 ϕ n ,ϕ n +1 = e u n − u n +1 ϕ n , ϕ n,x = 0 ,ϕ n,x = − e u n − u n ϕ n ,ϕ n,x = − e u n − u n ϕ n . (75)We give here x -integrals and n -integrals of system (74) found earlier in [42] J = e u n +2 − u n +3 + e u n +1 − u n +2 + e u n − u n +1 ; J = e u n +1 − u n +2 + u n +1 − u n +2 + e u n +1 − u n +2 + u n − u n +1 + e u n +1 − u n +2 + u n − u n +1 ; J = e u n − u n +1 + u n − u n +1 + u n − u n +1 and, respectively I = u n,x + u n,x + u n,x ; I = u n,xx + 2 u n,xx + u n,x u n,x + u n,x u n,x + u n,x u n,x ; I = u n,xxx + u n,x u n,xx + u n,xx u n,x + u n,xx u n,x + u n,x u n,x u n,x . Note that they can be readily derived also from the Lax pair (75) due to Lemma 5.2 and Lemma 5.1above.
By applying the manipulations similar to that fulfilled in the previous section by using the same change ofthe variables one can show that lattice (50) is reduced to the following finite field system which is of the form(3) u n +1 ,x = u n,x − e u n e u n + e u n +1 e u n ,u jn +1 ,x = u jn,x + e uj − n +1 e ujn +1 − e uj − n +1 e ujn − e ujn e uj +1 n + e ujn +1 e uj +1 n , < j < N − ,u Nn +1 ,x = u Nn,x + e uN − n +1 e uNn +1 − e uN − n +1 e uNn and its Lax pair is obtained from the system (57) by imposing additional conditions ϕ − n = 0 , ϕ N +1 n = 0 : ϕ n +1 = − e − u n +1 + u n (cid:0) ϕ n − ϕ n (cid:1) + ϕ n ,ϕ jn +1 = − e − u jn +1 + u jn (cid:0) ϕ j +1 n − ϕ jn (cid:1) + ϕ j +1 n , < j < N − ,ϕ Nn +1 = e − u Nn +1 + u Nn ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = e u j − n − u jn (cid:0) ϕ j − n − ϕ jn (cid:1) , < j < N − ,ϕ Nn,x = e u N − n − u Nn (cid:0) ϕ N − n − ϕ Nn (cid:1) . We will investigate in more detail the case when N = 1 . Then we obtain a system ( u n +1 ,x = u n,x − e u n − u n + e u n +1 − u n ,u n +1 ,x = u n,x + e u n +1 − u n +1 − e u n +1 − u n (76)and its Lax pair ( ϕ n +1 = − e − u n +1 + u n (cid:0) ϕ n − ϕ n (cid:1) + ϕ n ,ϕ n +1 = e − u n +1 + u n ϕ n , ( ϕ n,x = 0 ,ϕ n,x = e u n − u n (cid:0) ϕ n − ϕ n (cid:1) . (77) x -integrals and n -integrals of system (76) have the form J = (cid:16) e u n +1 − u n − (cid:17) (cid:16) e u n +1 − u n − (cid:17) ; J = (cid:16) e u n − u n +1 − (cid:17) (cid:16) e u n +1 − u n +2 − (cid:17) (78) nd, respectively I = u n,x + u n,x − e u n − u n ; I = u n,xx + u n,x u n,x − e u n − u n u n,x . (79)Let us briefly discuss the methods of searching the integrals. Note that for some cases the Lax pairgives a convenient tool for solving the problem, however this way is applied not always. Therefore in theother cases we use algebraic method. We explain both approaches with the example of the system (76). Forconstructing the x -integrals we use the method of characteristic algebras and we use the Lax pair for findingthe n -integrals.Assume that function H ( u n , u n , u n +1 , u n +1 , . . . ) is an x –integral for the system (76). Then according tothe definition the relation D x H ( u n , u n , u n +1 , u n +1 , . . . ) = 0 should hold. Due to the chain rule equation (5.2) implies K H = 0 , where K = u n,x ∂∂u n + u n,x ∂∂u n + u n +1 ,x ∂∂u n +1 + u n +1 ,x ∂∂u n +1 + u n − ,x ∂∂u n − + u n − ,x ∂∂u n − + . . . . We exclude the terms like u n ± i,x , u n ± i,x due to the system (76) and get K = u n,x ∂∂u n + u n,x ∂∂u n + (cid:16) u n,x − e u n − u n + e u n +1 − u n (cid:17) ∂∂u n +1 + (cid:16) u n,x + e u n +1 − u n +1 − e u n +1 − u n (cid:17) ∂∂u n +1 + . . . . Obviously equation K H = 0 is overdetermined since the coefficients of the equation depend on u n,x and u n,x while the solution H does not depend on these variables, in other words we are interested on the solution H of the equation (5.2) which solves in addition two more equations X H = 0 and X H = 0 , where X = ∂∂u n,x , X = ∂∂u n,x . More precisely we get a system of the first order linear partial differentialequations for one and the same unknown H K H = 0 ,X H = 0 ,X H = 0 . (80)Evidently operator K is represented as a linear combination of the vector fields K = u n,x Y + u n,x Y + W, where the coefficients of the operators Y , Y , W do not depend on the variables u n,x , u n,x : Y = ∂∂u n + ∂∂u n +1 + ∂∂u n +2 ,Y = ∂∂u n + ∂∂u n +1 ,W = (cid:16) e u n +1 − u n − e u n − u n (cid:17) ∂∂u n +1 + (cid:16) e u n +1 − u n +1 − e u n +1 − u n (cid:17) ∂∂u n +1 + (cid:16) e u n +2 − u n +1 − e u n +1 − u n +1 + e u n +1 − u n − e u n − u n (cid:17) ∂∂u n +2 . Thus (80) is reduced to the form Y H = 0 ,Y H = 0 ,W H = 0 . (81) t is checked that system (81) is closed, i.e. all of the commutators [ Y , Y ] , [ Y , W ] and [ Y , W ] are linearlyexpressed in terms of Y , Y , W such that: [ Y , Y ] = 0 , [ Y , W ] = W, [ Y , W ] = − W. Since the system contains three equations then in order to get two functionally independent solutions welook for a solution depending on five variables H = H (cid:0) u n , u n , u n +1 , u n +1 , u n +2 (cid:1) . Then system (81) turns into H u n + H u n +1 + H u n +2 = 0 ,H u n + H u n +1 = 0 , (cid:16) e u n +1 − u n − e u n − u n (cid:17) H u n +1 + (cid:16) e u n +1 − u n +1 − e u n +1 − u n (cid:17) H u n +1 + (cid:16) e u n +2 − u n +1 − e u n +1 − u n +1 + e u n +1 − u n − e u n − u n (cid:17) H u n +2 = 0 . (82)In order to solve the system we have to reduce it due to Jacobi method to a normal (triangular) form (cid:16) e u n +1 − u n − e u n − u n (cid:17) H u n − (cid:16) e u n +1 − u n +1 − e u n +1 − u n (cid:17) H u n +1 − (cid:16) e u n +2 − u n +1 − e u n +1 − u n +1 (cid:17) H u n +2 = 0 ,H u n + H u n +1 = 0 , (cid:16) e u n +1 − u n − e u n − u n (cid:17) H u n +1 + (cid:16) e u n +1 − u n +1 − e u n +1 − u n (cid:17) H u n +1 + (cid:16) e u n +2 − u n +1 − e u n +1 − u n +1 + e u n +1 − u n − e u n − u n (cid:17) H u n +2 = 0 . Then by applying the consecutive integration according to the Jacobi algorithm we get the result (78).To construct the complete set of the n -integrals for the system (76) we use the Lax pair. At first we makea change of the variables ϕ n = e − u n ψ n , ϕ n = e − u n ψ n in the system (77) and get a new form of the Lax pair Ψ n,x = V n Ψ n , Ψ n +1 = U n Ψ n , (83)where Ψ n = (cid:0) ψ n , ψ n (cid:1) T and V n = (cid:18) u n,x u n,x − e u n − u n (cid:19) , U n = (cid:18) e u n +1 − u n − e u n − u n (cid:19) . Now the Lax pair (83) satisfies all the requests of the Lemma 5.2. Therefore function r ( k ) N ,n is an n -integral.In this case we have r (2)21 ,n = u n,x + u n,x − e u n − u n ,r (3)21 ,n = u n,xx + 2 u n,xx + ( u n,x ) + u n,x u n,x + ( u n,x ) − e u n − u n u n,x + e ( u n − u n ) . After some slight simplifications we find integrals (79).
In equation (51) we put u n = c , u N +1 n = c N : u n +1 ,x = u n,x ( u n +1 ) c u n ,u jn +1 ,x = u jn,x ( u jn +1 ) u j − n +1 u j +1 n , < j < N − ,u Nn +1 ,x = u Nn,x ( u Nn +1 ) u N − n +1 c N . e will rewrite the Lax pair (58) with the above constraint for the field variables and additional constraintfor the eigenfunctions ϕ − n = 0 , ϕ N +1 n = 0 ϕ n +1 = − c u n (cid:0) ϕ n − ϕ n (cid:1) ,ϕ jn +1 = − u jn +1 u j +1 n (cid:0) ϕ jn − ϕ j +1 n (cid:1) , < j < N − ,ϕ Nn +1 = − u Nn +1 c N ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = − u jn,x u jn (cid:0) ϕ j − n − ϕ jn (cid:1) , < j < N − ,ϕ Nn,x = − u Nn,x u Nn (cid:0) ϕ N − n − ϕ Nn (cid:1) . We put N = 1 . Then we obtain an equation u n +1 ,x = u n,x (cid:0) u n +1 (cid:1) (84)having the Lax pair ( ϕ n +1 = − u n (cid:0) ϕ n − ϕ n (cid:1) ,ϕ n +1 = − u n +1 ϕ n , ( ϕ n,x = 0 ,ϕ n,x = − u n,x u n (cid:0) ϕ n − ϕ n (cid:1) .x -integrals and n -integrals of equation (84) have the form J = u n + 1 u n +1 and, respectively I = u n,xxx u n,x −
32 ( u n,xx ) ( u n,x ) . In the equation (52) we put u n = c , u N +1 n = c N and obtain a system u n +1 ,x = u n,x ( u n +1 − u n ) c − u n ,u jn +1 ,x = u jn,x ( u jn +1 − u j +1 n )( u j − n +1 − u jn ) , < j < N − ,u Nn +1 ,x = u Nn,x u Nn +1 − c N ( u N − n +1 − u Nn ) with the Lax pair ϕ n +1 = (cid:0) u n − c (cid:1) ϕ n + ϕ n ,ϕ jn +1 = (cid:16) u j +1 n − u jn +1 (cid:17) ϕ jn + ϕ j +1 n , < j < N − ,ϕ Nn +1 = (cid:0) c N − u Nn +1 (cid:1) ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = − u jn,x ϕ j − n , < j < N − ,ϕ Nn,x = − u Nn,x ϕ N − n . In the particular case N = 2 we get a system u n +1 ,x = u n,x ( u n +1 − u n ) c − u n ,u n +1 ,x = u n,x u n +1 − c ( u n +1 − u n ) (85)with the Lax pair ϕ n +1 = (cid:0) u n − c (cid:1) ϕ n + ϕ n ,ϕ n +1 = (cid:0) u n − u n +1 (cid:1) ϕ n + ϕ n ,ϕ n +1 = (cid:0) c − u n +1 (cid:1) ϕ n , ϕ n,x = 0 ,ϕ n,x = − u n,x ϕ n ,ϕ n,x = − u n,x ϕ n . (86) x -integrals and n -integrals of system (85) have the form J = (cid:0) u n − c (cid:1) (cid:0) u n +1 − c (cid:1) (cid:0) u n − u n +1 (cid:1) ; J = u n +1 u n +2 − c u n +2 − u n +1 u n +1 − c u n + u n u n +1 (87) nd, respectively I = u n,xxx u n,x + u n,xxx u n,x − u n,xx u n,x ! − u n,xx u n,xx u n,x u n,x − u n,xx u n,x ! ; I = u n,xxxx u n,x − u n,xx u n,xxx (cid:0) u n,x (cid:1) − u n,xx u n,xxx u n,x u n,x + 23 u n,xx u n,xxx (cid:0) u n,x (cid:1) + 13 u n,xx u n,xxx u n,x u n,x + 329 u n,xx u n,x ! + (cid:0) u n,xx (cid:1) u n,xx (cid:0) u n,x (cid:1) u n,x − u n,xx (cid:0) u n,xx (cid:1) u n,x (cid:0) u n,x (cid:1) − u n,xx u n,x ! . (88)Now we briefly explain how the integrals (87) and (88) were constructed. Let us begin with x -integrals.The Lax pair can be written as Ψ n,x = V n Ψ n , Ψ n +1 = U n Ψ n , (89)where V n = − u n,x − u n,x , U n = u n − c u n − u n +1
10 0 c − u n +1 . Since the Lax pair (89) satisfies the conditions of Lemma 5.1, we can use it for deriving the x -integrals of thesystem (85). To this end we evaluate the product P ( k ) n = U n + k U n + k − · · · U n and take its entry located at the right upper corner. It is easily checked that for k = 1 and k = 2 we gettrivial integrals J = 1 and J = c − c . For k = 3 and k = 4 we obtain p (3)13 ,n = (cid:0) u n +3 − c (cid:1) (cid:0) u n +2 − c (cid:1) + (cid:0) u n +2 − c (cid:1) (cid:0) u n +2 − c (cid:1) + (cid:0) u n +2 − c (cid:1) (cid:0) u n +1 − c (cid:1) ,p (4)13 ,n = (cid:0) u n +3 − c (cid:1) (cid:0) u n +2 − c (cid:1) (cid:0) u n +1 − c (cid:1) − (cid:0) u n +3 − c (cid:1) (cid:0) u n +1 − c (cid:1) (cid:0) u n +1 − c (cid:1) − (cid:0) u n +3 − u n +2 (cid:1) (cid:0) u n +2 − u n +1 (cid:1) (cid:0) u n +1 − c (cid:1) + ( c − c ) (cid:0) u n +1 − c (cid:1) (cid:0) u n − c (cid:1) . After some slight simplifications we reduce them into integrals J and J given in (87).The next step is to explain the algebraic method of looking for the n -integrals. Let us emphasize thatthis task is more difficult since the algebra L n has a complicated structure.It is convenient to introduce new notations by setting u n := u n and v n := u n . Thus we convert the system(85) to the form ( u n +1 ,x = F ,v n +1 ,x = F , (90)where F = u n,x ( u n +1 − v n ) c − u n , F = v n,x v n +1 − c ( u n +1 − v n ) . At first we evaluate the subalgebra L (1) n , generated by the vector fields Y , , Y , , X , and X , , where Y , = ∂∂u n + D − n (cid:18) ∂F n ∂u n +1 (cid:19) ∂∂u n,x + D − n (cid:18) ∂F n ∂u n +1 (cid:19) ∂∂v n,x + D − n ∂F n,x ∂u n +1 ! ∂∂u n,xx + D − n ∂F n,x ∂u n +1 ! ∂∂v n,xx + . . . ,Y , = ∂∂v n + D − n (cid:18) ∂F n ∂v n +1 (cid:19) ∂∂u n,x + D − n (cid:18) ∂F n ∂v n +1 (cid:19) ∂∂v n,x + D − n ∂F n,x ∂v n +1 ! ∂∂u n,xx + D − n ∂F n,x ∂v n +1 ! ∂∂v n,xx + . . . ,X , = ∂∂u n − , X , = ∂∂v n − . Explicit form of the operators Y , and Y , are Y , = ∂∂u n + u n,x u n − v n − ∂∂u n,x − v n,x u n − v n − ∂∂v n,x u n − v n − (cid:18) u n,xx + u n,x v n,x v n − c (cid:19) ∂∂u n,xx − u n − v n − v n,xx + 2 v n,x v n − c ! ∂∂v n,xx + . . . ,Y , = ∂∂v n + v n,x v n − c ∂∂v n,x + v n,xx v n − c ∂∂v n,xx + v n,xxx v n − c ∂∂v n,xxx . . . . It can be verified by analyzing the further terms in the series for Y , that it is decomposed into a sum Y , = Y + 1 u n − v n − Z , of two operators belonging the algebra L (1) n , since Z = ( u n − v n − ) [ X , , Y , ] . Here Y = ∂∂u n and thecoefficients in Z = u n,x ∂∂u n,x − v n,x ∂∂v n,x + (cid:18) u n,xx + u n,x v n,x v n − c (cid:19) ∂∂u n,xx − v n,xx + 2 v n,x v n − c ! ∂∂v n,xx + . . . do not depend on the variables u n − , v n − . Therefore operators Y , Z and Y , commute with X , and X , . Moreover the operators X , , X , , Y , Z and Y , constitute a basis in the algebra L (1) n . Thus algebra L (1) n is of dimension 5. However solution of the system X , I = 0 , X , I = 0 , Y I = 0 , Z I = 0 , Y , I = 0 does not produce any n –integral for the system (90), since the relation D n I = I is not satisfied.Therefore, we need to pass to the consideration of the second subalgebra L (2) n of the characteristic algebra L n generated by the operators X , , X , , X , , X , , Y , , Y , , Y , , Y , . In other words we add to L (1) n four extra operators X , = ∂∂u n − , X , = ∂∂v n − , Y , = D − n Y , D n , Y , = D − n Y , D n . Two last operators are of the form Y , = u n,x ( v n − − c )( u n − − v n − )( u n − − c ) ∂∂u n,x − v n,x u n − − v n − ∂∂v n,x + 1 u n − − v n − u n,xx ( v n − − c ) u n − − c + 2 u n,x ( v n − − c )( u n − − c )( u n − v n − ) + u n,x v n,x ( u n − c )( v n − − c )( v n − c ) ! ∂∂u n,xx − u n − − v n − v n,xx + 2 v n,x ( u n − c )( v n − c )( v n − − c ) + u n,x v n,x ( v n − − c )( u n − − c )( u n − v n − ) ! ∂∂v n,xx + . . . ,Y , = − u n,x u n − v n − ∂∂u n,x + v n,x ( u n − c )( u n − v n − )( v n − − c ) ∂∂v n,x − u n − v n − (cid:18) u n,xx + u n,x v n,x ( u n − c )( v n − c )( v n − − c ) (cid:19) ∂∂u n,xx + u n − c ( u n − v n − )( v n − − c ) v n,xx + 2 v n,x v n − c ! ∂∂v n,xx + . . . . It can be proved that these operators are represented as linear combinations of the operators Z = v n,x ∂∂v n,x − u n,x v n,x v n − c ∂∂u n,xx + v n,xx + 2 v n,x v n − c ! ∂∂v n,xx + . . . ,Y = u n,x ∂∂u n,x + u n,xx ∂∂u n,xx + u n,xxx ∂∂u n,xxx + . . . ,Z = 2 u n,x ∂∂u n,xx − u n,x v n,x ∂∂v n,xx + u n,x u n,xx + 3 u n,x v n,x v n − c ! ∂∂u n,xxx + . . . ,Z = 32 v n,x ∂∂v n,xx + (cid:18) u n,xx − u n,x v n,xx v n,x + 32 u n,x v n,x v n − c (cid:19) ∂∂u n,xxx + . . . , = 3 u n,x v n,x ∂∂u n,xxx − v n,x ∂∂v n,xxx + (14 v n,x u n,xx + 4 u n,x v n,xx ) ∂∂u n,xxxx + . . . which are simpler, since they commute with X , , X , , X , , X , .We prove that system of equations Y I = 0 , Y , I = 0 , Y I = 0 , Z I = 0 , Z I = 0 , Z I = 0 , Y I = 0 (91)is closed or in other words algebra L (2) n is of finite dimension. We look for the solution of (91) in the form I = I ( u n , v n , u n,x , v n,x , u n,xx , v n,xx , u n,xxx , v n,xxx , u n,xxxx ) . We reduce it to normal form and solve by Jacobi method. As a result we get (88), where one has to replace u n := u n and v n := u n . In equation (53) we put u n = c , u N +1 n = c N : u n +1 ,x = u n,x u n +1 ( u n +1 − u n ) u n ( c − u n ) ,u jn +1 ,x = u jn,x u jn +1 ( u jn +1 − u j +1 n ) u j +1 n ( u j − n +1 − u jn ) , < j < N − ,u Nn +1 ,x = u Nn,x u Nn +1 ( u Nn +1 − c N ) c N ( u N − n +1 − u Nn ) . We rewrite the Lax pair (60) with the above restrictions and by setting the conditions ϕ − n = 0 , ϕ N +1 n = 0 ϕ n +1 = (cid:16) − c u n (cid:17) ϕ n − c u n ϕ n ,ϕ jn +1 = (cid:18) − u jn +1 u j +1 n (cid:19) ϕ jn − u jn +1 u j +1 n ϕ j +1 n , < j < N − ,ϕ Nn +1 = (cid:16) − u Nn +1 c N (cid:17) ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = ϕ jn,x = u jn,x u jn (cid:0) ϕ j − n + ϕ jn (cid:1) , < j < N − ,ϕ Nn,x = u Nn,x u Nn (cid:0) ϕ N − n + ϕ Nn (cid:1) . We concentrate on a simplest case by taking N = 1 . Then we get an equation u n +1 ,x = u n,x u n +1 (cid:0) u n +1 − c (cid:1) c ( c − u n ) , (92)and its Lax pair ϕ n +1 = (cid:16) − c u n (cid:17) ϕ n − c u n ϕ n ,ϕ n +1 = (cid:16) − u n +1 c (cid:17) ϕ n , ( ϕ n,x = 0 ,ϕ n,x = u n,x u n (cid:0) ϕ n + ϕ n (cid:1) .x -integrals and n -integrals of equation (92) have the form J = u n u n +1 + c u n +1 + u n c and, respectively I = u n,xxx u n,x −
32 ( u n,xx ) ( u n,x ) . In equation (54) we put u n = c , u N +1 n = c N : u n +1 ,x = u n,x ( c − u n +1 )( u n +1 − u n ) ( c − u n )( u n − u n ) ,u jn +1 ,x = u jn,x ( u j − n +1 − u jn +1 )( u jn +1 − u j +1 n )( u j − n +1 − u jn )( u jn − u j +1 n ) , < j < N − ,u Nn +1 ,x = u Nn,x ( u N − n +1 − u Nn +1 )( u Nn +1 − c N )( u N − n +1 − u Nn ) ( u Nn − c N ) . e rewrite the Lax pair (61) with the constraint u n = c , u N +1 n = c N , ϕ − n = 0 , ϕ N +1 n = 0 ϕ n +1 = c − u n c − u n ϕ n + (cid:16) − c − u n c − u n (cid:17) ϕ n ,ϕ jn +1 = u jn +1 − u j +1 n u jn − u j +1 n ϕ jn + (cid:18) − u jn +1 − u j +1 n u jn − u j +1 n (cid:19) ϕ j +1 n , < j < N − ,ϕ Nn +1 = u Nn +1 − c N u Nn − c N ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = u jn,x u j − n − u jn (cid:0) ϕ j − n + ϕ jn (cid:1) , < j < N − ,ϕ Nn,x = u Nn,x u N − n − u Nn (cid:0) ϕ N − n + ϕ Nn (cid:1) , By putting N = 2 we arrive at a system u n +1 ,x = u n,x ( c − u n +1 )( u n +1 − u n ) ( c − u n )( u n − u n ) ,u n +1 ,x = u n,x ( u n +1 − u n +1 )( u n +1 − c N )( u n +1 − u n ) ( u n − c ) . (93)admitting the Lax pair ϕ n +1 = c − u n c − u n ϕ n + (cid:16) − c − u n c − u n (cid:17) ϕ n ,ϕ n +1 = u n +1 − u n u n − u n ϕ n + (cid:16) − u n +1 − u n u n − u n (cid:17) ϕ n ,ϕ n +1 = u n +1 − c u n − c ϕ n , ϕ n,x = 0 ,ϕ n,x = u n,x u n − u n (cid:0) ϕ n + ϕ n (cid:1) ,ϕ n,x = u n,x u n − u n (cid:0) ϕ n + ϕ n (cid:1) .x -integrals and n -integrals of system (93) have the form J = (cid:0) u n − c (cid:1) (cid:0) u n − c (cid:1) (cid:0) u n +1 − u n (cid:1)(cid:0) u n +1 − u n (cid:1) (cid:0) u n +1 − u n (cid:1) ; J = (cid:0) u n +2 − c (cid:1) (cid:0) u n − c (cid:1) (cid:0) u n +1 − u n +1 (cid:1) (cid:0) u n +2 − u n +1 (cid:1)(cid:0) u n +1 − c (cid:1) (cid:0) u n +2 − c (cid:1) (cid:0) u n +1 − u n (cid:1) (cid:0) u n +2 − u n +1 (cid:1) and, respectively I = u n,x u n,x ( u n − c ) ( u n − c ) ( u n − u n ) ; I = u n,xx u n,x − u n,x u n − c + u n,x (cid:0) u n − c (cid:1) ( u n − c ) ( u n − u n ) . In the equations (55), (62) we put u n = c , u N +1 n = c N , ϕ − n = 0 , ϕ N +1 n = 0 and get a system u n +1 ,x = u n,x sinh ( c − u n +1 ) sinh ( u n +1 − u n ) sinh( c − u n ) sinh( u n − u n ) ,u jn +1 ,x = u jn,x sinh ( u j − n +1 − u jn +1 ) sinh ( u jn +1 − u j +1 n ) sinh ( u j − n +1 − u jn ) sinh ( u jn − u j +1 n ) , < j < N − ,u Nn +1 ,x = u Nn,x sinh ( u N − n +1 − u Nn +1 ) sinh ( u Nn +1 − c N ) sinh ( u N − n +1 − u Nn ) sinh( u Nn − c N ) . with the Lax pair ϕ n +1 = e ( c − u n ) − e ( c − u n ) − ϕ n + (cid:18) − e ( c − u n ) − e ( c − u n ) − (cid:19) ϕ n ,ϕ jn +1 = e ( ujn +1 − uj +1 n ) − e ( ujn − uj +1 n ) − ϕ jn + (cid:18) − e ( ujn +1 − uj +1 n ) − e ( ujn − uj +1 n ) − (cid:19) ϕ j +1 n , < j < N − ,ϕ Nn +1 = e ( uNn +1 − cN ) − e ( uNn − cN ) − ϕ Nn , ϕ n,x = 0 ,ϕ jn,x = u jn,x e ( uj − n − ujn ) − (cid:0) ϕ j − n − ϕ jn (cid:1) , < j < N − ,ϕ Nn,x = u Nn,x e ( uN − n − uNn ) − (cid:0) ϕ N − n − ϕ Nn (cid:1) . e put N = 1 and obtain an equation u n +1 ,x = u n,x sinh (cid:0) c − u n +1 (cid:1) sinh (cid:0) u n +1 − c (cid:1) sinh ( c − u n ) sinh ( u n − c ) , (94)having the Lax pair ϕ n +1 = e ( c − u n ) − e ( c − u n ) − ϕ n + (cid:18) − e ( c − u n ) − e ( c − u n ) − (cid:19) ϕ n ,ϕ n +1 = e ( u n +1 − c ) − e ( u n − c ) − ϕ n , ϕ n,x = 0 ,ϕ n,x = u n,x e ( c − u n ) − (cid:0) ϕ n − ϕ n (cid:1) .x -integrals and n -integrals of equation (94) have the form J = 1sinh( c − c ) (cid:20) arctanh (cid:18) (cosh( c − c ) + 1) tanh( c + c − u n )sinh( c − c ) (cid:19) − arctanh (cid:18) (cosh( c − c ) + 1) tanh( c + c − u n +1 )sinh( c − c ) (cid:19)(cid:21) and, respectively I = u n,x sinh( c − u n ) sinh( c − u n ) . Here we present a list of Darboux integrable scalar equations obtained from the equations (51)–(55) byimposing degenerate boundary conditions: u n +1 ,x = u n,x ( u n +1 ) ;2) u n +1 ,x = u n,x c − u n +1 u n ;3) u n +1 ,x = u n,x ( u n +1 + c )( u n +1 + c ) u n ;4) u n +1 ,x = u n,x u n +1 ( u n +1 − c ) u n ( u n − c ) ;5) u n +1 ,x = u n,x sinh( u n +1 ) sinh( u n +1 − c )sinh( u n ) sinh( u n − c ) . Apparently they are new.
Conclusions
The paper proposes an algebraic method for the classification of integrable cases of differential-differenceequations of the form (1). It is shown that all known integrable equations of this form given in the listpresented in [17] admit reductions as finite systems of differential-difference equations integrable in sense ofDarboux. We used this property of the equation (1) as the basis for the classification criterion. It is convenientto formalize the property of Darboux integrability in terms of its characteristic Lie-Rinehart algebras. Thearticle gives a definition of these algebras for systems of differential-difference equations and investigates theirbasic properties that may be needed in solving the classification problem.
Funding
One of the authors A.R. Khakimova was supported in part by Young Russian Mathematics award.
References [1] Bogdanov, L.V. “Dunajski-Tod equation and reductions of the generalized dispersionless 2DTL hierar-chy.”
Phys. Lett. A
2] Bogdanov, L.V. and B.G. Konopelchenko. “On dispersionless BKP hierarchy and its reductions.”
J.Nonlinear Math. Phys.
12, suppl. 1 (2005): 64–73.[3] Calderbank, D.M.J. “Integrable background geometries.”
SIGMA
10, no. 034 (2014): 51 pp.[4] Calderbank, D.M.J. and B. Kruglikov. “Integrability via geometry: dispersionless differential equations inthree and four dimensions.”
Commun. Math. Phys. (2020). https://doi.org/10.1007/s00220-020-03913-y[5] Cl ´ e ry, F. and E.V. Ferapontov. “Dispersionless Hirota equations and the genus 3 hyper elliptic divisor.” Commun. Math. Phys.
Lecons sur la th´eorie g´en´erale des surfaces et les applications geometriques du calcul in-finitesimal . Paris: Gauthier-Villars, 1896.[7] Doubrov, B., E.V. Ferapontov, B. Kruglikov, and V.S. Novikov. “On integrability in Grassmann ge-ometries: integrable systems associated with fourfolds in Gr (3 , .” Proc. London Math. Soc. Gr (4 , .” Int. Math. Res. Not. IMRN
J. Geom. Phys.
37, no. 1-2 (2001): 63–93.[10] Ferapontov, E.V., I.T. Habibullin, M.N. Kuznetsova, and V.S. Novikov. “On a class of 2D integrablelattice equations.”
J. Math. Phys.
61, no. 7 (2020): 073505, 15 pp.[11] Ferapontov, E.V., L. Hadjikos, and K.R. Khusnutdinova. “Integrable equations of the dispersionlessHirota type and hypersurfaces in the Lagrangian Grassmannian.”
Int. Math. Res. Not. IMRN
J. Math. Phys.
45, no. 6 (2004): 2365–2377.[13] Ferapontov, E.V. and B.S. Kruglikov. “Dispersionless integrable systems in 3D and Einstein-Weyl ge-ometry”,
J. Differ. Geom.
97, no. 2 (2014): 215–254.[14] Ferapontov, E.V. and A. Moro. “Dispersive deformations of hydrodynamic reductions of 2D dispersionlessintegrable systems.”
J. Phys. A: Math. Theor.
42 (2009) 035211, 15 pp.[15] Ferapontov, E.V., A. Moro, and V. S. Novikov. “Integrable equations in 2+1 dimensions: deformationsof dispersionless limits.”
J. Phys. A: Math. Theor.
42 (2009) 18 pp.[16] Ferapontov, E.V. and A.V. Odesskii. “Integrable Lagrangians and modular forms.”
J. Geom. Phys.
Int. Math. Res. Not. IMRN
Integration of Classical Series An, Bn, Cn, of Exponential Systems .Krasnoyarsk: Krasnoyarsk State Pedagogical University Press, 2001.[19] Gibbons, J. and Y. Kodama. “A method for solving the dispersionless KP hierarchy and its exactsolutions. II.”
Phys. Lett. A
Phys. Lett. A
Phys. Lett. A
258 (1999): 263–270.[21] Goursat, M.E. “Recherches sur quelques e´quations aux de´rive´es partielles du second ordre.”
Annales dela Facult? des sciences de Toulouse: Math?matiques
1, no. 2 (1899): 31–78, 439–463.[22] Habibullin, I. “Characteristic Lie rings, finitely-generated modules and integrability conditions for (2+1)-dimensional lattices.”
Phys. Scr.
87, no. 6 (2013): 065005.[23] Habibullin, I.T. and M.N. Kuznetsova. “A classification algorithm for integrable two-dimensional latticesvia Lie-Rinehart algebras.”
Theoret. and Math. Phys.
J. Phys. A: Math. Theor.
53, no. 39 (2020): 395203 , 25 pp.[25] Habibullin, I., A. Pekcan, and N. Zheltukhina. “On Some Algebraic Properties of Semi-Discrete Hyper-bolic Type Equations.”
Turkish J. Math.
32, no. 3 (2008): 277–292.[26] Habibullin, I.T. and M.N. Poptsova. “Classification of a Subclass of Two-Dimensional Lattices via Char-acteristic Lie Rings.”
SIGMA
13, no. 073 (2017): 26 pp.[27] Habibullin, I., N. Zheltukhina, and A. Pekcan.“On the classification of Darboux integrable chains.”
J.Math. Phys.
49, no. 10 (2008): 1–39.
28] Habibullin, I., N. Zheltukhina, and A. Pekcan. “Complete list of Darboux integrable chains of the form t ,x = t x + d ( t, t ) .” J. Math. Phys.
50, no. 10 (2009): 1–23.[29] Habibullin, I., N. Zheltukhina, and A. Sakieva. “On Darboux-integrable semi-discrete chains.”
J. Phys.A: Math. Theor.
43, no. 43 (2010): 434017, 14 pp.[30] Habibullin, I., K. Zheltukhin, and M. Yangubaeva. “Cartan matrices and integrable lattice Toda fieldequations”,
J. Phys. A: Math. Theor.
44, no. 46 (2011): 465202.[31] Hitchin, N.J. “Complex manifolds and Einstein’s equations, Twistor geometry and nonlinear systems.”
Lecture Notes in Math. τ -function of the universal whitham hierarchy, matrix models and topological fieldtheories.” Comm. Pure Appl. Math.
47, no. 4 (1994): 437–475.[33] Leznov, A.N. and M.V. Savel’ev.
Group Methods of Integration of Nonlinear Dynamical Systems . Nauka:Moscow (in Russian), 1985.[34] Mikhailov A.V., M.A. Olshanetsky, and A.M. Perelomov. “Two-dimensional generalized Toda lattice.”
Commun. Math. Phys.
79, no. 4 (1981): 473–488.[35] Odesskii, A.V. and V.V. Sokolov. “Integrable pseudopotentials related to generalized hypergeometricfunctions.”
Sel. Math. New Ser.
16, (2010): 145–172.[36] Pavlov, M.V. “New integrable (2+1)-equations of hydrodynamic type.”
Russian Math. Surv.
58, no. 2(2003): 386–387.[37] Pavlov, M.V. “Classifying integrable egoroff hydrodynamic chains.”
Theoret. and Math. Phys.
Gen. Relat. Gravit.
7, no. 1 (1976): 31–52.[39] Poptsova, M.N. and I.T. Habibullin. “Algebraic properties of quasilinear two-dimensional lattices con-nected with integrability.”
Ufa Math. J.
10, no. 3 (2018): 86–105.[40] Rinehart, G. “Differential forms for general commutative algebras.”
Trans. Amer. Math. Soc.
108 (1963):195–222.[41] Shabat, A.B. and R.I. Yamilov. “Exponential Systems of Type I and the Cartan Matrices.”
Preprint,Soviet Academy of Sciences, Bashkirian Branch, Ufa (in Russian) (1981): 22 pp.[42] Smirnov, S.V. “Darboux integrability of discrete two-dimensional Toda lattices.”
Theoret. and Math.Phys. u xy = e u + vu y , v xy = − e u + vv y ”, Ufa Math. J.
5, no. 3 (2013): 20–27.[44] Ward, R.S. “Einstein-Weyl spaces and SU (infinity) Toda fields.”
Class. Quantum Grav.
7, no. 4 (1990):95–98.[45] Zakharov, V.E. “Dispersionless limit of integrable systems in 2+1 dimensions.”
Singular Limits of Dis-persive Waves, NATO ASI series
Vestnik UGATU (inRussian)
9, no. 7(25) (2007): 83–89.[47] Zhiber, A.V. and V.V. Sokolov. “Exactly integrable hyperbolic equations of Liouville type.”
RussianMath. Surveys
56, no. 1 (2001): 61–101.56, no. 1 (2001): 61–101.