Next Article in Journal
Synchronized Cyclograms to Assess Inter-Limb Symmetry during Gait in Women with Anorexia and Bulimia: A Retrospective Study
Next Article in Special Issue
Study of Quantum Difference Coupled Impulsive System with Respect to Another Function
Previous Article in Journal
A New Method for Finding Lie Point Symmetries of First-Order Ordinary Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Classical Symmetry Analysis of a Class of Nonlinear Lattice Equations

1
Faculty of Teacher Education, Qiongtai Normal University, Haikou 571127, China
2
School of Educational Sciences, Harbin Normal University, Harbin 150080, China
*
Authors to whom correspondence should be addressed.
Symmetry 2023, 15(12), 2199; https://doi.org/10.3390/sym15122199
Submission received: 18 October 2023 / Revised: 13 November 2023 / Accepted: 12 December 2023 / Published: 14 December 2023
(This article belongs to the Special Issue Symmetries in Differential Equations and Application—2nd Edition)

Abstract

:
In this paper, a non-classical symmetry method for obtaining the symmetries of differential–difference equations is proposed. The non-classical symmetry method introduces an additional constraint known as the invariant surface condition, which is applied after the infinitesimal transformation. By solving the governing equations that satisfy this condition, we can obtain the corresponding reduced equation. This allows us to determine the non-classical symmetry of the differential–difference equation. This method avoids the complicated calculation involved in extending the infinitesimal generator and allows for a wider range of symmetry forms. As a result, it enables the derivation of a greater number of differential–difference equations. In this paper, two kinds of (2+1)-dimensional Toda-like lattice equations are taken as examples, and their corresponding symmetric and reduced equations are obtained using the non-classical symmetry method.

1. Introduction

A lattice system is a mathematical representation of a network or lattice of particles. It can be used to describe various phenomena, such as the vibrations of atoms and molecules in a matrix or the oscillations of particles connected by a spring [1]. Nonlinear lattice dynamics has a wide range of applications in physics, biology, statistical physics, acoustics, optics, condensed matter physics, and other fields [2,3]. The mathematical model of a lattice system often takes the form of nonlinear differential–difference equations [4].
The Lie symmetry method is not only an effective approach for obtaining exact solutions of partial differential equations but also for differential–difference equations [5,6]. The Lie symmetry analysis method has also been proposed for finding similarity reduction and exact solutions of nonlinear evolution equations [7]. In 1993, Levi [8] proposed the classical Lie symmetry method for differential–difference equations. The Toda equation is a classical model of differential–difference equation [9,10,11,12], whose symmetries have been studied [13] and the differential–difference Lie symmetry method was applied to solve a class of Toda-like lattice equations [14]. In 2022, a survey of the connection between orthogonal polynomials, Toda lattices and related lattices, and Painlevé equations (discrete and continuous) was given [15].
In 1969, Bluman and Cole [16] generalized the classical symmetry method of Lie to a non-classical method. They proposed the non-classical symmetry method for the first time and obtained a new exact solution of the one-dimensional heat conduction equation using this new method. For a partial differential equation, we know that the given system of partial differential equations is invariant under symmetry in the usual symmetry method. In the non-classical method [17,18], not only is the given system invariant, but the invariant surface conditions and their differential consequences are also invariant under the corresponding symmetry. Moreover, the over-determined system of partial differential equations is nonlinear in the non-classical method. This implies that by using non-classical symmetry analysis, one could obtain new symmetries that go beyond those obtained by using the classical symmetry method. Bluman and Tian [19] utilized the non-classical symmetry method to solve the nonlinear Kompaneets equation. Xin et al. [20] obtained nonlocal symmetries and solutions for the (2+1) dimension integrable Burgers equation. ARRIGO et al. [21] discussed nonclassical symmetry solutions and the methods between Bluman–Cole and Clarkson–Kruskal. Subhankar and Raja Sekhar [22] studied nonclassical symmetry analysis and conservation laws of a one-dimensional macroscopic production model and the evolution of nonlinear waves.
We consider the invariance of a given differential equation
Δ = Δ ( x , t , u , u ( 1 ) , , u ( n ) ) = 0 ,
together with its invariant surface condition
Δ 0 = ξ ( x , t , u ) u x + η ( x , t , u ) u t ϕ ( x , t , u ) = 0 ,
where u = u ( x , t ) . Equation (2) describes the solution surface, which remains invariant under a one parameter group of transformations with the infinitesimal generator
Γ = i = 1 p ξ i ( x , u ) x i + ϕ ( x , u ) u .
Γ ( m ) is the m-order prolongation of Γ
Γ ( m ) = Γ + i = 1 p M ϕ [ x i ] M ( x , u ( m ) ) u x i ,
where
ϕ [ x i ] M ( x , u ) = D M ( ϕ [ x i ] i = 1 p ξ i u x i ) + i = 1 p ξ i u x i , M ,
with M = ( m 1 , , m i ) , 1 i p .
It is assumed that Equation (1) is invariant under the action of infinitesimal generators (3) and invariant surface conditions (2), if and only if the system of governing equations is satisfied
Γ ( n ) Δ | { Δ = 0 , Δ 0 = 0 } = 0 ,
In this way, the non-classical symmetry of Equation (1) can be obtained.
The question that naturally arises is whether the non-classical method can be extended to discuss discrete equations. To the best of our knowledge, this question has not been reported so far. Based upon the discrete Lie symmetry method [23], we extend the infinitesimal generator and its prolongation to the discrete performance to analyze the non-classical symmetries of differential–difference equations in this Letter.
In Section 2, the non-classical symmetry method for differential–difference equations is introduced. In Section 3 and Section 4, the non-classical symmetry method of differential–difference equations is applied to study two types of (2+1)-dimensional Toda-like lattice equations. The determinant equations of the equation are derived, and the corresponding reduced equations are obtained by solving the determinant equations. Conclusions are finally presented in Section 5.

2. Non-Classical Symmetry of Differential–Difference Equations

In this paper, we will introduce the non-classical symmetric method using an example of a differential–difference equation of order m that contains both continuous variables x = ( x 1 , x 2 , , x m ) and a discrete variable n. The equation is formulated as follows:
Δ = Δ n ( m ) Δ ( x , n , u ( n + k ) | k = a b , u x i ( n + k ) | k = a i b i , u x i x j ( n + k ) | k = a i j b i j ) = 0 ,
where all a , b , a i , b i , a i j , b i j are finite non-negative integers, m = Max { b , b i , b i j } . The infinitesimal generators acting on Equation (7) are as follows:
Γ = i = 1 p ξ i ( x , u ( n ) ) x i + ϕ ( x , u ( n ) , n ) u ( n ) .
Γ ( m ) is the m-order prolongation of Γ , which is formed as
Γ ( m ) = Γ + i = 1 p M ϕ [ x i ] M ( x , u m ( n ) , n ) u x i ( n ) ,
where
ϕ [ x i ] M ( x , u ( n ) , n ) = D M ( ϕ [ x i ] ( n ) i = 1 p ξ i u x i ( n ) ) + i = 1 p ξ i u x i , M ( n ) ,
with M = ( m 1 , , m i ) , 1 i p .
To better understand the non-classical symmetry of differential–difference equations, consider a differential–difference equation involving x 1 = x and x 2 = t as an example.
Δ Δ ( x , t , u ( n ) , u ( n + 1 ) , , u x ( n ) , u x ( n + 1 ) , , u t ( n ) , u t ( n + 1 ) , , u x t ( n ) , ) = 0 ,
Then the infinitesimal transformation acting on this equation is
x * = x + ξ ( x , t , u ( n ) ) ε + O ( ε 2 ) , t * = t + η ( x , t , u ( n ) ) ε + O ( ε 2 ) , u * = u + ϕ ( x , t , u ( n ) ) ε + O ( ε 2 ) .
The infinitesimal transformation is derived from the following equation
u ( x + ε ξ , t + ε η ) = u ( x , t , u ( n ) ) + ϕ ( x , t , u ( n ) ) ε + O ( ε 2 ) .
The Equation (13) is about taking the derivative ε and expanding it; thus, it can be written as follows
Δ 0 = ξ ( x , t , u ( n ) ) u x ( n ) + η ( x , t , u ( n ) ) u t ( n ) ϕ ( x , t , u ( n ) , n ) = 0 .
Equation (14) is called the invariant surface condition of differential–difference equation.
According to the classical Lie symmetry method, the differential–difference Equation (11) is assumed to be invariant under the infinitesimal transformations Δ 0 . The non-classical symmetry method adds invariant surface conditions on this basis, which is also invariant under infinitesimal transformations; that is, the governing equations are satisfied for Δ and Δ 0 at the same time:
Γ ( 2 ) Δ | { Δ = 0 , Δ 0 = 0 } = 0 ,
where the infinitesimal generator is
Γ = ξ ( x , t , u ( n ) ) x + η ( x , t , u ( n ) ) t + ϕ ( x , t , u ( n ) , n ) u ( n ) .
For the convenience of later calculation, the first order and second order continuation are given respectively
Γ ( 1 ) = Γ + i = n L n + M ϕ [ x ] ( x , t , u ( i ) , i ) u x ( i ) + i = n L n + M ϕ [ t ] ( x , t , u ( i ) , i ) u t ( i ) , Γ ( 2 ) = Γ ( 1 ) + i = n L n + M ϕ [ x x ] ( x , t , u ( i ) , i ) u x x ( i ) + i = n L n + M ϕ [ x t ] ( i ) ( x , t , u ( i ) , i ) u x t ( i ) + i = n L n + M ϕ [ t t ] ( x , t , u ( i ) , i ) u t t ( i ) ,
where n L i n + M , L , M are non-negative integers, for short ξ ( x , t , u ( n ) ) = ξ , η ( x , t , u ( n ) ) = η , ϕ ( x , t , u ( i ) , i ) = ϕ ( i ) ; hence,
ϕ [ x ] ( i ) = D x ( ϕ ( i ) ξ u x ( i ) η u t ( i ) ) + ξ u x x ( i ) + η u x t ( i ) ,
ϕ [ t ] ( i ) = D t ( ϕ ( i ) ξ u x ( i ) η u t ( i ) ) + ξ u x t ( i ) + η u t t ( i ) ,
ϕ [ x x ] ( i ) = D x x ( ϕ ( i ) ξ u x ( i ) η u t ( i ) ) + ξ u x x x ( i ) + η u x x t ( i ) ,
ϕ [ x t ] ( i ) = D x t ( ϕ ( i ) ξ u x ( i ) η u t ( i ) ) + ξ u x x t ( i ) + η u t t x ( i ) ,
ϕ [ t t ] ( i ) = D t t ( ϕ ( i ) ξ u x ( i ) η u t ( i ) ) + ξ u x t t ( i ) + η u t t t ( i ) .
We substituted (18)–(22) into the governing Equation (17) by solving the determining Equation (15); that is, the coefficient equations of the derivative of u to x , t , the expression about ξ , η , ϕ can be obtained and its corresponding characteristic equation can also be obtained under the integral of the invariant surface conditions
d x ξ ( x , t , u ( n ) ) = d t η ( x , t , u ( n ) ) = d u n ϕ ( x , t , u ( n ) , n ) .
By solving this characteristic equation, the invariants of Equation (1) can be found. In this article, we only need to consider two cases under the invariant surface condition (8): η 1 and η 0 , ξ 1 .
The invariant surface condition Δ 0 introduced a new correlation equation for the derivative of Equation (11). We substitute the invariant surface condition Δ 0 into the governing system (15), replacing the related terms and further reducing the governing system (15) to obtain the determining system of sum ξ , η and ϕ . The system of determining equations obtained by the non-classical symmetric method is usually nonlinear, while solving the system of determining equations gives a reduced equation with the same solution as the control equation. Finally, the group invariant solutions of Equation (11) can be obtained by solving the reduced equations. This is the non-classical symmetric approach to differential–difference equations. In simple terms, the non-classical symmetry method seeks the invariance of the original equation under the invariant surface condition.
Summarizing the above, two theorems are obtained.
Theorem 1.
Suppose σ ( U ) is a symmetry of the following differential–difference equations,
Δ ( n , t , δ U ) = 0 ,
then U is the group invariant solution of the equation corresponding to the invariant group of σ ( U ) if and only if
Δ ( n , t , δ U ) = 0 , σ ( U ) = 0 .
where δ is a difference operator and ∂ is a differetial operator.
Theorem 2.
Suppose (24) is a system of differential–difference equations defined on X × t × U with a maximum rank. G is the group of local transformations acting on X × t × U , and if any of the generators of G are true when
Γ ( n ) ( Δ ) | { Δ = 0 , Δ 0 = 0 } = 0 ,
then G is the invariant group of this system of equations, where n X .
Theorem 2 can be called the Lie criterion for non-classical Lie symmetry.

3. Non-Classical Symmetry of (2+1)-Dimensional Toda Lattice Equations

The (2+1)-dimensional Toda equation is a completely integrable differential–difference equation with a Lax pair, an infinitely multiple conservation law, a Backlund transform, a soliton solution, and the general characteristics of integrable models. The study of the Toda equation is of great significance in the study of fully integrable nonlinear systems. The (2+1)-dimensional Toda lattice Equation [8] has the following form
Δ = u x t ( n ) e u ( n 1 ) u ( n ) + e u ( n ) u ( n + 1 ) ,
where u ( n ) is the function of x , t .
The invariant surface condition is
Δ 0 = ξ ( x , t , u ( n ) ) u x ( n ) + η ( x , t , u ( n ) ) u t ( n ) ϕ ( x , t , u ( n ) , n ) .
As a brief note ξ ( x , t , u ( n ) ) = ξ , η ( x , t , u ( n ) ) = η , ϕ ( x , t , u ( n ) , n ) = ϕ ( n ) , the system of equations determining the non-classical symmetry of Equation (27) is obtained from the following governing Equation (15).
The infinitesimal generator Γ is
Γ = ξ x + η t + ϕ ( n ) u ( n ) ,
Its second-order prolongation Γ ( 2 ) is
Γ ( 2 ) = ϕ ( n 1 ) u ( n 1 ) + ϕ ( n ) u ( n ) + ϕ ( n + 1 ) u ( n + 1 ) + ϕ x ( n 1 ) u x ( n 1 ) + ϕ x ( n ) u x ( n ) + ϕ x ( n + 1 ) u x ( n + 1 ) + ϕ t ( n 1 ) u t ( n 1 ) + ϕ t ( n ) u t ( n ) + ϕ t ( n + 1 ) u t ( n + 1 ) + ϕ x x ( n 1 ) u x x ( n 1 ) + ϕ x x ( n ) u x x ( n ) + ϕ x x ( n + 1 ) u x x ( n + 1 ) + ϕ x t ( n 1 ) u x t ( n 1 ) + ϕ x t ( n ) u x t ( n ) + ϕ x t ( n + 1 ) u x t ( n + 1 ) + ϕ t t ( n 1 ) u t t ( n 1 ) + ϕ t t ( n ) u t t ( n ) + ϕ x t ( n + 1 ) u t t ( n + 1 ) ,
where
ϕ x ( n 1 ) u x ( n 1 ) = ϕ x ( n ) u x ( n ) = ϕ x ( n + 1 ) u x ( n + 1 ) = 0 ,
ϕ t ( n 1 ) u t ( n 1 ) = ϕ t ( n ) u t ( n ) = ϕ t ( n + 1 ) u t ( n + 1 ) = 0 ,
ϕ x x ( n 1 ) u x x ( n 1 ) = ϕ x x ( n ) u x x ( n ) = ϕ x x ( n + 1 ) u x x ( n + 1 ) = 0 ,
ϕ t t ( n 1 ) u t t ( n 1 ) = ϕ t t ( n ) u t t ( n ) = ϕ t t ( n + 1 ) u t t ( n + 1 ) = 0 ,
ϕ x t ( n 1 ) u x t ( n 1 ) = ϕ x t ( n + 1 ) u x t ( n + 1 ) = 0 ,
ϕ [ x t ] ( n ) = ϕ x t ( n ) u x t ( n ) = D x t ( ϕ ( n ) ξ u x ( n ) η u t ( n ) ) + ξ u x x t ( n ) + η u x t t ( n ) = ϕ x t + ϕ u x u t + ( ϕ u u u t ) u x + ϕ u u x t ( ξ x t + u t ξ x u ) u t ξ x u t t ξ u u x t u t ξ u u x u t t ( ξ u t + ξ u u u t ) u x u t ( η x t + η x u u t ) u x 2 η u u x u x t ( η u t + η u u u t ) u x 2 ( ξ t + ξ u u t ) u x t ( η t + η u u t ) u x x η x u x t .
Substituting (31)–(36) into Γ ( 2 ) (30), the control system (15) will take the following form:
Γ ( 2 ) Δ | { Δ = 0 , Δ 0 = 0 } = ϕ ( n 1 ) e u ( n 1 ) u ( n ) + ϕ ( n ) ( e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ) + ϕ ( n + 1 ) e u ( n ) u ( n + 1 ) + ϕ x t + ϕ x u ( ϕ ( n ) ξ u x ) + ϕ u t u x ξ x u ( ϕ ( n ) ξ u x ) 2 ξ u t u x ( ϕ ( n ) ξ u x ) ξ u u u x ( ϕ ( n ) ξ u x ) 2 ξ t ( e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ) + ϕ u u u x ( ϕ ( n ) ξ u x ) ξ u u x [ ϕ t ξ t u x ξ ( e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ) ] + ϕ u ( e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ) ξ x t ( ϕ ( n ) ξ u x ) ξ x ( ϕ t ξ t u x ξ ) ( e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ) 2 ξ u ( e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ) ( ϕ ( n ) ξ u x ) = 0 ,
According to the non-classical symmetry method, the non-classical symmetries of the (2+1) dimensional Toda-like lattice Equation (27) only need to be discussed in two cases; one is η 1 and the other case is η 0 , ξ 1 :
Case 1. When η 1 , the system of equations of ξ and ϕ is determined by
ϕ u t + ϕ ϕ u u ϕ ξ x t ϕ ξ u x ξ x ϕ t ϕ ξ u u + ϕ ξ x ξ u u = 0 , ξ t ξ u + ξ ξ u t ξ ξ x u ξ ϕ u u = 0 , ϕ u + ξ ξ x ϕ ξ u ξ t = 0 , ϕ ( n + 1 ) ϕ = 0 , ϕ ϕ ( n 1 ) = 0 , ξ ξ u u = 0 , ξ ξ u = 0 .
Using the symbolic computing software Maple to solve Equation (38), the solutions are
ξ = A ( x ) , ϕ = B ( t ) ,
where A ( x ) is an arbitrary function about x , B ( t ) is an arbitrary function of t .
Accordingly, the non-classical symmetry of Equation (27) is given by
Y 1 = A ( x ) x + t + B ( t ) u ( n ) .
Choose A ( x ) = a x 2 + b x + c , B ( t ) = k t + d , where k , a , b , c and d are arbitrary constants as an example. Substituting the non-classical symmetry Y 1 to the invariant surface condition (28), the following system of equations, also called symmetry reduction, is obtained:
u t ( n ) = k t + d ( a x 2 + b x + c ) u x ( n ) , u x t ( n ) = e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) .
By simplifying Equation (41), a reduced equation of the same dimension as Equation (27) can be obtained
( 2 a x + b ) u x ( n ) + ( a x 2 + b x + c ) u x x ( n ) + e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) = 0 .
If we choose a = 1 , b = 0 and c = 0 , Equation (42) will be simplified to
2 x u x ( n ) + x 2 u x x ( n ) + e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) = 0 .
Case 2. When η 0   a n d ξ 1 , the determining system of equations of ϕ is
ϕ u + ϕ ( n 1 ) + ϕ = 0 , ϕ ( n + 1 ) ϕ ϕ u = 0 , ϕ u x + ϕ ϕ u u = 0 , ϕ x t + ϕ ϕ u x = 0 ,
Using the symbolic computing software Maple to solve the Equation (44), we can obtain
ϕ = C ( t ) ,
where C ( t ) is an arbitrary function of t .
Accordingly, the non-classical symmetry of Equation (27) is given by
Y 2 = x + C ( t ) u ( n ) .
If we choose C ( t ) = p t + q , with arbitrary constants p and q as an example, and substitute the non-classical symmetry Y 2 to the invariant surface condition (28), we obtain the following system of equations:
u x ( n ) = p t + q , u x t ( n ) = e u ( n 1 ) u ( n ) e u ( n ) u ( n + 1 ) ,
By simplifying the Equation (47), a reduced equation of the same dimension as Equation (17) can be obtained
p e u ( n 1 ) u ( n ) + e u ( n ) u ( n + 1 ) = 0 .
With the help of Maple, the reduced Equation (48) could be solved under the initial condition u ( 0 ) = 0 , u ( 1 ) = 1 and p = 0 ; that is,
u ( n ) = n .
By using the non-classical symmetry method, the two symmetries Y 1 and Y 2 obtained are about the change of the functions A ( x ) , B ( t ) and C ( t ) , respectively. According to the different functions selected, the corresponding reduced equations of different forms can be obtained.
The above solving process can be summarized as the following theorem.
Theorem 3.
(2+1)-dimensional Toda Equation (27) has two non-classical symmetries, as follows
Y 1 = A ( x ) x + t + B ( t ) u ( n ) , η 1 ,
Y 2 = x + C ( t ) u ( n ) , η 0 , ξ 1 .
The two symmetries Y 1 and Y 2 are changed according to A ( x ) , B ( t ) and C ( t ) individually. Depending on which is selected between A ( x ) , B ( t ) and C ( t ) , different reduced equations can be obtained.
In [23], four Lie point symmetries were given
X 1 = t , X 2 = t u ( n ) , X 3 = u ( n ) , X 4 = t t + 2 n u ( n ) .
It could be easily seen that X i ( i = 1 , 2 , 4 ) was included by Y 1 .

4. Non-Classical Symmetry of Another Type of (2+1)-Dimensional Toda Lattice Equations

The other (2+1)-dimensional Toda lattice equation has the following form:
Δ = u x t ( n ) [ u t ( n ) + 1 ] [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] = 0 .
The invariant surface condition is
Δ 0 = ξ ( x , t , u ( n ) ) u x ( n ) + η ( x , t , u ( n ) ) u t ( n ) ϕ ( x , t , u ( n ) , n ) = 0 .
The system of equations determining the non-classical symmetry of Equation (50) is obtained from the following governing Equation (15). Where the infinitesimal generator Γ and its second-order prolongation are
Γ = ξ ( x , t , u ( n ) ) x + η ( x , t , u ( n ) ) t + ϕ ( x , t , u ( n ) , n ) u ( n ) ,
Γ ( 2 ) = ϕ ( n 1 ) u ( n 1 ) + ϕ ( n ) u ( n ) + ϕ ( n + 1 ) u ( n + 1 ) + ϕ [ t ] ( n ) u t ( n ) + ϕ [ x t ] ( n ) u x t ( n ) ,
where ξ ( x , t , u ( n ) ) = ξ , η ( x , t , u ( n ) ) = η , ϕ ( x , t , u ( n ) , n ) = ϕ ( n ) for short. Hence,
ϕ [ t ] ( n ) = D t ( ϕ ( n ) ξ u x ( n ) η u t ( n ) ) + ξ u x t ( n ) + η u t t ( n ) ,
ϕ [ x t ] ( n ) = D x t ( ϕ ( n ) ξ u x ( n ) η u t ( n ) ) + ξ u x x t ( n ) + η u t t x ( n ) .
The non-classical symmetries of the (2+1) dimensional Toda-like lattice Equation (50) in two cases, η 1 and η 0 , ξ 1 , are discussed as follows:
Case 3. When η 1 , the system of determining equations about ξ and ϕ ( n ) , is
ϕ x t ( n ) + ϕ ( n ) ϕ u x ( n ) ϕ ( n 1 ) + 2 ϕ ( n ) ϕ ( n + 1 ) ϕ ( n ) ϕ ( n 1 ) ϕ ( n ) ϕ ( n + 1 ) + 2 ϕ 2 ( n ) = 0 , ϕ u t ( n ) ξ x t ( n ) ϕ ( n ) ξ u x + ϕ ( n ) ϕ u u ( n ) ξ ϕ u x ( n ) + ξ ϕ ( n 1 ) 2 ξ ϕ ( n ) + ξ ϕ ( n + 1 ) = 0 , ϕ t ( n ) + ξ x + ϕ ( n ) ξ x ϕ u ( n ) ϕ ( n ) ϕ u ( n ) = 0 , ξ ξ u x ξ ϕ u u ( n ) ξ u t ϕ ( n ) ξ u u = 0 , ξ t ϕ ( n ) ξ u 2 ξ u + ξ ξ x = 0 , ξ t + ϕ ( n ) ξ u = 0 , ξ ξ u u = 0 , ξ ξ u = 0 .
Using the symbolic computing software Maple to solve the Equation (55), we will obtain
ξ ( x , t , u ( n ) ) = 0 , ϕ ( x , t , u ( n ) ) = a ,
ξ ( x , t , u ( n ) ) = 0 , ϕ ( x , t , u ( n ) ) = D ( x ) ,
ξ ( x , t , u ( n ) ) = b , ϕ ( x , t , u ( n ) ) = E ( x ) + t + u ( n ) F ( x ) t ,
where a , b are arbitrary constants, and D ( x ) , E ( x ) and F ( x ) are arbitrary functions of x.
Accordingly, the non-classical symmetry of Equation (50) are represented as
Y 1 = t + a u ( n ) , Y 2 = t + D ( x ) u ( n ) , Y 3 = b x + t + E ( x ) + t + u ( n ) F ( x ) t u ( n ) .
Substituting a = 1 2 into the non-classical symmetry Y 1 and applying it to the invariant surface condition (51) for example, we obtain the following system of equations
u t ( n ) = 1 2 , u x t ( n ) = [ u t ( n ) + 1 ] [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] .
By simplifying the Equation (60), a reduced equation of the same dimension as Equation (50) can be obtained
u ( n 1 ) 2 u ( n ) + u ( n + 1 ) = 0 .
With the help of Maple, the reduced Equation (61) could be solved under the initial condition u ( 0 ) = 1 , u ( 1 ) = 2 , that is
u ( n ) = n + 1 .
If we choose to substitute D ( x ) = x into the non-classical symmetry Y 2 to the invariant surface condition (51) as an example, the following system of equations are obtained:
u t ( n ) = x , u x t ( n ) = [ u t ( n ) + 1 ] [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] .
By simplifying the Equation (63), a reduced equation of the same dimension as Equation (50) can be obtained:
( x + 1 ) [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] 1 = 0 .
Substituting b = 1 and E ( x ) = 2 x , F ( x ) = x into the non-classical symmetry Y 2 and applying it to the invariant surface condition (51), the following system of equations is as follows:
u t ( n ) = u x ( n ) + 2 x + t + u ( n ) x t , u x t ( n ) = [ u t ( n ) + 1 ] [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] .
By simplifying the Equation (65), a reduced equation of the same dimension as Equation (50) can be obtained
( x t ) 2 u x x ( n ) + 3 t + u ( n ) [ ( x t ) 2 u x ( n ) + 3 x 2 3 x t + ( x t ) u ( n ) ] [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] = 0 .
Case 4. When η 0 , ξ 1 , the system of determining equations about ϕ ( n ) is
ϕ u x ( n ) + ϕ ( n ) ϕ u u ( n ) ϕ ( n 1 ) + 2 ϕ ( n ) ϕ ( n + 1 ) = 0 , ϕ x t ( n ) + ϕ ( n ) ϕ u t ( n ) ϕ ( n 1 ) + 2 ϕ ( n ) ϕ ( n + 1 ) = 0 , ϕ u ( n ) ϕ t ( n ) = 0 .
Using the symbolic computing software Maple to solve the Equation (67), then
ϕ = G ( x ) + H ( t + u ( n ) ) ,
where G ( x ) is an arbitrary function of x, and H ( t + u ( n ) ) is an arbitrary function of ( t + u ( n ) ) .
Accordingly, the non-classical symmetry of Equation (50) is expressed as
Y 4 = x + [ G ( x ) + H ( t + u ) ] u ( n ) .
If we choose G ( x ) = x and H ( t + u ) = t + u ( n ) for Y 4 in the invariant surface condition (51), we obtain
u x ( n ) = x + t + u ( n ) , u x t ( n ) = [ u t ( n ) + 1 ] [ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] .
By simplifying the Equation (70), a reduced equation of the same dimension as Equation (50) can be obtained
[ u ( n 1 ) 2 u ( n ) + u ( n + 1 ) ] 1 = 0 .
With the help of Maple, the reduced Equation (71) could be solved under the initial condition u ( 0 ) = 1 and u ( 1 ) = 1 ; that is,
u ( n ) = 1 2 n 2 n 2 + 1 = 1 2 ( n 1 2 ) 2 + 7 8
Figure 1 shows the image of u ( n ) (72).
The above solving process can be summarized as the following theorem:
Theorem 4.
The (2+1)-dimensional Toda equation (50) has four non-classical symmetries, as follows
Y 1 = t + a u ( n ) , Y 2 = t + D ( x ) u ( n ) , Y 3 = b x + t + E ( x ) + t + u ( n ) F ( x ) t u ( n ) , Y 4 = x + [ G ( x ) + H ( t + u ( n ) ) ] u ( n ) .
The four symmetries Y 1 , Y 2 , Y 3 and Y 4 are changed according to a , b , D ( x ) , E ( x ) , F ( x ) , G ( x ) and H ( t + u ( n ) ) individually. Depending on which function is selected, different reduced equations can be obtained.
The main steps for solving lattice equations by using the differential–difference non-classical symmetry method are as follows:
Step 1: The non-classical Lie symmetry algorithm is used to calculate the Lie transform group, which can be described by its infinitesimal generator, which preserves the form of the equation on its extension space.
Step 2: Extend the Lie group and its infinitesimal generators to include the space of independent variables, function variables, and finite-order derivatives of function variables, in order to obtain the extended expressions.
Step 3: Assuming that the differential–difference equation is invariant under infinitesimal transformations, the invariant surface condition remains unchanged under infinitesimal transformations. As a result, the system of governing equations is obtained.
Step 4: Further reduce the system of governing equations to obtain a system of deterministic equations and a definite system of equations, with respect to infinitesimal generators.
Step 5: The non-classical symmetry of the differential–difference equation can be obtained by solving the system of equations using Maple 2020 software.
Step 6: Depending on the specific circumstances, various parameter equations can be chosen to derive low-dimensional reduction equations that yield the same solutions as the differential–difference equation.

5. Conclusions

In this paper, we combine the classical Lie-symmetry method of differential–difference equations proposed by Levi with the non-classical symmetry method of nonlinear partial differential equations proposed by Bluman. Our aim is to extend the non-classical symmetry method to differential–difference equations in order to get more symmetries of differential–difference equations. This method incorporates the invariant surface condition (2) into the classical Lie symmetry method in order to derive the governing equations of the studied differential–difference equation. By solving the governing equations under the η 1 or η 0 , ξ 1 , we can obtain its determination systems and then work out the non-classical symmetries of the equation and the corresponding symmetry reduction. The non-classical symmetries are only applicable to solutions that satisfy the invariant surface conditions. This implies that nonclassical symmetry only leaves a small sub-manifold of the solution manifold invariant. Finally, the effectiveness of this method is illustrated by considering examples of two types of (2+1)-dimensional Toda-like lattice equations. It can be seen that this method not only simplifies the calculation steps with the help of the invariant surface condition as an additional equation but also yields many new functional symmetries (40) and (46) of Toda-like lattice Equation (27), which cannot be obtained by the classical Lie symmetry method [23]. According to the selected arbitrary functions A ( x ) , B ( t ) , C ( t ) of the infinitesimal generators, one can obtain the corresponding non-classical symmetries and the corresponding reduced equation of the (2+1)-dimensional Toda-like lattice equations.

Author Contributions

Methodology, W.L.; software, Y.C. and W.L.; investigation, W.L. and K.J.; writing, W.L. and Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Hainan Provincial Natural Science Foundation of China grant number 122QN332.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Toda, M. Theory of Nonlinear Lattices; Springer: Berlin/Heidelberg, Germany, 1989; pp. 159–176. [Google Scholar]
  2. Zhang, J.; Meng, J.; Zheng, C. Folded solitary waves and foldons in the (2+1) -dimensional breaking soliton Equation. Chaos Soliton Fractals 2004, 20, 523–527. [Google Scholar] [CrossRef]
  3. Elkenany, B.; Alfrnwani, O.A. Behavior of lattice dynamics, electronic, optical, and mechanical properties of pentanary GaxIn1-xPySbzAs1-y-z alloy lattice matched to GaSb substrate under temperature. J. Phys. Chem. Solids 2023, 175, 111192. [Google Scholar] [CrossRef]
  4. Antoniades, C.; Christofides, P.D. Feedback control of nonlinear differential–difference equation systems. Chem. Eng. Sci. 1999, 54, 5677–5709. [Google Scholar] [CrossRef]
  5. Olver, P. Applications of Lie Groups to Differential Equations; Springer: Berlin/Heidelberg, Germany, 1986; pp. 242–379. [Google Scholar]
  6. Bakkyaraj, T.; Thomas, R. Lie symmetry analysis and exact solution of (2+1)-dimensional nonlinear time-fractional differential–difference equations. Pramana J. Phys. 2022, 96, 225. [Google Scholar] [CrossRef]
  7. Sahoo, S.; Ray, S.S. Lie symmetry analysis and exact solutions of (3+1) dimensional Yu–Toda–Sasa–Fukuyama equation in mathematical physics. Comput. Math. Appl. 2017, 73, 253–260. [Google Scholar] [CrossRef]
  8. Levi, D.; Winternitz, P. Symmetries and conditional symmetries of differential–difference equations. J. Math. Phys. 1993, 34, 3713–3730. [Google Scholar] [CrossRef]
  9. Duarte, V.N. A class of solutions of the two-dimensional Toda lattice equation. Phys. Lett. A 2021, 385, 126979. [Google Scholar] [CrossRef]
  10. Wattis, J.A.D.; Gordoa, P.R.; Pickering, A. Behaviour of the extended Toda lattice. Commun. Nonlinear Sci. Numer. Simul. 2015, 28, 138–150. [Google Scholar] [CrossRef]
  11. Guseinov, G.S.H. A class of complex solutions to the finite Toda lattice. Math. Comput. Model. 2013, 57, 1190–1202. [Google Scholar] [CrossRef]
  12. Bracciali, C.F.; Pérez, T.E. Bivariate orthogonal polynomials, 2D Toda lattices and Lax-type pairs. Appl. Math. Comput. 2017, 309, 142–155. [Google Scholar] [CrossRef]
  13. Qian, X.; Lou, S. Symmetries of the (2+1)-dimensional discrete Toda equation. Acta Phys. Sin. 1996, 45, 721–728. [Google Scholar] [CrossRef]
  14. Shen, S.; Pan, Z.; Zhang, J. Symmetries of a (2+1)-Dimensional Toda-like Lattice. Commun. Theor. Phys. 2004, 42, 805–806. [Google Scholar]
  15. Van Assche, W. Orthogonal polynomials, Toda lattices and Painlevé equations. Phys. D Nonlinear Phenom. 2022, 434, 133214. [Google Scholar] [CrossRef]
  16. Bluman, G.; Cole, J. The general similarity solution of the heat equation. J. Appl. Math. Mech. 1969, 18, 1025–1042. [Google Scholar]
  17. Arrigo, D.J.; Beckham, J.R. Nonclassical symmetries of evolutionary partial differential equations and compatibility. J. Math. Anal. Appl. 2004, 289, 55–65. [Google Scholar] [CrossRef]
  18. Murata, S. Nonclassical symmetry analysis for hyperbolic partial differential equation. Commun. Nonlinear Sci. Numer. Simul. 2008, 13, 1472–1474. [Google Scholar] [CrossRef]
  19. Bluman, G.; Tian, S.; Yang, Z. Nonclassical analysis of the nonlinear Kompaneets equation. J. Eng. Math. 2014, 84, 87–97. [Google Scholar] [CrossRef]
  20. Xin, X.; Jin, M.; Yang, J.; Xia, Y. Nonlocal symmetries and solutions of the (2+1) dimension integrable Burgers equation. Appl. Math. Lett. 2024, 148, 108867. [Google Scholar] [CrossRef]
  21. Arrigo, D.J.; Broadbridge, P.; Hill, J.M. Nonclassical symmetry solutions and the methods of Bluman–Cole and Clarkson–Kruskal. J. Math. Phys. 1993, 34, 4692–4703. [Google Scholar] [CrossRef]
  22. Subhankar, S.I.L.; Sekhar, T.R. Nonclassical symmetry analysis, conservation laws of one-dimensional macroscopic production model and evolution of nonlinear waves. J. Math. Anal. Appl. 2021, 497, 124847. [Google Scholar]
  23. Levi, D.; Winternitz, P.; Yamilov, R.I. Lie point symmetries of differential–difference equations. J. Phys. A Math. Theor. 2010, 43, 292002. [Google Scholar] [CrossRef]
Figure 1. The plot of u ( n ) (72).
Figure 1. The plot of u ( n ) (72).
Symmetry 15 02199 g001
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

Li, W.; Chen, Y.; Jiang, K. Non-Classical Symmetry Analysis of a Class of Nonlinear Lattice Equations. Symmetry 2023, 15, 2199. https://doi.org/10.3390/sym15122199

AMA Style

Li W, Chen Y, Jiang K. Non-Classical Symmetry Analysis of a Class of Nonlinear Lattice Equations. Symmetry. 2023; 15(12):2199. https://doi.org/10.3390/sym15122199

Chicago/Turabian Style

Li, Wenting, Yueting Chen, and Kun Jiang. 2023. "Non-Classical Symmetry Analysis of a Class of Nonlinear Lattice Equations" Symmetry 15, no. 12: 2199. https://doi.org/10.3390/sym15122199

APA Style

Li, W., Chen, Y., & Jiang, K. (2023). Non-Classical Symmetry Analysis of a Class of Nonlinear Lattice Equations. Symmetry, 15(12), 2199. https://doi.org/10.3390/sym15122199

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