Next Article in Journal
Effective Improved NSGA-II Algorithm for Multi-Objective Integrated Process Planning and Scheduling
Next Article in Special Issue
On Global Solutions of Hyperbolic Equations with Positive Coefficients at Nonlocal Potentials
Previous Article in Journal
Optimization Models for the Vehicle Routing Problem under Disruptions
Previous Article in Special Issue
Regularization and Inverse Spectral Problems for Differential Operators with Distribution Coefficients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Miura-Type Transformations for Integrable Lattices in 3D

by
Ismagil T. Habibullin
*,
Aigul R. Khakimova
and
Alfya U. Sakieva
Ufa Federal Research Center, Department of Mathematical Physics, Institute of Mathematics, Russian Academy of Sciences, 450008 Ufa, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(16), 3522; https://doi.org/10.3390/math11163522
Submission received: 28 July 2023 / Revised: 11 August 2023 / Accepted: 12 August 2023 / Published: 15 August 2023

Abstract

:
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 n + 1 , x j = f ( u n , x j , u n j + 1 , u n j , u n + 1 j , u n + 1 j 1 ) up to point transformations. An efficient 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 efficiency 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 finite-field 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.

1. 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:
u n + 1 , x j = f ( u n , x j , u n j + 1 , u n j , u n + 1 j , u n + 1 j 1 ) ,
where the unknown function u = u n j ( x ) depends on the continuous x and the discrete n and j. The subscript x denotes the derivative u x = d d x 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]).
(L1)
u n + 1 , x j = u n , x j + e u n + 1 j u n j + 1 e u n + 1 j 1 u n j ,
(L2)
u n + 1 , x j = u n , x j + e u n + 1 j 1 u n + 1 j e u n + 1 j 1 u n j e u n j u n j + 1 + e u n + 1 j u n j + 1 ,
(L3)
u n + 1 , x j = u n , x j ( u n + 1 j ) 2 u n + 1 j 1 u n j + 1 ,
(L4)
u n + 1 , x j = u n , x j u n + 1 j u n j + 1 u n + 1 j 1 u n j ,
(L5)
u n + 1 , x j = u n , x j u n + 1 j ( u n + 1 j u n j + 1 ) u n j + 1 ( u n + 1 j 1 u n j ) ,
(L6)
u n + 1 , x j = u n , x j ( u n + 1 j 1 u n + 1 j ) ( u n + 1 j u n j + 1 ) ( u n + 1 j 1 u n j ) ( u n j u n j + 1 ) ,
(L7)
u n + 1 , x j = u n , x j + e u n j + 1 + u n j + u n + 1 j u n + 1 j 1 ,
(L8)
u n + 1 , x j = u n , x j e u n j + 1 + e u n j + e u n + 1 j e u n + 1 j 1 ,
(L9)
u n + 1 , x j = u n , x j + ( u n + 1 j u n j ) ( u n j + 1 u n + 1 j u n j + u n + 1 j 1 ) .
Equation (L1) coincides up to a simple change in the variables u ¯ ( j , n ) = e u n j with Equation (N-7) in [7]. Equation (L8) in a slightly different form,
x 1 v n , m + 1 v n + 1 , m = ε e v n + 1 , m + 1 + e v n , m e v n , m + 1 e v n + 1 , m ,
where ε = a b a b , 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 n + 1 m + n 1 , x 1 = x ε .
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.

2. 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:
( D n 1 ) A = ( D j 1 ) B
or the form
D x A = ( D j 1 ) B ,
where A and B are functions depending on the unknown function u = u n j ( 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 n j ( x ) by introducing a potential, i.e., assuming
( D j 1 ) v n j = A .
Then, substituting instead of the function A its representation (4) into relation (2), we arrive at the equality
( D j 1 ) ( D n 1 ) v n j = ( D j 1 ) B ,
from which, applying the operator inverse to ( D j 1 ) , we find another important relation connecting the variables u n j and v n j :
( D n 1 ) v n j = B + C .
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 n j from the pair of relations (4) and (6). In such a case, we arrive at an equation for the new function v n j .
The second way to look for the desired substitution uses representation (3). Here, we introduce a new function ω n j due to the relation below:
D x ω n j = B .
Now, the equation yields a relation
D x A = D x D j 1 ω n j .
Integrating the latter, we find, in addition to (7), the following constraint that the variable ω n j should satisfy,
D j 1 ω n j = A + c o n s t .
To obtain the lattice equation for ω n j , it is necessary to get rid of the variable u n j 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).

3. 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.

3.1. Relation between Equations (L1) and (L9)

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

3.2. 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 n j = e v n j :
v n + 1 , x j = v n , x j e v n + 1 j v n j + 1 1 e v n + 1 j 1 v n j 1 .
Then, we represent the obtained equation in the form of the local conservation law:
( D n 1 ) ln v n , x j = ( D j 1 ) ln e v n + 1 j 1 v n j 1 .
Let us introduce a new dependent variable ω n j as
D j 1 ω n j = ln v n , x j
or
v n , x j = e ω n j + 1 ω n j .
Then, from (14) and (15) we obtain the equation
D j 1 D n 1 ω n j = D j 1 ln e v n + 1 j 1 v n j 1 ,
which implies the following:
D n 1 ω n j = ln e v n + 1 j 1 v n j 1 + C .
One can easily verify that the variable v n j is completely eliminated from Formulas (15) and (16) and their consequences in a rather simple way. Indeed, by differentiating (16) we obtain
ω n + 1 , x j ω n , x j = e v n + 1 j 1 v n j v n + 1 , x j 1 v n , x j e v n + 1 j 1 v n j 1 .
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
ω n + 1 , x j ω n , x j = 1 + e ω n j ω n + 1 j + C e ω n + 1 j ω n + 1 j 1 e ω n j + 1 ω n j .
Up to the linear transformation ω n j = u n j + C n + π i j , i 2 = 1 , Equation (18) coincides with Equation (L2).

3.3. 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
ln u n + 1 , x j ln u n , x j = 2 ln u n + 1 j ln u n + 1 j 1 ln u n j + 1 ,
which implies a local conservation law of the form
D n 1 ln u n , x j u n j = D j 1 ln u n + 1 j 1 u n j .
Having introduced a new variable h n j by setting
D j 1 h n j = ln u n , x j u n j ,
we deduce from (20) the equality
D j 1 D n 1 h n j = D j 1 ln u n + 1 j 1 u n j ,
which gives rise to
D n 1 h n j = ln u n + 1 j 1 u n j + C ,
where C is assumed to be constant. Now, we differentiate (23) with respect to x and find
h n + 1 , x j h n , x j = u n + 1 , x j 1 u n + 1 j 1 u n , x j u n j ,
which is easily simplified due to (21). As a result, we obtain the lattice
h n + 1 , x j h n , x j = e h n + 1 j h n + 1 j 1 e h n j + 1 h n j .
The latter coincides with the discrete Toda lattice (L1) up to the linear change in the independent discrete variables j = j ¯ , n = j ¯ + n ¯ .

3.4. 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
D n 1 ln u n , x j = 2 D n D n D j 1 D j ln u n j .
Note that the relation (24) does not define any local conservation law.
Assuming that ln u n j = D n 1 v n j , we obtain
u n j = e v n + 1 j v n j .
Then, combining (24) and (25), we arrive at
u n , x j = e 2 v n + 1 j v n + 1 j 1 v n j + 1
and by differentiating (25) we find a similar relation
u n , x j = v n + 1 , x j v n , x j e v n + 1 j v n j .
Comparing the last two representations for u n , x j , one can obtain the equation
v n + 1 , x j v n , x j = e v n j + 1 + v n j + v n + 1 j v n + 1 j 1 ,
which is nothing but Equation (L7).

3.5. 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
D n 1 ln u n , x j = D j 1 ln u n + 1 j 1 u n j .
We introduce a new unknown s n j by taking
D j 1 ln s n j = ln u n , x j .
Then, from (29) it follows that
D j 1 D n 1 ln s n j = D j 1 ln u n + 1 j 1 u n j .
By applying to the latter equation the operator D j 1 1 , we obtain the second part of the searched substitution:
D n 1 ln s n j = ln u n + 1 j 1 u n j + ln C , C = c o n s t .
Let us simplify (31) a little, after which, applying the operator D x to both its sides, we obtain a nonlinear lattice,
s n + 1 , x j s n j s n + 1 j s n , x j s n j 2 = C s n + 1 j s n + 1 j 1 s n j + 1 s n j ,
which is reduced to the discrete Toda lattice (L1) by the point transformation s n j = e u n j , x C x .

3.6. 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,
D n 1 ln u n , x j = D n D j ln u n j 1 u n j + D j 1 ln u n + 1 j 1 u n j .
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:
D n 1 ln u n , x j u n j 1 u n j = D j 1 ln u n + 1 j 1 u n j u n j 1 u n j .
Then, we introduce a new variable h n j by the rule
D j 1 ln h n j = ln u n , x j u n j 1 u n j .
Applying the reasoning above to (32) and (33), one can derive:
D n 1 ln h n j = ln u n + 1 j 1 u n j u n j 1 u n j + ln C .
Let us rewrite Miura transformations (33) and (34) in a more convenient form
u n , x j u n j 1 u n j = h n j + 1 h n j , C u n + 1 j 1 u n j u n j 1 u n j = h n + 1 j h n j ,
which implies that
1 C u n , x j u n + 1 j 1 u n j = h n j + 1 h n + 1 j .
Differentiating (34), we find a relation
h n + 1 , x j h n + 1 j h n , x j h n j = u n + 1 , x j 1 u n , x j u n + 1 j 1 u n j u n , x j 1 u n , x j u n j 1 u n j = = u n + 1 , x j 1 u n + 1 j 1 u n j C h n j + 1 h n + 1 j u n , x j 1 u n j 1 u n j + h n j + 1 h n j .
We simplify some of the terms in (37) due to lattice (L6) and obtain the following:
u n + 1 , x j 1 u n + 1 j 1 u n j u n , x j 1 u n j 1 u n j = u n , x j 1 u n j 1 u n j u n + 1 j 2 u n + 1 j 1 u n + 1 j 2 u n j 1 1 = = C h n j h n + 1 j 1 u n j 1 u n + 1 j 1 u n j 1 u n j = C h n j h n + 1 j 1 1 u n + 1 j 1 u n j u n j 1 u n j = C h n j h n + 1 j 1 1 1 C h n + 1 j h n j .
Finally, we obtain a lattice,
h n + 1 , x j h n + 1 j h n , x j h n j = h n j + 1 h n j C h n j + 1 h n + 1 j + C h n j h n + 1 j 1 h n + 1 j h n + 1 j 1 ,
which is related to lattice (L2) by a point transformation h n j = C n e p n j :
p n + 1 , x j p n , x j = e p n + 1 j 1 p n + 1 j e p n + 1 j 1 p n j e p n j p n j + 1 + e p n + 1 j p n j + 1 .

3.7. 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:
D n 1 u n , x j = D j 1 e u n + 1 j 1 e u n j .
Supposing that
u n , x j = D j 1 p n , x j ,
we introduce a new dependent variable p n j . Then, obviously, it follows from (38) and (39) that
D j 1 D n 1 p n , x j = D j 1 e u n + 1 j 1 e u n j .
From the latter, one can derive the second part of the Miura-type transformation:
D n 1 p n , x j = e u n + 1 j 1 e u n j + C 1 .
The expression u n j = D j 1 p n j + C 2 is an evident consequence of (39); therefore, we have
e u n j = e C 2 e p n j + 1 p n j .
Let us eliminate the variables u n j , u n + 1 j 1 from (41) due to (42) and obtain a lattice of the form below:
p n + 1 , x j p n , x j = e C 2 e p n + 1 j p n + 1 j 1 e p n j + 1 p n j + C 1 .
Changing the variables by setting p n j = q n j + C 1 x n , we arrive at the lattice
q n + 1 , x j q n , x j = e C 2 e q n + 1 j q n + 1 j 1 e q n j + 1 q n j ,
which converts to the discrete Toda lattice (L1) under replacement j = j ¯ , n = j ¯ + n ¯ , x = e C 2 x ¯ .

4. 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]
ψ n + 1 j = v n j v n + 1 j ψ n j + 1 + ψ n j , ψ n , x j = v n j v n j 1 ψ n j ψ n j 1
is compatible if and only if function v = v n j ( x ) is a solution to lattice (L9).

4.1. 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:
δ 2 d d y v n + 1 j v n j = = v n + 1 j v n j v n j + 1 2 v n j + v n j 1 + v n + 1 j 1 v n j 1 v n + 1 j v n j .
Let us assume that the unknown function v n j ( x ) is of the form v n j ( 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:
v n + 1 j v n j = u ( t + δ 2 , z , y ) u ( t , z , y ) = δ 2 u t ( t , z , y ) + O ( δ 4 ) , δ 0 ,
v n + 1 j 1 v n j 1 = u ( t + δ 2 , z δ , y ) u ( t , z δ , y ) = = u ( t , z , y ) + δ 2 u t ( t , z , y ) δ u z ( t , z , y ) + 1 2 δ 2 u z z ( t , z , y ) u ( t , z , y ) + δ u z ( t , z , y ) 1 2 δ 2 u z z ( t , z , y ) + + O ( δ 3 ) = = δ 2 u t ( t , z , y ) + O ( δ 3 ) , δ 0 ,
v n j + 1 2 v n j + v n j 1 = u ( t , z + δ , y ) 2 u ( t , z , y ) + u ( t , z δ , y ) = = u ( t , z , y ) + δ u z ( t , z , y ) + 1 2 δ 2 u z z ( t , z , y ) 2 u ( t , z , y ) + + u ( t , z , y ) δ u z ( t , z , y ) + 1 2 δ 2 u z z ( z , t , y ) + = = δ 2 u z z ( t , z , y ) + O ( δ 3 ) , δ 0 .
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:
u t y = u t u z z ,
which is integrable as well [25].

4.2. 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 n j ( x ) depends on a small parameter ε > 0 as follows u n j ( x ) = ω j ( x , t ) with t = ε n . Under the condition ε 0 , lattice (L9) tends to an equation of the Toda lattice type [23]
ω x t j + ω t j ω j + 1 2 ω j + ω j 1 = 0
with a Lax pair of the form
φ t j = ω t j φ j + 1 + φ j , φ x j = ω j ω j 1 φ j φ j 1
which is easily obtained from (43) under the relation ψ n j ( 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
ω ¯ x j ω x j = ω ¯ j ω j ω j + 1 + ω j + ω ¯ j ω ¯ j 1 , ω t j + 1 + ω ¯ t j 1 ω ¯ j ω j + ω ¯ t j ω j + 1 ω ¯ j + 1 + ω t j ω j 1 ω ¯ j 1 = 0
to another solution ω ¯ j of the same lattice.
We note that eigenfunctions φ j and φ ¯ j of the system (48) and of the system
φ ¯ t j = ω ¯ t j φ ¯ j + 1 + φ ¯ j , φ ¯ x j = ω ¯ j ω ¯ j 1 φ ¯ j φ ¯ j 1
are connected with one another as follows:
φ ¯ j = ω j ω ¯ j φ j + 1 + φ j .
Iterations of the Bäcklund transformation (49) determine two nonlinear equations of the form
( i ) w n + 1 , x j = w n , x j + ( w n + 1 j w n j ) ( w n j + 1 + w n + 1 j + w n j w n + 1 j 1 ) , ( i i ) ω n , t j + 1 + ω n + 1 , t j 1 ω n + 1 j ω n j + + ω n + 1 , t j ω n j + 1 ω n + 1 j + 1 + ω n , t j ω n j 1 ω n + 1 j 1 = 0 .
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
φ n + 1 j = ω n j ω n + 1 j φ j + 1 + φ n j , φ n , t j = ω n , t j φ n j + 1 + φ n j .
Note that integrable lattices of this type have not been previously considered in the literature.

4.3. Darboux Integrable Reductions of Lattice (L9)

Lattice (L9) admits the following cutting-off boundary condition
v n j ( x ) | j = j 0 = C , C = c o n s t
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):
v n 1 = C 0 , v n + 1 , x j v n , x j = v n + 1 j v n j v n j + 1 + v n j + v n + 1 j v n + 1 j 1 , 0 j N 1 , v n N = C 1 ,
admitting the Lax representation of the form
ψ n + 1 j = v n j v n + 1 j ψ n j + 1 + ψ n j , 0 j N 1 , ψ n + 1 N = ψ n N ,
ψ n , x 0 = v n 0 C 0 ψ n 0 , ψ n , x j = v n j v n j 1 ψ n j ψ n j 1 , 1 j N
deduced from (43) under the constraints ψ n 1 = ψ n N + 1 = 0 , v n 1 = 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 of the system. Similarly, a function J depending on dynamical variables is an x-integral if it solves the equation D x J = 0 by means of the system. Integrals of the form I = I ( x ) and J = J ( n ) depending on x and n only are called trivial. We are only interested in nontrivial integrals. Let us give examples of integrals.
Example 1.
By setting N = 1 in (53), we obtain a scalar equation for v : = v 0
v n + 1 , x v n , x = v n + 1 v n v n + 1 + v n C 0 C 1 ,
which admits n-integral
I = v x ( v C 0 ) 2 ( v C 1 ) ( C 0 C 1 )
and x-integral
J = v n v n + 2 v n + 1 v n + 3 v n v n + 3 v n + 1 v n + 2 .
Example 2.
For the choice N = 2 , we arrive at the system for u n = v n 0 , v n = v n 1
u n + 1 , x u n , x = u n + 1 u n v n + u n + u n + 1 C 0 , v n + 1 , x v n , x = v n + 1 v n C 1 + v n + v n + 1 u n + 1 .
The obtained system admits a complete set of n-integrals, i.e., a set of two independent n-integrals
I 1 = u n , x + v n , x u n C 0 2 u n C 1 C 0 v n C 1 v n 2 , I 2 = v n , x x 2 v n C 1 v n , x C 0 v n u n , x v n C 0 u n C 1 u n v n C 0 + C 1
and a complete set of x-integrals
J 1 = τ n 0 τ n 1 + τ n + 1 0 τ n + 1 1 + τ n 0 τ n + 1 1 ( τ n + 1 0 τ n + 1 1 + τ n + 2 0 τ n + 2 1 + τ n + 1 0 τ n + 2 1 ) τ n 0 τ n + 2 1 τ n + 1 0 τ n 1 + τ n + 2 0 τ n + 1 1 + τ n + 2 0 τ n 1 , J 2 = τ n + 3 0 τ n 1 τ n + 1 0 τ n + 1 1 + τ n + 2 0 τ n + 2 1 + τ n + 1 0 τ n + 2 1 τ n + 2 0 τ n + 1 1 τ n + 1 0 τ n 1 + τ n + 2 0 τ n + 1 1 + τ n + 2 0 τ n 1 + τ n + 3 0 τ n + 2 1 + τ n + 3 0 τ n + 1 1 + τ n + 3 0 τ n 1 ,
where τ n 0 = u n + 1 u n , τ n 1 = 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].

4.4. 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
I 1 = φ 1 ( x ) , I 2 = φ 2 ( x ) , , I N = φ N ( x )
and the system of ordinary discrete equations
J 1 = ψ 1 ( n ) , J 2 = ψ 2 ( n ) , , J N = ψ N ( n ) ,
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 n j ( x ) j = 0 N 1 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
v n , x v n C 0 2 ε v n C 1 = φ ( x ) , ε = C 0 C 1 .
Since φ ( x ) is arbitrary, we can assume that φ ( x ) is represented as
φ ( x ) = z x ( z C 0 ) 2 ε ( z C 1 )
for function z = z ( x ) , which is now considered a new functional parameter instead of φ ( x ) . As a result, we obtain
v n , x v n C 0 2 ε v n C 1 = z x ( z C 0 ) 2 ε ( z C 1 ) .
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
s n , x s n 2 s n ε + 2 z C 0 = 0 ,
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:
g n , x = g n ε + 2 z C 0 1 .
In order to solve the corresponding homogeneous equation
g n , x = g n ε + 2 z C 0
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
p x = ε + 2 ( z C 0 ) .
Then, the homogeneous equation
g n , x = g n p x
is easily integrated into g n ( x ) = C e 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
g n = K ( n ) ω ω x so   that s n ( x ) = ω x K ( n ) ω .
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
z ( x ) = C 0 + 1 2 ( p x ε ) , e p ( x ) = ω x
and find
z ( x ) = 1 2 ω x x ω x + C 0 + C 1 2 .
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 ω x x ω x + ω x K ( n ) ω + C 0 + C 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 n 1 = C 0 , v n 1 = 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):
v n j ( x ) = C 0 for j 1 , ω x x 2 ω x + ω x K ( n ) ω + C 0 C 1 2 for j = 0 , C 1 for j 1 .
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.

5. 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.

Author Contributions

Conceptualization, I.T.H. and A.R.K.; methodology, I.T.H.; software, A.U.S.; validation, I.T.H., A.R.K. and A.U.S.; formal analysis, A.U.S.; investigation, I.T.H., A.R.K. and A.U.S.; writing—original draft preparation, I.T.H.; writing—review and editing, A.R.K.; visualization, A.U.S.; supervision, I.T.H.; project administration, A.R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The work of A.R. Khakimova is supported in part by the Young Russian Mathematics award.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bogdanov, L.V.; Konopelchenko, B.G. Analytic-bilinear approach to integrable hierarchies. I. Generalized KP hierarchy. J. Math. Phys. 1998, 39, 4683–4700. [Google Scholar] [CrossRef] [Green Version]
  2. Bogdanov, L.V.; Konopelchenko, B.G. Analytic-bilinear approach to integrable hierarchies. II. Multicomponent KP and 2D Toda lattice hierarchies. J. Math. Phys. 1998, 39, 4701–4728. [Google Scholar] [CrossRef] [Green Version]
  3. Konopelchenko, B.G.; Schief, W.K. Menelaus’ theorem, Clifford configurations and inversive geometry of the Schwarzian KP hierarchy. J. Phys. A Math. Gen. 2002, 35, 6125–6144. [Google Scholar] [CrossRef] [Green Version]
  4. Hirota, R. Nonlinear partial difference equations, II. Discrete-time Toda equations. J. Phys. Soc. Jpn. 1977, 43, 2074–2078. [Google Scholar] [CrossRef]
  5. Hirota, R. Discrete analogue of a generalized Toda equation. J. Phys. Soc. Jpn. 1981, 50, 3785–3791. [Google Scholar] [CrossRef]
  6. Miwa, T. On Hirota’s difference equation. Proc. Jpn. Acad. Ser. A 1982, 58, 9–12. [Google Scholar] [CrossRef]
  7. Date, E.; Jimbo, M.; Miwa, T. Method for generating discrete soliton equation. II. J. Phys. Soc. Jpn. 1982, 51, 4125–4131. [Google Scholar] [CrossRef]
  8. Adler, V.E.; Startsev, S.Y. Discrete analogues of the Liouville equation. Theoret. Math. Phys. 1999, 121, 1484–1495. [Google Scholar] [CrossRef] [Green Version]
  9. Shamsutdinov, M.A.; Khabibullin, I.T.; Kharisov, A.T.; Tankeyev, A.P. Dynamics of magnetic kinks in exchange-coupled ferromagnetic layers. Phys. Met. Metallogr. 2009, 108, 327–340. [Google Scholar] [CrossRef]
  10. Shamsutdinov, M.A.; Lomakina, I.Y.; Nazarov, V.N.; Kharisov, A.T.; Shamsutdinov, D.M. Ferro and Antiferromagnetodynamics. Nonlinear Oscillations, Waves and Solitons; Nauka: Moscow, Russia, 2009. (In Russian) [Google Scholar]
  11. Kunin, I.A. Elastic Media with Microstructure I: One-Dimensional Models; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1982. [Google Scholar]
  12. Aero, E.L.; Bulygin, A.N. The nonlinear theory of localized waves in a complex crystalline lattice as a discrete continual model. Comput. Contin. Mech. 2008, 1, 14–30. [Google Scholar] [CrossRef] [Green Version]
  13. Yang, X.; Wu, L.; Zhang, H. A space-time spectral order sinc-collocation method for the fourth-order nonlocal heat model arising in viscoelasticity. Appl. Math. Comput. 2023, 457, 128192. [Google Scholar] [CrossRef]
  14. Jiang, X.; Wang, J.; Wang, W.; Zhang, H. A predictor-corrector compact difference scheme for a nonlinear fractional differential equation. Fractal Fract. 2023, 7, 521. [Google Scholar] [CrossRef]
  15. Kuznetsova, M.N.; Habibullin, I.T.; Khakimova, A.R. On the problem of classifying integrable chain with three independent variables. Theoret. Math. Phys. 2003, 215, 242–268. [Google Scholar] [CrossRef]
  16. Habibullin, I.T.; Khakimova, A.R. On the classification of nonlinear integrable three-dimensional chains by means of characteristic Lie algebras. arXiv 2023, arXiv:2306.14585. [Google Scholar]
  17. Ferapontov, E.V.; Novikov, V.S.; Roustemoglou, I. On the classification of discrete Hirota-type equations in 3D. Int. Math. Res. Not. 2015, 2015, 4933–4974. [Google Scholar] [CrossRef] [Green Version]
  18. Adler, V.E. The tangential map and associated integrable equations. J. Phys. A Math. Theor. 2009, 42, 332004. [Google Scholar] [CrossRef]
  19. Habibullin, I. Characteristic Lie rings, finitely-generated modules and integrability conditions for (2+1)-dimensional lattices. Phys. Scr. 2013, 87, 065005. [Google Scholar] [CrossRef]
  20. Habibullin, I.; Poptsova, M. Classification of a subclass of two-dimensional lattices via characteristic Lie rings. SIGMA 2017, 13, 073. [Google Scholar] [CrossRef]
  21. Habibullin, I.T.; Khakimova, A.R. Characteristic Lie algebras of integrable differential-difference equations in 3D. J. Phys. A Math. Theor. 2021, 54, 295202. [Google Scholar] [CrossRef]
  22. Nijhoff, F.W.; Capel, H.W.; Wiersma, G.L.; Quispel, G.R.W. Bäcklund transformations and three-dimensional lattice equations. Phys. Lett. A 1984, 105, 267–272. [Google Scholar] [CrossRef]
  23. Shabat, A.B.; Yamilov, R.I. To a transformation theory of two-dimensional integrable systems. Phys. Lett. A 1997, 227, 15–23. [Google Scholar] [CrossRef]
  24. Levi, D.; Petrera, M.; Scimiterna, C.; Yamilov, R. On Miura transformations and Volterra-type equations associated with the Adler-Bobenko-Suris equations. SIGMA 2008, 4, 077. [Google Scholar] [CrossRef] [Green Version]
  25. Ferapontov, E.V.; Habibullin, I.T.; Kuznetsova, M.N.; Novikov, V.S. On a class of 2D integrable lattice equations. J. Math. Phys. 2020, 61, 073505. [Google Scholar] [CrossRef]
  26. Zhiber, A.V.; Kuznetsova, M.N. Integrals and characteristic Lie rings of semi-discrete systems of equations. Ufa Math. J. 2021, 13, 22–32. [Google Scholar] [CrossRef]
  27. Habibullin, I.T.; Khakimova, A.R. Integrals and characteristic algebras for systems of discrete equations on a quadrilateral graph. Theoret. Math. Phys. 2022, 213, 1589–1612. [Google Scholar] [CrossRef]
  28. Smirnov, S.V. Darboux integrability of discrete two-dimensional Toda lattices. Theoret. Math. Phys. 2015, 182, 189–210. [Google Scholar] [CrossRef] [Green Version]
  29. Leznov, A.N.; Savel’ev, M.V. Group Methods for the Integration of Nonlinear Dynamical Systems; Nauka: Moscow, Russia, 1985. (In Russian) [Google Scholar]
  30. Kostrigina, O.S.; Zhiber, A.V. Darboux-integrable two-component nonlinear hyperbolic systems of equations. J. Math. Phys. 2011, 52, 033503. [Google Scholar] [CrossRef]
  31. Kuznetsova, M.N. Construction of localized particular solutions of chains with three independent variables. Theoret. Math. Phys. 2023, 216, 1158–1167. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Habibullin, I.T.; Khakimova, A.R.; Sakieva, A.U. Miura-Type Transformations for Integrable Lattices in 3D. Mathematics 2023, 11, 3522. https://doi.org/10.3390/math11163522

AMA Style

Habibullin IT, Khakimova AR, Sakieva AU. Miura-Type Transformations for Integrable Lattices in 3D. Mathematics. 2023; 11(16):3522. https://doi.org/10.3390/math11163522

Chicago/Turabian Style

Habibullin, Ismagil T., Aigul R. Khakimova, and Alfya U. Sakieva. 2023. "Miura-Type Transformations for Integrable Lattices in 3D" Mathematics 11, no. 16: 3522. https://doi.org/10.3390/math11163522

APA Style

Habibullin, I. T., Khakimova, A. R., & Sakieva, A. U. (2023). Miura-Type Transformations for Integrable Lattices in 3D. Mathematics, 11(16), 3522. https://doi.org/10.3390/math11163522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop