Next Article in Journal
Lump, Breather, Ma-Breather, Kuznetsov–Ma-Breather, Periodic Cross-Kink and Multi-Waves Soliton Solutions for Benney–Luke Equation
Next Article in Special Issue
Exact and Data-Driven Lump Wave Solutions for the (3+1)-Dimensional Hirota–Satsuma–Ito-like Equation
Previous Article in Journal
Combinatorial Identities Concerning Binomial Quotients
Previous Article in Special Issue
Some Novel Fusion and Fission Phenomena for an Extended (2+1)-Dimensional Shallow Water Wave Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Symmetries of Integrable Quadrilateral Equations

1
School of Information Science and Engineering, Shandong Agricultural University, Taian 271018, China
2
Department of Mathematics, Shanghai University, Shanghai 200444, China
3
Newtouch Center for Mathematics, Shanghai University, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(6), 744; https://doi.org/10.3390/sym16060744
Submission received: 15 May 2024 / Revised: 2 June 2024 / Accepted: 11 June 2024 / Published: 14 June 2024
(This article belongs to the Special Issue Symmetry in Integrable Systems: Topics and Advances)

Abstract

:
In the paper, we describe a method for deriving generalized symmetries for a generic discrete quadrilateral equation that allows a Lax pair. Its symmetry can be interpreted as a flow along the tangent direction of its solution evolving with a Lie group parameter t. Starting from the spectral problem of the quadrilateral equation and assuming the eigenfunction evolves with the parameter t, one can obtain a differential-difference equation hierarchy, of which the flows are proved to be commuting symmetries of the quadrilateral equation. We prove this result by using the zero-curvature representations of these flows. As an example, we apply this method to derive symmetries for the lattice potential Korteweg–de Vries equation.

1. Introduction

In the study of integrable partial difference equations (P Δ Es), the existence of an infinite hierarchy of symmetries remains a fundamental integrability criterion. The concept of symmetries of integrable P Δ Es appeared as early as in 1991 in [1] as similarity constraints for integrable lattices. There have been many advancements in the research of symmetries of P Δ Es [2,3,4,5,6,7,8,9,10,11,12,13,14]. Various methods have been developed to construct these symmetries, such as the method via the spectral problem [2,3,4,6], using the recursion operator [11], using master symmetry [13], the Gardner method [12] and the Lax–Darboux scheme [9,10].
Among them, the Lax–Darboux scheme is an algebraic structure that encompasses integrable partial differential equations, differential-difference equations and P Δ Es within a unified framework. In particular, within this framework, the Bianchi permutability condition for the Darboux transformations leads to a P Δ E, and the compatibility of the Lax pairs and the Darboux transformations lead to an infinite hierarchy of differential-difference equations which provide commuting symmetries of the P Δ E [9].
In the study of symmetries of continuous and differential-difference equations with Lax pairs, the zero-curvature representations of isospectral and non-isospectral flows have proven to be a powerful tool in deriving symmetries and their Lie algebras for the isospectral equation hierarchy as well as non-isospectral equation hierarchy, e.g., [15,16,17,18,19,20,21,22,23,24,25,26]. In this paper, we aim to extend this method to Lax-integrable P Δ Es, with a particular focus on isospectral flows. Specifically, we will begin by deriving the isospectral flows and their zero-curvature representations associated with either the m-part or n-part of the Lax pair of a given P Δ E. Subsequently, we will use these representations, along with the zero-curvature representation of the P Δ E, to prove that these isospectral flows are symmetries of the P Δ E. The proof follows a similar approach to that used in the cases of the continuous and differential-difference equations involving Lax pairs. The method can be viewed as a more general version of the Lax–Darboux scheme.
This paper is organized as follows. Section 2 contains the notations and relevant definitions essential for the subsequent discussions. Section 3 elaborates on the method for deriving symmetries for Lax-integrable P Δ Es by zero-curvature representations. As a specific application of this method, in Section 4, we derive symmetries for the lattice potential Korteweg–de Vries (lpKdV) equation. Finally, Section 5 is devoted to concluding remarks.

2. Notations and Definitions

In this paper, our focus is solely on the following type of quadrilateral equation:
Q ( u , u ˜ , u ^ , u ˜ ^ ; p , q ) = 0 ,
where u = u ( n , m ) is a function defined on Z × Z ,
u ˜ = E n u = u ( n + 1 , m ) , u ^ = E m u = u ( n , m + 1 ) , u ˜ ^ = u ( n + 1 , m + 1 ) ,
E n and E m , respectively, serve as shift operators in the n and m directions, while p and q are spacing parameters of direction n and m, respectively.
We introduce some notations that will be used in deriving isospectral hierarchy. Assume that u = u ( n , m ) vanishes rapidly as | n | or | m | . Let V be a linear space consisting of all functions f, where each f is a function of u and its shifts, and f u = 0 = 0 . Each f is C differentiable with respect to n and m, and C -Gâteaux differentiable with respect to u and its shifts. Here, the Gâteaux (or Fr e ´ chet) derivative of f V (or f as an operator living on V ) in direction g V is defined as
f [ g ] = d d ϵ f ( u + ϵ g ) | ϵ = 0 .
By means of the Gâteaux derivative, one can define the Lie product for any f , g V as
[ [ f , g ] ] = f [ g ] g [ f ] .
We define a Laurent matrix polynomials space Q 2 composed of all 2 × 2 matrices M = ( q i j ( r , u ) ) 2 × 2 , where all the { q i j } are Laurent polynomials of r. In addition, we denote Q 2 ( 0 ) = { M Q 2 M u = 0 = 0 } . We also introduce operators Δ n ( ± ) = E n ± 1 , whose inverse operators can be denoted as (e.g., [4])
( Δ n ( ) ) 1 = k = 0 E n k , ( Δ n ( + ) ) 1 = k = 0 ( 1 ) k E n k .
Lie symmetries of Equation (1) are obtained by requiring the infinitesimal invariant condition
pr X Q | Q = 0 = 0 .
Here, X is the infinitesimal generator; pr X is the prolongation of the infinitesimal generator X . They are defined as follows
X = ξ u ,
pr X = ξ u + ξ ˜ u ˜ + ξ ^ u ^ + ξ ˜ ^ u ˜ ^ ,
where ξ is a function of u and its shifts u ( n + i , m + j ) , i , j Z . If ξ exclusively depends on u, then the symmetries are called point symmetries; otherwise, they are referred to as generalized symmetries. In this paper, we only focus on discussing the generalized symmetries of Equation (1). In terms of the Gâteaux derivative, to meet the condition (4a) is to find ξ such that
Q [ ξ ] | Q = 0 = 0 .
Corresponding to the infinitesimal generator (4b), a group transformation between solutions of Equation (1), namely, g t : u ( n , m ) u ( n , m ; t ) , can be derived in principle by solving the initial value problem
d u ( t ) d t = ξ ( u ( t ) , u ˜ ( t ) , u ^ ( t ) , ) , with u ( 0 ) = u ,
where t R is the Lie group parameter, u = u ( n , m ) and u ( t ) u ( n , m ; t ) . If u is a solution of Equation (1) and ξ is its symmetry, then u ( t ) determined by the initial problem (6) provides a solution for Equation (1) as well. Thus, (6) can be understood as that u ( t ) evolves with the parameter t along the tangent direction ξ in the solution manifold of Equation (1) (see [27]). In this context, if u ( t ) is a solution of (1), i.e.,
Q ( u ( t ) , u ˜ ( t ) , u ^ ( t ) , u ˜ ^ ( t ) ; p , q ) = 0 ,
we always have
d d t Q ( u ( t ) ) = 0 ,
where we drop off u ˜ ( t ) , u ^ ( t ) and u ˜ ^ ( t ) for short. If we assume d d t and shift operators commute, i.e., d d t E n = E n d d t , d d t E m = E m d d t , then (8) gives rise to (4a) (i.e., (5)) provided the initial problem (6) holds.

3. General Description of the Method

In this section, we elaborate on the method for deriving symmetries of P Δ E (1).

3.1. Isospectral Flows { K ( l ) }

Suppose (1) has a Lax pair:
ϕ ˜ = M ( u , u ˜ ; p , r ) ϕ ,
ϕ ^ = N ( u , u ^ ; q , r ) ϕ ,
where M and N are matrices in Q 2 , ϕ = ( ϕ 1 , ϕ 2 ) T and r stands for a spectral parameter. In other words, Equation (1) functions as the compatibility condition for (9), i.e.,
M ^ N = N ˜ M ,
if and only if u solves (1).
To construct symmetries for Equation (1), we introduce a parameter t into u. In this context, the linear problem (9a) is written as
ϕ ˜ = M ( u ( t ) , u ˜ ( t ) ; p , r ) ϕ ,
and we assume that ϕ evolves with t as
ϕ t = U ϕ ,
where U Q 2 is a matrix function of E n j u ( t ) and p and r. The compatibility condition between (11a) and (11b) is
M t = U ˜ M M U .
Next, we present some conditions (cf. [28]), under which we can construct an isospectral hierarchy
u t l = K ( l ) , l = 0 , 1 , 2 , .
Note that u can be considered as u ( t 0 , t 1 , t 2 , ) , while in the following, we denote it by u for short and without causing any confusion.
Condition 1.
Assume that the linear problems (11a) and (11b) satisfy the following conditions:
(1) 
The matrix equation
M [ X ] = U ˜ M M U
possesses a unique couple of non-zero solutions X V and U Q 2 satisfying U u = 0 = U 0 , where U 0 is a given matrix independent of u and meeting commutating relation U 0 M | u = 0 = M | u = 0 U 0 ;
(2) 
for any given Y 0 V , there exist solutions X V and U Q 2 ( 0 ) satisfying
M [ X r α Y ] = U ˜ M M U ,
where α is some integer related to (11a).
In fact, in light of the first assumption in Condition 1, there exists a unique K ( 0 ) V and a unique U ( 0 ) Q 2 , such that
M [ K ( 0 ) ] = U ˜ ( 0 ) M M U ( 0 ) ,
where U ( 0 ) u = 0 = U 0 . Thus, the first equation is given by
u t 0 = K ( 0 ) .
The expression (16) is commonly referred to as the zero-curvature representation of Equation (17) or the flow K ( 0 ) . Once K ( 0 ) is determined, by the second assumption in Condition 1, there exists a unique field K ( 1 ) V and a unique matrix V ( 1 ) Q 2 ( 0 ) , such that
M [ K ( 1 ) r α K ( 0 ) ] = V ˜ ( 1 ) M M V ( 1 ) .
Repeating this procedure, one can determine a unique sequence of fields { K ( l ) } V and a unique sequence of matrices { V ( l ) } Q 2 ( 0 ) , for l = 2 , 3 , , such that
M [ K ( l ) r α K ( l 1 ) ] = V ˜ ( l ) M M V ( l ) , l = 2 , 3 , .
From Equations (16), (18) and (19), we can derive
M [ K ( l ) ] = U ˜ ( l ) M M U ( l ) , l = 0 , 1 , 2 , ,
where
U ( l ) = j = 0 l V ( j ) r α ( l j ) Q 2 , V ( 0 ) = U ( 0 ) , U ( l ) u = 0 = U 0 r α l .
It means that we obtain the following isospectral hierarchy
u t l = K ( l ) , l = 0 , 1 , 2 , ,
with (20) as their zero-curvature representations. The Lax pair of the above hierarchy now is composed of (11a) and
ϕ t l = U ( l ) ϕ .
Note that the map Y X determined by (15) defines a recursion operator L such that X = L Y . This means the above hierarchy can be expressed alternatively as
u t l = K ( l ) = L l K ( 0 ) , l = 0 , 1 , 2 , .
In addition, it follows from Condition 1 that the following proposition is valid.
Proposition 1.
If U 0 = 0 in Condition 1, then the matrix Equation (14) has only zero solutions X = 0 and U = 0 .
Utilizing the identity [29]
M [ [ [ f , g ] ] ] = ( M [ f ] ) g ( M [ g ] ) f , f , g V ,
from the zero-curvature representations (20), we can derive that the isospectral flows { K ( l ) } satisfy (cf. [22,25])
M [ [ [ K ( l 1 ) , K ( l 2 ) ] ] ] = ( E n U ( l 1 ) , U ( l 2 ) ) M M U ( l 1 ) , U ( l 2 ) ,
where
U ( l 1 ) , U ( l 2 ) = U ( l 1 ) [ K ( l 2 ) ] U ( l 2 ) [ K ( l 1 ) ] + [ U ( l 1 ) , U ( l 2 ) ] ,
and [ U ( l 1 ) , U ( l 2 ) ] = U ( l 1 ) U ( l 2 ) U ( l 2 ) U ( l 1 ) . Meanwhile, from (20b), one can obtain
U ( l 1 ) , U ( l 2 ) | u = 0 = 0 .
Consequently, considering (25) and (26), and applying Proposition 1, we can obtain the following.
Proposition 2.
The isospectral flows { K ( l ) } satisfy
[ [ K ( l 1 ) , K ( l 2 ) ] ] = 0 .

3.2. Symmetries

We have constructed the isospectral hierarchy (21) associated with (11a) and (11b). Next, we will prove that under certain conditions these isospectral flows { K ( l ) } provide commuting symmetries of Equation (1).
We assume u ( t ) is a solution of Equation (1), i.e., (7) holds, so the Lax pair of (7) is then composed of (11a) and
ϕ ^ = N ( u ( t ) , u ^ ( t ) ; q , r ) ϕ .
Note that by u ( t ) , we denote a more general u ( t 0 , t 1 , t 2 , ) . Now, we introduce a second condition.
Condition 2.
For the matrices N in (28) and U ( l ) in (22), we assume
( N t l + N U ( l ) ) N 1 | u = 0 = U ^ ( l ) | u = 0 , l = 0 , 1 , 2 , .
From the compatibility of (11a) and (28), we have (cf. (10))
M ^ = N ˜ M N 1 .
Taking derivative w.r.t. t l yields
( M ^ ) t l = N ˜ t l M N 1 + N ˜ M t l N 1 M ^ N t l N 1 .
For the terms on the right-hand side, by using (30) and (12), we obtain
N ˜ t l M N 1 = N ˜ t l N ˜ 1 N ˜ M N 1 = N ˜ t l N ˜ 1 M ^ , N ˜ M t l N 1 = N ˜ ( U ˜ ( l ) M M U ( l ) ) N 1 = N ˜ U ˜ ( l ) N ˜ 1 N ˜ M N 1 N ˜ M N 1 N U ( l ) N 1 = N ˜ U ˜ ( l ) N ˜ 1 M ^ M ^ N U ( l ) N 1 .
Thus, we have
( M ^ ) t l = M ^ [ ( u ^ ) t l ] = { E n [ ( N t l + N U ( l ) ) N 1 ] } M ^ M ^ [ ( N t l + N U ( l ) ) N 1 ] ,
where M ^ [ ( u ^ ) t l ] represents the Gâteaux derivative of M ^ in the direction ( u ^ ) t l with respect to u ^ . Here, we remark that we have admitted the commutative relations E n E m u = E m E n u and E n d d t l = d d t l E n , while the relation
( u ^ ) t l = ( u t l ) ^
remains open. We will see that once (33) holds, then the flows { K ( l ) } are symmetries of Equation (1). (For a continuous equation
u t = S ( u , u x , u x x , ) ,
the flow K ( l ) ( u ) defined by u t l = K ( l ) ( u ) is a symmetry of the above equation if ( u t ) t l = ( u t l ) t . For the discrete Equation (1), if we formally write it as
u ^ = P ( E n j u ) ,
then the commutative relation (33) can be viewed as a counterpart of ( u t ) t l = ( u t l ) t in the discrete case). Let us proceed. By recalling (20), we can have
M ^ [ K ^ ( l ) ] = U ^ ˜ ( l ) M ^ M ^ U ^ ( l ) , l = 0 , 1 , 2 , ,
where M ^ [ K ^ ( l ) ] represents the Gâteaux derivative of M ^ in direction K ^ ( l ) with respect to u ^ . Now, comparing (32) and (34), noticing the condition (29) and following Proposition 1, we immediately obtain
( u ^ ) t l = K ^ ( l ) , l = 0 , 1 , 2 ,
and
( N t l + N U ( l ) ) N 1 = U ^ ( l ) , l = 0 , 1 , 2 , .
Let us focus on the first result (35). It immediately indicates the commuting relation (33) since K ^ ( l ) = ( u t l ) ^ . Consequently,
Q [ K ( l ) ] | Q = 0 = 0 ,
i.e., K ( l ) is a symmetry of Equation (1).
We can summarize the above results in the following theorem.
Theorem 1.
Assuming that the PΔE (1) has Lax pair (9), matrix M satisfies Condition 1, and Condition 2 holds for N and U ( l ) , then the flows { K ( l ) } in (21) are commuting symmetries of the PΔE (1).

3.3. Isospectral Flows { H ( l ) }

One can also start from (28) and in this turn, assuming ϕ evolves with the parameter t, as (we will see that with certain assumptions, there is W = U (mod Equation (7)))
ϕ t = W ϕ ,
where W Q 2 . Then, a second isospectral hierarchy can be derived by assuming the following conditions.
Condition 3.
Assume that the linear problems (28) and (37) satisfy the following conditions:
(1) 
The matrix equation
N [ X ] = W ^ N N W
possesses a unique couple of non-zero solutions X V and W Q 2 satisfying W u = 0 = W 0 , where W 0 is a matrix independent of u and meeting W 0 N u = 0 = N u = 0 W 0 ;
(2) 
for any given Y 0 V , there exist solutions X V and W Q 2 ( 0 ) satisfying
N [ X r β Y ] = W ^ N N W ,
where β is some integer related to (28).
With the above assumption, a second isospectral hierarchy
u t l = H ( l ) , l = 0 , 1 , 2 , ,
with the zero-curvature representations
N [ H ( l ) ] = W ^ ( l ) N N W ( l ) , W ( l ) u = 0 = W 0 r β l , l = 0 , 1 , 2 , ,
can be derived. At the same time, we can obtain the following two propositions.
Proposition 3.
If W 0 = 0 in Condition 3, then the matrix Equation (38) has only zero solutions X = 0 and W = 0 .
Proposition 4.
The isospectral flows { H ( l ) } satisfy
[ [ H ( l 1 ) , H ( l 2 ) ] ] = 0 .
Now, we have two sets of flows, namely, { K ( l ) } and { H ( l ) } . Apparently, K ( l ) is a function of E n j u and H ( l ) is a function of E m j u . However, we can unify them in light of Condition 1, 2, 3 and a further assumption:
U ( l ) | u = 0 = W ( l ) | u = 0 , l = 0 , 1 , 2 , .
To show the uniformity, we recall the results (36) together with (33), which yields
N t l = N [ u t l ] = N [ K ( l ) ] = U ^ ( l ) N N U ( l ) , l = 0 , 1 , 2 , .
Comparing the above equation with (41), in light of Proposition 3 and the assumption (43), we immediately obtain (modulo Equation (1))
K ( l ) = H ( l ) , l = 0 , 1 , 2 , .
Now that the two sets are the same, we conclude that the flows { H ( l ) } in (40) are also commuting symmetries of the P Δ E (1).

4. Symmetries for the lpKdV Equation

In this section, we employ the method described in Section 3 to derive symmetries for the lpKdV Equation.
The lpKdV Equation is given by [30,31]
Q = ( p q + u ^ u ˜ ) ( p + q + u u ˜ ^ ) p 2 + q 2 = 0 ,
where p and q serves as the spacing parameters in n-direction and m-direction, respectively. It is known as the H1 equation in the Adler–Bobenko–Suris (ABS) classification [32]. It also appears as the nonlinear superposition formula (Bianchi identity) of the solutions of the (potential) KdV Equation [33]. The lpKdV Equation has the following Lax pair [34,35]:
ϕ ˜ = M ϕ = v v 2 p 2 r 1 v ϕ ,
ϕ ^ = N ϕ = v 1 v 1 2 q 2 r 1 v 1 ϕ ,
where ϕ = ( ϕ 1 , ϕ 2 ) T , r is the spectral parameter, and v = u ˜ u p , v 1 = u ^ u q . Note that the Lax pair (46) exhibits a symmetric property in the sense of switching ( p ,   ˜ ) and ( q ,   ^ ) .
To obtain a differential-difference hierarchy associated with the spectral problem (46a), we introduce the parameter t and assume ϕ evolves with t as
ϕ t = U ϕ = A B C D ϕ , U Q 2 ,
where U satisfies
U u = 0 = U 0 = 0 r 1 0
and r is indepandent of t. In the method we described in Section 3, the flows of the differential-difference hierarchy determined by (46a) and (47) will provide symmetries for the lpKdV Equation (45). The compatibility of (46a) and (47), i.e., ( E n ϕ ) t = E n ( ϕ t ) , yields the zero-curvature Equation (12). With the above matrices M and U, the zero-curvature Equation (12) gives rise to
v t = v ( E n 1 ) A + B ˜ ( r p 2 + v 2 ) C ,
2 v v t = ( r p 2 + v 2 ) A ˜ + v ( E n 1 ) B ( r p 2 + v 2 ) D ,
0 = v ( E n 1 ) C + D ˜ A ,
v t = ( r p 2 + v 2 ) C ˜ + v ( E n 1 ) D B .
From (49), we can derive
( r + p 2 ) ( E n 1 ) ( A + D ) = 0 ,
which implies that ( E n 1 ) ( A + D ) is zero, i.e.,
D = A + a ,
where a is a constant. However, considering the boundary condition (48), we need to take a = 0 and thus we have
D = A .
The subsequent analysis indicates that the linear problems (46a) and (47) satisfies Condition 1, and simultaneously, we derive the isospectral hierarchy associated with these linear problems. We begin by substituting the matrix M from (46a) and the matrix U from (47) into Equation (15), with α = 1 , which yields (noting that the Gâteaux derivative of M in (15) is taken with respect to u not v).
Δ n ( ) ( X r Y ) = r Δ n ( ) C + L 1 C ,
A = ( Δ n ( + ) ) 1 v Δ n ( ) C ,
B = r C + ( Δ n ( + ) ) 1 [ 2 v Δ n ( ) ( Δ n ( + ) ) 1 v Δ n ( ) + ( v 2 p 2 ) Δ n ( + ) ] C ,
where
L 1 = [ 2 ( Δ n ( + ) ) 1 1 ] v Δ n ( ) ( Δ n ( + ) ) 1 v Δ n ( ) + [ 1 ( Δ n ( + ) ) 1 ] ( v 2 p 2 ) Δ n ( + ) v 2 + p 2 .
To derive the initial flow K ( 0 ) , one should take Y = 0 in (53a). Then, by comparing the coefficients of the same powers of r in (53a), we can have
C = 1 ,
X = ( Δ n ( + ) ) 1 ( v 2 p 2 ) .
Substituting (54) into (53b) and (53c) leads to
A = 0 ,
B = r p 2 + 2 ( Δ n ( + ) ) 1 v 2 .
Hence, there exists a unique K ( 0 ) V and a unique V ( 0 ) Q 2 , given by
K ( 0 ) = ( Δ n ( + ) ) 1 ( v 2 p 2 ) ,
V ( 0 ) = 0 r p 2 + 2 ( Δ n ( + ) ) 1 v 2 1 0 ,
which satisfies the following equation
M [ K ( 0 ) ] = V ˜ ( 0 ) M M V ( 0 ) .
Next, to derive the flow K ( 1 ) , we take Y = K ( 0 ) in (53a). Then, by comparing the coefficients of the same powers of r in (53a), we obtain
C = K ( 0 ) ,
X = L K ( 0 ) ,
where
L = ( Δ n ( ) ) 1 L 1 .
Substituting (60) into (53b) and (53c), we obtain
A = ( Δ n ( + ) ) 1 v Δ n ( ) K ( 0 ) ,
B = r K ( 0 ) + ( Δ n ( + ) ) 1 [ 2 v Δ n ( ) ( Δ n ( + ) ) 1 v Δ n ( ) + ( v 2 p 2 ) Δ n ( + ) ] K ( 0 ) .
Hence, there exists a unique K ( 1 ) V and a unique V ( 1 ) Q 2 ( 0 ) , given by
K ( 1 ) = L K ( 0 ) ,
V ( 1 ) = A B C A ,
where A, B and C are given by (63a), (63b) and (60), respectively, which satisfy the following equation
M [ K ( 1 ) r K ( 0 ) ] = V ˜ ( 1 ) M M V ( 1 ) .
Repeating the above procedure from K ( 0 ) to K ( 1 ) allows us to obtain higher order flows K ( l ) = L l K ( 0 ) V and V ( l ) Q 2 ( 0 ) , l = 2 , 3 , , which satisfy
M [ K ( l ) r K ( l 1 ) ] = V ˜ ( l ) M M V ( l ) , l = 2 , 3 , .
From Equations (59), (66) and (67), we can derive
M [ K ( l ) ] = U ˜ ( l ) M M U ( l ) , l = 0 , 1 , 2 , ,
where
U ( l ) = j = 0 l V ( j ) r l j Q 2 , U ( l ) u = 0 = U 0 r l .
Up to this point, we have obtained the following isospectral hierarchy
u t l = K ( l ) = L l K ( 0 ) , l = 0 , 1 , 2 , ,
where K ( 0 ) and L are given by (58a) and (62), with Equation (68) as their zero-curvature representations. The first equation in this hierarchy is
u t 0 = K ( 0 ) = ( Δ n ( + ) ) 1 ( v 2 p 2 ) ,
i.e.,
u ˜ t 0 + u t 0 = ( u ˜ u ) 2 2 p ( u ˜ u ) .
There also exists an isospectral hierarchy associated with the spectral problem (46b). Due to the symmetric property, this isospectral hierarchy and the corresponding zero-curvature representations can be, respectively, obtained from (69) and (68) by replacing ( p ,   ˜ ) with ( q ,   ^ ) . Their expressions are
u t l = H ( l ) = L ¯ H ( 0 ) , l = 0 , 1 , 2 , ,
N [ H ( l ) ] = W ^ ( l ) N N W ( l ) , l = 0 , 1 , 2 , ,
where W ( l ) can be obtained from U ( l ) in (68b) by replacing ( p ,   ˜ ) with ( q ,   ^ ) , and we have
W ( l ) | u = 0 = U 0 r l ,
H ( 0 ) = ( Δ m ( + ) ) 1 ( v 1 2 q 2 ) ,
L ¯ = ( Δ m ( ) ) 1 { [ 2 ( Δ m ( + ) ) 1 1 ] v 1 Δ m ( ) ( Δ m ( + ) ) 1 v 1 Δ m ( ) + [ 1 ( Δ m ( + ) ) 1 ] ( v 1 2 q 2 ) Δ m ( + ) v 1 2 + q 2 } .
Through straightforward calculations, we can obtain
( N t l + N U ( l ) ) N 1 | u = 0 = q r 1 q U 0 r l q r 1 q 1 = U 0 r l = U ^ ( l ) | u = 0 ,
U ( l ) | u = 0 = W ( l ) | u = 0 = U 0 r l .
Then, by Theorem 1 and the result in Section 3.3, we conclude the following:
  • Both flows { K ( l ) } and { H ( l ) } are symmetries of the lpKdV Equation (45).
  • The symmetries commute:
    [ [ K ( l 1 ) , K ( l 2 ) ] ] = 0 , [ [ H ( l 1 ) , H ( l 2 ) ] ] = 0 , l 1 , l 2 = 0 , 1 , 2 , .
  • { K ( l ) } and { H ( l ) } are same up to modulo the lpKdV Equation (45).

5. Concluding Remarks

In this paper, we have described an approach to deriving symmetries for Lax-integrable P Δ Es. The primary technique employed in this paper involves the zero-curvature representations of flows, which have been used in deriving symmetries for continuous and differential-difference Lax-integrable equations. In references [3,4,6], similar results have already been obtained for several specific equations. The difference is that they proved this point in the space of the spectral data while we employed zero-curvature representations of flows. In addition, our method can be viewed as a more general version of the Lax–Darboux scheme (cf. [9]).
For continuous and differential-difference integrable systems, non-isospectral flows can also be used to construct symmetries, and in particular, some non-isospectral flows act as master symmetries to generate more symmetries. Non-isospectral flows are usually characterized by containing explicitly independent variables. In the case of P Δ Es, it has been pointed out that all the ABS equations have symmetries that depend explicitly on n (or m) [13], and are closely related to master symmetries [13]. We can also explore non-isospectral symmetries and master symmetries, provided that we can drive non-isospectral flows in our approach. However, for the lpKdV Equation, we are currently unable to derive non-isospectral flows from its Lax pair (46) (unless we change the asymptotic matrix U ( 0 ) ). This method we described in this paper can also be applied to other quadrilateral equations such as ABS equations [32] and equations from the Lax–Darboux scheme (e.g., [35]) and multi-component P Δ Es, which we will explore in the future research. In addition, Toda-type lattice equations have been shown to have symmetries [13]. These equations were classified by Adler [36,37]. Their integrability is understood as 2[1,1] type ABS-coupled equations (see [38]) by eliminating one component. It would be interesting to make clear how their symmetries are related to those of the ABS equations.
In this paper, the symmetry ξ of the P Δ E (1) is interpreted as a flow along the tangent direction of solution u ( t ) , which is in the solution manifold of the P Δ E. This approach relies on introducing t as a parameter. Such an idea has been used in describing Hamiltonian structures as well as symmetries for the continuous and differential-difference integrable equations (e.g., [19,27,39,40]). One can introduce a parameter as a continuous independent variable by implementing continuum limits. For the lpKdV Equation (45) that we considered in the previous section, in its continuum limits, the leading term gives rise to a differential-difference equation (see Equation (5.13) in [34] or Equation (4.6) in [41]) which is nothing but the first equation u t 0 = K ( 0 ) written in the form (71). There is also a way to introduce high order parameters { t l } and obtain a hierarchy u t l = S ( l ) from the continuum limits of the lpKdV Equation, see [42]. The obtained flows { S ( l ) } should be symmetries of the lpKdV Equation. In principle, it is understood the lpKdV Equation itself includes its infinitely many commuting symmetries via its continuum limits. One can also consider the spacing parameter p (or q) as a parameter. In that case, one can have non-autonomous (or non-isospectral) equations that involve explicit independent variables (see Chapter 10.4 of [34] and cf. [13]), which can also provide symmetries. It would be interesting to understand the connections of these symmetries obtained from different ways and also to understand how these symmetries are related to the solutions and conservation laws of the considered discrete integrable equations.

Author Contributions

Conceptualization, D.-j.Z.; Funding acquisition, D.-j.Z.; Investigation, J.C., J.L. and D.-j.Z.; Methodology, J.C., J.L. and D.-j.Z.; Project administration, D.-j.Z.; Supervision, D.-j.Z.; Writing–original draft, J.C. and J.L.; Writing–review and editing, D.-j.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the NSF of China (No. 12271334).

Data Availability Statement

No new data are created.

Acknowledgments

This project is supported by the NSF of China (No. 12271334).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nijhoff, F.W.; Papageorgiou, V.G. Similarity reductions of integrable lattices and discrete analogues of the Painlevé II equation. Phys. Lett. A 1991, 153, 337–344. [Google Scholar] [CrossRef]
  2. Fordy, A.; Xenitidis, P. Symmetries of ZN graded discrete integrable systems. J. Phys. A Math. Theor. 2020, 53, 235201. [Google Scholar] [CrossRef]
  3. Levi, D.; Martina, L. Integrable hierarchies of nonlinear difference-difference equations and symmetries. J. Phys. A Math. Gen. 2001, 34, 10357–10368. [Google Scholar] [CrossRef]
  4. Levi, D.; Petrera, M. Continuous symmetries of the lattice potential KdV equation. J. Phys. A Math. Theor. 2007, 40, 4141–4159. [Google Scholar] [CrossRef]
  5. Garifullin, R.N.; Gudkova, E.V.; Habibullin, I.T. Method for searching higher symmetries for quad-graph equations. J. Phys. A Math. Theor. 2011, 44, 325202. [Google Scholar] [CrossRef]
  6. Levi, D.; Petrera, M.; Scimiterna, C. The lattice Schwarzian KdV equation and its symmetries. J. Phys. A Math. Theor. 2007, 40, 12753–12761. [Google Scholar] [CrossRef]
  7. Levi, D.; Yamilov, R.I. The generalized symmetry method for discrete equations. J. Phys. A Math. Theor. 2009, 42, 454012. [Google Scholar] [CrossRef]
  8. Levi, D.; Winternitz, P.; Yamilov, R.I. Continuous Symmetries and Integrability of Discrete Equations; AMS: Providence, RI, USA, 2023. [Google Scholar]
  9. Mikhailov, A.V. Formal Diagonalisation of the Lax-Darboux Scheme and Conservation Laws of Integrable Partial Differential, Differential Difference and Partial Difference. DIS a Follow-Up Meeting, 8–12 July 2013 (Isaac Newton Institute for Mathematical Sciences). Available online: https://www.newton.ac.uk/event/disw05/ (accessed on 1 December 2023).
  10. Mikhailov, A.V. Formal diagonalisation of Lax-Darboux schemes. Model. Anal. Inform. Syst. 2015, 22, 795–817. [Google Scholar] [CrossRef]
  11. Mikhailov, A.V.; Wang, J.P.; Xenitidis, P. Recursion operators, conservation laws and integrability conditions for difference equations. Theor. Math. Phys. 2011, 167, 421–443. [Google Scholar] [CrossRef]
  12. Rasin, A.G. Infinitely many symmetries and conservation laws for quad-graph equations via the Gardner method. J. Phys. A Math. Theor. 2010, 43, 235201. [Google Scholar] [CrossRef]
  13. Rasin, O.G.; Hydon, P.E. Symmetries of integrable difference equations on the quad-graph. Stud. Appl. Math. 2007, 119, 253–269. [Google Scholar] [CrossRef]
  14. Tongas, A.; Tsoubelis, D.; Xenitidis, P. Affine linear and D4 symmetric lattice equations: Symmetry analysis and reductions. J. Phys. A Math. Theor. 2007, 40, 13353–13384. [Google Scholar] [CrossRef]
  15. Carpentier, S.; Mikhailov, A.V.; Wang, J.P. Rational recursion operators for integrable differential-difference equations. Commun. Math. Phys. 2019, 370, 807–851. [Google Scholar] [CrossRef]
  16. Chen, D.Y.; Zhang, D.J. Lie algebraic structures of (1+1)-dimensional Lax integrable systems. J. Math. Phys. 1996, 37, 5524–5538. [Google Scholar] [CrossRef]
  17. Chen, D.Y.; Zhang, H.W. Lie algebraic structure for the AKNS system. J. Phys. A Math. Gen. 1991, 24, 377–383. [Google Scholar] [CrossRef]
  18. Chen, D.Y.; Xin, H.W.; Zhang, D.J. Lie algebraic structures of some (1+2)-dimensional Lax integrable systems. Chaos Solitons Fractals 2003, 15, 761–770. [Google Scholar] [CrossRef]
  19. Fu, W.; Huang, L.; Tamizhmani, K.M.; Zhang, D.J. Integrability properties of the differential-difference Kadomtsev-Petviashvili hierarchy and continuum limits. Nonlinearity 2013, 26, 3197–3229. [Google Scholar] [CrossRef]
  20. Liu, J.; Zhang, D.J.; Zhao, X.H. Symmetries of the DΔmKP hierarchy and their continuum limits. Stud. Appl. Math. 2024, 152, 404–430. [Google Scholar] [CrossRef]
  21. Ma, W.X. K symmetries and τ symmetries of evolution equations and their Lie algebras. J. Phys. A: Math. Gen. 1990, 23, 2707–2716. [Google Scholar] [CrossRef]
  22. Ma, W.X.; Fuchssteiner, B. Algebraic structure of discrete zero curvature equations and master symmetries of discrete evolution equations. J. Math. Phys. 1999, 40, 2400–2418. [Google Scholar] [CrossRef]
  23. Zhang, D.J.; Chen, S.T. Symmetries for the Ablowitz-Ladik hierarchy: Part I. Four-potential case. Stud. Appl. Math. 2010, 125, 393–418. [Google Scholar] [CrossRef]
  24. Zhang, D.J.; Chen, S.T. Symmetries for the Ablowitz-Ladik hierarchy: Part II. Integrable discrete nonlinear Schrödinger equations and discrete AKNS hierarchy. Stud. Appl. Math. 2010, 125, 419–443. [Google Scholar] [CrossRef]
  25. Zhang, D.J.; Ning, T.K.; Bi, J.B.; Chen, D.Y. New symmetries for the Ablowitz-Ladik hierarchies. Phys. Lett. A 2006, 359, 458–466. [Google Scholar] [CrossRef]
  26. Zhang, J.B.; Chen, Q. The τ-symmetries and Lie algebra structure of the Blaszak-Marciniak lattice equation. Math. Method Appl. Sci. 2024, 47, 8160–8170. [Google Scholar] [CrossRef]
  27. Fuchssteiner, B. Hamiltonian structure and integrability. Math. Sci. Eng. 1991, 185, 211–256. [Google Scholar]
  28. Zhang, D.J.; Chen, D.Y. Hamiltonian structure of discrete soliton systems. J. Phys. A Math. Gen. 2002, 35, 7225–7241. [Google Scholar] [CrossRef]
  29. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  30. Nijhoff, F.W.; Capel, H.W. The discrete Korteweg-de Vries equation. Acta Appl. Math. 1995, 39, 133–158. [Google Scholar] [CrossRef]
  31. Nijhoff, F.W.; Quispel, G.R.W.; Capel, H.W. Direct linearization of nonlinear difference-difference equations. Phys. Lett. A 1983, 97, 125–128. [Google Scholar] [CrossRef]
  32. Adler, V.E.; Bobenko, A.I.; Suris, Y.B. Classification of integrable equationson quad-graphs. The consistency approach. Commun. Math. Phys. 2003, 233, 513–543. [Google Scholar] [CrossRef]
  33. Wahlquist, H.D.; Estabrook, F.B. Bäcklund transformation for solutions of the Korteweg-de Vries equation. Phys. Rev. Lett. 1973, 31, 1386–1390. [Google Scholar] [CrossRef]
  34. Hietarinta, J.; Joshi, N.; Nijhoff, F.W. Discrete Systems and Integrability; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
  35. Cao, C.W.; Zhang, G.Y. Lax pairs for discrete integrable equations via Darboux transformations. Chin. Phys. Lett. 2012, 29, 050202. [Google Scholar] [CrossRef]
  36. Adler, V.E. On the structure of the Bäcklund transformations for the relativistic lattices. J. Nonlinear Math. Phys. 2000, 7, 34–56. [Google Scholar] [CrossRef]
  37. Adler, V.E. Discrete equations on planar graphs. J. Phys. A Math. Gen. 2001, 34, 10453–10460. [Google Scholar] [CrossRef]
  38. Zhang, D.D.; van der Kamp, P.H.; Zhang, D.J. Multi-component extension of CAC systems. Symmetry Integr. Geom. Meth. Appl. 2020, 16, 060. [Google Scholar] [CrossRef]
  39. Magri, F. A simple model of the integrable Hamiltonian equation. J. Math. Phys. 1978, 19, 1156–1162. [Google Scholar] [CrossRef]
  40. Miwa, T.; Jimbo, M.; Date, E. Solitons: Differential Equations, Symmetries and Infinite Dimensional Algebras; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  41. Levi, D.; Rodriguez, M.A. Yamilov’s theorem for differential and difference equations. Ufa Math. J. 2021, 13, 152–159. [Google Scholar] [CrossRef]
  42. Wiersma, G.L.; Capel, H.W. Lattice equations, hierarchies and Hamiltonian structures. Phys. A 1987, 142, 199–244. [Google Scholar] [CrossRef]
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

Cheng, J.; Liu, J.; Zhang, D.-j. On Symmetries of Integrable Quadrilateral Equations. Symmetry 2024, 16, 744. https://doi.org/10.3390/sym16060744

AMA Style

Cheng J, Liu J, Zhang D-j. On Symmetries of Integrable Quadrilateral Equations. Symmetry. 2024; 16(6):744. https://doi.org/10.3390/sym16060744

Chicago/Turabian Style

Cheng, Junwei, Jin Liu, and Da-jun Zhang. 2024. "On Symmetries of Integrable Quadrilateral Equations" Symmetry 16, no. 6: 744. https://doi.org/10.3390/sym16060744

APA Style

Cheng, J., Liu, J., & Zhang, D. -j. (2024). On Symmetries of Integrable Quadrilateral Equations. Symmetry, 16(6), 744. https://doi.org/10.3390/sym16060744

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