Next Article in Journal
Reconstructed Interpolating Differential Operator Method with Arbitrary Order of Accuracy for the Hyperbolic Equation
Next Article in Special Issue
Hopf Bifurcation Analysis of a Diffusive Nutrient–Phytoplankton Model with Time Delay
Previous Article in Journal
Survey of Eight Modern Methods of Hamiltonian Mechanics
Previous Article in Special Issue
Degenerated and Competing Dirichlet Problems with Weights and Convection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Darboux Transformation and N-Soliton Solutions of Gerdjikov–Ivanov Equation on a Time–Space Scale

College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(4), 294; https://doi.org/10.3390/axioms10040294
Submission received: 12 September 2021 / Revised: 31 October 2021 / Accepted: 3 November 2021 / Published: 5 November 2021
(This article belongs to the Special Issue Nonlinear Dynamical Systems with Applications)

Abstract

:
The Gerdjikov–Ivanov (GI) equation is one type of derivative nonlinear Schrödinger equation used widely in quantum field theory, nonlinear optics, weakly nonlinear dispersion water waves and other fields. In this paper, the coupled GI equation on a time–space scale is deduced from Lax pairs and the zero curvature equation on a time–space scale, which can be reduced to the classical and the semi-discrete GI equation by considering different time–space scales. Furthermore, the Darboux transformation (DT) of the GI equation on a time–space scale is constructed via a gauge transformation. Finally, N-soliton solutions of the GI equation are given through applying its DT, which are expressed by the Cayley exponential function. At the same time, one-solition solutions are obtained on three different time–space scales ( X = R , X = C and X = K p ).

1. Introduction

There are some practical problems that cannot be solved accurately by using only continuous or discrete analysis. In order to unify continuous and discrete analysis, a time scale was initiated by Stefan Hilger in 1988, which is an arbitrary nonempty closed subset of the real numbers [1,2,3]. In recent years, extensive research about time scales has been conducted, particularly in stability, oscillation and initial-boundary value problems [4,5,6,7,8]. In addition, time scale dynamic equations have wide application prospects in many areas, such as population dynamic models [9], epidemic models [10,11] and models of the financial consumption process [12,13].
Toda’s lattice, Hirota’s network and nonlinear Schrödinger dynamic equations were derived on a time–space scale by extending an Ablowitz–Ladik hierarchy of integrable dynamic systems on a time–space scale [14]. This extension facilitates a variety of modeling applications of Ablowitz–Ladik hierarchies, including optics and chaos in dispersion numerical schemes [15]. The formulas for solutions of boundary value problem of Burgers equation and heat equation were derived on a time–space scale by using the Cole–Hopf transformation. These formulas may be used to study the wave motion on a time–space scale. Sine-Gordon equation was obtained on a time–space scale and its solution expressed by the Cayley exponential function was given [16,17,18]. However, the development of time–space scales is relatively slow in nonlinear dynamical systems compared to other fields.
There are important applications regarding the derivative nonlinear Schrödinger (DNLS) equation in many fields [19]. In particular, in situations where higher order nonlinear effects need to be restored, a family of DNLS equations was investigated [20]. There are three famous DNLS equations, which are the DNLS I equation [21,22], DNLS II equation [23,24] and DNLS III equation [25]. The forms of these three equations are as follows
i q t + q x x + i q 2 q * x = 0 ,
i q t + q x x + i q q * q x = 0 ,
i q t + q x x + 1 2 q 3 q * 2 i q 2 q x * = 0 ,
where q * represents the complex conjugate of q. They can be transformed into each other by a gauge transformation [26]. Specifically, the last equation is also known as the Gerdjikov–Ivanov (GI) equation, which was discovered by Gerdjikov and Ivanov [27]. In recent years, several useful methods have been proposed for obtaining solutions of the GI equation, such as the Darboux transformation (DT) [28,29], algebra-geometric solution [30,31,32,33], Wronskian type solution [29,34] and Hamiltonian structures [35,36].
The advantage of DT is that new solutions can be obtained successively through iteration. The explicit soliton-like solution of the GI equation was obtained by its DT [26]. The explicit N-fold DT with multiparameters for the GI equation was constructed with the help of a gauge transformation [28]. The dark soliton, bright soliton, breather solution and periodic solution are given explicitly from different seed solutions. In this paper, the coupled GI equation on a time–space scale is deduced by the Lax matrix equation extended on a time–space scale. This extension will provide a wider range of nonlinear integrable dynamic models and promote solutions to practical problems.
This paper is organized as follows. In Section 2, the coupled GI equation on a time–space scale is obtained, which can be reduced to the classical and the semi-discrete GI equation. In Section 3, N-fold DT and N-soliton solutions of the GI equation on a time–space scale are constructed with the help of a gauge transformation. In particular, one-soliton solutions of the GI equation on three different time–space scales are obtained from seed solution. The last section is our conclusions.

2. GI Equation on a Time–Space Scale

For constructing the GI equation on a time–space scale, jump operators, graininess functions and the derivative are introduced as follows [1,2,3].
Definition 1.
For t , x T × X , backward jump operators are defined as
σ : T T , ρ : X X ,
σ t = sup s T : s < t , ρ x = sup y X : y < x .
For x X , the forward jump operator β x : X X is defined as β x = ρ 1 x = inf y X : y > x .
Definition 2.
The derivative associated with t (time) and x (space) variables is defined as
t f t , x = lim p μ t f t , x f σ t , x p ,
x f t , x = lim q ν x f t , x f ρ t , x q ,
where the graininess functions μ : T [ 0 , + ) , ν : X [ 0 , + ) are defined as
μ t = t σ t , ν x = x ρ x .
Note that,
f σ ( t , x ) : = f ( σ ( t ) , x ) = f t , x μ t t f t , x ,
f ρ ( t , x ) : = f ( t , ρ ( x ) ) = f t , x ν x x f t , x .
Definition 3.
The Cayley exponential function on a time scale is defined by
e α x , x 0 : = exp ι 0 x ζ μ ( s ) ( α ( s ) ) Δ s , e α ( x ) : = e α ( x , 0 ) ,
where α = α ( x ) is a given rd-continuous regressive function and
ζ h ( z ) : = 1 h log 1 + 1 2 z h 1 1 2 z h , h > 0 , ζ 0 ( z ) : = z .
When X = R and X = h Z , the Cayley exponential function becomes
e α ( x ) = e 0 x α ( s ) d s a n d e α ( x ) = 1 + 1 2 α h 1 1 2 α h x h ,
respectively.
Lemma 1.
Take T × X = R × R . The backward jump operators
σ t = sup , t = t , ρ x = sup , x = x ,
and the graininess functions
μ t = t σ t = 0 , ν x = x ρ x = 0 .
Lemma 2.
Take T × X = R × Z . The backward jump operators
σ t = sup , t = t , ρ x = sup { x 1 , x 2 , } = x 1 ,
and the graininess functions
μ t = t σ t = 0 , ν x = x ρ x = 1 .
Lemma 3.
When X = R , X = Z and X = K p , the derivative becomes
x f ( x ) = f x ( x ) , x f ( x ) = f ( x ) f ( x ) a n d x f ( x ) = f ( x ) f p 1 x 1 p 1 x ,
respectively.
In what follows, based on Lax pairs of DNLS equation from the generalized Kaup–Newell spectrum problem [32], a ∇-dynamical system is introduced
x ψ t , x = U t , x ψ t , x , t ψ t , x = V t , x ψ t , x ,
where
U = i λ 2 1 2 i q r λ q λ r i λ 2 + 1 2 i q r , V = A t , x B t , x C t , x A t , x ,
with ψ = ψ 1 t , x ψ 2 t , x , q and r are potential functions, and λ is a spectral parameter.
According to the compatibility condition x t ψ = t x ψ and derivative product rules [15], the zero curvature equation on a time–space scale is obtained
t U x V + U σ V V ρ U = 0 .
Then, substituting Equation (12) into Equation (13), we find
i C + C ρ λ 2 r σ A + r A ρ + t r λ 1 2 i ( q r ) σ C 1 2 i q r C ρ + x C = 0 , i B + B ρ λ 2 q σ A + q A ρ + t q λ 1 2 i ( q r ) σ B 1 2 i q r B ρ x B = 0 , i A A ρ λ 2 + q σ C r B ρ λ 1 2 i ( q r ) σ A + 1 2 i q r A ρ 1 2 i t ( q r ) x A = 0 , i A A ρ λ 2 + λ r σ B q C ρ 1 2 i ( q r ) σ A + 1 2 i q r A ρ 1 2 i t ( q r ) + x A = 0 .
Take A, B and C as quaternary polynomials of λ ,
A = j = 0 4 a j λ j , B = j = 0 4 b j λ j , C = j = 0 4 c j λ j .
Then, by substituting Equation (15) into Equation (14), these relations are obtained
a 4 = 2 i , a 1 = a 3 = b 0 = b 2 = b 4 = c 0 = c 2 = c 4 = 0 , b 3 = b 3 ρ + 2 q σ + q , c 3 = c 3 ρ + 2 r σ + r = 0 , b 1 = b 1 ρ + i q a 2 ρ + i q σ a 2 1 2 q r b 3 ρ 1 2 ( q r ) σ b 3 + i x b 3 , c 1 = c 1 ρ + i r a 2 ρ + i r σ a 2 1 2 q r c 3 ρ 1 2 ( q r ) σ c 3 i x c 3 , a 2 = x 1 ( 1 2 q σ c 1 1 2 r σ b 1 1 2 r b 1 ρ + 1 2 q c 1 ρ ) , a 0 = x 1 ( 1 2 i q r σ a 0 ρ + 1 4 q r r σ b 1 ρ 1 2 i r σ b 1 x 1 2 i q σ r σ a 0 + 1 4 q r σ b 1 r σ + 1 2 i q r a 0 ρ 1 4 q 2 r c 1 ρ + 1 2 i r σ q a 0 1 4 q r σ q c 1 1 2 i c 1 x q ) ,
and evolution equations on a time–space scale are obtained
t q = q a 0 ρ + q σ a 0 + 1 2 i q r b 1 ρ + 1 2 i q r σ b 1 + x b 1 ,
t r = r a 0 ρ r σ a 0 1 2 i q r c 1 ρ 1 2 i q r σ c 1 + x c 1 .
According to Equations (5) and (6), Equation (16) is reduced to
b 3 = 2 2 ν ( x ) x 1 q + q σ ,
c 3 = 2 2 ν ( x ) x 1 r + r σ ,
b 1 = 2 i 2 ν ( x ) x 1 m 1 a 2 + 1 2 2 ν ( x ) x 1 m 4 q + q σ ,
c 1 = 2 i 2 ν ( x ) x 1 m 2 a 2 + 1 2 2 ν ( x ) x 1 m 3 r + r σ ,
x a 0 = 1 2 i m 5 a 0 + 1 2 i r σ m 4 m 1 q m 4 m 2 a 2 + 1 8 r σ m 4 2 q + q σ 1 8 q m 4 m 3 r + r σ ,
x a 2 = 1 2 i m 1 m 2 m 2 m 1 2 ν ( x ) x a 2 + 1 2 m 1 m 3 r + r σ + 1 2 m 2 m 4 q + q σ ,
with
m 1 = q σ + q 1 ν ( x ) x 2 ν ( x ) x 1 , m 2 = r σ + r 1 ν ( x ) x 2 ν ( x ) x 1 , m 3 = ( q r ) σ + ( q r ) 1 ν ( x ) x + 2 i x 2 ν ( x ) x 1 , m 4 = ( q r ) σ + ( q r ) 1 ν ( x ) x 2 i x 2 ν ( x ) x 1 , m 5 = q r q r σ 1 ν ( x ) x + q r σ q σ r σ .
Then, the coupled GI equation on a time–space scale is obtained
t q = q 1 ν ( x ) x a 0 + q σ a 0 + 1 2 i q r 1 ν ( x ) x b 1 + 1 2 i ( q r ) σ b 1 + x b 1 , t r = r 1 ν ( x ) x a 0 r σ a 0 1 2 i q r 1 ν ( x ) x c 1 1 2 i ( q r ) σ c 1 + x c 1 ,
where a 0 , b 1 , c 1 are defined by Equations (21)–(23), respectively.
In the following, two special kinds of equations are given as follows.
Case I: Taking T × X = R × R , we find μ t = 0 , ν x = 0 .
Equations (21)–(23) are reduced to
b 1 = i q x , c 1 = i r x , a 0 = 1 2 r q x q r x + 1 4 i q 2 r 2 .
Then, Equation (25) is reduced to the coupled GI equation
i q t + q x x + i q 2 r x + 1 2 q 3 r 2 = 0 , i r t r x x + i r 2 q x 1 2 q 2 r 3 = 0 .
When r = q * , the classical GI equation is obtained
i q t + q x x + 1 2 q 3 q * 2 i q 2 q x * = 0 .
Case II: Taking T × X = R × Z , we find μ t = 0 , ν x = 1 .
f σ ( x , t ) = f ( x , t ) , f ρ ( x , t ) = E f ( x , t ) = f ( x , t ) ( 1 E ) f ( x , t ) ,
where E is the shift operator. Then, Equations (19)–(24) are reduced to
b 3 = 4 ( 1 + E ) 1 q ,
c 3 = 4 ( 1 + E ) 1 r ,
a 2 = ( 1 E ) 1 q r 2 + r m 7 q ,
b 1 = 2 i ( 1 + E ) 1 q ( 1 E ) 1 q r 2 + r m 7 q + ( 1 + E ) 1 m 7 q ,
c 1 = 2 i ( 1 + E ) 1 r ( 1 E ) 1 q r 2 + r m 7 q + ( 1 + E ) 1 m 6 q ,
a 0 = 1 2 i ( 1 E ) 1 r m 7 q q m 7 r ( 1 E ) 1 q r 2 + r m 7 q + 1 4 ( 1 E ) 1 r m 7 2 q q m 7 m 6 r ,
with
m 6 = q r + 2 i ( 1 E ) ( 1 + E ) 1 , m 7 = q r 2 i ( 1 E ) ( 1 + E ) 1 .
Therefore, the semi-discrete coupled GI equation is obtained
q t = q ( 1 + E ) a 0 + 1 2 i q r ( 1 + E ) b 1 + ( 1 E ) b 1 , r t = r ( 1 + E ) a 0 1 2 i q r ( 1 + E ) c 1 + ( 1 E ) c 1 ,
where a 0 , b 1 , and c 1 are defined by Equations (32)–(34), respectively.

3. DT of GI Equation on a Time–Space Scale

In this section, we construct a DT for GI equation and give its N-soliton solutions on a time–space scale.

3.1. Construction of DT on a Time–Space Scale

First, it can be shown by long calculations that Equation (12) is transformed to
U = i λ 2 σ 3 + λ Q 1 2 i Q 2 σ 3 , V = 2 i σ 3 λ 4 + B 3 λ 3 + a 2 σ 3 λ 2 + B 1 λ + a 0 σ 3 ,
with σ 3 is a Pauli matrix where Q = 0 q q * 0 , B 1 = 0 b 1 c 1 0 , B 3 = 0 b 3 c 3 0 , a j ( j = 0 , 2 ) , b j , c j ( j = 1 , 3 ) are defined by Equations (19)–(24), respectively.
Then, the ∇-dynamical system Equation (11) is transformed into
x ψ [ 1 ] = U [ 1 ] ψ [ 1 ] , t ψ [ 1 ] = V [ 1 ] ψ [ 1 ] ,
under a gauge transformation
ψ 1 = T 1 ψ .
Substituting Equation (38) into Equation (37), we find
U [ 1 ] T [ 1 ] = x T [ 1 ] + T [ 1 ] ρ U ,
V [ 1 ] T [ 1 ] = t T [ 1 ] + T [ 1 ] σ V ,
where
U 1 = i λ 2 σ 3 + λ Q 1 1 2 i Q 1 2 σ 3 , V 1 = 2 i σ 3 λ 4 + B 3 1 λ 3 + a 2 1 σ 3 λ 2 + B 1 1 λ + a 0 1 σ 3 ,
with Q 1 = 0 q 1 q 1 * 0 , B 1 1 = 0 b 1 1 c 1 1 0 , B 3 1 = 0 b 3 1 c 3 1 0 .
Assume
T 1 = T 0 + T 1 λ
where T 0 = a b c d , T 1 = a 11 b 12 c 21 d 22 .
Substituting Equation (42) into Equation (39) and comparing the coefficients in the terms of the same powers λ j ( j = 0 , , 5 ) on both sides of equation, we find
c 21 = b 12 = 0 , a 11 = d 22 = 1 , q 1 = q + i b + i b ρ , q 1 * = q * + i c + i c ρ .
Setting S = T 0 = s 11 s 12 s 21 s 22 , we obtain
T 1 = λ I S ,
q 1 = q i s 12 i s 12 ρ .
Substituting Equation (44) into Equation (39), we obtain
x S = 1 2 i S ρ Q 2 σ 3 + 1 2 i S Q 2 σ 3 + Q S 2 S ρ Q S + i S ρ S 2 σ 3 + i S 3 σ 3 .
Assume
S = H Λ H 1
with Λ = λ 0 0 λ * is an eigenvalue matrix, H = ψ 1 ψ 2 * ψ 2 ψ 1 * is a fundamental solution matrix and satisfies
x H = i σ 3 H Λ 2 + Q H Λ 1 2 i Q 2 σ 3 H , t H = 2 i σ 3 H Λ 4 + B 3 Λ 3 + a 2 σ 3 H Λ 2 + B 1 H Λ + a 0 σ 3 H .
It is easy to obtain
x S = x H Λ H 1 = 1 2 i S ρ Q 2 σ 3 + 1 2 i S Q 2 σ 3 + Q S 2 S ρ Q S + i S ρ S 2 σ 3 + i S 3 σ 3 ,
which means that Equation (47) yields Equation (46). From T 1 t + T 1 σ V = V 1 T 1 , we find
t S + λ I S σ 2 i σ 3 λ 4 + B 3 λ 3 + a 2 σ 3 λ 2 + B 1 λ + a 0 σ 3 = 2 i σ 3 λ 4 + B 3 1 λ 3 + a 2 1 σ 3 λ 2 + B 1 1 λ + a 0 1 σ 3 λ I S .
Comparing the coefficients in terms of the same powers λ j ( j = 0 , , 5 ) on both sides of Equation (50), we obtain
λ 0 : t S a 0 S σ σ 3 = a 0 1 σ 3 S , λ 1 : a 0 σ 3 S σ B 1 = B 1 1 S + a 0 1 σ 3 , λ 2 : B 1 a 2 S σ σ 3 = B 1 1 a 2 1 σ 3 S , λ 3 : a 2 σ 3 S σ B 3 = a 2 1 σ 3 B 3 1 S , λ 4 : B 3 + 2 i S σ σ 3 = B 3 1 + 2 i σ 3 S , λ 5 : 2 i σ 3 = 2 i σ 3 .
Then, the gauge transformations Equations (44) and (45) are proven to be DT of the GI equation on a time–space scale.

3.2. Soliton Solutions of the GI Equation on a Time–Space Scale

Soliton solutions of the GI equation on a time–space scale are constructed by applying its DT. First, Equation (11) is transformed to
x ψ 0 = U 0 ψ 0 = i λ 2 + 1 2 i q 0 q 0 * λ q 0 λ q 0 * i λ 2 1 2 i q 0 q 0 * ψ 0 , t ψ 0 = V 0 ψ 0 = 2 i λ 4 + a 2 0 λ 2 + a 0 0 b 3 0 λ 3 + b 1 0 λ c 3 0 λ 3 + c 1 0 λ 2 i λ 4 a 2 0 λ 2 a 0 0 ψ 0 ,
where ψ [ 0 ] = ψ 1 [ 0 ] ψ 2 [ 0 ] .
Let us set the spectral parameter λ = λ 1 . A one-fold DT of the GI equation on a time–space scale is constructed
ψ 1 = T 1 ψ 0 = λ I S 0 ψ 0 = λ s 11 0 s 12 0 s 21 0 λ s 22 0 ψ 0 , q 1 = q 0 i s 12 0 i s 12 0 ρ = q 0 i λ 1 λ 1 * ψ 1 [ 0 ] ψ 2 * [ 0 ] Δ 0 i λ 1 λ 1 * ψ 1 ρ [ 0 ] ψ 2 * ρ [ 0 ] Δ 0 ρ ,
where
S = 1 Δ 0 λ 1 ψ 1 [ 0 ] 2 + λ 1 * ψ 2 [ 0 ] 2 λ 1 * λ 1 ψ 1 [ 0 ] ψ 2 * [ 0 ] λ 1 * λ 1 ψ 1 * [ 0 ] ψ 2 [ 0 ] λ 1 * ψ 1 [ 0 ] 2 λ 1 ψ 2 [ 0 ] 2 ,
with Δ 0 = ψ 1 [ 0 ] 2 ψ 2 [ 0 ] 2 .
Under the DT (53), the -dynamical system (52) is transformed into
x ψ 1 = U 1 ψ 1 = i λ 2 + 1 2 i q 1 q 1 * λ q 1 λ q 1 * i λ 2 1 2 i q 1 q 1 * ψ 1 , t ψ 1 = V 1 ψ 1 = 2 i λ 4 + a 2 1 λ 2 + a 0 1 b 3 1 λ 3 + b 1 1 λ c 3 1 λ 3 + c 1 1 λ 2 i λ 4 a 2 1 λ 2 a 0 1 ψ 1 .
In what follows, taking the “seed solution” q 0 = 0 , we obtain eigenvectors ψ [ 0 ] of Equation (52) with λ = λ 1
ψ [ 0 ] = ψ 1 [ 0 ] ψ 2 [ 0 ] = e i λ 1 2 ( x , 0 ) e 2 i λ 1 4 ( t , 0 ) e i λ 1 2 ( x , 0 ) e 2 i λ 1 4 ( t , 0 ) ,
ψ ρ [ 0 ] = ψ 1 ρ [ 0 ] ψ 2 ρ [ 0 ] = 1 i λ 1 2 ν ( x ) e i λ 1 2 ( x , 0 ) e 2 i λ 1 4 ( t , 0 ) 1 + i λ 1 2 ν ( x ) e i λ 1 2 ( x , 0 ) e 2 i λ 1 4 ( t , 0 ) ,
where e ± i λ 1 2 ( x , 0 ) and e ± 2 i λ 1 4 ( t , 0 ) are Cayley exponential functions [18]. Then, a one-soliton solution of the GI equation on a time–space scale is obtained
q [ 1 ] = i λ 1 λ 1 * E 3 E 1 + E 2 + i λ 1 λ 1 * 1 i λ 1 * 2 ν ( x ) E 3 1 + i λ 1 * 2 ν ( x ) E 1 + 1 i λ 1 * 2 ν ( x ) E 2 ,
where
E 1 = e i λ 1 2 λ 1 * 2 ( x , 0 ) e 2 i λ 1 4 λ 1 * 4 ( t , 0 ) , E 2 = e i λ 1 2 λ 1 * 2 ( x , 0 ) e 2 i λ 1 4 λ 1 * 4 ( t , 0 ) , E 3 = e i λ 1 2 + λ 1 * 2 ( x , 0 ) e 2 i λ 1 4 + λ 1 * 4 ( t , 0 ) .
Similarly, we take the spectral parameter λ = λ 2 . A two-fold DT of the GI equation on a time–space scale is constructed
ψ 2 = T 2 ψ 1 = λ I S 1 ψ 1 = λ s 11 1 s 12 1 s 21 1 λ s 22 1 ψ 1 = T 2 T 1 ψ 0 , q 2 = q 1 i s 12 1 i s 12 1 ρ = q 1 i λ 2 λ 2 * ψ 1 [ 1 ] ψ 2 * [ 1 ] Δ 1 i λ 2 λ 2 * ψ 1 ρ [ 1 ] ψ 2 * ρ [ 1 ] Δ 1 ρ ,
where
S [ 1 ] = 1 Δ 1 λ 2 ψ 1 [ 1 ] 2 + λ 2 * ψ 2 [ 1 ] 2 λ 2 * λ 2 ψ 1 [ 1 ] ψ 2 [ 1 ] * λ 2 * λ 2 ψ 1 [ 1 ] * ψ 2 [ 1 ] λ 2 * ψ 1 [ 1 ] 2 λ 2 ψ 2 [ 1 ] 2 ,
with Δ 1 = ψ 1 [ 1 ] 2 ψ 2 [ 1 ] 2 .
When the spectral parameter λ = λ N , N-fold DT is constructed as follows
ψ N = T N ψ N 1 = λ I S N 1 ψ N 1 = λ s 11 N 1 s 12 N 1 s 21 N 1 λ s 22 N 1 ψ N 1 = T N T 3 T 2 T 1 ψ 0 , q N = q N 1 i s 12 N 1 i s 12 ρ N 1 = q 0 + i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 + i j = 1 N λ j λ j * ψ 1 ρ [ j 1 ] ψ 2 * ρ [ j 1 ] ψ 1 ρ [ j 1 ] 2 + ψ 2 ρ [ j 1 ] 2 .
An N-soliton solution of the GI equation on a time–space scale is obtained
q N = i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 + i j = 1 N λ j λ j * ψ 1 ρ [ j 1 ] ψ 2 * ρ [ j 1 ] ψ 1 ρ [ j 1 ] 2 + ψ 2 ρ [ j 1 ] 2 .
In what follows, N-fold DT and N-soliton solutions of the GI equation on three special time–space scales are obtained as follows.
Case I: Taking T × X = R × R , we obtain an N-fold DT of the classical GI equation
ψ N = T N ψ N 1 = λ I S N 1 ψ N 1 = λ s 11 N 1 s 12 N 1 s 21 N 1 λ s 22 N 1 ψ N 1 = T N T 3 T 2 T 1 ψ 0 , q N = q N 1 2 i s 12 N 1 = q 0 + 2 i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 .
When N = 1 , q 0 = 0 and the spectral parameter λ 1 = α 1 + i η 1 , we obtain a one-soliton solution of Equation (27)
q [ 1 ] = 2 η 1 e 2 i Y 1 sech 2 X 1 ,
where
X 1 = 4 α 1 η 1 x + 16 α 1 3 η 1 α 1 η 1 3 t , Y 1 = 2 α 1 2 η 1 2 x 4 α 1 2 6 α 1 2 η 1 2 t .
The profile of the one-soliton in Figure 1.
When N = 2 , q 0 = 0 and the spectral parameter λ 2 = α 2 + i η 2 , we obtain a two-soliton solution of Equation (27)
q [ 2 ] = 2 η 1 e 2 i Y 1 sech 2 X 1 4 η 2 M 1 M 2 e 2 i Y 2 α 1 M 1 sech 2 X 1 e 2 X 2 2 i Y 1 M 1 2 e 2 X 2 + M 2 2 e 2 X 2 + M 4 + M 5 + M 6 + 4 η 2 i η 1 M 2 sech 2 X 1 e 2 X 2 2 i Y 1 M 3 M 1 2 e 2 X 2 + M 2 2 e 2 X 2 + M 4 + M 5 + M 6 ,
where
M 1 = α 2 α 1 tanh 2 X 1 + η 2 η 1 i , M 2 = α 2 α 1 + η 2 η 1 tanh 2 X 1 i , M 3 = i α 1 η 1 sech 2 2 X 1 e 4 i Y 1 + 2 i Y 2 , M 4 = 2 i sech 2 X 1 sinh 2 i Y 1 2 i Y 2 η 1 α 2 + α 1 η 2 , M 5 = 2 sech 2 X 1 cosh 2 i Y 1 2 i Y 2 η 1 η 2 η 1 2 + α 1 α 2 α 1 2 , M 6 = sech 2 2 X 1 η 1 2 e 2 X 2 + α 1 2 e 2 X 2 , X 2 = 4 α 2 η 2 x + 16 α 2 3 η 2 α 2 η 2 3 t , Y 2 = 2 α 2 2 η 2 2 x 4 α 2 2 6 α 2 2 η 2 2 t .
Case II: Taking T × X = R × C , we find
μ ( t ) = 0 , ν ( x ) = 1 3 m + 1 , x L , 0 , x C \ L ,
where C is a Cantor set. L contains left discrete elements of C ,
L = k = 1 m a k 3 k + 1 3 m + 1 : m N , a k { 0 , 2 } , 1 k m .
Then, an N-fold DT of the GI equation is constructed
ψ [ N ] = T [ N ] ψ [ N 1 ] = ( λ I S [ N 1 ] ) ψ [ N 1 ] = λ s 11 [ N 1 ] s 12 [ N 1 ] s 21 [ N 1 ] λ s 22 [ N 1 ] ψ [ N 1 ] = T [ N ] T [ 3 ] T [ 2 ] T [ 1 ] ψ [ 0 ] , q [ N ] = q [ 0 ] + i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 + i j = 1 N λ j λ j * ψ 1 ρ [ j 1 ] ψ 2 * ρ [ j 1 ] ψ 1 ρ [ j 1 ] 2 + ψ 2 ρ [ j 1 ] 2 , x L , t R , q [ 0 ] + 2 i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 , x C \ L , t R .
According to Definition 3, we have
e ± 2 i λ 1 4 ( x , 0 ) = 1 ± i λ 1 4 3 m + 1 1 i λ 1 4 3 m + 1 x 3 m + 1 , e ± i λ 1 2 ( x , 0 ) = 1 ± i λ 1 2 2 × 3 m + 1 1 i λ 1 2 2 × 3 m + 1 x 3 m + 1 .
When N = 1 , q 0 = 0 and the spectral parameter λ 1 = α 1 + i η 1 , a one-soliton solution is obtained
q [ 1 ] = 1 N 1 3 m + 1 i α 1 2 2 M 7 i η 1 2 2 α 1 η 1 2 M 7 i α 1 2 i η 1 2 N 2 , x L , t R , 2 η 1 e 2 i Y 1 sech 2 X 1 , x C \ L , t R ,
where
N 1 = E 1 λ 1 = α 1 + i η 1 + E 2 λ 1 = α 1 + i η 1 , N 2 = E 1 λ 1 = α 1 + i η 1 E 2 λ 1 = α 1 + i η 1 , M 7 = 2 η 1 1 + i α 1 2 3 m + 1 2 i η 1 2 2 i α 1 2 3 m + 1 x e 24 i α 1 2 η 1 2 4 i α 1 4 t .
Case III: Taking T × X = R × K p , we find
μ t = 0 , ν ( x ) = ( 1 p 1 ) x , x = p k p Z , 0 , x = 0 ,
where p > 1 , p Z = p k : k Z and K p = p Z { 0 } .
Then, an N-fold DT is constructed
ψ [ N ] = T [ N ] ψ [ N 1 ] = ( λ I S [ N 1 ] ) ψ [ N 1 ] = λ s 11 [ N 1 ] s 12 [ N 1 ] s 21 [ N 1 ] λ s 22 [ N 1 ] ψ [ N 1 ] = T [ N ] T [ 3 ] T [ 2 ] T [ 1 ] ψ [ 0 ] , q [ N ] = q [ 0 ] + i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 + i j = 1 N λ j λ j * ψ 1 ρ [ j 1 ] ψ 2 * ρ [ j 1 ] ψ 1 ρ [ j 1 ] 2 + ψ 2 ρ [ j 1 ] 2 , x p Z , t R , q [ 0 ] + 2 i j = 1 N λ j λ j * ψ 1 [ j 1 ] ψ 2 * [ j 1 ] ψ 1 [ j 1 ] 2 + ψ 2 [ j 1 ] 2 , x = 0 , t R .
According to
a b f ( x ) x = 1 p 1 x = a b x f ( x ) , a b f ( ρ ( x ) ) x = ( p 1 ) x = a b x f ( x ) ,
we have
e ± i λ 1 2 ( x , 0 ) = e 1 p 1 x = 0 p k ± i λ 1 2 , e ± i λ 1 2 ( ρ ( x ) , 0 ) = e ( p 1 ) x = 0 p k ± i λ 1 2 .
When N = 1 , q 0 = 0 and the spectral parameter λ 1 = α 1 + i η 1 , a one-soliton solution is obtained
q [ 1 ] = η 1 e 16 α 1 η 1 3 16 α 1 3 η 1 4 i α 1 4 + 24 i α 1 2 η 1 2 t M 8 , x p Z , t R , 2 η 1 e 2 i Y 1 sech 2 X 1 , x = 0 , t R ,
where
M 8 = e 1 p 1 x = 0 p k 2 i α 1 2 + 2 i η 1 2 x sech 1 p 1 x = 0 p k 4 α 1 η 1 x + e ( p 1 ) k = 0 p k 2 i α 1 2 + 2 i η 1 2 x sech p 1 x = 0 p k 4 α 1 η 1 x .

4. Conclusions

In this paper, the coupled GI equation on a time–space scale was obtained by extending the Lax matrix equation on a time–space scale, which can be reduced to the classical GI equation. In particular, the semi-discrete GI equation was given by providing parallel computations for the discrete and continuous case. The standard DT of the GI equation was extended on a time–space scale. On this basis, its N-soliton solutions on a time–space scale were obtained, which were expressed using Cayley exponential functions.
The extension provides a wider range of nonlinear integrable dynamic models and promotes the study of nonlinear dynamic systems. By taking the “seed solution” q = 0 and λ = α + i β , one-solition solutions of the GI equation were obtained on three different time–space scales ( X = R , X = C and X = K p ). In one case, the exact solution (64) and its dynamic figure were obtained when x R . In the other cases, when x C \ L and x = 0 , exact solutions (68) and (71) were obtained and were similar to Equation (64). However, when x L and x p Z , the structures of solutions (68) and (71) were more complicated and their values were different from those of Equation (64) at those discontinuity points.
Due to the limitations of the computer, it was difficult to obtain their dynamic figures at this stage. Furthermore, there is another well-known equation, the Eckhaus equation, which possesses a very similar structure. The Eckhaus equation is also integrable and has soliton-like solutions expressed in terms of the hyperbolic functions [37,38]. Therefore, we will find the most effective way to reduce structures of solutions (68) and (71) on C and K p , and study the Eckhaus equation on a time–space scale in our future work.

Author Contributions

Conceptualization, H.D., Y.Z. and X.H.; methodology, H.D.; software, Y.Z.; validation, H.D., Y.Z. and Y.F.; formal analysis, M.L.; investigation, Y.Z. and M.L.; writing—original draft preparation, Y.Z. and X.H.; project administration, H.D.. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 11975143, 12105161, 61602188), Natural Science Foundation of Shandong Province (Grant No. ZR2019QD018), CAS Key Laboratory of Science and Technology on Operational Oceanography (Grant No. OOST2021-05), Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (Grant No. 2017RCJJ068, 2017RCJJ069).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their thanks to the editors and the reviewers for their kind comments to improve our paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hilger, S. Analysis on Measure Chains-A Unified Approach to Continuous and Discrete Calculus. Results Math. 1990, 18, 18–56. [Google Scholar] [CrossRef]
  2. Bohner, M.; Peterson, A. Advances in Dynamic Equations on a Time Scale; Birkhauser: Boston, MA, USA, 2003. [Google Scholar]
  3. Bohner, M.; Peterson, A. Dynamic Equations on a Time Scale: An Introduction with Applications; Springer Science and Business Media: Berlin, Germany, 2012. [Google Scholar]
  4. Meng, H.; Wang, L. Multiple periodic solutions in shifts δ ± for an impulsive functional dynamic equation on time scales. Adv. Differ. Equ. 2014, 152, 1–52. [Google Scholar]
  5. Agarwal, R.P.; Bohner, M.; Peterson, A. Dynamic equations on time scales: A survey. J. Comput. Appl. Math. 2002, 141, 1–26. [Google Scholar] [CrossRef] [Green Version]
  6. Saker, S.H. Oscillation of nonlinear dynamic equations on time scales. Appl. Math. Comput. 2004, 148, 81–91. [Google Scholar] [CrossRef]
  7. Kosmatov, N. Multi-point boundary value problems on time scales at resonance. J. Math. Anal. Appl. 2006, 323, 253–266. [Google Scholar] [CrossRef] [Green Version]
  8. Anderson, D.R. Eigenvalue intervals for a two-point boundary value problem on a measure chain. J. Comput. Appl. Math. 2002, 141, 57–64. [Google Scholar] [CrossRef] [Green Version]
  9. Christiansen, F.B.; Fenchel, T.M. Theories of Populations in Biological Communities, volume 20 of Lecture Notes in Ecological Studies; Springer: Berlin, Germany, 1997. [Google Scholar]
  10. Andreasen, V. Multiple Time Scales in the Dynamics of Infectious Diseases; Springer: Berlin/Heidelberg, Germany, 1978. [Google Scholar]
  11. Bowman, C.; Gumel, A.B.; Dricssche, P. A mathematical model for assessing control strategies against West Nile virus. Bull. Math. Biol. 2005, 67, 1107–1133. [Google Scholar] [CrossRef] [PubMed]
  12. Atici, F.M.; Biles, D.C.; Lebedinsky, A. An application of time scales to economics. Math. Comput. Model. 2006, 43, 718–726. [Google Scholar] [CrossRef]
  13. Dryl, M.; Torres, D. A General Delta-Nabla Calculus of Variations on Time Scales with Application to Economics. Int. J. Dyn. Syst. Differ. 2014, 5, 42–71. [Google Scholar] [CrossRef]
  14. Hovhannisyan, G. Ablowitz–Ladik hierarchy of integrable equations on a time–space scale. J. Appl. Phys. 2014, 55, 102701. [Google Scholar] [CrossRef] [Green Version]
  15. Hovhannisyan, G.; Bonecutter, L.; Mizer, A. On Burgers equation on a time–space scale. Adv. Differ. Equ. 2015, 289, 1–19. [Google Scholar] [CrossRef] [Green Version]
  16. Cieslinski, J.L.; Nikiciuk, T.; Wakiewicz, K. The sine-Gordon equation on time scales. J. Math. Anal. Appl. 2015, 423, 1219–1230. [Google Scholar] [CrossRef]
  17. Cieslinski, J.L. Pseudospherical surfaces on time scales: A geometric definition and the spectral approach. J. Phys. A-Math. Theor. 2007, 40, 12525–12538. [Google Scholar] [CrossRef] [Green Version]
  18. Cieslinski, J.L. New definitions of exponential, hyperbolic and trigonometric functions on time scales. J. Math. Anal. Appl. 2012, 388, 8–22. [Google Scholar] [CrossRef] [Green Version]
  19. Clarkson, P.A.; Tuszynski, J.A. Exact solutions of the multidimensional derivative nonlinear Schrodinger equation for many-body systems of criticality. J. Phys. A 1990, 23, 4269. [Google Scholar] [CrossRef]
  20. Hisakado, M.; Wadati, M. Integrable MultiComponent Hybrid Nonlinear Schrödinger Equations. J. Phys. Soc. Jap. 1995, 64, 408–413. [Google Scholar] [CrossRef]
  21. Kaup, D.J.; Newell, A.C. An exact solution for a derivative nonlinear Schrödinger equation. J. Math. Phys. 1978, 19, 798–801. [Google Scholar] [CrossRef]
  22. Xu, S.; He, J.; Wang, L. The Darboux transformation of the derivative nonlinear Schrödinger equation. J. Phys. A 2011, 44, 6629–6636. [Google Scholar] [CrossRef] [Green Version]
  23. Chen, H.H.; Lee, Y.C.; Liu, C.S. Integrability of Nonlinear Hamiltonian Systems by Inverse Scattering Method. Phys. Scr. 2007, 20, 490. [Google Scholar] [CrossRef]
  24. Ren, H.Y.; Li, D.S. Self-adjointness and Conservation Laws of Chen-Lee-Liu Equation. J. Pingdingshan Univ. 2013, 28, 14–17. [Google Scholar]
  25. Gerdjikov, V.S.; Ivanov, M.I. The quadratic bundle of general form and the nonlinear evolution equations. Bulg. J. Phys. 1983, 10, 35. [Google Scholar]
  26. Fan, E.G. Darboux transformation and soliton-like solutions for the Gerdjikov–Ivanov equation. J. Phys. A 2000, 33, 6925. [Google Scholar] [CrossRef]
  27. Gerdjikov, V.S.; Ivanov, M.I. A quadratic pencil of general type and nonlinear evolution equations. II. Hierarchies of Hamiltonian structures. Bulg. J. Phys. 1983, 10, 130–143. [Google Scholar]
  28. Pei, L.; Li, B. The Darboux transformation of the Gerdjikov–Ivanov equation from non-zero seed. In Proceedings of the 2011 International Conference on Consumer Electronics, Communications and Networks (CECNet), Xianning, China, 16–18 April 2011; pp. 5320–5323. [Google Scholar]
  29. Guo, L.; Zhang, Y.; Xu, S. The higher order Rogu’e Wave solutions of the Gerdjikov–Ivanov equation. Phys. Scr. 2014, 89, 240. [Google Scholar] [CrossRef]
  30. Dai, H.H.; Fan, E.G. Variable separation and algebro-geometric solutions of the Gerdjikov–Ivanov equation. Chaos Soliton. Fract. 2004, 22, 93–101. [Google Scholar] [CrossRef]
  31. Arshed, S. Two Reliable Techniques for the Soliton Solutions of Perturbed Gerdjikov–Ivanov Equation. Optik 2018, 33, 6925. [Google Scholar] [CrossRef]
  32. Nematollah, K.; Hossein, J. Analytical solutions of the Gerdjikov–Ivanov equation by using exp(ϕ(ξ))-expansion method. Optik 2017, 139, 72–76. [Google Scholar]
  33. Cao, C.; Geng, X.; Wang, H. Algebro-geometric solution of the 2+1 dimensional Burgers equation with a discrete variable. J. Math. Phys. 2002, 43, 621–643. [Google Scholar] [CrossRef]
  34. Saburo, K.; Tetsuya, K. Solutions of a derivative nonlinear schrödinger hierarchy and its similarity reduction. Glasgow Math. J. 2005, 47, 99–107. [Google Scholar] [CrossRef] [Green Version]
  35. Fan, E.G. Integrable evolution systems based on Gerdjikov–Ivanov equations, bi-Hamiltonian structure, finite-dimensional integrable systems and N-fold Darboux transformation. J. Math. Phys. 2000, 41, 7769–7782. [Google Scholar] [CrossRef]
  36. Fan, E.G. Bi-Hamiltonian Structure and Liouville Integrability for a Gerdjikov–Ivanov Equation Hierarchy. Chinese Phys. Lett. 2001, 18, 1. [Google Scholar]
  37. Cherniha, R. Galilean-invariant Nonlinear PDEs and their Exact Solutions. J. Nonlinear Math. Phys. 1995, 2, 374–383. [Google Scholar] [CrossRef] [Green Version]
  38. Calogero, F.; Lillo, S.D. The Eckhaus PDE t + ψxx + 2(|ψ|2)xψ + |ψ|4ψ = 0. Inverse Probl. 1999, 3, 633. [Google Scholar] [CrossRef]
Figure 1. One-soliton solution (64) with α 1 = 0.7 , η 1 = 0.6 .
Figure 1. One-soliton solution (64) with α 1 = 0.7 , η 1 = 0.6 .
Axioms 10 00294 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dong, H.; Huang, X.; Zhang, Y.; Liu, M.; Fang, Y. The Darboux Transformation and N-Soliton Solutions of Gerdjikov–Ivanov Equation on a Time–Space Scale. Axioms 2021, 10, 294. https://doi.org/10.3390/axioms10040294

AMA Style

Dong H, Huang X, Zhang Y, Liu M, Fang Y. The Darboux Transformation and N-Soliton Solutions of Gerdjikov–Ivanov Equation on a Time–Space Scale. Axioms. 2021; 10(4):294. https://doi.org/10.3390/axioms10040294

Chicago/Turabian Style

Dong, Huanhe, Xiaoqian Huang, Yong Zhang, Mingshuo Liu, and Yong Fang. 2021. "The Darboux Transformation and N-Soliton Solutions of Gerdjikov–Ivanov Equation on a Time–Space Scale" Axioms 10, no. 4: 294. https://doi.org/10.3390/axioms10040294

APA Style

Dong, H., Huang, X., Zhang, Y., Liu, M., & Fang, Y. (2021). The Darboux Transformation and N-Soliton Solutions of Gerdjikov–Ivanov Equation on a Time–Space Scale. Axioms, 10(4), 294. https://doi.org/10.3390/axioms10040294

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