Next Article in Journal
New Self-Adaptive Inertial-like Proximal Point Methods for the Split Common Null Point Problem
Next Article in Special Issue
A New Method for Blow-Up to Scale-Invariant Damped Wave Equations with Derivatives and Combined Nonlinear Terms
Previous Article in Journal
Method of Wavelet-Decomposition to Research Cosmic Ray Variations: Application in Space Weather
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Relativistic Toda Lattice Hierarchy, Discrete Generalized (m,2Nm)-Fold Darboux Transformation and Diverse Exact Solutions

1
School of Applied Science, Beijing Information Science and Technology University, Beijing 100192, China
2
Department of Mathematics and Information Technology, The Education University of Hong Kong, New Territories, Hong Kong, China
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(12), 2315; https://doi.org/10.3390/sym13122315
Submission received: 11 October 2021 / Revised: 14 November 2021 / Accepted: 19 November 2021 / Published: 3 December 2021

Abstract

:
This paper investigates a relativistic Toda lattice system with an arbitrary parameter that is a very remarkable generalization of the usual Toda lattice system, which may describe the motions of particles in lattices. Firstly, we study some integrable properties for this system such as Hamiltonian structures, Liouville integrability and conservation laws. Secondly, we construct a discrete generalized ( m , 2 N m ) -fold Darboux transformation based on its known Lax pair. Thirdly, we obtain some exact solutions including soliton, rational and semi-rational solutions with arbitrary controllable parameters and hybrid solutions by using the resulting Darboux transformation. Finally, in order to understand the properties of such solutions, we investigate the limit states of the diverse exact solutions by using graphic and asymptotic analysis. In particular, we discuss the asymptotic states of rational solutions and exponential-and-rational hybrid solutions graphically for the first time, which might be useful for understanding the motions of particles in lattices. Numerical simulations are used to discuss the dynamics of some soliton solutions. The results and properties provided in this paper may enrich the understanding of nonlinear lattice dynamics.

1. Introduction

In recent years, the study of nonlinear lattice equations (NLEs), viewed as spatially discrete counterparts of nonlinear partial differential equations, has attracted increasing interest. The NLEs may well model some physical phenomena in various fields such as nonlinear optics, lattice dynamics, electric circuits, population dynamics and so on [1,2,3,4,5,6,7,8,9]. The Toda lattice (TL) system is the first example of Lax integrable NLEs, which can describe a lattice of particles interacting with their nearest neighbors via forces exponentially dependent on distances [1,2,3], as shown in Figure 1. In the picture, the wavy lines represent springs, and small spheres represent particles, which are connected graphically.
TL system is also used in studying nonlinear waves in nonlinear lattice dynamics [3]. In order to better investigate the physical phenomena described by the TL system [10], proposed a new Hamiltonian function H = n = 1 N e α p n 1 + α p n α 2 + n = 1 N e x n x n 1 α p n for which its Hamilton’s motion equation is the following relativistic Toda lattice (RTL) equation (see also Equation (8.2.28) in Ref. [10]):
x n , t = α 1 ( 1 e α p n ) α e x n x n 1 α p n , p n , t = e x n + 1 x n α p n + 1 e x n x n 1 α p n ,
where x n = x ( n , t ) and p n = p ( n , t ) are the real functions of discrete variable n and time variable t. When α 0 , Equation (1) reduces to the famous TL system as a very remarkable generalization of TL system; thus, the study of Equation (1) is of great significance for understanding the physical phenomena described by the TL system. Equation (1) can also be considered as a one-parameter perturbation of TL system, and in a certain physical interpretation, the small parameter α has the meaning of the inverse speed of light [10,11]. The author has used the abbreviated name RTL ( α ) for Equation (1). For convenience, hereafter we still use RTL ( α ) for Equation (1). According to Ref. [10], the Lax pair of Equation (1) admits the following:
E ϕ n = U n ϕ n = z 2 e α p n 1 α z e x n α z e x n + α p n 0 ϕ n ,
ϕ n , t = V n ϕ n = 0 z e x n α p n z e x n 1 z 2 1 α + α e α p n x n 1 + x n ϕ n ,
where z is the spectral parameter independent of t, E is the shift operator defined by E f ( n , t ) = f ( n + 1 , t ) , E 1 f ( n , t ) = f ( n 1 , t ) and ϕ n = ( φ n , ψ n ) T is an eigenfunction vector. Equation (1) can be recovered from the integrability condition U n , t = ( E V n ) U n U n V n between the spatial part (2) and time evolution part (3), and they have been verified to be correct. It should be noted here that we have made some small changes to the Lax pair, (2) and (3), for the convenience of later discussion. We divide every element of matrix U n in the spatial part (2) by 1 z , which has no effect on the compatibility condition of Equation (1) because each term in the compatibility condition has U n . The Lax pairs (2) and (3) are inconvenient for constructing Darboux transformations (DTs); thus, we need to simplify them before performing this. Here, we take q n = e α p n , r n = e x n , λ = z , then the Lax pair (2) and (3) can be transformed into the following.
E ϕ n = U n ϕ n = λ 2 q n 1 α λ r n α λ q n r n 0 ϕ n ,
ϕ n , t = V n ϕ n = 0 λ r n q n λ r n 1 λ 2 1 α + α r n q n r n 1 ϕ n .
At this time, Equation (1) is equivalent to the following equation.
q n , t = α ( q n r n 1 r n + 1 q n + 1 r n 2 ) q n + 1 r n r n 1 , r n , t = r n ( α 2 r n q n r n 1 + r n 1 ) α q n r n 1 .
Below, we first investigate Equation (6), then we use transformations p n = ln q n α and x n = ln r n to study Equation (1) when q n > 0 , r n > 0 . There are some known research results on the RTL system such as Lax pair [12,13,14,15], Hamiltonian structures [12,16,17], DT and exact solutions [17,18,19,20,21] and so on. In Ref. [11], Suris proposed some discrete time Toda systems including the RTL ( α ) system (see System (8.2) in Ref. [11]). Although Equation (1) has the same abbreviation as system (8.2) in the literature [11], by careful comparison, we found that they are different from each other. Here, we need to point out that the RTL ( α ) system that we study here is different from the RTL equation in Refs. [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. To our knowledge, some integrable properties, various exact solutions via the discrete generalized ( m , 2 N m ) -fold DT and asymptotic analysis for Equation (6) or Equation (1) have not been studied before.
In soliton theory, searching for a new Lax integrable lattice hierarchy is an important subject [27,28] from which we not only find some new higher-order Lax integrable lattice systems but also so that we can investigate some relevant integrable properties such as Lax pairs, infinitely many conservation laws, Hamiltonian structure and Liouville integrability by means of the Tu scheme [29,30,31]. Furthermore, explicit exact solutions to NLEs play an essential role in understanding complex natural phenomena; thus, searching for exact solutions, particularly soliton solutions, is an important subject in the study of NLEs. Recently, a discrete generalized ( m , 2 N m ) -fold DT based on Lax pair was proposed to solve Lax integrable NLEs [17,32]. Compared with the usual N-fold DT, which only provides soliton solutions, this generalized method is a useful technique that not only can present soliton solutions but also provide rational and semi-rational solutions and hybrid solutions of exponential and rational functions [17,32].
Therefore, the aim of this paper can be summarized in the following two points: (i) constructing the lattice hierarchy associated with Equation (6) and then studying its relevant integrable properties including Hamiltonian structure, Liouville integrability and infinite conservation laws; (ii) establishing the discrete generalized ( m , 2 N m ) -fold DT for Equation (6) to provide various exact solutions such as discrete soliton solutions, rational and semi-rational solutions and hybrid solutions and then carrying out the asymptotic analysis technique in order to analyze the limit states of such exact solutions. It should be noted here that the properties mentioned above are studied for the first time for Equation (6) or (1).
The structure of this paper is as follows. Section 2 constructs relativistic Toda lattice hierarchies associated with Equation (6) and studies some integrable properties including Hamiltonian structures and Liouville integrability via the Tu scheme [29,30,31]. Section 3 constructs the discrete generalized ( m , 2 N m ) -fold DT of Equation (6). Various exact solutions such as multi-soliton solutions, rational and semi-rational solutions and soliton-and-rational hybrid solutions for Equation (6) are given by using the special cases of the resulting generalized DT in Section 4. Some conclusions are provided in Section 5.

2. The Integrable Properties of Equation (6)

In this section, we will study the integrable properties of Equation (6) including its lattice hierarchy, Hamiltonian structures, Liouville integrability and infinite conservation laws.

2.1. A Discrete Integrable Lattice Hierarchy Associated with Equation (6)

First, we will construct the lattice hierarchy associated with Equation (6) from the spectral problem (4) by using the Tu scheme technique [29,30,31]. Let the following be the case.
P n = A n B n C n A n = j = 0 A n ( j ) λ 2 j j = 0 B n ( j ) λ 2 j 1 j = 0 C n ( j ) λ 2 j 1 j = 0 A n ( j ) λ 2 j .
By solving the following stationary discrete zero-curvature equation P n + 1 U n U n P n = 0 , we have the following recursion relations:
q n ( A n + 1 ( j + 1 ) A n ( j + 1 ) ) A n + 1 ( j ) + A n ( j ) α q n r n B n + 1 ( j ) α r n C n ( j ) = 0 , α r n ( A n + 1 ( j + 1 ) + A n ( j + 1 ) ) q n B n ( j + 1 ) + B n ( j ) = 0 , α q n r n ( A n + 1 ( j + 1 ) + A n ( j + 1 ) ) C n + 1 ( j ) + q n C n + 1 ( j + 1 ) = 0 , α r n C n + 1 ( j ) + α q n r n B n ( j ) = 0 ,
where A n ( j ) , B n ( j ) and C n ( j ) are functions of q n , r n . Let A n ( 0 ) = 1 2 α be the initial values, the recursion relations (7) will determine the rest A n ( j ) , B n ( j ) , C n ( j ) uniquely and the first few coefficients are listed below.
B n ( 0 ) = r n q n , C n ( 0 ) = 1 r n 1 , A n ( 1 ) = α r n q n r n 1 + 1 α , B n ( 1 ) = 1 q n ( α 2 r n + 1 q n + 1 + α 2 r n 2 q n r n 1 + r n q n + 2 r n ) , C n ( 1 ) = 1 q n 1 r n 1 + α 2 r n q n r n 1 2 + α 2 q n 1 r n 2 2 r n 1 , A n ( 2 ) = α q n ( 2 r n r n 1 r n r n 1 q n 1 r n q n r n 1 α 2 r n 2 q n r n 1 2 α 2 r n q n 1 r n 2 α 2 r n + 1 q n + 1 r n 1 ) , B n ( 2 ) = 1 q n 2 ( α 2 r n + 1 q n + 1 + α 2 r n 2 q n r n 1 + r n q n ) + 2 r n q n 2 + α r n q n [ α q n + 1 ( 2 r n + 1 r n r n + 1 r n q n r n + 1 q n + 1 r n α 2 r n + 1 2 q n + 1 r n 2 α 2 r n + 1 q n r n 1 α 2 r n + 2 q n + 2 r n ) + α q n ( 2 r n r n 1 r n r n 1 q n 1 r n q n r n 1 α 2 r n 2 q n r n 1 2 α 2 r n q n 1 r n 2 α 2 r n + 1 q n + 1 r n 1 ) ] , C n ( 2 ) = 1 q n 1 ( 1 q n 1 r n 1 + α 2 r n q n r n 1 2 + α 2 q n 1 r n 2 2 r n 1 ) 1 r n 1 [ α q n ( 2 r n r n 1 r n r n 1 q n 1 r n q n r n 1 α 2 r n 2 q n r n 1 2 α 2 r n q n 1 r n 2 α 2 r n + 1 q n + 1 r n 1 ) + α q n 1 ( 2 r n 1 r n 2 r n 1 r n 2 q n 2 r n 1 q n 1 r n 2 α 2 r n 1 2 q n 1 r n 2 2 α 2 r n 1 q n 2 r n 3 α 2 r n q n r n 2 ) ] . . . .
For any integer m 0 , we define the following.
Γ n ( m ) = λ 2 m + 2 P n = j = 0 m A n ( j ) λ 2 j + 2 m + 2 j = 0 m B n ( j ) λ 2 j + 2 m + 1 j = 0 m C n ( j ) λ 2 j + 2 m + 1 j = 0 m A n ( j ) λ 2 j + 2 m + 2 , m 0 .
By means of the coefficients in (8), we obtain the following.
E Γ n ( m ) U n U n Γ n ( m ) = q n ( A n + 1 ( m + 1 ) A n ( m + 1 ) ) λ B n ( m ) λ C n + 1 ( m ) 0 .
From (4), we know that the following is the case.
U n , t = λ 2 q n , t α λ r n , t α λ ( q n r n , t r n q n , t ) r n 2 0 .
In order to render (9) and (10) compatible, we need to modify (9); thus, we here take a modification matrix:
Δ n ( m ) = 0 0 0 A n ( m + 1 ) ,
and define V n ( m ) = P n ( m ) + δ n such that the following is obtained.
E V n ( m ) U n U n V n ( m ) = q n ( A n + 1 ( m + 1 ) A n ( m + 1 ) ) λ B n ( m ) + α λ r n A n ( m + 1 ) λ C n + 1 ( m ) + α λ q n r n A n + 1 ( m + 1 ) 0 .
Assume that the time evolution of ϕ n meets equations ϕ n , t m = V n ( m ) ϕ n for which its compatibility condition with (4) is as follows:
U n , t m = ( E V n ( m ) ) U n U n V n ( m ) , m 0 ,
which yields the following discrete integrable lattice hierarchy.
q n , t m = q n ( A n + 1 ( m + 1 ) A n ( m + 1 ) ) , r n , t m = B n ( m ) + α r n A n ( m + 1 ) α .
(i)
When m = 0 , Equation (14) reduces to (6) with t 0 = t as follows.
q n , t 0 = q n ( A n + 1 ( 1 ) A n ( 1 ) ) = α ( q n r n + 1 r n + 1 q n + 1 r n 2 ) q n + 1 r n r n 1 , r n , t 0 = B n ( 0 ) + α r n A n ( 1 ) α = r n ( α 2 r n q n r n 1 + r n 1 ) α q n r n 1 .
At this time, its time part of the Lax pair is as follows.
ϕ n , t 0 = V n ( 0 ) ϕ n = λ 2 2 α λ r n q n λ r n 1 λ 2 2 α 1 α + α r n r n 1 q n ϕ n .
(ii)
When m = 1 , Equation (14) reduces to a new higher-order RTL ( α ) system.
q n , t 1 = q n ( A n + 1 ( 2 ) A n ( 2 ) ) = q n [ α q n ( 2 r n r n 1 r n r n 1 q n 1 r n q n r n 1 α 2 r n 2 q n r n 1 2 α 2 r n q n 1 r n 2 α 2 r n + 1 q n + 1 r n 1 ) α q n + 1 ( 2 r n + 1 r n r n + 1 r n q n r n + 1 q n + 1 r n α 2 r n + 1 2 q n + 1 r n 2 α 2 r n + 1 q n r n 1 α 2 r n + 2 q n + 2 r n ) ] , r n , t 1 = B n ( 1 ) + α r n A n ( 2 ) α = 1 α q n ( α 2 r n + 1 q n + 1 + α 2 r n 2 q n r n 1 + r n q n ) + 2 r n α q n + r n [ α q n ( 2 r n r n 1 r n r n 1 q n 1 r n q n r n 1 α 2 r n 2 q n r n 1 2 α 2 r n q n 1 r n 2 α 2 r n + 1 q n + 1 r n 1 ) ] ,
Moreover, its time part of the Lax pair is as follows:
ϕ n , t 1 = V n ( 1 ) ϕ n = λ 4 A n ( 0 ) + λ 2 A n ( 1 ) λ 3 B n ( 0 ) + λ B n ( 1 ) λ 3 C n ( 0 ) + λ C n ( 1 ) λ 4 A n ( 0 ) λ 2 A n ( 1 ) 1 α A n ( 1 ) ϕ n V 11 ( n ) V 12 ( n ) V 21 ( n ) V 22 ( n ) ϕ n ,
with the following case.
V 11 ( n ) = λ 4 2 α α λ 2 r n q n r n 1 + λ 2 α , V 12 ( n ) = λ 3 r n q n λ q n ( α 2 r n + 1 q n + 1 + α 2 r n 2 q n r n 1 + r n q n + 2 r n ) , V 21 ( n ) = λ 3 r n 1 + λ ( 1 q n 1 r n 1 + α 2 r n q n r n 1 2 + α 2 q n 1 r n 2 2 r n 1 ) , V 22 ( n ) = λ 4 2 α + α λ 2 r n q n r n 1 λ 2 α α q n ( 2 r n r n 1 r n r n 1 q n 1 r n q n r n 1 α 2 r n 2 q n r n 1 2 α 2 r n q n 1 r n 2 α 2 r n + 1 q n + 1 r n 1 ) .
Remark 1.
It should be pointed out that the new higher-order RTL ( α ) system (17) has more nonlinear terms than Equation (6) and may have some novel properties that are worthy of further study. The main purpose of this paper is to study system (6); thus, we do not specifically study system (17) here.

2.2. Discrete Hamiltonian Structures

Next, we will construct the Hamiltonian structures of Equation (14) by using the trace identity technique [29,30,31]. Let < U , V > = t r ( U V ) , where U and V are arbitrary square matrices, and t r means the trace of a matrix [29,30,31].
By defining the following:
V n = P n U n 1 = A n B n C n A n 0 r n α λ q n 1 α λ r n λ 2 q n 1 α 2 λ 2 q n = B n α λ r n r n A n α λ q n + ( λ 2 q n 1 ) B n α 2 λ 2 q n A n α λ r n r n C n α λ q n ( λ 2 q n 1 ) A n α 2 λ 2 q n ,
we obtain the following case.
< V n , U n λ > = q n B n α r n + B n α λ 2 r n , < V n , U n q n > = A n q n + B n α λ q n r n , < V n , U n r n > = 2 A n r n + λ q n B n α r n 2 B n α λ r n 2 .
Substituting them into the trace identity [29,30,31]:
δ δ u n Z < V n , U n λ > = ( λ ε λ λ ε ) < V n , U n u i > , i = 1 , 2 ,
results in the following case.
δ δ u n Z ( q n B n α r n + B n α λ 2 r n ) = ( λ ε λ λ ε ) A n q n + B n α λ q n r n 2 A n r n + λ q n B n α r n 2 B n α λ r n 2 .
Comparing the coefficients of λ 2 m 3 on both sides of Equation (22), we obtain the following.
δ δ q n δ δ r n n Z ( q n B n ( m + 1 ) + B n ( m ) α r n ) = ( ε 2 m 2 ) A n ( m + 1 ) q n + B n ( m ) α q n r n q n B n ( m + 1 ) α r n 2 B n ( m ) α r n 2 2 A n ( m + 1 ) r n .
Setting m = 0 in the above equation, we can fix ε = 0 . Let H n ( m ) = n Z q n B n ( m + 1 ) + B n ( m ) 2 α ( m 1 ) r n , thus the following is obtained.
δ H n ( m ) δ u = A n ( m + 1 ) q n + B n ( m ) α q n r n q n B n ( m + 1 ) α r n 2 B n ( m ) α r n 2 2 A n ( m + 1 ) r n f n ( m + 1 ) g n ( m + 1 ) .
From the previous recursion relations, we can obtain A n ( m ) = ( E 1 ) 1 r n g n ( m ) , B n ( m ) = α r n q n f n ( m ) + α r n ( E 1 ) 1 r n + 1 g n + 1 ( m ) q n , C n ( m ) = α q n 1 f n 1 ( m ) + α ( E 1 ) 1 r n g n ( m ) r n 1 ; thus, Equation (14) can be rewritten as the following Hamiltonian form:
U t m = q n , t m r n , t m = J δ H n ( m ) δ u = J f n ( m + 1 ) g n ( m + 1 ) ,
with the following case.
J = 0 q n r n q n r n 0 .
Taking η = η 11 η 12 η 21 η 22 to satisfy δ H n ( m ) δ u = η δ H n ( m 1 ) δ u and by still using the recursion relations (7), we obtain the following.
η 11 = 1 q n [ ( E 1 ) 1 ( α 2 r n + 1 E r n α 2 r n q n 1 E 1 q n r n 1 ) + 1 ] , η 12 = 1 q n ( E 1 ) 1 ( α 2 r n + 1 ( E 1 ) 1 r n + 2 E 2 r n q n + 1 α 2 r n ( E 1 ) 1 r n q n r n 1 + r n q n ) + 1 q n 2 ( E 1 ) 1 r n + 1 E η 21 = α 2 q n 1 E 1 q n r n 1 + α 2 r n + 1 E r n 2 , η 22 = α 2 ( E 1 ) 1 r n q n r n 1 + α 2 r n + 1 ( E 1 ) 1 r n + 2 E 2 q n + 1 r n 2 + 1 q n .
Hence, Equation (14) can be rewritten as the recursive Hamiltonian form [29,30,31] below.
U t m = J δ H n ( m ) δ u = J f n ( m + 1 ) g n ( m + 1 ) = J η f n ( m ) g n ( m ) = = J η m f n ( 0 ) g n ( 0 ) = J η m + 1 1 α q n 0 .
In this manner, the lattice hierarchy (14) can be written into the above Hamiltonian form (27) successfully. It can be verified that J and J η are skew-symmetric operators, and the Hamiltonian functions H n ( m ) ( m 0 ) described by Equation (27) are pairwise involutory concerning the Poisson bracket. Here, we will not provide a detailed derivation of their involution, and the readers can refer to Refs. [29,30,31] for the detailed process. Therefore, we can say that the lattice hierarchy (14) or the Hamiltonian form (27) has Liouville integrability.

2.3. An Infinite Number of Conservation Laws of Equation (6)

Existing infinite conservation laws in an NLE play an important role in further verification of its integrability [33,34,35]. Therefore, the study of infinitely many conservation laws for NLE has important theoretical significance and practical values. In this subsection, we will derive the infinitely many conservation laws of Equation (6).
Setting θ n = ψ n φ n , from the space part of Lax pair (4) of Equation (6), we have the following.
φ n + 1 φ n = λ 2 q n 1 + α λ r n θ n , ψ n + 1 ψ n = α λ q n r n θ n .
Let the above two expressions be divided by each other, and we obtain the following.
( λ 2 q n 1 + α λ r n θ n ) θ n + 1 + α λ q n r n = 0 .
Inserting θ n = j = 0 n θ n ( j ) λ j into (29) and collecting the coefficients of same powers of λ yield the following.
θ n ( 0 ) = 0 , θ n ( 1 ) = α r n 1 , θ n ( 2 ) = 0 , θ n ( 3 ) = α q n 1 ( 1 r n 1 + α 2 r n 2 ) , θ n ( 4 ) = 0 , θ n ( 5 ) = α q n 1 2 r n 1 α 3 q n 1 2 r n 2 α 3 r n 1 q n 1 ( 1 r n 2 r n 1 q n 1 + α 2 r n 2 2 q n 1 + 1 r n 2 r n 1 q n 2 + α 2 r n 3 r n 1 q n 2 ) , , θ n ( 2 j ) = 0 , θ n ( 2 j + 1 ) = θ n ( 2 j + 1 ) α r n 1 i = 0 2 j θ n 1 ( i ) θ n ( 2 j i ) q n 1 , j 3 .
Combining the time part (16) of the Lax pair, we can obtain the following.
[ ln ( λ 2 q n + α λ r n θ n 1 ) ] t = ( E 1 ) ( λ 2 2 α λ r n q n θ n ) .
Substitution of (30) into (31) and comparing the same powers of λ on both sides of (31) can provide an infinite number of conservation laws for Equation (6). The first three conservation laws are listed below:
( T k ) t = ( E 1 ) ( X k ) , ( k = 1 , 2 , 3 ) ,
with the following being the case:
T 1 = ln q n , X 1 = α r n q n r n 1 , T 2 = 1 q n ( α 2 r n r n 1 1 ) , X 2 = α r n q n 1 q n ( 1 r n 1 + α 2 r n 2 ) , T 3 = α 2 r n q n 1 q n ( 1 r n 1 + α 2 r n 2 ) 1 2 q n 2 ( α 4 r n 2 r n 1 2 + 2 α 2 r n r n 1 + 1 ) , X 3 = α r n q n 1 2 q n r n 1 + α 3 r n q n 1 2 q n r n 2 + α 3 r n 1 r n q n 1 q n ( 1 r n 2 r n 1 q n 1 + α 2 r n 2 2 q n 1 + 1 r n 2 r n 1 q n 2 + α 2 r n 3 r n 1 q n 2 ) ,
where T k denotes the conserved densities, and X k denotes associated fluxes. The above three conservation laws usually represent energy conservation, momentum conservation and Hamiltonian conservation, respectively [33]. Existing infinitely many conservation laws means that Equation (6) is a discrete integrable system.

3. Discrete Generalized ( m , 2 N m ) -Fold DT

In the previous section, we have discussed the integrability of Equation (6). In this part, we will construct the discrete generalized ( m , 2 N m ) -fold DT of Equation (6). For this purpose, we consider the following gauge transformation:
ϕ ˜ n = T n ϕ n ,
which can transform the Lax pair (4) and (16) into the following forms.
ϕ ˜ n + 1 = U ˜ n ϕ ˜ n = T n + 1 U n T n 1 ϕ ˜ n , ϕ ˜ n , t = V ˜ n ( 0 ) ϕ ˜ n = ( T n , t + V n ( 0 ) T n ) T n 1 ϕ ˜ n .
In order to construct DT, we require that U ˜ n , V ˜ n ( 0 ) have the same forms as U n , V n ( 0 ) in addition to replacing the new potential functions q ˜ n , r ˜ n with the old potential functions q n , r n . In order to achieve this special purpose, we need to construct a particular Darboux matrix T n in the following form:
T n = a n b n c n d n = 1 + j = 1 N a n ( 2 N 2 j ) λ 2 j j = 1 N b n ( 2 j 1 ) λ 2 j 1 j = 1 N c n ( 2 j 1 ) λ 2 j 1 λ 2 N + j = 0 N 1 d n ( 2 N 2 j ) λ 2 j ,
in which the number N is a positive integer representing the order of DT, and a n ( 2 N 2 j ) , b n ( 2 j 1 ) , c n ( 2 j 1 ) and d n ( 2 N 2 j ) ( j = 1 , 2 , N ) are functions of the variables n , t that are determined later. According to (34), we can obtain the following transformations between the new potentials and old ones as follows.
q ˜ n = a n + 1 ( 0 ) q n a n ( 0 ) , r ˜ n = α r n + b n ( 2 N 1 ) α d n ( 2 N ) .
In order to provide the new solutions of Equation (6), the unknown functions a n ( 0 ) , b n ( 2 N 1 ) , d n ( 2 N ) must be determined.
Next, we will determine their concrete expressions by constructing them directly. According to Darboux matrix T n , we know that det T n is a 4 N -th order polynomial of λ , which should have 4 N roots. Let us assume that ϕ i , n ( λ i ) = ( φ i , n ( λ i , n ) , ψ i , n ( λ i ) ) T ( φ i , n , ψ i , n ) T are n solutions of Lax pair (4) and (16) for n different spectral parameters λ i ( k = 1 , 2 , , n ) . Moreover, we also want them to be n roots of det T n . It should be noted that this can easily meet the conditions. In fact, we only need to solve Lax pairs at n known roots of det   T n . Moreover, for these roots, we must guarantee that T ( λ i ) ϕ i , n ( λ i ) = 0 ( i = 1 , 2 , , m ) , which will only produce 2 m algebraic equations such that we cannot determine the 4 N unknown functions a n ( 2 N 2 j ) , b n ( 2 j 1 ) , c n ( 2 j 1 ) ,   d n ( 2 N 2 j ) when m < N .
In order to derive a linear algebraic system of 4 N equations with 4 N unknown functions a n ( 2 N 2 j ) , b n ( 2 j 1 ) , c n ( 2 j 1 ) ,   d n ( 2 N 2 j ) , for every λ i , we need to expand the following:
T ( λ i + ε ) ϕ i , n ( λ i + ε ) = K = 0 N 1 j = 0 k T ( j ) ( λ i ) ϕ i , n ( k j ) ( λ i ) ε k ,
where ε is a small parameter, φ n ( k ) ( λ i ) = 1 k ! k λ i k φ n ( λ i ) is given by the Taylor series expansion of ϕ i , n ( λ i + ε ) = ϕ i , n ( 0 ) ( λ i ) + ϕ i , n ( 1 ) ( λ i ) ε + ϕ i , n ( 2 ) ( λ i ) ε 2 + ϕ i , n ( 3 ) ( λ i ) ε 3 + around ε = 0 , T n ( j ) is given by a binomial expansion of T ( λ i + ε ) = T n ( 0 ) + T n ( 1 ) ε + + T n ( m i ) ε m i . Thus, the formula lim ε 0 T n ( λ i + ε ) ϕ n ( λ i + ε ) ε k i = 0 ( i = 1 , 2 , , m , k i = 0 , 1 , , v i , 2 N = m + i = 1 m v i ) can produce 4 N algebraic equations for 4 N unknown functions a n ( 2 N 2 j ) , b n ( 2 j 1 ) , c n ( 2 j 1 ) ,   d n ( 2 N 2 j ) . In other words, we have the following:
T ( 0 ) ( λ i ) ϕ i , n ( 0 ) ( λ i ) = 0 , T ( 0 ) ( λ i ) ϕ i , n ( 1 ) ( λ i ) + T ( 1 ) ( λ i ) ϕ i , n ( 0 ) ( λ i ) = 0 , T n ( 0 ) ( λ i ) ϕ i , n ( 2 ) ( λ i ) + T n ( 1 ) ( λ i ) ϕ i , n ( 1 ) ( λ i ) + T n ( 2 ) ( λ i ) ϕ i , n ( 0 ) ( λ i ) = 0 , , j = 0 v i T ( j ) ( λ i ) ϕ i , n ( v i j ) ( λ i ) = 0 , .
from which the determinant of the coefficients for system (4) is nonzero when the m spectral parameters λ i are suitably chosen so that a n ( 2 N 2 j ) , b n ( 2 j 1 ) , c n ( 2 j 1 ) ,   d n ( 2 N 2 j ) in the Darboux matrix T n are uniquely determined by (38). Thus, we can provide the undetermined functions a n ( 0 ) , b n ( 2 N 1 ) , d n ( 2 N ) in (36) as follows:
a n ( 0 ) = Δ a n ( 0 ) Δ 1 , n , b n ( 2 N 1 ) = Δ b n ( 2 N 1 ) Δ 1 , n , d n ( 2 N ) = Δ d n ( 2 N ) Δ 2 , n
with
Δ 1 = ( Δ 1 ( 1 ) , Δ 1 ( 2 ) , . . . , Δ 1 ( m ) ) T ,   Δ 2 = ( Δ 2 ( 1 ) , Δ 2 ( 2 ) , . . . , Δ 2 ( m ) ) T ,   Δ 1 ( i ) = ( Δ 1 , j , s ( i ) ) 2 ( v i + 1 ) × 2 N ,   Δ 2 ( i ) = ( Δ 2 , j , s ( i ) ) 2 ( v i + 1 ) × 2 N , in which Δ 1 , j , s ( i ) , Δ 2 , j , s ( i ) ( 1 j 2 ( v i + 1 ) , 1 s 2 N , i = 1 , 2 , , m ) given as follows:
Δ 1 , j , s ( i ) = k = 0 j 1 C 2 N 2 s + 2 k λ i 2 N 2 s k + 2 φ i , n ( j 1 k ) for l + i = 1 l 1 v i j l + i = 1 l ( 1 l m ) , 1 s N , k = 0 j 1 C 4 N 2 s + 1 k λ i 4 N 2 s k + 1 ψ i , n ( j 1 k ) for l + i = 1 l 1 v i j l + i = 1 l ( 1 l m ) , N + 1 s 2 N , Δ 2 , j , s ( i ) = k = 0 j 1 C 2 N 2 s + 1 k λ i 2 N 2 s k + 1 φ i , n ( j 1 k ) for l + i = 1 l 1 v i j l + i = 1 l ( 1 l m ) , 1 s N , k = 0 j 1 C 4 N 2 s + 2 k λ i 4 N 2 s k + 2 ψ i , n ( j 1 k ) for l + i = 1 l 1 v i j l + i = 1 l ( 1 l m ) , N + 1 s 2 N ,
where Δ a n ( 0 ) and b n ( 2 N 1 ) are given from the determinant Δ 1 by replacing their first and last columns by the column vector ( f 1 ( 1 ) , f 2 ( 1 ) , . . . , f ( v 1 + 1 ) ( 1 ) , . . . , f 1 ( i ) , f 2 ( i ) , . . . , f ( v i + 1 ) ( i ) , . . . , f 1 ( m ) , f 2 ( m ) , . . . , f ( v m + 1 ) ( m ) ) with f j ( i ) = φ i , n ( j 1 ) ( 1 j ( v i + 1 ) , 1 i m ) , respectively, while Δ d n ( 2 N ) is obtained from the determinant Δ 2 by replacing the last columns by the column vector ( r 1 ( 1 ) , r 2 ( 1 ) , . . . , r ( v 1 + 1 ) ( 1 ) , . . . , r 1 ( i ) , r 2 ( i ) , . . . , r ( v i + 1 ) ( i ) , . . . , r 1 ( m ) , r 2 ( m ) , . . . , r ( v m + 1 ) ( m ) ) with r j ( i ) = λ i 2 N ψ i , n ( j 1 ) ( 1 j ( v i + 1 ) , 1 i m ) .
Different from the previous usual N-fold DT using N spectral parameters, we call the transformations (33) and (36) using m spectral parameters as the generalized ( m , 2 N m ) -fold DT. Here, the following notation is used: m in the name of the generalized ( m , 2 N m ) -fold DT stands for the number of spectral parameters used, while 2 N m represents the sum of the used Taylor expansion order number of the eigenfunction ϕ i , n .
Remark 2.
When m = 2 N and m i = 0 , the transformations (33) and (36) reduce to the discrete generalized ( 2 N , 0 ) -fold DT which includes the discrete 2 N -fold DT if we do not use the Taylor series expansion for every ϕ i , n ( λ i ) . If m = 1 and m i = 2 N 1 , the transformations (33) and (36) reduce to the discrete generalized ( 1 , 2 N 1 ) -fold DT, which can provide higher-order rational and semi-rational solutions of Equation (6). If m = 2 and m i = 2 N 2 , the transformations (33) and (36) reduce to the discrete generalized ( 2 , 2 N 2 ) -fold DT, which are hybrid solutions of soliton and rational or semi-rational solutions. If 2 < m < 2 N , the transformations (33) and (36) can also provide some other new discrete mixed solutions that will not be discussed here.

4. Explicit Exact Solutions and Their Asymptotic Analysis

In this section, we will apply the special cases of the discrete generalized ( m , 2 N m ) -fold DT to provide some explicit exact solutions including soliton, rational and semi-rational solutions and their hybrid solutions. Then, we will analyze the limit states of such solutions by asymptotic analysis technique.

4.1. Multi-Soliton Solutions and Dynamics

When m = 2 N , the discrete generalized ( m , 2 N m ) -fold DT reduces to the discrete generalized ( 2 N , 0 ) -fold DT. If we do not make a Taylor series expansion for every eigenfunction ϕ i , n ( λ i ) ( i = 1 , 2 , . . . , 2 N ) , at this time the discrete generalized ( 2 N , 0 ) -fold DT is just the usual 2 N -fold DT, which can provide multi-soliton solutions. Substitution of the initial seed solutions q n = 1 + α 2 , r n = 1 into (4) and (5) yields one basic solution with λ = λ i ( i = 1 , 2 , . . . , 2 N ) as follows:
ϕ i , n = φ i , n ψ i , n = C i 1 τ i 1 n e ρ i 1 t + C i 2 τ i 2 n e ρ i 2 t α λ i ( α 2 + 1 ) C i 1 τ i 1 n 1 e ρ i 1 t α λ i ( α 2 + 1 ) C i 2 τ i 2 n 1 e ρ i 2 t
with
τ i 1 = 1 2 ( α 2 λ i 2 + λ i 2 1 + α 4 λ i 4 4 α 4 λ i 2 + 2 α 2 λ i 4 6 α 2 λ i 2 + λ i 4 2 λ i 2 + 1 ) , τ i 2 = 1 2 ( α 2 λ i 2 + λ i 2 1 α 4 λ i 4 4 α 4 λ i 2 + 2 α 2 λ i 4 6 α 2 λ i 2 + λ i 4 2 λ i 2 + 1 ) , ρ i 1 = 1 α ( α 2 + 1 ) ( α 2 λ i 2 λ i 2 + 1 + α 4 λ i 4 4 α 4 λ i 2 + 2 α 2 λ i 4 6 α 2 λ i 2 + λ i 4 2 λ i 2 + 1 ) , ρ i 2 = 1 α ( α 2 + 1 ) ( α 2 λ i 2 λ i 2 + 1 α 4 λ i 4 4 α 4 λ i 2 + 2 α 2 λ i 4 6 α 2 λ i 2 + λ i 4 2 λ i 2 + 1 ) .
According to (36), we can obtain exact 2 N -soliton solutions of Equation (6). By direct calculation, we find that even soliton solutions of higher orders can degenerate into odd soliton solutions of lower order when we take one of λ i equal to 1 + α α 2 + 1 , 1 α α 2 + 1 , 1 + α α 2 + 1 , 1 α α 2 + 1 . Hereafter, in the degenerate case, we uniformly choose λ = 1 + α α 2 + 1 . Next, as an example, we discuss the two cases: N = 1 , 2 .
Case (1) When N = 1 , λ = λ i ( i = 1 , 2 ) based on the 2-fold DT from (36), we can provide the two-fold exact solutions as follows:
q ˜ n = ( 1 + α 2 ) a n + 1 ( 0 ) a n ( 0 ) , r ˜ n = α + b n ( 1 ) α d n ( 2 ) ,
where a n ( 0 ) = Δ a n ( 0 ) Δ 1 , b n ( 1 ) = Δ b n ( 1 ) Δ 1 and d n ( 2 ) = Δ d n ( 2 ) Δ 2 , in which the following is the case.
Δ 1 = λ 1 2 λ 1 δ 1 , n λ 2 2 λ 2 δ 2 , n , Δ 2 = λ 1 δ 1 , n λ 2 δ 2 , n , Δ a n ( 0 ) = 1 λ 1 δ 1 , n 1 λ 2 δ 2 , n , Δ b n ( 1 ) = λ 1 2 1 λ 2 2 1 , Δ d n ( 2 ) = λ 1 λ 1 2 δ 1 , n λ 2 λ 2 2 δ 2 , n .
For the sake of analysis, the analytical expressions of solution (36) can be rewritten as follows:
q ˜ n = ( 1 + α 2 ) [ λ 2 cosh ( ξ 1 + X 1 ) cosh ξ 2 λ 1 cosh ξ 1 cosh ( ξ 2 + X 2 ) ] [ λ 2 cosh ( ξ 1 X 1 ) cosh ξ 2 λ 1 cosh ξ 1 cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ξ 1 cosh ( ξ 2 + X 2 ) λ 1 cosh ( ξ 1 + X 1 ) cosh ξ 2 ] [ λ 2 cosh ξ 1 cosh ( ξ 2 X 2 ) λ 1 cosh ( ξ 1 X 1 ) cosh ξ 2 ] , r ˜ n = R 1 α 1 + α 2 λ 1 2 λ 2 2 [ λ 2 cosh ( ξ 1 X 1 ) cosh ξ 2 λ 1 cosh ξ 1 cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ξ 1 cosh ( ξ 2 X 2 ) λ 1 cosh ( ξ 1 X 1 ) cosh ξ 2 ] ,
in which the following is the case.
R 1 = [ α 1 + α 2 λ 1 2 λ 2 cosh ( ξ 1 X 1 ) cosh ξ 2 α 1 + α 2 λ 1 λ 2 2 cosh ξ 1 cosh ( ξ 2 X 2 ) + ( λ 2 2 λ 1 2 ) cosh ξ 1 cosh ξ 2 ] [ λ 1 cosh ( ξ 1 X 1 ) cosh ξ 2 λ 2 cosh ξ 1 cosh ( ξ 2 X 2 ) ] , ξ i = 1 2 ( ln τ i 1 ln τ i 2 ) n + 1 2 ( ρ i 1 ρ i 2 ) t + 1 2 ( ln C i 1 ln C i 2 ) , X i = 1 2 ( ln τ i 1 ln τ i 2 ) , i = 1 , 2 .
The parameters are suitably chosen so that the solutions (41) may be one-soliton or two-soliton solutions, for which their evolution plots are indicated in Figure 2 and Figure 3. When both λ 1 and λ 2 are not equal to 1 + α α 2 + 1 , the solutions (41) are two-soliton ones. It is worth mentioning that when one of the λ 1 , λ 2 is 1 + α α 2 + 1 , without loss of generality, we set λ 1 = 1 + α α 2 + 1 = A , and the solutions (41) reduce to one-soliton ones of the following form:
q ˜ n = ( 1 + α 2 ) [ λ 2 cosh ξ 2 A cosh ( ξ 2 + X 2 ) ] [ λ 2 cosh ξ 2 A cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ( ξ 2 + X 2 ) A cosh ξ 2 ] [ λ 2 cosh ( ξ 2 X 2 ) A cosh ξ 2 ] , r ˜ n = [ α 1 + α 2 λ 2 A 2 cosh ξ 2 α 1 + α 2 λ 2 2 A cosh ( ξ 2 X 2 ) + ( λ 2 2 A 2 ) cosh ξ 2 ] [ A cosh ξ 2 λ 2 cosh ( ξ 2 X 2 ) ] α λ 2 2 A 2 1 + α 2 [ λ 2 cosh ξ 2 A cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ( ξ 2 X 2 ) A cosh ξ 2 ] ,
for which their corresponding evolution plots are shown in Figure 2. Figure 2a1,a2 show dark one-soliton structures of the component q ˜ n under nonzero seed background (i.e., q ˜ n = 1.6 ). Figure 2b1,b2 show the kink-shaped one-soliton propagation structure for solution r ˜ n . We know that p ˜ n = ln | q ˜ n | α , x ˜ n = ln | r ˜ n | are the solutions of Equation (1), and Figure 2c1,c2 show dark one-soliton propagation structures on nonzero background. Figure 2d1,d2 show the kink-shaped one-soliton structures for solution x ˜ n . From Figure 2, we can clearly observe that one soliton retains the same amplitude and shape during its propagation.
When neither of the two parameters λ 1 and λ 2 are 1 + α α 2 + 1 , the solutions (41) are two-soliton ones. In order to further study the elastic interaction properties of two-soliton solutions q ˜ n , r ˜ n in (41), we perform an asymptotic analysis to investigate their limit states, which are listed as follows.
Before collision ( t ) , the following cases are observed:
(i)
If ξ 1 is fixed , ξ 2 + , we can calculate the initial limit states of the first soliton as follows;
q ˜ n q n 1 = ( 1 + α 2 ) [ λ 2 cosh ( ξ 1 + X 1 ) λ 1 e X 2 cosh ξ 1 ] [ λ 2 cosh ( ξ 1 X 1 ) λ 1 e X 2 cosh ξ 1 ] [ λ 2 e X 2 cosh ξ 1 λ 1 cosh ( ξ 1 + X 1 ) ] [ λ 2 e X 2 cosh ξ 1 λ 1 cosh ( ξ 1 X 1 ) ] , r ˜ n r n 1 = [ α 1 + α 2 λ 1 2 λ 2 cosh ( ξ 1 X 1 ) α 1 + α 2 λ 1 λ 2 2 e X 2 cosh ξ 1 + ( λ 2 2 λ 1 2 ) cosh ξ 1 ] [ λ 1 cosh ( ξ 1 X 1 ) λ 2 cosh ξ 1 e X 2 ) ] α 1 + α 2 λ 1 2 λ 2 2 [ λ 2 cosh ( ξ 1 X 1 ) λ 1 cosh ξ 1 e X 2 ] [ λ 2 cosh ξ 1 e X 2 λ 1 cosh ( ξ 1 X 1 ) ] .
(ii)
If ξ 2 is fixed , ξ 1 + , we can calculate the initial limit states of the second soliton as follows.
q ˜ n q n 2 = ( 1 + α 2 ) [ λ 2 e X 1 cosh ξ 2 λ 1 cosh ( ξ 2 + X 2 ) ] [ λ 2 e X 1 cosh ξ 2 λ 1 cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ( ξ 2 + X 2 ) λ 1 e X 1 cosh ξ 2 ] [ λ 2 cosh ( ξ 2 X 2 ) λ 1 e X 1 cosh ξ 2 ] , r ˜ n r n 2 = [ α 1 + α 2 λ 1 2 λ 2 e X 1 cosh ξ 2 α 1 + α 2 λ 1 λ 2 2 cosh ( ξ 2 X 2 ) + ( λ 2 2 λ 1 2 ) cosh ξ 2 ] [ λ 1 e X 1 cosh ξ 2 λ 2 cosh ( ξ 2 X 2 ) ] α 1 + α 2 λ 1 2 λ 2 2 [ λ 2 e X 1 cosh ξ 2 λ 1 cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ( ξ 2 X 2 ) λ 1 e X 1 cosh ξ 2 ] .
After collision ( t + ) , we observe the following:
(iii)
If ξ 1 is fixed , ξ 2 , we can calculate the end limit states of the first soliton as follows;
q ˜ n q n 1 = ( 1 + α 2 ) [ λ 2 cosh ( ξ 1 + X 1 ) λ 1 e X 2 cosh ξ 1 ] [ λ 2 cosh ( ξ 1 X 1 ) λ 1 e X 2 cosh ξ 1 ] [ λ 2 e X 2 cosh ξ 1 λ 1 cosh ( ξ 1 + X 1 ) ] [ λ 2 e X 2 cosh ξ 1 λ 1 cosh ( ξ 1 X 1 ) ] , r ˜ n r n 1 = [ α 1 + α 2 λ 1 2 λ 2 cosh ( ξ 1 X 1 ) α 1 + α 2 λ 1 λ 2 2 e X 2 cosh ξ 1 + ( λ 2 2 λ 1 2 ) cosh ξ 1 ] [ λ 1 cosh ( ξ 1 X 1 ) λ 2 cosh ξ 1 e X 2 ) ] α 1 + α 2 λ 1 2 λ 2 2 [ λ 2 cosh ( ξ 1 X 1 ) λ 1 cosh ξ 1 e X 2 ] [ λ 2 cosh ξ 1 e X 2 λ 1 cosh ( ξ 1 X 1 ) ] .
(iv)
If ξ 2 is fixed , ξ 1 , we can calculate the end limit states of the second soliton as follows:
q ˜ n q n 2 = ( 1 + α 2 ) [ λ 2 e X 1 cosh ξ 2 λ 1 cosh ( ξ 2 + X 2 ) ] [ λ 2 e X 1 cosh ξ 2 λ 1 cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ( ξ 2 + X 2 ) λ 1 e X 1 cosh ξ 2 ] [ λ 2 cosh ( ξ 2 X 2 ) λ 1 e X 1 cosh ξ 2 ] , r ˜ n r n 2 = [ α 1 + α 2 λ 1 2 λ 2 e X 1 cosh ξ 2 α 1 + α 2 λ 1 λ 2 2 cosh ( ξ 2 X 2 ) + ( λ 2 2 λ 1 2 ) cosh ξ 2 ] [ λ 1 e X 1 cosh ξ 2 λ 2 cosh ( ξ 2 X 2 ) ] α 1 + α 2 λ 1 2 λ 2 2 [ λ 2 e X 1 cosh ξ 2 λ 1 cosh ( ξ 2 X 2 ) ] [ λ 2 cosh ( ξ 2 X 2 ) λ 1 e X 1 cosh ξ 2 ] ,
where q n 1 , q n 2 , r n 1 , r n 2 stand for the asymptotic state expressions of q ˜ n and r ˜ n , the ‘−’ sign indicates the limit states before collision and the ‘+’ sign denotes the limit states after collision.
According to the above analysis, we know that the interactions between two solitons for q ˜ n and r ˜ n are elastic. When the parameters are given by λ 1 = 1 4 , λ 2 = 3 , α = 1 , C 11 = C 12 = C 21 = 1 , C 22 = 1 2 , the corresponding evolution structure plots are shown in Figure 3. Figure 3a1,a2 demonstrate the head-on elastic interactions between one dark soliton and one anti-bell-shaped soliton of the component q ˜ n on nonzero constant background for q ˜ n . Figure 3b1,b2 display the head-on inelastic interactions between two kink solitons for solution r ˜ n . Figure 3c1,c2 demonstrate the head-on elastic interactions between one dark soliton and one anti-bell-shaped soliton of the component p ˜ n . Figure 3d1,d2 exhibit the head-on elastic interactions between two kink solitons for solution r ˜ n . From Figure 3, we can clearly observe that the amplitudes and shapes of two solitons do not change before and after the interactions, which is consistent with the results of our asymptotic analysis.
Next, we implement numerical simulations in order to display the dynamical behaviors of the previous one-soliton and two-soliton solutions of Equation (6) by using the finite difference method [36]. Figure 4 and Figure 5 exhibit exact one-soliton and two-soliton solutions (42) and (41) of Equation (6), with numerical solutions without any noises and the perturbed numerical solutions with small noises 2 % and 10 % . Figure 4a1,a2,b1,b2 and Figure 5a1,a2,b1,b2 show that time evolutions of one-soliton solutions (42) and two-soliton solutions (41) without any noise are almost consistent with one of the corresponding exact one-soliton and two-soliton solutions. We can clearly observe that the numerical solutions almost recover the analytical exact solutions. In other words, these solutions have stable evolutions without any noise, which also indicates the accuracy of our numerical scheme. The numerical results in Figure 4c1,c2,d1,d2 and Figure 5c1,c2,d1,d2 are the numerical solutions obtained by adding 2 % and 10 % noises to both the initial exact solutions q ˜ n and r ˜ n , respectively. For 2 % noise, they almost evolve as before without any noise, while for 10 % noise, the evolutions show obvious fluctuation. That is to say that the exact solutions (42) and (41) have almost steady evolutions against a small noise in a relatively short time.
Case (2) When N = 2 , λ = λ i ( i = 1 , 2 , 3 , 4 ) , from (36), the 4-fold DT will provide the exact solutions as follows:
q ˜ n = ( 1 + α 2 ) a n + 1 ( 0 ) a n ( 0 ) , r ˜ n = α + b n ( 3 ) α d n ( 4 ) ,
where a n ( 0 ) = Δ a n ( 0 ) Δ 1 , b n ( 3 ) = Δ b n ( 3 ) Δ 1 and d n ( 4 ) = Δ d n ( 4 ) Δ 2 , in which the following is the case.
Δ 1 , n = λ 1 4 λ 1 2 λ 1 3 δ 1 , n λ 1 δ 1 , n λ 2 4 λ 2 2 λ 2 3 δ 2 , n λ 2 δ 2 , n λ 3 4 λ 3 2 λ 3 3 δ 3 , n λ 3 δ 3 , n λ 4 4 λ 4 2 λ 4 3 δ 4 , n λ 4 δ 4 , n , Δ 2 , n = λ 1 3 λ 1 λ 1 2 δ 1 , n δ 1 , n λ 2 3 λ 2 λ 2 2 δ 2 , n δ 2 , n λ 3 3 λ 3 λ 3 2 δ 3 , n δ 3 , n λ 4 3 λ 4 λ 4 2 δ 4 , n δ 4 , n , Δ a n ( 0 ) = 1 λ 1 2 λ 1 3 δ 1 , n λ 1 δ 1 , n 1 λ 2 2 λ 2 3 δ 2 , n λ 2 δ 2 , n 1 λ 3 2 λ 3 3 δ 3 , n λ 3 δ 3 , n 1 λ 4 2 λ 4 3 δ 4 , n λ 4 δ 4 , n , Δ b n ( 3 ) = λ 1 4 λ 1 2 λ 1 3 δ 1 , n 1 λ 2 4 λ 2 2 λ 2 3 δ 2 , n 1 λ 3 4 λ 3 2 λ 3 3 δ 3 , n 1 λ 4 4 λ 4 2 λ 4 3 δ 4 , n 1 , Δ d n ( 4 ) = λ 1 3 λ 1 λ 1 2 δ 1 , n λ 1 4 δ 1 , n λ 2 3 λ 2 λ 2 2 δ 2 , n λ 2 4 δ 2 , n λ 3 3 λ 3 λ 3 2 δ 3 , n λ 3 4 δ 3 , n λ 4 3 λ 4 λ 4 2 δ 4 , n λ 4 4 δ 4 , n .
when the appropriate parameters are selected, the solutions (43) may be the three-soliton or four-soliton solutions for which its corresponding evolution plots are shown in Figure 6. If one of λ i ( i = 1 , 2 , 3 , 4 ) is equal to 1 + α α 2 + 1 , the solutions reduce to three-soliton ones. We here choose α = 3 4 ; that is to say that one component of λ is 8 5 . Without loss of generality, we permit λ 1 = 8 5 . When the parameters are λ 1 = 8 5 , λ 2 = 2 5 , λ 3 = 3 , λ 4 = 1 4 , C 11 = C 12 = C 21 = C 22 = C 31 = C 32 = C 41 = 1 and C 42 = 1 2 , the corresponding evolution plots are displayed in Figure 6a1,b1,a2,b2. Figure 6a1,a2 demonstrate the elastic interactions between one bell-shaped soliton and two dark solitons of the component q ˜ n on nonzero constant background. Figure 6b1,b2 display the collisions among three anti-kink solitons of the component r ˜ n . When none of the parameters λ i ( i = 1 , 2 , 3 , 4 ) are equal to 1 + α α 2 + 1 , the solutions (43) are four-soliton solutions. By choosing the parameters λ 1 = 2 , λ 2 = 3 , λ 3 = 5 2 , λ 4 = 1 4 , α = 3 4 , C 11 = C 12 = C 21 = C 22 = C 31 = C 32 = C 41 = 1 and C 42 = 1 2 , we plot the corresponding evolution in Figure 6c1,d1,c2,d2. Figure 6c1,c2 demonstrate the elastic interactions between one bell-shaped soliton and three dark solitons of the component q ˜ n on nonzero constant background. Figure 6d1,d2 display the collision interactions among four anti-kink solitons of the component r ˜ n . It should be noted here that for the component r ˜ n , no matter how we choose the parameters, there will always be one anti-kink soliton with high amplitude, which makes the other three solitons look shorter. In fact, they represent four anti-kink solitons if we look through local magnification. Similarly, similarly to one-soliton or two-soliton solutions, we can draw the components p ˜ n and x ˜ n of Equation (1), and here we omit their figures.

4.2. Rational and Semi-Rational Solutions and Their Mathematical Characteristics

When m = 1 , the discrete generalized ( m , 2 N m ) -fold DT reduces to the discrete generalized ( 1 , 2 N 1 ) -fold DT, which can investigate some rational or semi-rational solutions of Equation (6). In order to obtain a more general solution, we rewrite the solution (39) as follows:
ϕ i , n = φ i , n ψ i , n = C i 1 τ i 1 n e ρ i 1 t + ζ ( ε ) + C i 2 τ i 2 n e ρ i 2 t ζ ( ε ) α λ i ( α 2 + 1 ) C i 1 τ i 1 n 1 e ρ i 1 t + ζ ( ε ) α λ i ( α 2 + 1 ) C i 2 τ i 2 n 1 e ρ i 2 t ζ ( ε ) ,
where ζ ( ε ) = α 4 λ i 4 4 α 4 λ i 2 + 2 α 2 λ i 4 6 α 2 λ i 2 + λ i 4 2 λ i 2 + 1 j = 0 2 N e j ε j , e j , ( j = 1 , 2 , . . . 2 N ) are real parameters. ε is an artificially introduced small parameter, while τ i 1 , τ i 2 , ρ i 1 and ρ i 2 are the same as ones defined in the previous subsection. We expand the vector function ϕ 1 , n with λ = λ 1 + ε in (44) as two Taylor series around ε = 0 given by the following.
ϕ 1 , n ( ε ) = ϕ 1 , n ( 0 ) + ϕ 1 , n ( 1 ) ε + ϕ 1 , n ( 2 ) ε 2 + ϕ 1 , n ( 3 ) ε 3 + ϕ 1 , n ( 4 ) ε 4 + ϕ 1 , n ( 5 ) ε 5 + .
In fact, by choosing different parameters C 11 and C 12 , we can obtain different rational or semi-rational solutions, and we here enumerate two types of Taylor series expansions of the vector function ϕ 1 , n by choosing two sets of different parameters C 11 and C 12 :
  • The first kind of expansion. By setting α = 4 3 , C 11 = C 12 = 1 , one can obtain ϕ 1 , n ( i ) as follows.
    ϕ 1 , n ( 0 ) = φ 1 , n ( 0 ) ψ 1 , n ( 0 ) = 2 ( 4 ) n e 27 25 t 10 3 ( 4 ) n e 27 25 t , ϕ 1 , n ( 1 ) = φ 1 , n ( 1 ) ψ 1 , n ( 1 ) , ϕ 1 , n ( 2 ) = φ 1 , n ( 2 ) ψ 1 , n ( 2 ) .
    Since the expressions of ϕ 1 , n ( 1 ) and ϕ 1 , n ( 2 ) are relatively long, we list them in the Appendix A. The rest ( φ 1 , n ( j ) , ψ 1 , n ( j ) ) T ( j = 3 , 4 , 5 , ) are omitted here.
  • The second kind of expansion. By setting α = 4 3 , C 11 = C 12 = 1 ε , one can obtain the new expansions, and we here only list the following first two expansions:
    ϕ 1 , n ( 0 ) = φ 1 , n ( 0 ) ψ 1 , n ( 0 ) = 2 15 4 n e 27 25 t ( 25 n 27 t + 200 e 0 ) 2 9 4 n e 27 25 t ( 25 n 27 t + 200 e 0 25 ) , ϕ 1 , n ( 1 ) = φ 1 , n ( 1 ) ψ 1 , n ( 1 ) ,
    where the following is the case.
    φ 1 , n ( 1 ) = 2 162 , 000 4 n e 27 25 t ( 62 , 500 n 3 202 , 500 n 2 t + 1 , 500 , 000 n 2 e 0 + 218 , 700 n t 2 3 , 240 , 000 n t e 0 + 12 , 000 , 000 n e 0 2 78 , 732 t 3 + 1 , 749 , 600 t 2 e 0 12 , 960 , 000 t e 0 2 + 32 , 000 , 000 e 0 3 + 150 , 000 n 2 + 202 , 500 n t + 1 , 200 , 000 n e 0 393 , 660 t 2 + 2 , 916 , 000 t e 0 23 , 125 n 204 , 525 t + 1 , 515 , 000 e 0 + 2 , 160 , 000 e 1 ) , ψ 1 , n ( 1 ) = 2 97 , 200 4 n e 27 25 t ( 62 , 500 n 3 202 , 500 n 2 t + 1 , 500 , 000 n 2 e 0 + 218 , 700 n t 2 3 , 240 , 000 n t e 0 + 12 , 000 , 000 n e 0 2 78 , 732 t 3 + 1 , 749 , 600 t 2 e 0 12 , 960 , 000 t e 0 2 + 32 , 000 , 000 e 0 3 37 , 500 n 2 + 607 , 500 n t 1 , 800 , 000 n e 0 612 , 360 t 2 + 6 , 156 , 000 t e 0 12 , 000 , 000 e 0 2 + 14 , 375 n 771 , 525 t + 3 , 015 , 000 e 0 + 2 , 160 , 000 e 1 39 , 375 ) .
In what follows, for rational or semi-rational solutions, we also discuss two cases: N = 1 , 2 .
Case (1) When N = 1 , using the first kind expansion, the discrete generalized ( 1 , 1 ) -fold DT yields the first-order rational solutions of Equation (6) as the follows:
q ˜ n = ( 1 + α 2 ) a n + 1 ( 0 ) a n ( 0 ) , r ˜ n = α + b n ( 1 ) α d n ( 2 ) ,
where a n ( 0 ) = Δ a n ( 0 ) Δ 1 , b n ( 1 ) = Δ b n ( 1 ) Δ 1 and d n ( 2 ) = Δ d n ( 2 ) Δ 2 in which the following is the case.
Δ 1 = λ 1 2 φ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) , Δ 2 = λ 1 φ 1 , n ( 0 ) ψ 1 , n ( 0 ) λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) ψ 1 , n ( 1 ) , Δ a n ( 0 ) = φ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) φ 1 , n ( 1 ) λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) , Δ b n ( 1 ) = λ 1 2 φ 1 , n ( 0 ) φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) φ 1 , n ( 1 ) , Δ d n ( 2 ) = λ 1 φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 0 ) λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 1 ) 2 λ 1 ψ 1 , n ( 0 ) .
Direct calculation results in the simplification expressions of solutions (48) below.
q ˜ n = 25 9 + 50 , 000 9 50 n 54 t + 400 e 0 + 45 50 n 54 t + 400 e 0 45 , r ˜ n = 1 + 5300 n 5724 t + 42 , 400 e 0 3520 4050 n 4374 t + 32 , 400 e 0 3645 .
We can see that q ˜ n possesses singularities at two paralleled straight lines 50 n 54 t + 400 e 0 + 45 = 0 and 50 n 54 t + 400 e 0 45 = 0 , while r ˜ n has singularity at one straight line 4050 n 4374 t + 32 , 400 e 0 3645 = 0 from which we can move these lines with singularities parallel to any position by changing the parameter e 0 . Moreover, we can conclude that q ˜ n 25 9 , r ˜ 1 1 as n ± , t ± . What we need to explain here is that, in the first-order rational solutions, we have an arbitrary parameter e 0 that can control the position of the rational solutions; that is to say that we can move the first-order rational solutions to any position we need. Through the transformations p n = ln | q n | α , x n = ln | r n | , we can provide the solutions of Equation (1) as follows.
p ˜ n = 3 4 ln | 25 9 + 50 , 000 9 50 n 54 t + 400 e 0 + 45 50 n 54 t + 400 e 0 45 | , x ˜ n = ln | 1 + 5300 n 5724 t + 42 , 400 e 0 3520 4050 n 4374 t + 32 , 400 e 0 3645 | .
It should be noted that if we permit λ 1 9 5 (e.g., λ 1 = 3 ) and expand ϕ 1 , n in (39) around λ 1 = 3 , then by using the above similar process we can obtain the hybrid solutions of polynomial functions and exponential functions. Compared with the above rational solutions, we here call this kind of hybrid solutions the semi-rational solutions. From (48), the generalized ( 1 , 1 ) -fold DT yields the first-order semi-rational solutions as follows:
q ˜ n = Q 1 Q 2 , r ˜ n = R 1 R 2 ,
where the following is the case:
Q 1 = 350 [ 140 , 625 e ξ 1 + ξ 2 48 14 ( 584 , 200 e 0 + 682 , 752 e 1 + 14 , 859 n 78 , 867 t + 27 14 ) e ξ 2 48 14 ( 584 , 200 e 0 682 , 752 e 1 14 , 859 n + 78 , 867 t + 27 14 ) e ξ 2 + 2 ( 6 , 132 , 672 n 2 65 , 100 , 672 n t + 482 , 227 , 200 n e 0 + 563 , 576 , 832 n e 1 + 172 , 767 , 168 t 2 2 , 559 , 513 , 600 t e 0 2 , 991 , 292 , 416 t e 1 + 9 , 479 , 680 , 000 e 0 2 + 22 , 157 , 721 , 600 e 0 e 1 + 12 , 947 , 816 , 448 e 1 2 158 , 769 ) + 9 e ξ 1 ξ 2 ] Q 2 = 9 [ 75 ( 84 + 23 14 ) e ξ 1 + 75 ( 84 + 23 14 ) e ξ 1 + 56 ( 234 n 1242 t + 9200 e 0 + 10 , 752 e 1 225 ) ] [ 75 ( 84 + 23 14 ) e ξ 2 + 75 ( 84 + 23 14 ) e ξ 2 + 56 ( 234 n 1242 t + 9200 e 0 + 10 , 752 e 1 + 225 ) ] , R 1 = 56 ( 9200 e 0 + 10 , 752 e 2 + 234 n 1242 t + 9 ) + 3 ( 127 14 84 ) e ξ 1 3 ( 127 14 + 84 ) e ξ 1 , R 2 = 504 ( 234 n 1242 t + 9200 e 0 + 10 , 752 e 1 225 ) 675 ( 84 + 23 14 ) e ξ 1 + 675 ( 84 + 23 14 ) e ξ 1
with
ξ 1 = n 2 ln 6576 1728 14 6576 + 1728 14 + 4 14 75 ( 27 t 200 e 0 ) , ξ 2 = n 2 ln 2192 576 14 2192 + 576 14 + 4 14 75 ( 27 t 200 e 0 ) .
With the aid of symbolic computation Maple, one can verify the solutions (48) and (49) by substituting them into Equation (6).
Case (2) When N = 2 , using the first kind expansion, the generalized ( 1 , 3 ) -fold DT yields the second-order rational solutions of Equation (6) as follows:
q ˜ n = ( 1 + α 2 ) a n + 1 ( 0 ) a n ( 0 ) = Q 1 Q 2 , r ˜ n = α + b n ( 3 ) α d n ( 4 ) = R 1 R 2 ,
where a n ( 0 ) = Δ a n ( 0 ) Δ 1 , b n ( 3 ) = Δ b n ( 3 ) Δ 1 and d n ( 4 ) = Δ d n ( 4 ) Δ 2 , in which the following is the case.
Δ 1 , n = λ 1 4 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) λ 1 4 φ 1 , n ( 1 ) + 4 λ 1 3 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 1 ) + 3 λ 1 2 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) λ 1 4 φ 1 , n ( 2 ) + 4 λ 1 3 φ 1 , n ( 1 ) + 6 λ 1 2 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 2 ) + 2 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 2 ) + 3 λ 1 2 ψ 1 , n ( 1 ) + 3 λ 1 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 2 ) + ψ 1 , n ( 1 ) λ 1 4 φ 1 , n ( 3 ) + 4 λ 1 3 φ 1 , n ( 2 ) + 6 λ 1 2 φ 1 , n ( 1 ) + 4 λ 1 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 3 ) + 2 λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 3 ψ 1 , n ( 3 ) + 3 λ 1 2 ψ 1 , n ( 2 ) + 3 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 3 ) + ψ 1 , n ( 2 ) , Δ 2 , n = λ 1 3 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 0 ) ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 1 ) + 3 λ 1 2 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 1 ) + 2 λ 1 ψ 1 , n ( 0 ) ψ 1 , n ( 1 ) λ 1 3 φ 1 , n ( 2 ) + 3 λ 1 2 φ 1 , n ( 1 ) + 3 λ 1 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 2 ψ 1 , n ( 2 ) + 2 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) ψ 1 , n ( 2 ) λ 1 3 φ 1 , n ( 3 ) + 3 λ 1 2 φ 1 , n ( 2 ) + 3 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 φ 1 , n ( 3 ) + φ 1 , n ( 2 ) λ 1 2 ψ 1 , n ( 3 ) + 2 λ 1 ψ 1 , n ( 2 ) + ψ 1 , n ( 1 ) ψ 1 , n ( 3 ) , Δ a n ( 0 ) = φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) φ 1 , n ( 1 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 1 ) + 3 λ 1 2 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) φ 1 , n ( 2 ) λ 1 2 φ 1 , n ( 2 ) + 2 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 2 ) + 3 λ 1 2 ψ 1 , n ( 1 ) + 3 λ 1 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 2 ) + ψ 1 , n ( 1 ) φ 1 , n ( 3 ) λ 1 2 φ 1 , n ( 3 ) + 2 λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 3 ψ 1 , n ( 3 ) + 3 λ 1 2 ψ 1 , n ( 2 ) + 3 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 3 ) + ψ 1 , n ( 2 ) , Δ b n ( 3 ) = λ 1 4 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 0 ) φ 1 , n ( 0 ) λ 1 4 φ 1 , n ( 1 ) + 4 λ 1 3 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 1 ) + 3 λ 1 2 ψ 1 , n ( 0 ) φ 1 , n ( 1 ) λ 1 4 φ 1 , n ( 2 ) + 4 λ 1 3 φ 1 , n ( 1 ) + 6 λ 1 2 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 2 ) + 2 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 2 ) + 3 λ 1 2 ψ 1 , n ( 1 ) + 3 λ 1 ψ 1 , n ( 0 ) φ 1 , n ( 2 ) λ 1 4 φ 1 , n ( 3 ) + 4 λ 1 3 φ 1 , n ( 2 ) + 6 λ 1 2 φ 1 , n ( 1 ) + 4 λ 1 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 3 ) + 2 λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 3 ψ 1 , n ( 3 ) + 3 λ 1 2 ψ 1 , n ( 2 ) + 3 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) φ 1 , n ( 3 ) , Δ d n ( 4 ) = λ 1 3 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 0 ) λ 1 4 ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 1 ) + 3 λ 1 2 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 1 ) + 2 λ 1 ψ 1 , n ( 0 ) λ 1 4 ψ 1 , n ( 1 ) 4 λ 1 3 ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 2 ) + 3 λ 1 2 φ 1 , n ( 1 ) + 3 λ 1 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 2 ψ 1 , n ( 2 ) + 2 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) λ 1 4 ψ 1 , n ( 2 ) 4 λ 1 3 ψ 1 , n ( 1 ) 6 λ 1 2 ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 3 ) + 3 λ 1 2 φ 1 , n ( 2 ) + 3 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 φ 1 , n ( 3 ) + φ 1 , n ( 2 ) λ 1 2 ψ 1 , n ( 3 ) + 2 λ 1 ψ 1 , n ( 2 ) + ψ 1 , n ( 1 ) λ 1 4 ψ 1 , n ( 3 ) 4 λ 1 3 ψ 1 , n ( 2 ) 6 λ 1 2 ψ 1 , n ( 1 ) 4 λ 1 ψ 1 , n ( 0 ) .
Since the simplified expressions are very complex, we list only their analytical expressions with e 0 = e 1 = e 2 = 0 of solutions q ˜ n = Q 1 Q 2 , r ˜ n = R 1 R 2 with the following:
Q 1 = 25 9 ( 16 ξ 6 + 240 ξ 5 2 , 125 , 000 n ξ 3 + 120 , 000 ξ 4 + 15 , 937 , 500 n ξ 2 907 , 500 ξ 3 + 56 , 445 , 312 , 500 n 2 2 , 711 , 718 , 750 n ξ + 2 , 278 , 125 ξ 2 9 , 357 , 421 , 875 n + 455 , 625 , 000 ξ ) ( 16 ξ 6 240 ξ 5 2 , 125 , 000 n ξ 3 + 120 , 000 ξ 4 15 , 937 , 500 n ξ 2 + 907 , 500 ξ 3 + 56 , 445 , 312 , 500 n 2 2 , 711 , 718 , 750 n ξ + 2 , 278 , 125 ξ 2 + 9 , 357 , 421 , 875 n 455 , 625 , 000 ξ ) , Q 2 = ( 16 ξ 6 2160 ξ 5 2 , 125 , 000 n ξ 3 143 , 437 , 500 n ξ 2 + 5 , 467 , 500 ξ 3 + 56 , 445 , 312 , 500 n 2 5 , 899 , 218 , 750 n ξ + 184 , 528 , 125 ξ 2 + 12 , 498 , 046 , 875 n ) ( 16 ξ 6 + 2160 ξ 5 2 , 125 , 000 n ξ 3 + 143 , 437 , 500 n ξ 2 5 , 467 , 500 ξ 3 + 56 , 445 , 312 , 500 n 2 5 , 899 , 218 , 750 n ξ + 184 , 528 , 125 ξ 2 12 , 498 , 046 , 875 n ) , R 1 = 10 , 000 ξ 6 150 , 000 ξ 5 + 75 , 000 , 000 ξ 4 + ( 1 , 328 , 125 , 000 n + 567 , 187 , 500 ) ξ 3 + ( 9 , 960 , 937 , 500 n + 1 , 423 , 828 , 125 ) ξ 2 + ( 1 , 694 , 824 , 218 , 750 n 284 , 765 , 625 , 000 ) ξ + 35 , 278 , 320 , 312 , 500 n 2 + 5 , 848 , 388 , 671 , 875 n , R 2 = 104 , 976 ξ 6 + 14 , 171 , 760 ξ 5 + ( 13 , 942 , 125 , 000 n 35 , 872 , 267 , 500 ) ξ 3 + ( 941 , 093 , 437 , 500 n + 1 , 210 , 689 , 028 , 125 ) ξ 2 38 , 704 , 774 , 218 , 750 n ξ + 370 , 337 , 695 , 312 , 500 n 2 81 , 999 , 685 , 546 , 875 n ,
in which ξ = 25 n 27 t . Next, we implement asymptotic analysis to study the second-order rational solutions q ˜ n and r ˜ n . Using ξ 1 = ξ + ( 459 100 + 1377 500 5 ) 1 3 t 1 3 , ξ 2 = ξ + ( 459 100 1377 500 5 ) 1 3 t 1 3 and c = ( 459 100 + 1377 500 5 ) 1 3 ( 459 100 1377 500 5 ) 1 3 > 0 , we obtain the solutions q ˜ n and r ˜ n with two different asymptotic states when | t | , which are listed as follows:
(i)
If ξ 1 = ξ + ( 459 100 + 1377 500 5 ) 1 3 t 1 3 is fixed, from ξ 2 = ξ 1 c t 1 3 we have ξ 2 when t ± . Calculating the limits of solutions q ˜ n and r ˜ n in (50) then yields the following.
q ˜ n q 1 ± = 25 9 + 50 , 000 9 ( 4 ξ 1 2 2025 ) , r ˜ n r 1 ± = 1 + 11 , 872 ξ 1 + 298 , 370 13 , 122 ξ 1 295 , 245 .
(ii)
If ξ 2 = ξ + ( 459 100 1377 500 5 ) 1 3 t 1 3 is fixed, from ξ 1 = ξ 2 + c t 1 3 we have ξ 1 ± when t ± , then calculating the limits of solutions q ˜ n and r ˜ n in (50) has the following asymptotic expressions in the following form.
q ˜ n q 2 ± = 25 9 + 50 , 000 9 ( 4 ξ 2 2 2025 ) , r ˜ n r 2 ± = 1 + 11 , 872 ξ 2 + 298 , 370 13 , 122 ξ 2 295 , 245 .
It can be observed that q 1 ± and q 2 ± possess singularities at four curves 2 ξ 1 ± 45 = 0 , 2 ξ 2 ± 45 = 0 , which also are the four center trajectories of solution q ˜ n , while r 1 ± and r 2 ± possess singularities at two curves 13 , 122 ξ 1 295 , 245 = 0 , 13 , 122 ξ 2 295 , 245 = 0 , which are also the two center trajectories of solution r ˜ n . In order to show the correctness of our analysis results, we draw the density plots of the rational solutions and the trajectory plots obtained after asymptotic analysis, respectively, which are shown in Figure 7. Through comparison, we find that the singularity of rational solutions is completely consistent with these trajectories, which also shows the correctness of our asymptotic analysis results of second-order rational solutions. In addition, from asymptotic expressions (51) and (52), we can clearly observe that the asymptotic expressions of second-order rational solutions are consistent with the first-order rational solutions. The main difference is that first-order rational solutions’ trajectories are straight lines, while higher-order rational solutions’ trajectories are curves.
Next, we will no longer calculate higher-order rational solutions, but we will summarize some mathematical characteristics of the higher-order rational solutions for Equation (6) listed in Table 1 and Table 2. In the two tables, the first column represents the order number of the rational solutions, the second and fourth columns represent the highest power of numerator (HPN) polynomials involved in solutions q n , r n , the third and fifth columns represent the highest power of denominator (HPD) polynomials involved in solutions q n , r n and the sixth and seventh columns mean the background levels of solutions q n , r n . In Table 1, these rational solutions are given by using the first kind of expansion that was represented in the previous calculation of rational solutions, while these rational solutions in Table 2 are obtained by using the second kind of expansion, which is omitted from the calculation. From Table 1, we can easily observe that for the rational solution q n of order j, HPN and HPD are both 2 j ( 2 j 1 ) , whereas for the rational solution r n of order j, HPN and HPD are j ( 2 j 1 ) . From Table 2, we know that HPN and HPD of the rational solutions q n , r n of order j are 2 j ( 2 j + 1 ) and j ( 2 j + 1 ) , respectively.
Remark 3.
It should be noted here that, unlike our previous research in [17,32], we have obtained rational and semi-rational solutions with arbitrary controllable parameters. By using these parameters, we can control the positions where the rational and semi-rational solutions appear, which is not available in our previous researches in Refs. [17,32]. The introduction of these controllable parameters is an important innovation of this paper.

4.3. Hybrid Solutions of Exponential Function and Rational Solutions

When m = 2 , the discrete generalized ( m , 2 N m ) -fold DT reduces to the discrete generalized ( 2 , 2 N 2 ) -fold DT, which can provide hybrid solutions of exponential function and rational solutions of Equation (6). We will also discuss two cases: N = 1 , 2 .
Case (1) When N = 1 , for λ 1 = 9 5 ( i.e., α = 4 3 ) , by using the second kind expansion for λ 2 9 5 ( e.g., λ 2 = 3 ) without using Taylor expansion, the discrete generalized ( 2 , 0 ) -fold DT will yield hybrid collision solutions of exponential function and first-order rational solutions of Equation (6) as follows:
q ˜ n = ( 1 + α 2 ) a n + 1 ( 0 ) a n ( 0 ) = Q 1 Q 2 , r ˜ n = α + b n ( 1 ) α d n ( 2 ) = R 1 R 2 ,
where a n ( 0 ) = Δ a n ( 0 ) Δ 1 , b n ( 1 ) = Δ b n ( 1 ) Δ 1 and d n ( 2 ) = Δ d n ( 2 ) Δ 2 , in which the following is the case.
Δ 1 = λ 1 2 φ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) λ 2 2 φ 2 , n λ 2 ψ 2 , n , Δ 2 = λ 1 φ 1 , n ( 0 ) ψ 1 , n ( 0 ) λ 2 φ 2 , n ψ 2 , n , Δ a n ( 0 ) = φ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) φ 2 , n λ 2 ψ 2 , n , Δ b n ( 1 ) = λ 1 2 φ 1 , n ( 0 ) φ 1 , n ( 0 ) λ 2 2 φ 2 , n φ 2 , n , Δ d n ( 2 ) = λ 1 φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 0 ) λ 2 φ 2 , n λ 2 2 ψ 2 , n .
Direct calculation provides the simplified analytic expressions of the solutions (53) as follows.
Q 1 = [ ( 400 9 ) n e 36 14 25 t ( 1 , 270 , 000 n 2 + 2 , 743 , 200 n t 20 , 320 , 000 n e 0 1 , 481 , 328 t 2 + 21 , 945 , 600 t e 0 81 , 280 , 000 e 0 2 + 1 , 562 , 500 ) 125 ( 2192 9 64 14 ) n e 72 14 25 t ( 400 e 0 + 5 + 15 14 + 50 n 54 t ) ( 400 e 0 5 + 15 14 + 50 n 54 t ) 500 ( 2192 9 + 64 14 ) n ( 50 n 5 + 15 14 + 54 t 400 e 0 ) ( 50 n + 5 + 15 14 + 54 t 400 e 0 ) ] , Q 2 = 9 [ ( 14 + 3 ) ( 50 n 54 t + 400 e 0 45 + 15 14 ) ( 12 8 3 14 ) n e 36 14 25 t + 2 ( 3 + 14 ) ( 50 n + 54 t 400 e 0 + 15 14 + 45 ) ( 12 + 8 3 14 ) n ] [ ( 3 + 14 ) ( 50 n 54 t + 400 e 0 + 15 14 + 45 ) ( 12 8 3 14 ) n e 36 14 25 t + 2 ( 14 + 3 ) ( 50 n + 54 t 400 e 0 45 + 15 14 ) ( 12 + 8 3 14 ) n ] , R 1 = ( 1 3 14 ) ( 15 14 + 400 e 0 + 50 n 54 t + 5 ) ( 12 8 3 14 ) n e 36 14 25 t 2 ( 1 + 3 14 ) ( 400 e 0 50 n + 54 t + 15 14 5 ) ( 12 + 8 3 14 ) n , R 2 = 27 ( 14 + 3 ) ( 50 n 54 t + 400 e 0 + 15 14 45 ) ( 12 8 3 14 ) n e 36 14 25 t 54 ( 14 3 ) ( 400 e 0 50 n + 54 t + 15 14 + 45 ) ( 12 + 8 3 14 ) n .
For convenience of analysis, we chose e 0 = 0 and let ξ 1 = 25 n 27 t , ξ 2 = n ln 9 2 14 9 + 2 14 + 36 14 25 t ; the asymptotic states for solutions (53) when t ± are given as follows:
Before collision t , we have the following:
(i)
If ξ 1 is fixed , ξ 2 , we have the following.
q ˜ n q n 1 = 25 9 + 12 [ 4 ξ 1 2 60 ( 14 + 2 2 ) ξ 1 + 25 ( 617 + 72 7 ) ] 27 [ 4 ξ 1 2 60 ( 14 + 2 2 ) ξ 1 + 225 ( 13 + 8 7 ) ] , r ˜ n r n 1 = 1 10 ( 4 2 3 7 ) [ 124 ξ 1 5 ( 159 14 + 372 2 + 424 ) ] 837 ( 3 2 2 7 ) [ 2 ξ 1 15 ( 14 + 2 2 + 3 ) ] .
(ii)
If ξ 2 is fixed , ξ 1 + , we have the following.
q ˜ n q n 2 = 25 9 32 [ 89 + 5 cosh ( ξ 2 ln 2 ) ] 375 [ 16 2 cosh ( ξ 2 ln 2 ) ] , r ˜ n r n 2 = 1 20 31 27 10 cosh [ 1 2 ln 10 ( 95 + 24 14 ) 31 ( 46 + 12 14 ) ] 20 31 27 10 sinh [ 1 2 ln 10 ( 95 + 24 14 ) 31 ( 46 + 12 14 ) ] coth [ 1 2 ξ 2 + 1 2 ln 23 + 6 14 10 ] .
After collision t + , we observe the following cases:
(iii)
If ξ 1 is fixed , ξ 2 + , we have the following.
q ˜ n q n 1 + = 25 9 + 4 [ 4 ξ 1 2 + 60 ( 14 2 2 ) ξ 1 + 25 ( 617 72 7 ) ] 9 [ 4 ξ 1 2 + 60 ( 14 2 2 ) ξ 1 + 225 ( 13 8 7 ) ] , r ˜ n r n 1 + = 1 10 ( 4 2 + 3 7 ) [ 124 ξ 1 + 5 ( 159 14 372 2 424 ) ] 837 ( 3 2 + 2 7 ) [ 2 ξ 1 15 ( 14 2 2 3 ) ] .
(iv)
If ξ 2 is fixed , ξ 1 , we have the following.
q ˜ n q n 2 + = 25 9 32 [ 89 + 5 cosh ( ξ 2 ln 2 ) ] 375 [ 16 2 cosh ( ξ 2 ln 2 ) ] , r ˜ n r n 2 + = 1 20 31 27 10 cosh [ 1 2 ln 10 ( 95 + 24 14 ) 31 ( 46 + 12 14 ) ] 20 31 27 10 sinh [ 1 2 ln 10 ( 95 + 24 14 ) 31 ( 46 + 12 14 ) ] coth [ 1 2 ξ 2 + 1 2 ln 23 + 6 14 10 ] .
Due to the relationship between exponential function and hyperbolic function, in the results of asymptotic analysis, we turn the exponential function into a hyperbolic function. From the above analysis, we can observe that hyperbolic form solutions and rational solutions possess singularities before and after collisions, as shown in Figure 8. Next, let us analyze q ˜ n and r ˜ n in (53), respectively:
  • For solution q ˜ n , before collision q ˜ n has four singular lines L 1 , L 2 , L 3 , L 4 , in which the two singular lines L 1 , L 2 are obtained by solving 4 ξ 1 2 60 ( 14 + 2 2 ) ξ 1 + 225 ( 13 + 8 7 ) = 0 , while the other two singular lines L 3 , L 4 can be obtained by solving cosh ( ξ 2 ln 2 ) 8 = 0 . As t , the solution q ˜ n possesses singularities in four lines L 1 , L 2 , L 3 and L 4 , which also are its four trajectories. After collision, q ˜ n has four singular lines L 1 + , L 2 + , L 3 + and L 4 + , in which the two singular lines L 1 + and L 2 + are obtained by solving 4 ξ 1 2 + 60 ( 14 + 2 2 ) ξ 1 + 225 ( 13 + 8 7 ) = 0 , while the other two singular lines L 3 + and L 4 + are still given by solving cosh ( ξ 2 ln 2 ) 8 = 0 . As t + , solution q ˜ n possesses singularities in four lines L 1 + , L 2 + , L 3 + , L 4 + which also are its four trajectories;
  • For the solution r ˜ n , before collision r ˜ n has two singular lines L 1 (i.e., 2 ξ 1 15 ( 14 + 2 2 + 3 ) = 0 ) and L 2 (i.e., ξ 2 + ln 23 + 6 14 10 = 0 ); that is to say that as t , the solution r ˜ n possesses singularities in two lines L 1 , L 2 , which also are its two trajectories. After collision, r ˜ n has two singular lines L 1 + (i.e., 2 ξ 1 15 ( 14 2 2 3 ) = 0 ) and L 2 + (i.e., ξ 2 + ln 23 + 6 14 10 = 0 ); in other words, as t + , the solution r ˜ n possesses singularities in two lines L 1 + , L 2 + , which also are its two trajectories;
  • Through the above discussion, we can observe that the hyperbolic solutions and rational solutions do not change their directions before and after collisions, and the positions of singular lines of hyperbolic solutions do not change, while the singular lines of rational solutions have changed their positions. In order to show the correctness of our asymptotic analysis results, we draw the three-dimensional plots of hybrid solutions and the trajectory plots after asymptotic analysis, respectively, as shown in Figure 8. By comparing Figure 8a1,a2 with Figure 8b1,b2, we find that the singularities in the three-dimensional plots are completely consistent with the trajectories in the two-dimensional picture, which also verifies our asymptotic analysis’ correctness of hybrid solutions;
  • Here, the authors would like to say the following: In Section 4.1, when the soliton solutions are discussed separately, they are nonsingular; however, in the hybrid solutions, although taking the same parameters, these hyperbolic soliton solutions do become singular, and the possible reason is that the rational solutions in the hybrid solutions result in their singularity. This new property is worthy of further discussion.
Case (2) When N = 2 , for λ 1 = 9 5 ( i.e., α = 4 3 ) , by using the first kind expansion, for λ 2 9 5 ( e.g., λ 2 = 3 ) without using Taylor expansion, the discrete generalized ( 2 , 1 ) -fold DT will provide hybrid collision solutions of standard one-soliton and second-order rational solutions of Equation (6) as follows:
q ˜ n = ( 1 + α 2 ) a n + 1 ( 0 ) a n ( 0 ) , r ˜ n = α + b n ( 3 ) α d n ( 4 ) ,
where a n ( 0 ) = Δ a n ( 0 ) Δ 1 , b n ( 3 ) = Δ b n ( 3 ) Δ 1 and d n ( 4 ) = Δ d n ( 4 ) Δ 2 , in which the following is the case.
Δ 1 , n = λ 1 4 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) λ 1 4 φ 1 , n ( 1 ) + 4 λ 1 3 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 1 ) + 3 λ 1 2 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) λ 1 4 φ 1 , n ( 2 ) + 4 λ 1 3 φ 1 , n ( 1 ) + 6 λ 1 2 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 2 ) + 2 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 2 ) + 3 λ 1 2 ψ 1 , n ( 1 ) + 3 λ 1 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 2 ) + ψ 1 , n ( 1 ) λ 2 4 φ 2 , n λ 2 2 φ 2 , n λ 2 3 ψ 2 , n λ 2 ψ 2 , n , Δ 2 , n = λ 1 3 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 0 ) ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 1 ) + 3 λ 1 2 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 1 ) + 2 λ 1 ψ 1 , n ( 0 ) ψ 1 , n ( 1 ) λ 1 3 φ 1 , n ( 2 ) + 3 λ 1 2 φ 1 , n ( 1 ) + 3 λ 1 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 2 ψ 1 , n ( 2 ) + 2 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) ψ 1 , n ( 2 ) λ 2 3 φ 2 , n λ 2 φ 2 , n λ 2 2 ψ 2 , n ψ 2 , n , Δ a n ( 0 ) = φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 0 ) φ 1 , n ( 1 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 1 ) + 3 λ 1 2 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) φ 1 , n ( 2 ) λ 1 2 φ 1 , n ( 2 ) + 2 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 2 ) + 3 λ 1 2 ψ 1 , n ( 1 ) + 3 λ 1 ψ 1 , n ( 0 ) λ 1 ψ 1 , n ( 2 ) + ψ 1 , n ( 1 ) φ 2 , n λ 2 2 φ 2 , n λ 2 3 ψ 2 , n λ 2 ψ 2 , n ,
Δ b n ( 3 ) = λ 1 4 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 0 ) φ 1 , n ( 0 ) λ 1 4 φ 1 , n ( 1 ) + 4 λ 1 3 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 1 ) + 2 λ 1 φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 1 ) + 3 λ 1 2 ψ 1 , n ( 0 ) φ 1 , n ( 1 ) λ 1 4 φ 1 , n ( 2 ) + 4 λ 1 3 φ 1 , n ( 1 ) + 6 λ 1 2 φ 1 , n ( 0 ) λ 1 2 φ 1 , n ( 2 ) + 2 λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 3 ψ 1 , n ( 2 ) + 3 λ 1 2 ψ 1 , n ( 1 ) + 3 λ 1 ψ 1 , n ( 0 ) φ 1 , n ( 2 ) λ 2 4 φ 2 , n λ 2 2 φ 2 , n λ 2 3 ψ 2 , n φ 2 , n , Δ d n ( 4 ) = λ 1 3 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 0 ) λ 1 4 ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 1 ) + 3 λ 1 2 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 1 ) + φ 1 , n ( 0 ) λ 1 2 ψ 1 , n ( 1 ) + 2 λ 1 ψ 1 , n ( 0 ) λ 1 4 ψ 1 , n ( 1 ) 4 λ 1 3 ψ 1 , n ( 0 ) λ 1 3 φ 1 , n ( 2 ) + 3 λ 1 2 φ 1 , n ( 1 ) + 3 λ 1 φ 1 , n ( 0 ) λ 1 φ 1 , n ( 2 ) + φ 1 , n ( 1 ) λ 1 2 ψ 1 , n ( 2 ) + 2 λ 1 ψ 1 , n ( 1 ) + ψ 1 , n ( 0 ) λ 1 4 ψ 1 , n ( 2 ) 4 λ 1 3 ψ 1 , n ( 1 ) 6 λ 1 2 ψ 1 , n ( 0 ) λ 2 3 φ 2 , n λ 2 φ 2 , n λ 2 2 ψ 2 , n λ 2 4 ψ 2 , n .
Direct calculation provides the simplified analytic expressions with e 0 = e 1 = 0 as q ˜ n = F 1 F 2 , r ˜ n = G 1 G 2 . The expressions of F 1 , F 2 , G 1 and G 2 are very complex; thus, we list them in the Appendix B.
For the sake of analysis, let ξ 1 = 25 n 27 t , ξ 2 = n ln 9 2 14 9 + 2 14 + 36 14 25 t be the asymptotic analysis for solutions (54) when | t | ± is given as follows:
Before collision t :
(i)
If ξ 2 is fixed , ξ 1 + , we have the following:
q ˜ n q n 1 = 25 9 + Q 1 Q 2 , r ˜ n r n 1 = 1 + R 1 R 2 ,
where the following is the case.
Q 1 = 4 46 , 875 [ 4 ξ 1 4 120 ξ 1 3 + 24 , 900 ξ 1 2 3125 ( 45 14 + 68 n ) ξ 1 + 1875 ( 1386 + 850 14 n ) ] , Q 2 = 1 3 , 906 , 250 , 000 [ 8 ( 14 + 3 ) ξ 1 3 900 ξ 1 2 100 ( 567 + 54 14 ) ξ 1 + 12 , 500 ( 17 14 + 51 ) n + 151 , 875 14 ] [ 8 ( 14 3 ) ξ 1 3 900 ξ 1 2 + 100 ( 567 54 14 ) ξ 1 + 12 , 500 ( 17 14 51 ) n 151 , 875 14 ] , R 1 = 32 ( 1056 14 3293 ) ξ 1 3 3 , 561 , 600 ξ 1 2 + 1200 ( 207 , 959 20 , 508 14 ) ξ 1 + 50 , 000 ( 17 , 952 14 55 , 981 ) n 671 , 332 , 500 14 , R 2 = 34 , 992 ( 14 3 ) ξ 1 3 + 3 , 936 , 600 ξ 1 2 + 11 , 809 , 800 ( 21 2 14 ) ξ 1 10 , 935 , 000 ( 598 23 14 ) n + 664 , 301 , 250 14 .
(ii)
If ξ 2 is fixed , ξ 1 + , we have the following.
q ˜ n q n 2 = 25 9 + 448 45 cosh ( ξ 2 ln 2 ) 207 , r ˜ n r n 2 = 1 2 953 , 611 2187 cosh [ 1 2 ln 971 , 111 + 49 , 050 14 953 , 611 ] 2 953 , 611 2187 sinh [ 1 2 ln 971 , 111 + 49 , 050 14 953 , 611 ] coth [ 1 2 ξ 2 + 1 2 ln 23 + 6 14 10 ] .
After collision t + :
(iii)
If ξ 1 is fixed , ξ 2 + , we have the following:
q ˜ n q n 1 + = 16 25 + Q 1 Q 2 , r ˜ n r n 1 + = 1 + R 1 R 2 ,
with the following case.
Q 1 = 1 45 , 875 [ 4 ξ 1 4 + 120 ξ 1 3 + 24 , 900 ξ 1 2 + 3125 ( 45 14 68 n ) ξ 1 + 1875 ( 1386 850 14 n ) ] , Q 2 = 1 15 , 625 , 000 , 000 [ 8 ( 14 3 ) ξ 1 3 + 900 ξ 1 2 + 100 ( 567 54 14 ) ξ 1 + 12 , 500 ( 17 14 51 ) n + 151 , 875 14 ] [ 8 ( 14 + 3 ) ξ 1 3 + 900 ξ 1 2 100 ( 567 + 54 14 ) ξ 1 + 12 , 500 ( 17 14 + 51 ) n + 151 , 875 14 ] , R 1 = 16 ( 1056 14 + 3293 ) ξ 1 3 1 , 780 , 800 ξ 1 2 + 600 ( 20 , 508 14 + 207 , 959 ) ξ 1 25 , 000 ( 17 , 592 14 + 55 , 981 ) n + 335 , 666 , 250 14 , R 2 = 17 , 496 ( 14 + 3 ) ξ 1 3 + 1 , 968 , 300 ξ 1 2 900 ( 118 , 089 14 + 1 , 240 , 029 ) ξ 1 + 2500 ( 185 , 886 14 + 557 , 685 ) n 332 , 150 , 625 14 .
(iv)
If ξ 2 is fixed , ξ 1 , we have the following.
q ˜ n q n 2 + = 25 9 + 448 45 cosh ( ξ 2 ln 2 ) 207 , r ˜ n r n 2 + = 1 2 953 , 611 2187 cosh [ 1 2 ln 971 , 111 + 49 , 050 14 953 , 611 ] 2 953 , 611 2187 sinh [ 1 2 ln 971 , 111 + 49 , 050 14 953 , 611 ] coth [ 1 2 ξ 2 + 1 2 ln 23 + 6 14 10 ] .
From the above analysis, we can observe that the standard solitons and rational solutions due to the asymptotic results of hybrid solutions are formally consistent with those of standard solitons and rational solutions alone in the previous two subsections, which also shows the correctness of our analysis from another aspect. We can also draw their structures by imitating the case of N = 1 , but the figures in this case are more complex and they do not look good, and we omit the discussion here.
Remark 4.
It should be emphasized here that in our previous research [17,32], we did not carry out the asymptotic analysis of hybrid soliton-and rational solutions. However, in this paper, we have overcome this problem and realized the asymptotic analysis of these hybrid solutions, which were not available in our previous research in [17,32]. Therefore, the asymptotic analysis of hybrid solutions also is another important innovation of this paper.

5. Conclusions

In this paper, we have investigated Equation (6), which may describe particle vibrations in lattices with an exponential interaction force. The main contributions of this paper are as follows: (i) Some integrable properties of Equation (6) have been investigated, such as lattice hierarchy, Hamiltonian structure and infinite conservation laws. (ii) The discrete ( m , 2 N m )-fold DT for Equation (6) has been constructed in detail for the first time. (iii) By applying the resulting discrete DT, we have obtained the soliton solutions, rational and semi-rational solutions and hybrid solutions of Equation (6), and their limit states were discussed by using the asymptotic analysis technique. Multi-soliton solutions and relevant numerical simulations are shown in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. (iv) Asymptotic analysis and mathematical characteristics of rational solutions are shown in Figure 7 and listed in Table 1 and Table 2, from which we find that the asymptotic analysis results of rational solutions are consistent with their exact solutions. (v) We have obtained the rational solutions (48) and semi-rational solutions (49) with arbitrary controllable parameters, which can control their positions. (vi) Asymptotic states of the exponential-and-rational hybrid solutions are shown in Figure 8, from which we can observe that the trajectories of our asymptotic analysis agree with hybrid solutions.
The properties of Equation (1) or (6) discussed above are reported for the first time. We hope that these results in this paper might provide new insights in understanding lattice dynamics.

Author Contributions

Formal analysis, M.-L.Q.; funding acquisition, X.-Y.W.; methodology, M.-L.Q. and X.-Y.W.; validation, M.-L.Q. and M.Y.; writing—Original draft, M.-L.Q.; writing—revision, X.-Y.W. and M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially supported by National Natural Science Foundation of China under Grant No. 12071042 and Beijing Natural Science Foundation under Grant No. 1202006. M.-L.Q. is supported by the Postgraduate Science and Technology Innovation Project of Beijing Information Science and Technology University under Grant No. 5112111017.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

φ 1 , n ( 1 ) = 4 n 450 e 27 25 t ( 625 n 2 1350 n t + 10 , 000 n e 0 + 729 t 2 10 , 800 t e 0 + 40 , 000 e 0 2 + 500 n + 1215 t ) , ψ 1 , n ( 1 ) = 4 n 270 e 27 25 t ( 625 n 2 1350 n t + 10 , 000 n e 0 + 729 t 2 10 , 800 t e 0 + 40 , 000 e 0 2 750 n + 2565 t 10 , 000 e 0 + 625 ) , φ 1 , n ( 2 ) = 4 n 2 , 430 , 000 e 27 25 t ( 390 , 625 n 4 1 , 687 , 500 n 3 t + 12 , 500 , 000 n 3 e 0 + 2 , 733 , 750 n 2 t 2 40 , 500 , 000 n 2 t e 0 + 150 , 000 , 000 n 2 e 0 2 1 , 968 , 300 n t 3 + 43 , 740 , 000 n t 2 e 0 324 , 000 , 000 n t e 0 2 + 800 , 000 , 000 n e 0 3 + 531 , 441 t 4 15 , 746 , 400 t 3 e 0 + 174 , 960 , 000 t 2 e 0 2 864 , 000 , 000 t e 0 3 + 1 , 600 , 000 , 000 e 0 4 + 1 , 875 , 000 n 3 + 506 , 250 n 2 t + 30 , 000 , 000 n 2 e 0 7 , 654 , 500 n t 2 + 40 , 500 , 000 n t e 0 + 120 , 000 , 000 n e 0 2 + 5 , 314 , 410 t 3 78 , 732 , 000 t 2 e 0 + 291 , 600 , 000 t e 0 2 + 171 , 875 n 2 843 , 750 n t + 33 , 250 , 000 n e 0 + 54 , 000 , 000 n e 1 + 9 , 950 , 850 t 2 81 , 810 , 000 t e 0 58 , 320 , 000 t e 1 + 303 , 000 , 000 e 0 2 + 432 , 000 , 000 e 0 e 1 750 , 000 n + 1 , 822 , 500 t ) , ψ 1 , n ( 2 ) = 4 n 1 , 458 , 000 e 27 25 t ( 390 , 625 n 4 1 , 687 , 500 n 3 t + 12 , 500 , 000 n 3 e 0 + 2 , 733 , 750 n 2 t 2 40 , 500 , 000 n 2 t e 0 + 150 , 000 , 000 n 2 e 0 2 1 , 968 , 300 n t 3 + 43 , 740 , 000 n t 2 e 0 324 , 000 , 000 n t e 0 2 + 800 , 000 , 000 n e 0 3 + 531 , 441 t 4 15 , 746 , 400 t 3 e 0 + 174 , 960 , 000 t 2 e 0 2 864 , 000 , 000 t e 0 3 + 1 , 600 , 000 , 000 e 0 4 + 312 , 500 n 3 + 5 , 568 , 750 n 2 t 7 , 500 , 000 n 2 e 0 13 , 122 , 000 n t 2 + 121 , 500 , 000 n t e 0 180 , 000 , 000 n e 0 2 + 7 , 282 , 710 t 3 122 , 472 , 000 t 2 e 0 + 615 , 600 , 000 t e 0 2 800 , 000 , 000 e 0 3 1 , 234 , 375 n 2 10 , 968 , 750 n t + 40 , 750 , 000 n e 0 + 54 , 000 , 000 n e 1 + 22 , 526 , 100 t 2 195 , 210 , 000 t e 0 58 , 320 , 000 t e 1 + 453 , 000 , 000 e 0 2 + 432 , 000 , 000 e 0 e 1 + 718 , 750 n + 12 , 555 , 000 t 45 , 750 , 000 e 0 54 , 000 , 000 e 1 187 , 500 ) .

Appendix B

F 1 = 1 1 , 125 , 000 , 000 e 72 14 25 t ( 2192 9 64 14 ) n ( 112 , 500 14 n 2 243 , 000 14 t n + 131 , 220 14 t 2 + 125 , 000 n 3 405 , 000 n 2 t + 437 , 400 n t 2 157 , 464 t 3 22 , 500 14 n + 24 , 300 t 14 37 , 500 n 2 + 81 , 000 n t 43 , 740 t 2 + 1125 14 + 385 , 000 n 186 , 300 t 47 , 250 ) ( 112 , 500 14 n 2 243 , 000 14 t n + 131 , 220 14 t 2 + 125 , 000 n 3 405 , 000 n 2 t + 437 , 400 n t 2 157 , 464 t 3 + 22 , 500 14 n 24 , 300 t 14 + 37 , 500 n 2 81 , 000 n t + 43 , 740 t 2 + 1125 14 + 385 , 000 n 186 , 300 t + 47 , 250 ) 1 17 , 578 , 125 , 000 e 36 14 25 t ( 400 9 ) n ( 992 , 187 , 500 , 000 n 6 6 , 429 , 375 , 000 , 000 n 5 t + 17 , 359 , 312 , 500 , 000 n 4 t 2 24 , 997 , 410 , 000 , 000 n 3 t 3 + 20 , 247 , 902 , 100 , 000 n 2 t 4 8 , 747 , 093 , 707 , 200 n t 5 + 1 , 574 , 476 , 867 , 296 t 6 5 , 110 , 703 , 125 , 000 n 4 + 25 , 721 , 550 , 000 , 000 n 3 t 47 , 571 , 077 , 250 , 000 n 2 t 2 + 38 , 500 , 735 , 320 , 000 n t 3 11 , 542 , 579 , 655 , 400 t 4 + 9 , 134 , 693 , 750 , 000 n 2 8 , 943 , 291 , 000 , 000 n t + 2 , 348 , 608 , 365 , 000 t 2 138 , 427 , 734 , 375 ) + 1 281 , 250 , 000 ( 2192 9 + 64 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 125 , 000 n 3 + 405 , 000 n 2 t 437 , 400 n t 2 + 157 , 464 t 3 22 , 500 14 n + 24 , 300 t 14 + 37 , 500 n 2 81 , 000 n t + 43 , 740 t 2 + 1125 14 385 , 000 n + 186 , 300 t + 47 , 250 ) ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 125 , 000 n 3 + 405 , 000 n 2 t 437 , 400 n t 2 + 157 , 464 t 3 + 22 , 500 14 n 24 , 300 t 14 37 , 500 n 2 + 81 , 000 n t 43 , 740 t 2 + 1125 14 385 , 000 n + 186 , 300 t 47 , 250 ) , F 2 = [ 1 125 , 000 ( 14 + 3 ) e 36 14 25 t ( 12 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 + 125 , 000 n 3 405 , 000 n 2 t + 437 , 400 n t 2 157 , 464 t 3 202 , 500 14 n + 218 , 700 t 14 337 , 500 n 2 + 729 , 000 n t 393 , 660 t 2 + 91 , 125 14 + 685 , 000 n 510 , 300 t 425 , 250 ) + 1 62 , 500 ( 3 + 14 ) ( 12 + 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 125 , 000 n 3 + 405 , 000 n 2 t 437 , 400 n t 2 + 157 , 464 t 3 202 , 500 14 n + 218 , 700 14 t + 337 , 500 n 2 729 , 000 n t + 393 , 660 t 2 + 91 , 125 14 685 , 000 n + 510 , 300 t + 425 , 250 ) ] [ 1 125 , 000 ( 3 + 14 ) e 36 14 25 t ( 12 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 + 125 , 000 n 3 405 , 000 n 2 t + 437 , 400 n t 2 157 , 464 t 3 + 202 , 500 14 n 218 , 700 14 t + 337 , 500 n 2 729 , 000 n t + 393 , 660 t 2 + 91 , 125 14 + 685 , 000 n 510 , 300 t + 425 , 250 ) + 1 62 , 500 ( 14 + 3 ) ( 12 + 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 125 , 000 n 3 + 405 , 000 n 2 t 437 , 400 n t 2 + 157 , 464 t 3 + 202 , 500 14 n 218 , 700 14 t 337 , 500 n 2 + 729 , 000 n t 393 , 660 t 2 + 91 , 125 14 685 , 000 n + 510 , 300 t 425 , 250 ) ] , G 1 = 25 ( 3 14 1 ) e 36 14 25 t ( 12 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 14 t n + 131 , 220 14 t 2 + 125 , 000 n 3 405 , 000 n 2 t + 437 , 400 n t 2 157 , 464 t 3 + 22 , 500 14 n 24 , 300 t 14 + 37 , 500 n 2 81 , 000 n t + 43 , 740 t 2 + 1125 14 + 385 , 000 n 186 , 300 t + 47 , 250 ) + 50 ( 1 + 3 14 ) ( 12 + 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 125 , 000 n 3 + 405 , 000 n 2 t 437 , 400 n t 2 + 157 , 464 t 3 + 22 , 500 14 n 24 , 300 14 t 37 , 500 n 2 + 81 , 000 n t 43 , 740 t 2 + 1125 14 385 , 000 n + 186 , 300 t 47 , 250 ) , G 2 = 2187 ( 14 + 3 ) e 36 14 25 t ( 12 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 14 t n + 131 , 220 14 t 2 + 125 , 000 n 3 405 , 000 n 2 t + 437 , 400 n t 2 157 , 464 t 3 202 , 500 14 n + 218 , 700 t 14 337 , 500 n 2 + 729 , 000 n t 393 , 660 t 2 + 91 , 125 14 + 685 , 000 n 510 , 300 t 425 , 250 ) + 4374 ( 3 + 14 ) ( 12 + 8 3 14 ) n ( 112 , 500 14 n 2 243 , 000 t 14 n + 131 , 220 14 t 2 125 , 000 n 3 + 405 , 000 n 2 t 437 , 400 n t 2 + 157 , 464 t 3 202 , 500 14 n + 218 , 700 t 14 + 337 , 500 n 2 729 , 000 n t + 393 , 660 t 2 + 91 , 125 14 685 , 000 n + 510 , 300 t + 425 , 250 ) .

References

  1. Toda, M. Vibration of a chain with nonlinear interaction. J. Phys. Soc. Jpn. 1967, 22, 431–436. [Google Scholar] [CrossRef]
  2. Toda, M. Waves in nonlinear lattice. Prog. Theor. Phys. Suppl. 1970, 45, 174–200. [Google Scholar] [CrossRef]
  3. Wadati, M. Transformation theories for nonlinear discrete systems. Prog. Theor. Phys. Suppl. 1976, 59, 36–63. [Google Scholar] [CrossRef] [Green Version]
  4. Hirota, R.; Suzuki, K. Studies on lattice solitons by using electrical networks. J. Phys. Soc. Jpn. 1970, 28, 1366–1367. [Google Scholar] [CrossRef]
  5. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  6. Matveev, V.B.; Salle, M.A. Darboux Transformations and Solitons; Springer: Berlin, Germany, 1991. [Google Scholar]
  7. Ma, W.X.; Xu, X.X. A modified Toda spectral problem and its hierarchy of bi-hamiltonian lattice equations. J. Phys. A 2004, 37, 1323–1336. [Google Scholar] [CrossRef]
  8. Ma, W.X.; You, Y.C. Rational solutions of the Toda lattice equation in Casoratian form. Chaos Solitons Fractals 2004, 22, 395–406. [Google Scholar] [CrossRef]
  9. Ma, W.X. A Darboux transformation for the Volterra lattice equation. Anal. Math. Phys. 2019, 9, 1711–1718. [Google Scholar] [CrossRef]
  10. Suris, Y.B. The Problem of Integrable Discretization: Hamiltonian Approach; Birkhäuser Verlag: Basel, Switzerland, 2003. [Google Scholar]
  11. Suris, Y.B. Discrete time Toda systems. J. Phys. A 2018, 51, 333001. [Google Scholar] [CrossRef] [Green Version]
  12. Ruijsenaars, S.N.M. Relativistic Toda systems. Commun. Math. Phys. 1990, 133, 217–247. [Google Scholar] [CrossRef] [Green Version]
  13. Bracciali, C.F.; Silva, J.S.; Sri Ranga, A. Extended Relativistic Toda lattice, L-orthogonal polynomials and associated Lax pair. Acta. Appl. Math. 2019, 164, 137–154. [Google Scholar] [CrossRef] [Green Version]
  14. Bruschi, M.; Ragnisco, O. Lax representation and complete integrability for the periodic relativistic Toda lattice. Phys. Lett. A 1989, 134, 365–370. [Google Scholar] [CrossRef]
  15. Suris, Y.B. On some integrable systems related to the Toda lattice. J. Phys. A 1997, 30, 2235–2249. [Google Scholar] [CrossRef] [Green Version]
  16. Yang, H.X.; Shen, D.; Zhu, L.L. A hierarchy of Hamiltonian lattice equations associated with the relativistic Toda type system. Phys. Lett. A 2009, 373, 2695–2703. [Google Scholar] [CrossRef]
  17. Qin, M.L.; Wen, X.Y.; Yuan, C.L. Integrability, multi-soliton and rational solutions, and dynamical analysis for a relativistic Toda lattice system with one perturbation parameter. Commun. Theor. Phys. 2021, 73, 065003. [Google Scholar] [CrossRef]
  18. Yang, H.X.; Xu, X.X.; Sun, Y.P.; Ding, H.Y. Integrable relativistic Toda type lattice hierarchies, associated coupling systems and the Darboux transformation. J. Phys. A 2006, 39, 3933–3947. [Google Scholar] [CrossRef]
  19. Fan, F.C.; Xu, Z.G.; Shi, S.Y. N-fold Darboux transformations and exact solutions of the combined Toda lattice and relativistic Toda lattice equation. Anal. Math. Phys. 2020, 10, 31. [Google Scholar] [CrossRef]
  20. Andrew, P.; Zhu, Z.N. Darboux-Bäcklund transformation and explicit solutions to a hybrid lattice of the relativistic Toda lattice and the modified Toda lattice. Phys. Lett. A 2014, 378, 1510–1513. [Google Scholar]
  21. Zhou, R.G.; Jiang, Q.Y. A Darboux transformation and an exact solution for the relativistic Toda lattice equation. J. Phys. A 2005, 38, 7735–7742. [Google Scholar] [CrossRef]
  22. Maruno, K.; Kajiwara, K.; Oikawa, M. Casorati determinant solution for the discrete-time relativistic Toda lattice equation. Phys. Lett. A 1998, 241, 335–343. [Google Scholar] [CrossRef]
  23. Darvishi, M.T.; Khani, F. New exact solutions of a relativistic Toda lattice system. Chin. Phys. Lett. 2012, 29, 094101. [Google Scholar] [CrossRef]
  24. Choudhury, A.G. Bäcklund transformation for the first flows of the relativistic Toda hierarchy and associated properties. Cent. Eur. J. Phys. 2013, 11, 49–58. [Google Scholar] [CrossRef] [Green Version]
  25. Hamedi-Nezhad, S.; Sofla, M.Z.; Kavitha, L.; Kumar, V.S. New rational solutions for relativistic discrete Toda lattice system. Commun. Theor. Phys. 2014, 62, 363–372. [Google Scholar] [CrossRef]
  26. Vineta, L.; Yu, G.F.; Zhang, Y.N. On an integrable system related to the relativistic Toda lattice- Bäcklund transformation and integrable discretization. J. Differ. Equ. Appl. 2015, 21, 403–417. [Google Scholar] [CrossRef] [Green Version]
  27. Ablowitz, M.J.; Musslimani, Z.H. Integrable discrete P T symmetric model. Phys. Rev. E 2014, 90, 032912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Ablowitz, M.J.; Musslimani, Z.H. Integrable Nonlocal Nonlinear Equations. Stud. Appl. Math. 2017, 139, 7–59. [Google Scholar] [CrossRef] [Green Version]
  29. Tu, G.Z. A trace identity and its applications to the theory of discrete integrable systems. J. Phys. A 1990, 23, 3903–3922. [Google Scholar]
  30. Tu, G.Z. On Liouville integrability of zero-cuvature equations and the Yang hierarchy. J. Phys. A 1989, 22, 2375–2392. [Google Scholar]
  31. Tu, G.Z. The trace identity, a powerful tool for constructing the Hamiltonian structure of integrable systems. J. Math. Phys. 1989, 30, 330–338. [Google Scholar] [CrossRef]
  32. Yuan, C.L.; Wen, X.Y. Integrability, discrete kink multi-soliton solutions on an inclined plane background and dynamics in the modified exponential Toda lattice equation. Nonlinear Dyn. 2021, 105, 643–669. [Google Scholar] [CrossRef]
  33. Wadati, M.; Sanuki, H.; Konno, K. Relationships among inverse method, Bäcklund transformation and an infinite number of conservation laws. Progr. Theoret. Phys. 1975, 53, 419–436. [Google Scholar] [CrossRef] [Green Version]
  34. Wadati, M.; Watanabe, M. Conservation laws of a Volterra system and nonlinear self-dual network equation. Prog. Theor. Phys. 1977, 57, 808–811. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, D.J.; Chen, D.Y. The conservation laws of some discrete soliton systems. Chaos Solitons Fractals 2002, 14, 573–579. [Google Scholar] [CrossRef]
  36. Trefethen, L.N. Spectral Methods in MATLAB; SIAM: Philadelphia, PA, USA, 2000. [Google Scholar]
Figure 1. A one-dimensional lattice of particles with fixed ends (see also the first figure in Ref. [3]).
Figure 1. A one-dimensional lattice of particles with fixed ends (see also the first figure in Ref. [3]).
Symmetry 13 02315 g001
Figure 2. (Color online) One-soliton structures via solutions (42) with parameters λ 1 = 8 5 , λ 2 = 2 , α = 3 4 and C 11 = C 12 = C 21 = 1 = C 22 = 1 . (a1d1) The profiles of one-soliton solutions are q ˜ n , r ˜ n , p ˜ n and x ˜ n . (a2d2) The propagation processes for q ˜ n , r ˜ n , p ˜ n and x ˜ n at t = 15 (dashdot line), t = 0 (longdash line) and t = 15 (solid line).
Figure 2. (Color online) One-soliton structures via solutions (42) with parameters λ 1 = 8 5 , λ 2 = 2 , α = 3 4 and C 11 = C 12 = C 21 = 1 = C 22 = 1 . (a1d1) The profiles of one-soliton solutions are q ˜ n , r ˜ n , p ˜ n and x ˜ n . (a2d2) The propagation processes for q ˜ n , r ˜ n , p ˜ n and x ˜ n at t = 15 (dashdot line), t = 0 (longdash line) and t = 15 (solid line).
Symmetry 13 02315 g002
Figure 3. (Color online) Two-soliton interaction structures via solutions (41) with parameters λ 1 = 1 4 , λ 2 = 3 , α = 1 , C 11 = C 12 = C 21 = 1 and C 22 = 1 2 . (a1d1) The profiles of two-soliton solutions q ˜ n , r ˜ n , p ˜ n and x ˜ n . (a2d2) The propagation processes for q ˜ n , r ˜ n , p ˜ n and x ˜ n at t = 10 (dashdot line), t = 0 (longdash line) and t = 10 (solid line).
Figure 3. (Color online) Two-soliton interaction structures via solutions (41) with parameters λ 1 = 1 4 , λ 2 = 3 , α = 1 , C 11 = C 12 = C 21 = 1 and C 22 = 1 2 . (a1d1) The profiles of two-soliton solutions q ˜ n , r ˜ n , p ˜ n and x ˜ n . (a2d2) The propagation processes for q ˜ n , r ˜ n , p ˜ n and x ˜ n at t = 10 (dashdot line), t = 0 (longdash line) and t = 10 (solid line).
Symmetry 13 02315 g003
Figure 4. (Color online) One-soliton solutions (42) with the same parameters as Figure 2. (a1,a2) Exact solutions. (b1,b2) Numerical solutions without any noise. (c1,c2) Numerical solutions with a 2% noise. (d1,d2) Numerical solutions with a 10% noise.
Figure 4. (Color online) One-soliton solutions (42) with the same parameters as Figure 2. (a1,a2) Exact solutions. (b1,b2) Numerical solutions without any noise. (c1,c2) Numerical solutions with a 2% noise. (d1,d2) Numerical solutions with a 10% noise.
Symmetry 13 02315 g004
Figure 5. (Color online) Two-soliton solutions (42) with the same parameters as Figure 3. (a1,a2) Exact solutions. (b1,b2) Numerical solutions without any noise. (c1,c2) Numerical solutions with a 2% noise. (d1,d2) Numerical solutions with a 10% noise.
Figure 5. (Color online) Two-soliton solutions (42) with the same parameters as Figure 3. (a1,a2) Exact solutions. (b1,b2) Numerical solutions without any noise. (c1,c2) Numerical solutions with a 2% noise. (d1,d2) Numerical solutions with a 10% noise.
Symmetry 13 02315 g005
Figure 6. (Color online) (a1,b1) The profiles of three-soliton solutions q ˜ n , r ˜ n ; (a2,b2) the propagation processes of three-soliton solutions at t = 5 (dashdot line), t = 0 (longdash line) and t = 5 (solid line); (c1,d1) the profiles of four-soliton solutions q ˜ n , r ˜ n ; (c2,d2) the propagation processes of four-soliton solutions at t = 5 (dashdot line), t = 0 (longdash line) and t = 5 (solid line).
Figure 6. (Color online) (a1,b1) The profiles of three-soliton solutions q ˜ n , r ˜ n ; (a2,b2) the propagation processes of three-soliton solutions at t = 5 (dashdot line), t = 0 (longdash line) and t = 5 (solid line); (c1,d1) the profiles of four-soliton solutions q ˜ n , r ˜ n ; (c2,d2) the propagation processes of four-soliton solutions at t = 5 (dashdot line), t = 0 (longdash line) and t = 5 (solid line).
Symmetry 13 02315 g006
Figure 7. (Color online) Second-order rational solutions: (a1) the density plot of solution q ˜ n in (50); (b1) the trajectory plot of solution q ˜ n by expressions (51) and (52); (a2) the density plot of solution r ˜ n in (50); (b2) the trajectory plot of solution r ˜ n by expressions (51) and (52).
Figure 7. (Color online) Second-order rational solutions: (a1) the density plot of solution q ˜ n in (50); (b1) the trajectory plot of solution q ˜ n by expressions (51) and (52); (a2) the density plot of solution r ˜ n in (50); (b2) the trajectory plot of solution r ˜ n by expressions (51) and (52).
Symmetry 13 02315 g007
Figure 8. (Color online) Hybrid solutions of hyperbolic and rational solutions: (a1) the three-dimensional structure of solution q ˜ n in (53); (b1) the trajectory plot of solution q ˜ n ; (a2) the three-dimensional structure of solution r ˜ n in (53); (b2) the trajectory plot of solution r ˜ n .
Figure 8. (Color online) Hybrid solutions of hyperbolic and rational solutions: (a1) the three-dimensional structure of solution q ˜ n in (53); (b1) the trajectory plot of solution q ˜ n ; (a2) the three-dimensional structure of solution r ˜ n in (53); (b2) the trajectory plot of solution r ˜ n .
Symmetry 13 02315 g008
Table 1. Main mathematical characteristics of the j-th rational order solutions q n , r n .
Table 1. Main mathematical characteristics of the j-th rational order solutions q n , r n .
jHPN of q n HPD of q n HPN of r n HPD of r n Background of q n Background of r n
12211 1 + α 2 1
2121266 1 + α 2 1
330301515 1 + α 2 1
j 2 j ( 2 j 1 ) 2 j ( 2 j 1 ) j ( 2 j 1 ) j ( 2 j 1 ) 1 + α 2 1
Table 2. Main mathematical characteristics of the j-th rational order solutions q n , r n .
Table 2. Main mathematical characteristics of the j-th rational order solutions q n , r n .
jHPN of q n HPD of q n HPN of r n HPD of r n Background of q n Background of r n
16633 1 + α 2 1
220201010 1 + α 2 1
342422121 1 + α 2 1
j 2 j ( 2 j + 1 ) 2 j ( 2 j + 1 ) j ( 2 j + 1 ) j ( 2 j + 1 ) 1 + α 2 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qin, M.-L.; Wen, X.-Y.; Yuen, M. A Relativistic Toda Lattice Hierarchy, Discrete Generalized (m,2Nm)-Fold Darboux Transformation and Diverse Exact Solutions. Symmetry 2021, 13, 2315. https://doi.org/10.3390/sym13122315

AMA Style

Qin M-L, Wen X-Y, Yuen M. A Relativistic Toda Lattice Hierarchy, Discrete Generalized (m,2Nm)-Fold Darboux Transformation and Diverse Exact Solutions. Symmetry. 2021; 13(12):2315. https://doi.org/10.3390/sym13122315

Chicago/Turabian Style

Qin, Meng-Li, Xiao-Yong Wen, and Manwai Yuen. 2021. "A Relativistic Toda Lattice Hierarchy, Discrete Generalized (m,2Nm)-Fold Darboux Transformation and Diverse Exact Solutions" Symmetry 13, no. 12: 2315. https://doi.org/10.3390/sym13122315

APA Style

Qin, M. -L., Wen, X. -Y., & Yuen, M. (2021). A Relativistic Toda Lattice Hierarchy, Discrete Generalized (m,2Nm)-Fold Darboux Transformation and Diverse Exact Solutions. Symmetry, 13(12), 2315. https://doi.org/10.3390/sym13122315

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