Miura-Type Transformations for Integrable Lattices in 3D Ismagil

: This article studies a class of integrable semi-discrete equations with one continuous and two discrete independent variables. At present, in the literature there are nine integrable equations of the form u jn + 1, x = f ( u jn , x , u j + 1 n , u jn , u jn + 1 , u j − 1 n + 1 ) up to point transformations. An efﬁcient method based on some relation that generalizes the notion of the local conservation law is proposed for searching for Miura-type transformations relating to semi-discrete equations in 3D. The efﬁciency of the method is illustrated with the equations from the list. For one of the equations, which is little studied, the continuum limit is calculated. For this equation, the problem of ﬁnite-ﬁeld reductions in the form of Darboux integrable systems of equations of a hyperbolic type is discussed. For reductions of small orders, N = 1 and N = 2, complete sets of characteristic integrals are presented. Note that the existence of characteristic integrals makes it possible to construct particular solutions to the original lattice. For the case N = 1, an explicit solution was found in this paper. A new semi-discrete equation is found that lies beyond the considered class. For this equation, the Lax pair is presented.


Introduction
Nonlinear integrable differential-difference equations with one continuous and two discrete variables are intermediate objects between the class of integrable nonlinear equations of dimension 1 + 2, such as the Toda lattice, the Kadomtsev-Petviashvili equation, the Davey-Stewartson equation, etc. (see [1][2][3]), and the class of fundamental discrete Hirota-Miwa-type equations (see [4][5][6][7]). They also have applications in the theory of Laplace transformations in the class of linear differential-difference equations of a hyperbolic type [8]. More precisely, the sequence of Laplace invariants satisfies some integrable differential-difference equations of the form (1), related by a point transformation to Equation (L7) from the list below. Such lattices can find application in studying the nonlinear dynamics of localized magnetic inhomogeneities in such magnetically ordered substances as ferromagnets, antiferromagnets with weak ferromagnetism, and magnetoelastic and magnetoelectric interactions [9,10]. Equations of the type under consideration have potential applications in the problem of describing dislocations in media with a microstructure, as well as in the nonlinear theory of elastic and nonelastic deformations accompanied with the deep reconstruction of an initially ideal lattice: the switching of interatomic bonds, the changing of the class of symmetry, the formation of new phases, singular defects and heterogeneities, and the fragmentation of the lattice (see, for instance, [11,12]). Lattices of the type under consideration can be used in the derivation of difference schemes for finding numerical solutions to nonlinear partial differential equations (see [13,14]).
Consider the class of integrable differential-difference equations with one continuous and two discrete independent variables of the following form: where the unknown function u = u j n (x) depends on the continuous x and the discrete n and j. The subscript x denotes the derivative u x = d dx u. Below, we give a list of known integrable models of the form (1), without pretending to be complete. The problem of the complete classification of equations of this class remains open (see, for example, [15,16]).
Equation (L1) coincides up to a simple change in the variablesū(j, n) = e u j n with Equation (N-7) in [7]. Equation (L8) in a slightly different form, where ε = a−b ab , was also obtained in [7] (see (N-2)). Indeed, it is easy to check that it goes into (L8) after a point change in variables of the form v n,m = u m+n−1 −n+1 , Equation (L7) was found, to the best of our knowledge, in [8] as an equation for the Laplace invariants of the linear differential-difference equations of the hyperbolic type with variable coefficients. Equations (L2)-(L6) were obtained in [17] by an original method based on a preliminary study of the dispersionless limits of lattices for integrability, followed by a return transition to the class of lattices (1) as a result of dispersion deformations. Equation (L9) was found in [18]. A distinctive feature of nonlinear lattices with three independent variables is that they allow reductions in the form of systems of differential-difference equations with complete sets of characteristic integrals. Such systems of equations are called Darboux integrable. As a rule, it is possible to find in an explicit form general solutions to Darboux integrable systems at least for the case of small orders. This circumstance allows one to find localized particular solutions for the lattices in 3D via the Darboux integrable reductions. In our articles [15,[19][20][21], an algorithm for classifying integrable equations of the form (1) is developed, where the presence of a hierarchy of the Darboux integrable reductions is taken as a classification criterion.
One can check that the equations from the list above are not related to each other by point changes in variables. However, the equations could be transformed from one to the other using more complex irreversible Miura-type transformations, which would be useful in finding explicit solutions to the equations (see, for instance, [22][23][24]). Note that the problem of the integrable classification of three-dimensional lattices is rather complicated and has not been solved so far, so the problem of finding new integrable equations by alternative methods, for example, using the Miura transformation, is topical. Therefore, the task of searching for such transformations is in demand. The study of this problem is one of the aims of the present work.
Let us briefly discuss the content of the article. In Section 2, we suggest an algorithm for searching differential or/and difference substitutions generated by a suitable representation of an equation that generalizes a local conservation law. In Section 3, we apply the algorithm to Equations (L1)-(L9). We establish relations between the equations. In Section 4, we give the Lax pair for lattice (L9) and calculate the continuum limit of two types for it. We also consider examples of the small-order Darboux integrable reductions admitted by the lattice. For the reduced systems, the complete sets of the characteristic integrals are constructed. An explicit particular solution to the lattice (L9) is presented, which is obtained via the reduced equation. A new example of an integrable equation with one continuous and two discrete variables (see (51)), which does not belong to the class of lattices (1), is found. The Lax pair for this model is presented.

Local Conservation Laws and Miura-Type Transformations
To find the Miura transformation that converts an equation of a given class (1) into some other equation of the same class, we use the following scheme. Let us first represent the equation as a local conservation law of the form below: where A and B are functions depending on the unknown function u = u j n (x), as well as on its derivatives with respect to x and shifts with respect to discrete variables j and n. Here, D n and D j denote operators for shifting a function argument in n and j, for example, D n y(n) = y(n + 1), D j y(j) = y(j + 1), and D x stands for the operator of the total differentiation with respect to x. If the representation (2) takes place, we define a new function v = v j n (x) by introducing a potential, i.e., assuming Then, substituting instead of the function A its representation (4) into relation (2), we arrive at the equality from which, applying the operator inverse to (D j − 1), we find another important relation connecting the variables u j n and v j n : Since we do concentrate only on autonomous transformations, we assume that C is a constant.
In the examples considered below, it is usually possible to completely exclude the variable u j n from the pair of relations (4) and (6). In such a case, we arrive at an equation for the new function v j n . The second way to look for the desired substitution uses representation (3). Here, we introduce a new function ω j n due to the relation below: Now, the equation yields a relation Integrating the latter, we find, in addition to (7), the following constraint that the variable ω j n should satisfy, To obtain the lattice equation for ω j n , it is necessary to get rid of the variable u j n by means of the relations (7) and (9).
As one can see from the examples below, this trick very often leads to success. In particular, it is in this way that we managed to find Miura-type transformations for each of the equations from the list under consideration.

Remark 1.
Note that some equations do not admit either of the representations (2) or (3) of the form of a conservation law. In this case, one can use a relation of a more general form (see Section 3.4).

Searching for the Transformations
In this section, we apply the algorithm discussed above to the lattices (L1)-(L9). It is interesting to note that the set of equations falls into two nonoverlapping subclasses. Most of the equations belong to the first subclass, consisting of different forms of the discrete version of 2D Toda lattice, namely (L1), (L3), (L4), (L7), (L8), and (L9). Only three equations fall into the second class: (L2), (L5), and (L6) (these equations were found in [17] and remain less studied). Equations from these two subclasses might be related to Miura-type transformations that are not of the form discussed in Section 2.

Relation between Equations (L1) and (L9)
Obviously, Equation (L1) can be rewritten as follows: According to the rule above, we introduce a new variable v j n by setting then Formulas (10) and (11) imply Let us differentiate (12) with respect to x and replace derivatives u j−1 n+1,x , u j n,x due to (11). As a result, we arrive at the equation that is reduced to (L9) by a point transformation

Relation between Equations (L2) and (L5)
Now, we derive a substitution generated by Equation (L5). To this end, we first simplify (L5) by the change in variables u Then, we represent the obtained equation in the form of the local conservation law: Let us introduce a new dependent variable ω j n as or Then, from (14) and (15) we obtain the equation which implies the following: One can easily verify that the variable v j n is completely eliminated from Formulas (15) and (16) and their consequences in a rather simple way. Indeed, by differentiating (16) we obtain The right-hand side of the formula is greatly simplified due to equalities (15) and (16), and, as a result, the equation reduces to the form Up to the linear transformation ω

Relation between Equations (L3) and (L1)
Let us turn to Equation (L3). Taking the logarithm of both sides of the equation, we obtain an identity which implies a local conservation law of the form Having introduced a new variable h j n by setting we deduce from (20) the equality which gives rise to where C is assumed to be constant. Now, we differentiate (23) with respect to x and find which is easily simplified due to (21). As a result, we obtain the lattice The latter coincides with the discrete Toda lattice (L1) up to the linear change in the independent discrete variables j = −j, n =j +n.

Relation between Equations (L3) and (L7)
We note that for Equation (L3) one can derive one more substitution. To this end, we use Equation (19), rewritten as Note that the relation (24) does not define any local conservation law. Assuming that ln u j n = (D n − 1)v j n , we obtain Then, combining (24) and (25), we arrive at and by differentiating (25) we find a similar relation Comparing the last two representations for u j n,x , one can obtain the equation which is nothing but Equation (L7).

Relation between Equations (L4) and (L1)
Now, we study Equation (L4) in the list. It can easily be represented as a local concervation law of the form We introduce a new unknown s j n by taking Then, from (29) it follows that By applying to the latter equation the operator D j − 1 −1 , we obtain the second part of the searched substitution: Let us simplify (31) a little, after which, applying the operator D x to both its sides, we obtain a nonlinear lattice,

Relation between Equations (L6) and (L2)
Let us look for the Miura-type transformations generated by Equation (L6). By taking the logarithm of both sides of (L6), we obtain a relation, Afterwards, by using the decomposition D n − D j = D n − 1 − D j − 1 , one can rewrite the relation as a local conservation law of the form below: Then, we introduce a new variable h j n by the rule Applying the reasoning above to (32) and (33), one can derive: Let us rewrite Miura transformations (33) and (34) in a more convenient form Differentiating (34), we find a relation We simplify some of the terms in (37) due to lattice (L6) and obtain the following: Finally, we obtain a lattice,

Relation between Equations (L8) and (L1)
Let us look for a substitution that would connect lattices (L8) and (L1). At first, we rewrite Equation (L8) in the form of a local conservation law: Supposing that u j n, we introduce a new dependent variable p j n . Then, obviously, it follows from (38) and (39) that From the latter, one can derive the second part of the Miura-type transformation: (D n − 1)p j n,x = e u j−1 n+1 − e u j n + C 1 .
The expression u j n = D j − 1 p j n + C 2 is an evident consequence of (39); therefore, we have e u j n = e C 2 e p j+1 n −p j n .
Let us eliminate the variables u j n , u j−1 n+1 from (41) due to (42) and obtain a lattice of the form below: Changing the variables by setting p j n = q j n + C 1 xn, we arrive at the lattice which converts to the discrete Toda lattice (L1) under replacement j = −j, n =j +n, x = e −C 2x .

Investigation of Lattice (L9)
In this section, we focus on a detailed study of lattice (L9): we calculate the continuum limit, study Darboux integrable reductions, construct explicit solutions, and show that it defines a Bäcklund transformation for an equation of the Toda lattice type. For reduced systems of small dimensions, characteristic integrals are found in both directions n and x. We also briefly discuss the method of constructing particular solutions of lattice (L9) by means of its reduction.
Let us begin with the Lax pair. It is checked straightforwardly that the system of linear equations [18] is compatible if and only if function v = v j n (x) is a solution to lattice (L9).

Computation of the Continuum Limit for Lattice (L9)
In this section, we will calculate the continuum limit of lattice (L9) under the unbounded refinement of the difference grid in both directions n and j. At first, we perform a scaling transformation y = δ 2 x in lattice (51) and represent it as follows: Let us assume that the unknown function v j n (x) is of the form v j n (x) = u(t, z, y), where t = δ 2 n, z = δj, and δ is a small positive parameter. Now, we estimate the difference derivatives with respect to the variables n and j appearing in (44) by virtue of the Taylor formula: Now, we substitute the found asymptotic representations of finite differences into Equation (44) and pass to the limit at δ → 0. As a result, we obtain the following dispersionless partial differential equation: which is integrable as well [25].

Intermediate Continuum Limit and Bäcklund Transformation-A New Integrable Lattice
Now, we study the continuum limit of lattice (L9), assuming that unknown function u j n (x) depends on a small parameter ε > 0 as follows u j n (x) = ω j (x, t) with t = εn. Under the condition ε → 0, lattice (L9) tends to an equation of the Toda lattice type [23] with a Lax pair of the form which is easily obtained from (43) under the relation ψ j n (x) = ϕ j (x, t), when ε → 0. A remarkable fact is that lattices (L9) and (47) have a close connection. Indeed, by comparing the Lax pairs one can conclude that lattice (L9) realizes a sequence of the Bäcklund transformations for lattice (47) (to the best of our knowledge, this observation is new). More precisely, any solution ω j of (47) is mapped by the Bäcklund transformation to another solutionω j of the same lattice.
We note that eigenfunctions ϕ j andφ j of the system (48) and of the system are connected with one another as follows: Iterations of the Bäcklund transformation (49) determine two nonlinear equations of the form The first lattice coincides with (L9) up to substitution x → −x; the second one does not belong to the class (1). It is also integrable, since it admits the Lax pair Note that integrable lattices of this type have not been previously considered in the literature.

Darboux Integrable Reductions of Lattice (L9)
Lattice (L9) admits the following cutting-off boundary condition preserving the integrability property. Imposing this kind of the condition on two points j = −1 and j = N, we obtain a finite field reduction of lattice (L9): admitting the Lax representation of the form deduced from (43) under the constraints ψ −1 n = ψ N+1 n = 0, v −1 n = C 0 , v N n = C 1 . Let us concentrate on the finite field system (53). It is worth noting that (53) admits complete sets of integrals in both directions n and x. Recall that a function I, depending on the dynamical variables, is an n-integral of the differential-difference system if it satisfies the relation D n I = I by virtue which admits n-integral and x-integral Example 2. For the choice N = 2, we arrive at the system for u n = v 0 n , v n = v 1 The obtained system admits a complete set of n-integrals, i.e., a set of two independent n-integrals and a complete set of x-integrals where τ 0 n = u n+1 − u n , τ 1 n = v n+1 − v n (about the complete set of integrals, see, for instance, [26,27]).
It can be proved that system (53) admits complete sets of x-integrals and n-integrals for any choice of N. For N ≥ 3, x-integrals have extremely large expressions. They are found by the method of characteristic Lie-Rinehart algebras (see [21]). As can be seen from the above examples, the n-integrals of the reduced system are much simpler. Note that the complete set of n-integrals to the system (53) for an arbitrary natural N ≥ 1 can be efficiently derived using the algorithm proposed in Lemma 2 of [21]. An alternative approach to obtain integrals is proposed in [28].

Explicit Solutions to Lattice (L9)
The existence of complete sets of integrals in both characteristic directions x and n makes it possible to completely separate the variables. Indeed, according to the definition, each solution v = (v 0 , v 1 , . . . , v N−1 ) of system (53) simultaneously satisfies two systems of equations, namely the system of ordinary differential equations and the system of ordinary discrete equations where ϕ 1 (x), ϕ 2 (x), . . . , ϕ N (x) are arbitrary functions of x and similarly ψ 1 (n), ψ 2 (n), . . . , ψ N (n) are arbitrary functions of the variable n.
In this case, the common solution v j n (x) of the systems provides a solution to the system (53). We remark that the solution of the reduced system (53) can obviously be prolonged to a particular solution of the original lattice (L9). However, for arbitrary N, system (60) is rather difficult to solve.
Here, we discuss a comparatively simple case, N = 1. Let us derive a formula describing the general solution to lattice (56). Due to the formula D n I = I, one can conclude that I = ϕ(x), i.e., function v n (x) is a solution to the equation Since ϕ(x) is arbitrary, we can assume that ϕ(x) is represented as for function z = z(x), which is now considered a new functional parameter instead of ϕ(x). As a result, we obtain Now, we have Equation (63), where z = z(x) is an arbitrary given function of x and v n = v n (x) is an unknown function. Let us make a change in the variables in Equation (63) by setting v n (x) = z(x) + s n (x). Then, we obtain which has a form of the Bernoulli equation. Let us replace s n (x) = 1 g n (x) and arrive at a linear inhomogeneous equation: In order to solve the corresponding homogeneous equation we again change the functional parameter. We pass from the functional parameter z = z(x) to a new parameter p(x) in such a way that Then, the homogeneous equation g n,x = −g n p x is easily integrated into g n (x) = Ce −p(x) . Now, supposing C = C(x), we return to Equation (65) and find C x = −e p(x) . We change the parameter in the equation again by setting e p(x) = ω x , where ω = ω(x) is a new parameter. As a result, we obtain C(x) = −ω(x) + K(n), where K(n) is a constant of integration that might depend on n. Therefore, we have To write down the final form of v n (x), we need an expression of z(x) in terms of ω(x). To this end, we exclude parameter p(x) from the relations Now, we substitute the found expressions into v n (x) = s n (x) + z(x) and obtain a general solution to lattice (56): v n (x) = 1 2 where ω = ω(x) is an arbitrary function of x and K(n) is an arbitrary function of n. By direct calculation, one can verify that function (66) is a solution to lattice (56) for any functions ω(x) and K(n). In a more general situation, to clarify the form of the solution, it is also necessary to use integrals in another direction.
Due to the fact that the termination conditions v −1 n = C 0 , v 1 n = C 1 are fully consistent with the dynamics, by virtue of the original three-dimensional lattice (L9), the solution of the reduced system (56) extends in an obvious way to the solution of lattice (L9): A problem of constructing explicit solutions via characteristic integrals is discussed in [29][30][31]. The papers convince that this approach is effective but needs further development, since the existing algorithms are technically complex.

Conclusions
In this article, integrable lattices with three independent variables, x, n, and j, are considered. A list of currently known examples is given. Differential and difference substitutions linking equations of the list with each other are presented. One of the lattices is studied in more detail. For that one, the continuum limits are computed. The degenerate cutoff condition is pointed out, which reduces the lattice to a Darboux integrable finite system of hyperbolic-type differential-difference equations when it is imposed on two different points, (x, n, j = 0) and (x, n, j = N). For the reductions of small orders, characteristic integrals are found. A method for constructing explicit particular solutions to the lattice via reductions is proposed.
A new integrable lattice is found that is not of the form (1). For this model, the Lax pair is presented.