N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem
Next Article in Journal
Investment Opportunities Identification Based on Macroeconomic Development, the Multiple Criteria Decision Approach
Next Article in Special Issue
On Some Initial and Initial Boundary Value Problems for Linear and Nonlinear Boussinesq Models
Previous Article in Journal
Directional Thermodynamic Formalism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem

College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(6), 826; https://doi.org/10.3390/sym11060826
Submission received: 6 June 2019 / Revised: 17 June 2019 / Accepted: 20 June 2019 / Published: 22 June 2019
(This article belongs to the Special Issue Ordinary and Partial Differential Equations: Theory and Applications)

Abstract

:
In this paper, a kind of nonlinear Schrödinger (NLS) equation, called an NLS-like equation, is Riemann–Hilbert investigated. We construct a 2 × 2 Lax pair associated with the NLS equation and combine the spectral analysis to formulate the Riemann–Hilbert (R–H) problem. Then, we mainly use the symmetry relationship of potential matrix Q to analyze the zeros of det P + and det P ; the N-soliton solutions of the NLS-like equation are expressed explicitly by a particular R–H problem with an unit jump matrix. In addition, the single-soliton solution and collisions of two solitons are analyzed, and the dynamic behaviors of the single-soliton solution and two-soliton solutions are shown graphically. Furthermore, on the basis of the R–H problem, the evolution equation of the R–H data with the perturbation is derived.

1. Introduction

It is known that soliton theory plays an important role in many fields. There are many methods to study soliton equations, of which the inverse scattering method [1,2,3] and the Riemann–Hilbert (R–H) method [4,5,6,7] are two important techniques. The former uses the nonlinear Fourier method [8], in which the calculation procedure is extremely complicated. Conversely, the latter can provide an equivalent and more direct method to solve integrable equations, especially for soliton solutions. Thus, this method has been constantly developed [9,10,11,12,13,14,15]. Furthermore, there are many techniques and transformations for finding exact solutions for soliton equations, such as the Darboux transformation method [16,17,18], the Bäcklund transformation method [19,20], the Hirota bilinear method [21,22,23], the homogeneous balance method [24,25], Frobenius integrable decompositions [26,27,28], and Wronskian technology [29,30]. These methods have greatly promoted the development of soliton theory. From the specific limit of the general soliton solution, lump solutions [31,32,33,34], periodic solutions [35,36], and complex solutions [37,38] can be obtained. In recent years, the initial value problem of integrable equations on the half-line and finite interval [39,40,41] have also been discussed by formulating an associated R–H problem.
As we all know, the soliton solutions of soliton equations with important physical backgrounds have been widely studied. Among them, the nonlinear Schrödinger (NLS) equation is a very significant integrable model in Mathematical physics, which describes water wave theory, nonlinear optics, plasma physics, and so on. It has the following form:
i q t + q x x ± 2 q | q | 2 = 0 .
On the basis of this, an NLS-like equation
i q t + q x x 2 q | q | 2 + 2 a ( | q | 2 ) x q = 0 , a R
is derived as follows. We consider a soliton equation which has the following Lax pair:
Φ x = M Φ , M = i λ 2 q r i λ 2 ,
Φ t = N 0 Φ , N 0 = A B C A ,
where Φ ( x , t , λ ) is a matrix function, A, B, C contain the spectral parameter λ and function q, r and its derivatives. The related stationary zero curvature equation is
N 0 x = [ M , N 0 ] .
Then the Equation (5) becomes
A x = q C B r , B x = 2 i λ 2 B 2 A q , C x = 2 i λ 2 C + 2 A r .
Let us take A, B, C as the six polynomial of λ ,
A = j = 0 6 a j λ j , B = j = 0 6 b j λ j , C = j = 0 6 c j λ j .
Therefore, the Equation (6) has following equivalence relation:
b 6 = c 6 = 0 , a 6 x = 0 , a j x = q c j b j r , ( j = 0 , 1 , 2 , 3 , 4 , 5 ) b j 2 = i 2 b j x + i a j q , c j 2 = i 2 c j x + i a j r .
We choose a 6 = α = const , and have
b 4 = i 2 b 6 x + i a 6 q = i α q , c 4 = i 2 c 6 x + i a 6 r = i α r , a 4 x = q c 4 b 4 r = i q α r i q α r = 0 .
The following equations can be obtained by setting a 4 = β = const :
b 2 = i 2 b 4 x + i a 4 q = i 2 ( i α q ) x + i β q = α 2 q x + i β q , c 2 = i 2 c 4 x + i a 4 r = i 2 ( i α r ) x + i β r = α 2 r x + i β r , a 2 x = q c 2 b 2 r = q ( α 2 r x + i β r ) ( α 2 q x + i β q ) r = α 2 ( q r x + q x r ) = α 2 ( q r ) x .
In addition, we can get a 2 = α 2 q r + a 20 ( a 20 = const ) and
b 0 = i 2 b 2 x + i a 2 q = i 4 α ( q x x + 2 q 2 r ) 1 2 β q x + i q a 20 , c 0 = i 2 c 2 x + i a 2 r = i 4 α ( r x x + 2 q r 2 ) + 1 2 β r x + i r a 20 , a 0 x = q c 0 b 0 r = i 4 α ( q r x r q x ) x + 1 2 β ( q r ) x .
Similarly, a 0 = i 4 α ( q r x r q x ) + 1 2 β q r + a 00 , ( a 00 = const ) can be obtained. We also choose b 5 = c 5 = 0 , through the same steps as above and have a 5 = 0 , b 3 = b 1 = 0 , c 3 = c 1 = 0 . Thus,
A = a 6 λ 6 + a 5 λ 5 + a 4 λ 4 + a 3 λ 3 + a 2 λ 2 + a 1 λ 1 + a 0 λ 0 = i 4 α ( q r x r q x ) + 1 2 β q r + a 00 + ( α 2 q r + a 20 ) λ 2 + β λ 4 + α λ 6 , B = b 6 λ 6 + b 5 λ 5 + b 4 λ 4 + b 3 λ 3 + b 2 λ 2 + b 1 λ 1 + b 0 λ 0 = i 4 α ( q x x + 2 q 2 r ) 1 2 β q x + i q a 20 + ( α 2 q x + i β q ) λ 2 + ( i α q ) λ 4 , C = c 6 λ 6 + c 5 λ 5 + c 4 λ 4 + c 3 λ 3 + c 2 λ 2 + c 1 λ 1 + c 0 λ 0 = i 4 α ( r x x + 2 q r 2 ) + 1 2 β r x + i r a 20 + ( α 2 r x + i β r ) λ 2 + ( i α r ) λ 4 .
The matrix spectral problem can be obtained by taking α = a 00 = a 20 = 0 , β = 2 i , r = q ¯ ,
M = i λ 2 q q ¯ i λ 2 , N 0 = 2 i λ 4 + i | q | 2 i q x + 2 λ 2 q i q ¯ x 2 λ 2 q ¯ 2 i λ 4 i | q | 2 .
A direct calculation shows q t i q x x 2 i q | q | 2 = 0 . Notice that the term 2 i q | q | 2 is independent of λ , thus, it can become positive as we assume
N 1 = 2 i λ 4 i | q | 2 i q x + 2 λ 2 q i q ¯ x 2 λ 2 q ¯ 2 i λ 4 + i | q | 2 ,
and using the same method as above, we have
N = 2 i λ 4 i | q | 2 + A 1 i q x + 2 λ 2 q + B 1 i q ¯ x 2 λ 2 q ¯ + C 1 2 i λ 4 + i | q | 2 A 1 .
The corresponding zero curvature equation is M t N x + [ M , N ] = 0 , and we get
q t i q x x B 1 x + 2 i q | q | 2 2 A 1 q 2 i λ 2 B 1 = 0 , q ¯ t i q ¯ x x C 1 x + 2 i q ¯ | q | 2 2 A 1 q ¯ + 2 i λ 2 C 1 = 0 .
Taking B 1 = C 1 = 0 , A 1 = i a ( | q | 2 ) x , the 2 × 2 Lax pair can be obtained:
Φ x = M Φ , M = i λ 2 q q ¯ i λ 2 ,
Φ t = N Φ , N = 2 i λ 4 + i a | q | x 2 i | q | 2 i q x + 2 λ 2 q i q ¯ x 2 λ 2 q ¯ 2 i λ 4 i a | q | x 2 + i | q | 2 ,
where the symbol represents complex conjugation. In our analysis, we assume the potential q is smooth enough and decays rapidly to zero when x ± . Furthermore, at any later time t, we look for solution q ( x , t ) with the initial condition q ( x , 0 ) . When setting a = 0 , Equation (2) becomes a classical nonlinear Schrödinger Equation (1).
In this paper, we study the perturbation theory of the NLS-like equation. Obviously, small perturbations of integrable conditions can be regarded as perturbations of integrable models. Our formalism is in view of the R–H problem related to the NLS-like equation. The main advantage of the proposed method is its algebraic property, which is different from the method using the Gel’Fand Leviaon integral equation [42]. The R–H problem has many applications in dealing with disturbed soliton dynamics [43]. Modern versions of perturbation theory for the R–H problem have been published in a series of papers [44,45,46,47]. The direct perturbation theory is another form of soliton perturbation theory, which develops on the basis of the perturbation solution expansion into the square eigenfunction of the linearized soliton equation [48].
The main structure of this article is as follows. In Section 2, we give the Lax pair of the NLS-like equation. Then, the properties of the equivalent space matrix spectral problem of matrix eigenfunctions are analyzed, and the R–H problem related to the newly introduced space matrix spectral problem is formulated. In Section 3, through the special reductive R–H problem, in which the jump matrix is a unit matrix, the explicit expressions of the N-soliton solutions are obtained. In addition, the single soliton solution and collisions of the two-soliton solutions are analyzed. The perturbation theory based on Section 2 and the evolution equation of R–H data with perturbation are given in Section 4. Finally, the paper is summarized and further questions are given in Section 5.

2. The Riemann–Hilbert Problem

In what follows, we set a = 1 , and constructed a R–H formulation for Equation (2) with scattering and inverse scattering methods.
Let T = ( t 1 , t 2 ) T be a solution of Lax pair (8) and (9) and the following relation can be obtained by defining μ = t 1 / t 2 :
( ln t 2 ) x = i λ 2 q ¯ μ ,
( ln t 2 ) t = ( 2 i λ 4 i | q | x 2 + i | q | 2 ) + ( i q ¯ x 2 λ 2 q ¯ ) μ .
By taking the derivative of the right-hand side of these two equations with respect to t and x, we get
( i λ 2 q ¯ μ ) t = [ ( 2 i λ 4 i | q | x 2 + i | q | 2 ) + ( i q ¯ x 2 λ 2 q ¯ ) μ ] x .
Throughout this work, we consider Φ in (8)–(9) to be a fundamental matrix. In addition, q 0 as x ± and we get Φ e i λ 2 σ 3 x 2 i λ 4 σ 3 t , where σ 3 is a diagonal constant matrix,
σ 3 = diag ( 1 , 1 ) .
It is convenient to introduce a new matrix spectral function J = J ( x , t ; λ ) , which can be defined as
Φ = J e i λ 2 σ 3 x 2 i λ 4 σ 3 t .
Hence, the variable x and t in the new matrix J are independent at infinity. By inserting Equation (14) into (8)–(9), the original Lax pair (8)–(9) can be rewritten as
J x = i λ 2 [ σ 3 , J ] + Q J ,
J t = 2 i λ 4 [ σ 3 , J ] + V J ,
with
Q = 0 q q ¯ 0 , V = 2 λ 2 Q + i | q | x 2 i | q | 2 i q x i q ¯ x i | q | x 2 + i | q | 2 .
and [ σ 3 , J ] = σ 3 J J σ 3 . From Equation (17), it can be seen that t r ( Q ) = t r ( V ) = 0 , and
Q = Q , V = V ,
where the superscript represents the Hermitian of a matrix. In the scattering process, we start from the x-part of the Lax pair and regard t as a parameter and omit it.
Let us first introduce the Jost solutions, which have the following asymptotic properties:
J ± ( x , λ ) I , x ± ,
where I is a 2 × 2 identity matrix, and the subscript of J ± indicates that the boundary conditions are at ± , respectively. By using the large-x asymptotic condition (19), the x part of Equation (15) can be transformed into the Volterra integral equation of J ±
J ( x , λ ) = I + x e i λ 2 σ 3 ( y x ) Q ( y ) J ( y , λ ) e i λ 2 σ 3 ( x y ) d y ,
J + ( x , λ ) = I x e i λ 2 σ 3 ( y x ) Q ( y ) J ( y , λ ) e i λ 2 σ 3 ( x y ) d y .
Through the direct analysis of Equations (20) and (21), because of the structure of the potential Q in Equation (17), it can be seen that the first column of J contains only the exponential factor e i λ 2 ( x y ) , as λ C + = { λ | arg λ ( 0 , π 2 ) ( π , 3 π 2 ) } , since y < x , e i λ 2 ( x y ) decays. In addition, the second column of J + includes only exponential factor e i λ 2 ( y x ) , as λ C + , since y > x , e i λ 2 ( y x ) also decays. Thus, we believe these two columns can be analytic for λ C + and continuous for λ C + R i R . By using a similar analysis, the second column of J as well as the first column of J + can also be analytic for λ C = { λ | arg λ ( π 2 , π ) ( 3 π 2 , 2 π ) } and continuous for λ C R i R .
Remark 1.
The implementation of the inverse scattering transform of the NLS-like equation is different to the implementation of the derivative NLS equation. The differences are shown below:
(i) One needs to see the continuity of λ 2 differently;
(ii) It is necessary to distinguish between the upper half plane and the lower half plane of λ 2 .
Able’s formula tells us that
det Φ ( x ) = det Φ ( x 0 ) e x 0 x t r A ( ζ ) d ζ ,
and applying this identity to the Equation (8) and using relation (14), we see that no matter what the value of x is, det J ( x , λ ) is a constant. Then, we obtain det J ± ( x , λ ) = 1 by utilizing the boundary conditions (19). Thus, denoting
E ( x , λ ) = e i λ 2 σ 3 x ,
and
φ J E , Ψ J + E .
In fact, two solutions of the Equation (8), φ ( x , λ ) and Ψ ( x , λ ) , are linearly related. Their relationship can be stated as follows:
φ ( x , λ ) = Ψ ( x , λ ) S ( λ ) , λ R i R .
That is,
J E = J + E S ( λ ) , λ R i R ,
where,
S ( λ ) = s 11 s 12 s 21 s 22 ,
which is called the scattering matrix. In view of Equation (26) and det J ± ( x , λ ) = 1 , we notice that the scattering matrix S ( λ ) satisfies
det ( S ( λ ) ) = 1 .
Thus, defining ( φ , Ψ ) as a collection of columns, which can be read as
φ = ( φ 1 , φ 2 ) , Ψ = ( ψ 1 , ψ 2 ) .
In addition, when the Jost solutions are
P + = ( φ 1 , ψ 2 ) e i λ 2 σ 3 x = J + E S + E 1 = J E S E 1 = J H 1 + J + H 2 ,
which are analytic in λ C + , where the matrices S ± ,
S + ( λ ) = s 11 0 s 21 1 , S ( λ ) = 1 s 21 * 0 s 22 * ,
and H 1 = diag ( 1 , 0 ) , H 2 = diag ( 0 , 1 ) . The matrices S ± provide a factorization of the scattering matrix S S = S + . In addition, when the Jost solutions are
( ψ 1 , φ 2 ) e i λ 2 σ 3 x = J + H 1 + J H 2 ,
which are analytic in λ C . Furthermore, using the Volterra integral Equations (20)–(21), the asymptotic properties of these analytic functions at large- λ can be obtained:
P + ( x , λ ) I , λ C + ,
and
( ψ 1 , φ 2 ) e i λ 2 σ 3 x I , λ C .
In order to get the corresponding analysis of P + in C , the adjoint equation of Equation (15) is introduced:
K x = i λ 2 [ σ 3 , K ] K Q .
As a result of the relationship
( J J 1 ) x = 0 = J x J 1 + J ( J 1 ) x ,
and the scattering Equation (15), we have
( J 1 ) x = i λ 2 [ σ 3 , J 1 ] J 1 Q ,
so that J ± 1 satisfies the adjoint scattering Equation (34). Using the same technique as above, the collection of rows φ 1 and Ψ 1 can be expressed as
φ 1 = φ ˜ 1 φ ˜ 2 , Ψ 1 = ψ ˜ 1 ψ ˜ 2 .
It can be seen that the adjoint Jost solutions are
P = e i λ 2 σ 3 x φ ˜ 1 ψ ˜ 2 = E T + E 1 J + 1 = E T E 1 J 1 = H 1 J 1 + H 2 J + 1 ,
which are analytic in λ C , where
T + ( λ ) = s 11 * s 21 * 0 1 , T ( λ ) = 1 0 s 21 s 22 .
The matrices T ± also provide a factorization of the scattering matrix: T + S = T . In addition, when the adjoint Jost solutions are
e i λ 2 σ 3 x ψ ˜ 1 φ ˜ 2 = H 1 J + 1 + H 2 J 1 ,
which are analytic in λ C + . Taking the Volterra integral equation again, we have
P ( x , λ ) I , λ C
and
e i λ 2 σ 3 x ψ ˜ 1 φ ˜ 2 I , λ C + .
The analytical properties of the above Jost solutions are summarized below:
φ = ( φ 1 + , φ 2 ) , Ψ = ( ψ 1 , ψ 2 + ) ,
φ 1 = ( φ ˜ 1 , φ ˜ 2 + ) T , Ψ 1 = ( ψ ˜ 1 + , ψ ˜ 2 ) T .
Here, the superscript “±” represents that the basic quantity is analyzed in C ± . Obviously, we know the analytic properties of the Jost solutions, hence the analytic properties of the scattering matrix S ( λ ) can be easily analyzed. Because of the relation
S = Ψ 1 φ = ψ ˜ 1 + ψ ˜ 2 ( φ 1 + , φ 2 ) ,
and
S 1 = φ 1 Ψ = φ ˜ 1 φ ˜ 2 + ( ψ 1 , ψ 2 + ) ,
the analysis structure of the scattering matrices S and S 1 can be obtained, which is expressed as
S = s 11 + s 12 s 21 s 22 , S 1 = s ˜ 11 s ˜ 12 s ˜ 21 s ˜ 22 + .
According to the relationship between S 1 and S, the following equations are obtained:
s ˜ 11 = s 22 , s ˜ 22 = s 11 , s ˜ 12 = s 12 , s ˜ 21 = s 21 .
So far, the matrix functions P + ( x , λ ) and P ( x , λ ) are constructed, which are analytic in C + and C , respectively. By defining
G + ( x , λ ) = lim μ C + , μ λ P + ( x , μ ) , ( G ) 1 ( x , λ ) = lim μ C , μ λ P ( x , μ ) , λ R i R ,
the two matrix functions G + and G are related to each other by using Equations (26), (30), (38), and (48):
G + ( x , λ ) ( G ) 1 = G ( x , λ ) ,
where
G = E S + T + E 1 = E S T E 1 = E 1 s ˜ 12 s 21 1 E 1 .
Equation (49) accurately gives the R–H problem of a correlation matrix. From Equations (32) and (40), the asymptotic properties of the above R–H problem shows
P ± ( x , λ ) I , λ C ± ,
and the canonical normalization condition
G ± ( x , λ ) I , λ R i R .
We know that a key step for solving soliton solutions is to calculate the potential matrix Q through P ± ( x , λ ) . In view of P + being the solution of the scattering problem (15), we expand the P + at large- λ as
P + ( x , λ ) = I + λ 1 P 1 + ( x ) + λ 2 P 1 + ( x ) + o ( λ 3 ) , λ ,
and by taking Equation (53) into Equation (15) and comparing the term of o ( 1 ) , we have
Q = i [ σ 3 , P 1 + ] = i 0 2 ( P 1 + ) 12 2 ( P 1 + ) 21 0 .
Therefore, the reconstructed solution q can be represented by P + as
q = 2 i ( P 1 + ) 12 .
At this point, the inverse scattering process has been completed.
Similarly, we obtain
diag ( P 1 + ) x = diag ( Q P 1 + ) .
Through the large-x asymptotic of P + ( x , λ ) from Equation (19) and Equations (55)–(56), we find the full matrix P 1 + ( x , λ ) can be expressed as
P 1 + ( x ) = 1 2 i x | q ( y ) | 2 d y q ( x ) q ¯ ( x ) x | q ( y ) | 2 d y .
By the same method, we can get the asymptotic expansion of P .
It is well known that soliton solutions for the R–H problem with zeros can be obtained by transforming them into a problem without zeros. As long as the det P ± is specified at the zeros in C ± , and the structure of the ker P ± at these zeros can be determined, then the uniqueness of the solution for each of the associated R–H problems defined by the Equations (49)–(50) is established. From the definitions of Equations (30) and (38) as well as the scattering relation (26), we get
det P + ( x , λ ) = s ˜ 22 ( x , λ ) = s 11 ( x , λ ) ,
det P ( x , λ ) = s 22 ( x , λ ) = s ˜ 11 ( x , λ ) .
Let N be an arbitrary nature, we assume that s ˜ 22 has N zeros { λ k C + , 1 K N } and s 22 has N zeros { λ ^ k C , 1 K N } . To get N- Solitons, we suppose that all zeros λ k and λ ^ k are simple zeros. In this case, each of ker P + ( λ k ) , 1 K N , which includes only a single column vector v k ; each of ker P ( λ ^ k ) , 1 K N , which includes only a single row vector v ^ k . That is,
P + ( λ k ) v k = 0 , v ^ k P ( λ ^ k ) = 0 , 1 k N .
The potential matrix Q possesses a symmetry property (18), which yields a symmetry property in the scattering matrix and Jost functions. In addition, we notice that the scattering Equation (15) has the Hermitian property, then by utilizing the anti-Hermitian property of the first equation in (18), we have
J x = i λ ¯ 2 [ σ 3 , J ] J Q .
Therefore, Equation (61) shows that J ± ( x , λ ¯ ) satisfies the adjoint scattering Equation (34). From Equation (35), we know that J ± 1 ( x , λ ) also satisfies the adjoint equation. Therefore, J ± ( x , λ ¯ ) and J ± 1 ( x , λ ) must be linearly related to each other. That is, J ± ( x , λ ¯ ) = A J ± 1 ( x , λ ) , where A is x-independent. We can get A = 1 by utilizing the boundary conditions (19) of Jost solutions J ± . That is,
J ± ( x , λ ¯ ) = J ± 1 ( x , λ ) .
By utilizing this involution property and definitions Equation (30) as well as (38) for P ± , we see that the analytic solutions P ± also possess the involution property:
( P + ) ( λ ¯ ) = P ( λ ) .
In addition, from the scattering relationship (26) between J + and J , the involution property is also suitable for S:
S ( λ ¯ ) = S 1 ( λ ) .
Considering the zeros of the scattering coefficients s ˜ 22 ( λ ) and s 22 ( λ ) are λ k and λ ^ k , respectively; the involution relation from the involution property (64) shows
λ ^ k = λ ¯ k .
In order to get the symmetric properties of the eigenvectors v ^ k and v k , we use the Hermitian of the equation P + ( λ k ) v k = 0 , and take the involution properties Equations (63) and (65). Thus, we have
v k P ( λ ^ k ) = 0 .
Then, we compare it with the equation v ^ k P ( λ ^ k ) = 0 , and see that eigenvectors ( v k , v k ^ ) have the involution property
v ^ k = v k .
To obtain soliton solutions in the R–H problems above, we set G = I . When we set s 21 = s ˜ 12 = 0 , this means that the reflection does not exist in the scattering problem. By factoring a rational matrix Γ ( λ ) , the solution of the non-regular R–H problem with zeros is
P + ( λ ) = P ˜ + ( λ ) Γ ( λ ) , P ( λ ) = Γ 1 ( λ ) P ˜ ( λ ) .
P ˜ ± is the solution to the following regularized R–H problem:
P ˜ ( λ ) P ˜ + ( λ ) = Γ ( λ ) G ( λ ) Γ 1 , λ R i R ,
and P ˜ ± I as λ .
The rational matrix functions Γ and Γ 1 , which are defined as
Γ ( λ ) = I + k , l = 1 N v k ( M 1 ) k l v ^ l λ λ ^ l , Γ 1 ( λ ) = I k , l = 1 N v k ( M 1 ) k l v ^ l λ λ k ,
where
M k l = v ^ k v l λ ^ k λ l , 1 k , l N ,
and det Γ ( λ ) = Π l = 1 N λ λ l λ λ ^ l . Γ ( λ ) and Γ 1 ( λ ) have the same zeros as P + ( λ ) and P ( λ ) , respectively, as well as the null spaces:
Γ ( λ k ) v k = 0 , v ^ k Γ 1 ( λ ^ k ) = 0 .
Because the zeros λ k and λ ^ k are contants, i.e., they do not rely on spatial variable x and time variable t, it is easy to determine the temporal and spatial evolution of vectors v k ( x , t ) and v ^ k ( x , t ) ( 1 k N ) in ker P ± . We calculate the derivatives of x on both sides of equation P + v k = 0 . By utilizing Equation (15), we have
P + ( x , λ k ) d v k d x + i λ k 2 σ 3 v k = 0 , 1 k N .
Thus, we can draw a conclusion that for 1 k N , the vector d v k d x + i λ k 2 σ 3 v k must be in the kernel of P + ( x , λ k ) and it must be a scalar function of the vector. We set the constant vanishes and have
d v k d x = i λ k 2 σ 3 v k , 1 k N .
In a similar way, the time dependency of v k can be obtained by utilizing (16):
d v k d t = 2 i λ k 4 σ 3 v k , 1 k N .
Finally, we can get
v k ( x , t ) = e i λ k 2 σ 3 x 2 i λ k 4 σ 3 t v k 0 , 1 k N ,
v ^ k ( x , t ) = v ^ k 0 e i λ ¯ k 2 σ 3 x + 2 i λ ¯ k 4 σ 3 t , 1 k N ,
where v k 0 ( 1 k N ) is an arbitrary constant number column vector, and v ^ k 0 ( 1 k N ) is an arbitrary constant row vector.

3. The N-Soliton Solutions and Their Danamics

By using the relation P + in Equations (53) and (68) as well as the reconstructed potential Q in Equation (54), we get the N-soliton solutions of Equation (2), which are expressed as
P 1 + = k , l = 1 N v k ( M 1 ) k l v ^ l
and
q ( x , t ) = 2 i ( P 1 + ) 12 = 2 i k , l = 1 N v k ( M 1 ) k l v ^ l 12 .
Here, vectors v k are given by Equation (74), v ^ k = v k . In addition, matrix M has been defined in Equation (71). Without loss of generality, let v l 0 = ( c l , 1 ) T . Therefore, the solution q can be expressed explicitly as
q ( x , t ) = 2 i k , l = 1 N c k e θ k θ ¯ l ( M 1 ) k l ,
where the N × N matrix M is
M k l = 1 λ ¯ k λ l e ( θ ¯ k + θ l ) + c ¯ k c l e θ ¯ k + θ l ,
where θ k = i λ k 2 x 2 i λ k 4 t .
In order to simplify the calculation, solution q can also be expressed by matrix determinants [48]:
q ( x , t ) = 2 i det F det M ,
where
F = 0 e θ ¯ 1 e θ ¯ N c 1 e θ 1 M 11 M 1 N c N e θ N M N 1 M N N .

3.1. Single-Soliton Solutions

When N = 1 , the solution q ( x , t ) is
q ( x , t ) = 2 i ( λ ¯ 1 λ 1 ) c 1 e θ 1 θ ¯ 1 e ( θ 1 + θ ¯ 1 ) + | c 1 | 2 e θ 1 + θ ¯ 1 .
We assume that
λ 1 = ζ + i τ , c 1 = e 4 ζ τ x 0 + i σ 0 ,
here, ζ is real part of λ 1 , and τ is imaginary part of λ 1 . In addition, x 0 , σ 0 are real parameters; thus, the solution of Equation (79) can be expressed as
q ( x , t ) = 2 τ sech { 4 ζ τ [ x + 4 ( ζ 2 τ ζ τ 2 ) t ] x 0 } exp 2 i ( ζ 2 τ 2 ) x 4 i ( ζ 4 + τ 4 6 ζ 2 τ 2 ) t + i σ 0 } .
Notice that the shape of the amplitude function | q ( x , t ) | is a hyperbolic secant, and its peak amplitude is 2 τ and the velocity is 4 ( ζ 2 τ ζ τ 2 ) . We can see the soliton’s peak amplitude only relies on τ , consequently, the peak cannot change after the soliton collisions. The phase of the solution depends not only linearly on spatial x but also on t. In addition, the parameter x 0 represents the initial position of the solitary wave, and σ 0 represents the phase of the solitary wave. We choose the appropriate parameter and give the evolution characteristics of single soliton solutions in Figure 1.
Remark 2.
Because of the difference in Lax pairs, the corresponding R–H problem of the spatial matrix spectral problem is also different. We can clearly see through the solution (85) that the NLS-like equation is different from the single soliton solution of the derivative NLS equation, which also makes the pulse width and velocity of the NLS-like equation corresponding to the derivative NLS equation different.

3.2. Two-Soliton Solutions

When N = 2 , the two-soliton solutions can also be written out explicitly; however, it is quite complicated. For simplicity, we take N = 2 into Equation (81) and get
q ( x , t ) = 2 i c 1 e θ 1 θ ¯ 1 M 22 + c 2 e θ 2 θ ¯ 1 M 12 | M | × c 1 e θ 1 θ ¯ 2 M 21 c 2 e θ 2 θ ¯ 2 M 11 | M | ,
where
M 11 = e ( θ 1 + θ ¯ 1 ) + | c 1 | 2 e θ 1 + θ ¯ 1 λ ¯ 1 λ 1 , M 12 = e ( θ 2 + θ ¯ 1 ) + c ¯ 1 c 2 e θ 2 + θ ¯ 1 λ ¯ 1 λ 2 , M 21 = e ( θ 1 + θ ¯ 2 ) + c ¯ 2 c 1 e θ 1 + θ ¯ 2 λ ¯ 2 λ 1 , M 22 = e ( θ 2 + θ ¯ 2 ) + | c 2 | 2 e θ 2 + θ ¯ 2 λ ¯ 2 λ 2 .
Then, we show the collision of the two-soliton in Figure 2 and Figure 3. One for the case of ζ 1 ζ 2 and the other one for the case of ζ 1 = ζ 2 . Here,
λ k = ζ k + i τ k , k = 1 , 2 ,
where ζ k is real part of λ k and τ k is imaginary part of λ k .
Case I. We set λ 1 = 1 + 0.75 i , λ 2 = 0.65 + 0.55 i .
In this case, we assume that ζ 1 > ζ 2 , which indicates that soliton-2 is on the left side of soliton-1 as t as well as moves fast. After collision, their position and phase will be scattered. Through asymptotic analysis, we can explain this change.
When t , by simple calculation, the asymptotic state of Equation (86) can be expressed as
q ( x , t ) 2 i ( λ ¯ 1 λ 1 ) c 1 e θ 1 θ ¯ 1 e ( θ 1 + θ ¯ 1 ) + | c 1 | 2 e θ 1 + θ ¯ 1 , t ,
where c 1 = c 1 λ 1 λ 2 λ 1 λ ¯ 2 and the asymptotic solution has the same peak amplitude 2 τ 1 and velocity 4 ( ζ 1 2 τ 1 ζ 1 τ 1 2 ) as Equation (86).
When t + , the asymptotic state of Equation (86) also can be expressed as
q ( x , t ) 2 i ( λ ¯ 1 λ 1 ) c 1 + e θ 1 θ ¯ 1 e ( θ 1 + θ ¯ 1 ) + | c 1 + | 2 e θ 1 + θ ¯ 1 , t + ,
where c 1 + = c 1 λ 1 λ ¯ 2 λ 1 λ 2 . The asymptotic solution also has same peak amplitude 2 τ 1 and velocity 4 ( ζ 1 2 τ 1 ζ 1 τ 1 2 ) as Equation (86). As Figure 2 shows, the solution does not change their velocity and shape after collision, but the initial positions and phases of solitary waves have shifted.
The position shift is
Δ x 01 = 1 4 ζ 1 τ 1 ( ln | c 1 + | ln | c 1 | ) = 1 2 ζ 1 τ 1 ln λ 1 λ 2 λ 1 λ ¯ 2 .
It is easy to see Δ x 01 < 0 , because λ k C + . Then, the phase shift is
Δ σ 01 = arg ( c 1 + ) arg ( c 1 ) = 2 arg λ 1 λ 2 λ 1 λ ¯ 2 .
Similarly, as t ± , the asymptotic solutions have the same peak amplitude 2 τ 2 and velocity 4 ( ζ 2 2 τ 2 ζ 2 τ 2 2 ) with single solitons, where c 2 = c 1 λ 2 λ ¯ 1 λ 2 λ 1 , c 2 + = c 1 λ 2 λ 1 λ 2 λ ¯ 1 . Therefore, the second soliton position shift is
Δ x 02 = 1 2 ζ 2 τ 2 ln λ 2 λ 1 λ 2 λ ¯ 1 > 0 ,
and the phase shift is
Δ σ 02 = 2 arg λ 2 λ 1 λ 2 λ ¯ 1 .
From Equations (90) and (92), we get
Δ x 02 Δ x 01 = ζ 1 τ 1 ζ 2 τ 2 .
This shows that the position deviation of each soliton is inversely proportional to its amplitude.
Case II. We set λ 1 = 0.8 + 0.55 i , λ 2 = 0.8 + 0.45 i .
In the second case, we assume that ζ 1 = ζ 2 . As Figure 3 shows, the two solitons have the same velocity, and the amplitude function | q ( x , t ) | has periodic oscillations with time.

4. Evolution of the R–H Data in the Perturbed NLS-Like Equation

In this part, a perturbed NLS-like equation
i q t + q x x 2 q | q | 2 + 2 a ( | q | 2 ) x q = R ( q )
is analyzed, where Equation (95) is called a nearly integrable system. Here, R is a perturbation term, and ε 1 . We give the symbol δ δ t to the perturbations in order to distinguish the integrable and perturbation contributions. Consequently,
i δ Q δ t = R , R = 0 R R ¯ 0 .
Generally, a perturbation leads to a slow evolution of the R–H data. In fact, a perturbation causes the variational δ Q of the potential of the spectral Equation (15), which leads to the variation of the Jost solutions,
δ J ± x = i λ 2 [ σ 3 , δ J ± ] + δ Q J ± + Q δ J ± .
By solving the equation, we have
δ J ± = J ± E ( ± x d x E 1 J ± 1 δ Q J ± E ) E 1 .
Hence, by utilizing the Equations (26), (30), and (38), we get a variation of a scattering matrix:
δ S δ t = i ε S + d x E 1 ( P + ) 1 R P + E S 1 = i ε T + 1 d x E 1 P R ( P ) 1 E T .
Here, S ± and T ± are the matrices defined in Section 2. Notice that the analytic solution P ± is naturally installed in the equation. We can thus denote that
Y + ( a , b ) = a b d x E 1 ( P + ) 1 R P + E , Y ( a , b ) = a b d x E 1 P R ( P ) 1 E , Y ± ( λ ) Y ± ( , ) .
Then,
δ S δ t = i ε S + Y + ( λ ) S 1 = i ε T + 1 Y ( λ ) T .
Matrix Y ± is interrelated through matrix G into the R–H problem (49):
Y ( λ ) = G Y + ( λ ) G 1 .
From Equations (30) and (38), the variations of the analytic solutions show
δ P + δ t = i ε P + E H + E 1 , δ P δ t = i ε E H E 1 P ,
where H + are the evolution functionals, which are defined as
H + = Y + ( λ ) H 1 Y + ( x , ) , H = H 1 Y ( λ ) Y ( x , ) .
From Equation (63), we get the connection between H + and H , H = H + , λ R i R . They contain all the basic information about a perturbation. The additional terms obtained from the perturbation evolution equation P ± are defined by H ± ,
( P + ) t = 2 i k 4 [ σ 3 , P + ] + V P + i ε P + E H + E 1 , ( P ) t = 2 i k 4 [ σ 3 , P ] P V + i ε E H E 1 P .
In addition, the evolution equation for the matrix G of the R–H problem also has the form
G t = 2 i k 4 [ σ 3 , G ] i ε ( G H + H G ) .
Indeed, the equation provides the evolution of the continuous R–H data.
Next, we derive the perturbation induced evolution equation for the discrete R–H data, i.e., for the zero λ k and the eigenvector v k . Vectors v K = ( v K 1 , v K 2 ) T are constant without perturbation. Under perturbation, vectors v K = ( v K 1 , v K 2 ) T have slow t dependence. Let us start with the equation
P + ( λ k ) v k = P + ( λ k ) e ( i λ k 2 x 2 i d t λ k 4 ) σ 3 v P = 0
which is unaffected by a perturbation. Here, the integral in the exponential takes account of the time dependence of the zero λ k caused by the possible perturbation. Taking the total derivative in t, we get
( P + ( λ ) e i λ 2 σ 3 x 2 i d t λ k 4 σ 3 ) t + ( P + ( λ ) e i λ 2 σ 3 x 2 i d t λ k 4 σ 3 ) λ λ t λ = λ k + P + ( λ k ) e i λ k 2 σ 3 x 2 i d t λ k 4 σ 3 ( v P ) t = 0 .
The first term with ( P + ) t is given by Equation (98), which includes the evolution functional H + . Note that the evolution functional H + ( λ ) is defined by Y + in Equation (97) which depend on ( P + ) 1 , conversely. Thus, the function H + , which has the simple pole in λ k , is meromorphic in C + , where P + has zero,
H + ( λ ) = H + ( reg ) ( λ ) + 1 λ λ k Res [ H + ( λ ) , λ k ] ,
where H + ( Res ) represents the regular part of H + in the point λ k . The perturbation evolution of the vector v K is given by
( v K ) t = i ε e 2 i d t λ k 4 σ 3 H + ( reg ) ( λ k ) e 2 i d t λ k 4 σ 3 v K .
In order to derive the evolution equation of λ k , we start with the equation det P + ( λ k ) = 0 . Taking the total derivative in t, we get
[ ( det P + ( λ ) ) t ] λ = λ k + [ ( det P + ( λ ) ) λ ] λ = λ k ( λ k ) t = 0 .
From Equations (68) and (70), we have
det P + ( λ ) = λ λ k λ λ ¯ k det P ˜ + ( λ ) ,
where det P ˜ + ( λ ) 0 . From Equations (73)–(74), we get
v k = exp ( θ k σ 3 ) v k 0 = exp ( θ k + a k + i σ 0 k ) e θ k .
Here, θ k = i λ k 2 x 2 i λ k 4 t , and defined exp ( a k + i σ 0 k ) = ( v k 0 ) 1 / ( v k 0 ) 2 , with a k and σ 0 k are real constants. Denote [ det P + ( λ ) ] t = i ε tr H + det P + ( λ ) , we finally get a simple evolution equation for spectral λ k , a k , σ 0 k ,
( λ k ) t = i ε Res [ Y + 22 ( λ ) , λ k ] ,
and
( a k + i σ 0 k ) t = Y 22 ( reg ) ( λ k ) exp 4 i 0 t d t λ k 4 a k i σ 0 k Y 12 ( reg ) ( λ k ) .
Perturbation-induced evolution of R–H data is determined by Equations (99), (102), and (103). Notice that these equations are accurate because we have not mentioned any small part of ε . In addition, these equations cannot be applied directly, because Y ± goes in and depends on the unknown solutions P ± for the spectral problem of the perturbed potential Q.

5. Conclusions

In this paper, an NLS-like equation associated with a 2 × 2 Lax pair is studied. We start from the spectral analysis of the Lax pair of Equation (2). By using the R–H method, when the scattering coefficients vanish, the regularization condition at the infinity on the real axis can be used to solve the corresponding R–H problem. When the jump matrix G is the unit matrix, the N-soliton solutions of the integrable equation can be obtained. The R–H method is a very useful tool, especially for soliton solutions. As we all know, the R–H approach has been widely used to solve many integrable equations, for example, the generalized Sasa–Satsuma equation [14], the general coupled nonlinear Schrödinger equation [49], and the Harry–Dym equation [50]. Furthermore, based on the R–H problem, the evolution functional is derived and the R–H data in the perturbed NLS-like equation is obtained. Perturbation theory has many applications, such as propagation of arbitrarily polarized optical pulses in optical fibers, multicomponent soliton equations, and soliton pulses in various Bose–Einstein condensations.
It is very meaningful to study the exact solutions and other types of integrable equations, and to analyze perturbation theory based on R–H problems. Further research on how to apply the R–H problem to the generalized integrable correspondence equation combined with perturbation theory will be one of our future considerations.

Author Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Funding

This research was funded by Natural Science Foundation of Shandong Province (Grant No. ZR2019QD018), National Natural Science Foundation of China (Grant No. 61602188), Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (Grant No. 2017RCJJ068, 2017RCJJ069).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ablowitz, M.J.; Kaup, D.J.; Newell, A.C.; Segur, H. The inverse scattering Transform-Fourier analysis for nonlinear problems. Stud. Appl. Math. 2015, 53, 249–315. [Google Scholar] [CrossRef]
  2. Ji, J.L.; Zhu, Z.N. Soliton solutions of an integrable nonlocal modified Korteweg-de Vries equation through inverse scattering transform. J. Math. Anal. Appl. 2017, 453, 973–984. [Google Scholar] [CrossRef]
  3. Maeda, M.; Sasaki, H.; Segawa, E.; Suzuki, A.; Suzuki, K. Scattering and inverse scattering for nonlinear quantum walks. Discrete Contin. Dyn. Syst. 2018, 38, 18. [Google Scholar] [CrossRef]
  4. Shchesnovich, V.S. The soliton perturbation theory based on the Riemann-Hilbert spectral problem. Chaos Solitons Fractals 1995, 5, 2121–2133. [Google Scholar] [CrossRef]
  5. Kaup, D.J.; Yang, J.K. The inverse scattering transform and squared eigenfunctions for a degenerate 3 × 3 operator. Inverse Probl. 2009, 25, 105010–105021. [Google Scholar] [CrossRef]
  6. Guo, B.L.; Ling, L.M. Riemann-Hilbert approach and N-soliton formula for coupled derivative Schrödinger equation. J. Math. Phys. 2012, 53, 73506. [Google Scholar] [CrossRef]
  7. Wang, D.S.; Ma, Y.Q.; Li, X.G. Prolongation structures and matter-wave solitons in F = 1 spinor Bose-Einstein condensate with time-dependent atomic scattering lengths in an expulsive harmonic potential. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 3556–3569. [Google Scholar] [CrossRef]
  8. Huang, L.; Lenells, J. Nonlinear Fourier transforms for the sine-Gordon equation in the quarter plane. J. Differ. Equ. 2017, 264, 3445–3499. [Google Scholar] [CrossRef]
  9. Zakharov, V.E.; Shabat, A.B. Integration of nonlinear equations of mathematical physics by the method of inverse scattering. II. Funct. Anal. Appl. 1979, 13, 166–174. [Google Scholar] [CrossRef]
  10. Novikov, S.; Manakov, S.V.; Pitaevskii, L.P.; Zakharov, V.E. Theory of Solitons: The Inverse Scattering Method; Springer Science and Business Media: New York, NY, USA, 1984. [Google Scholar]
  11. Deift, P.; Zhou, X. A steepest descent method for oscillatory Riemann-Hilbert problems. Bull. Aust. Math. Soc. 1991, 26, 295–368. [Google Scholar] [CrossRef]
  12. Shchesnovich, V.S.; Yang, J.K. General soliton matrices in the Riemann-Hilbert problem for integrable nonlinear equations. J. Math. Phys. 2003, 44, 4604–4639. [Google Scholar] [CrossRef]
  13. Webb, J.R.L.; Infante, G. Positive solutions of nonlocal boundary value problems: A unified approach. J. Lond. Math. Soc. 2006, 74, 673–693. [Google Scholar] [CrossRef]
  14. Geng, X.G.; Wu, J.P. Riemann-Hilbert approach and N-soliton solutions for a generalized Sasa-Satsuma equation. Wave Motion 2016, 60, 62–72. [Google Scholar] [CrossRef]
  15. Ma, W.X. Riemann-Hilbert problems and N-soliton solutions for a coupled mKdV system. J. Geom. Phys. 2018, 132, 45–54. [Google Scholar] [CrossRef]
  16. Xu, X.X.; Sun, Y.P. An integrable coupling hierarchy of Dirac integrable hierarchy, its Liouville integrability and Darboux transformation. J. Nonlinear Sci. Appl. 2017, 10, 3328–3343. [Google Scholar] [CrossRef]
  17. Chen, J.C.; Ma, Z.Y.; Hu, Y.H. Nonlocal symmetry, Darboux transformation and soliton-cnoidal wave interaction solution for the shallow water wave equation. J. Math. Anal. Appl. 2018, 460, 987–1003. [Google Scholar] [CrossRef]
  18. Zhang, N.; Xia, T.C.; Jin, Q.Y. N-Fold Darboux transformation of the discrete Ragnisco-Tu system. Adv. Differ. Equ. 2018, 2018, 302. [Google Scholar] [CrossRef]
  19. Zhu, S.D.; Song, J.F. Residual symmetries, nth Bäcklund transformation and interaction solutions for (2+1)-dimensional generalized Broer-Kaup equations. Appl. Math. Lett. 2018, 83, 33–39. [Google Scholar] [CrossRef]
  20. Chen, H.H. General derivation of Bäcklund transformations from inverse scattering problems. Phys. Rev. Lett. 1974, 33, 925–928. [Google Scholar] [CrossRef]
  21. Wazwaz, A.M. Multiple-soliton solutions for the KP equation by Hirota’s bilinear method and by the tanh-coth method. Appl. Math. Comput. 2007, 190, 633–640. [Google Scholar] [CrossRef]
  22. Hu, X.B.; Wu, Y.T. Application of the Hirota bilinear formalism to a new integrable differential-difference equation. Phys. Lett. A 1998, 246, 523–529. [Google Scholar] [CrossRef]
  23. Ma, W.X.; Zhou, Y. Lump solutions to nonlinear partial differential equations via Hirota bilinear forms. J. Differ. Equ. 2018, 264, 2633–2659. [Google Scholar] [CrossRef]
  24. Fan, E.G.; Zhang, H.Q. The homogeneous balance method for solving nonlinear soliton equations. Acta Phys. Sin. 1998, 47, 353–362. [Google Scholar]
  25. Wang, M.L.; Zhou, Y.B.; Li, Z.B. Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics. Phys. Lett. A 1996, 216, 67–75. [Google Scholar] [CrossRef]
  26. Ma, W.X.; Wu, H.Y.; He, J.S. Partial differential equations possessing Frobenius integrable decompositions. Phys. Lett. A 2007, 364, 29–31. [Google Scholar] [CrossRef]
  27. Gao, L.; Ma, W.X.; Xu, W. Frobenius integrable decompositions for ninth-order partial differential equations of specific polynomial type. Appl. Math. Comput. 2010, 216, 2728–2733. [Google Scholar] [CrossRef]
  28. Fang, Y.; Dong, H.H.; Hou, Y.J.; Kong, Y. Frobenius integrable decompositions Of nonlinear evolution equations with modified term. Appl. Math. Comput. 2014, 226, 435–440. [Google Scholar] [CrossRef]
  29. Xu, T.; Tian, B. An extension of the Wronskian technique for the multicomponent Wronskian solution to the vector nonlinear Schrödinger equation. J. Math. Phys. 2010, 51, 359. [Google Scholar] [CrossRef]
  30. Corona-Corona, G. A Wronskian of Jost solutions. J. Math. Phys. 2004, 45, 4282–4287. [Google Scholar] [CrossRef]
  31. Yang, J.Y.; Ma, W.X.; Qin, Z.Y. Lump and lump-soliton solutions to the (2+1)-dimensional Ito equation. Anal. Math. Phys. 2018, 8, 427–436. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Dong, H.; Zhang, X.; Yang, H. Rational solutions and lump solutions to the generalized (3+1)-dimensional Shallow Water-like equation. Comput. Math. Appl. 2017, 73, 246–252. [Google Scholar] [CrossRef]
  33. Yong, X.L.; Ma, W.X.; Huang, Y.H.; Liu, Y. Lump solutions to the Kadomtsev-Petviashvili I equation with a self-consistent source. Comput. Math. Appl. 2018, 75, 2633–2659. [Google Scholar] [CrossRef]
  34. Wang, H. Lump and interaction solutions to the (2+1)-dimensional Burgers equation. Appl. Math. Lett. 2018, 85, 27–34. [Google Scholar] [CrossRef]
  35. Zeng, Z.J. Periodic solutions of a discrete time non-autonomous ratio-dependent predator-prey system with control. Commun. Korean Math. Soc. 2007, 22, 465–474. [Google Scholar] [CrossRef]
  36. Rostworowski, A. Higher order perturbations of Anti-de Sitter space and time-periodic solutions of vacuum Einstein equations. Phys. Rev. D 2017, 95, 16. [Google Scholar] [CrossRef]
  37. Savulescu, J. Wicked problems, complex solutions, and the cost of trust. J. Med. Ethics 2018, 44, 147–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Datsko, B.; Gafiychuk, V. Complex spatio-temporal solutions in fractional reaction-diffusion systems near a bifurcation point. Fract. Calc. Appl. Anal. 2018, 21, 237–253. [Google Scholar] [CrossRef] [Green Version]
  39. Zhu, Q.Z.; Xu, J.; Fan, E.G. The Riemann-Hilbert problem and long-time asymptotics for the Kundu-Eckhaus equation with decaying initial value. Appl. Math. Lett. 2017, 76, 81–89. [Google Scholar] [CrossRef]
  40. Tian, S.F. Initial-boundary value problems for the general coupled nonlinear Schrödinger equation on the interval via the Fokas method. J. Differ. Equ. 2017, 262, 506–558. [Google Scholar] [CrossRef]
  41. Zhang, N.; Xia, T.C.; Fan, E.G. A Riemann-Hilbert Approach to the Chen-Lee-Liu Equation on the Half Line. Acta. Math. Appl. Engl. Ser. 2018, 34, 493–515. [Google Scholar] [CrossRef]
  42. Song, G.H.; Shin, S.Y. Design of corrugated waveguide filters by the Gel’fand-Levitan-Marchenko inverse-scattering method. J. Opt. Soc. Am. A 1985, 2, 1905–1915. [Google Scholar] [CrossRef]
  43. Kivshar, Y.S.; Zhang, F.; Takeno, S. Multistable nonlinear surface modes. Phys. D 1998, 119, 125–133. [Google Scholar] [CrossRef]
  44. McLaughlin, D.W.; Scott, A.C. Perturbation analysis of fluxon dynamics. Phys. Rev. A 1978, 18, 1652. [Google Scholar] [CrossRef]
  45. Kivshar, Y.S.; Malomed, B.A. Dynamics of solitons in nearly integrable systems. Rev. Mod. Phys. 1989, 61, 763. [Google Scholar] [CrossRef]
  46. Shchesnovich, V.S. Perturbation theory for nearly integrable multicomponent nonlinear PDEs. J. Math. Phys. 2002, 43, 1460–1486. [Google Scholar] [CrossRef] [Green Version]
  47. Doktorov, E.V.; Wang, J.D.; Yang, J.K. Perturbation theory for bright spinor Bose-Einstein condensate solitons. Phys. Rev. A 2008, 77, 043617. [Google Scholar] [CrossRef] [Green Version]
  48. Kaup, D.J. Integrable systems and squared eigenfunctions. Theor. Math. Phys. 2009, 159, 806–818. [Google Scholar] [CrossRef]
  49. Wang, D.S.; Zhang, D.J.; Yang, J.K. Integrable properties of the general coupled nonlinear Schrödinger equations. J. Math. Phys. 2010, 51, 133–148. [Google Scholar] [CrossRef]
  50. Xiao, Y.; Fan, E.G. A Riemann-Hilbert approach to the Harry-Dym equation on the Line. Chin. Ann. Math. Ser. B 2014, 37, 373–384. [Google Scholar] [CrossRef]
Figure 1. Modulus of the soliton q ( x , t ) in Equation (85) with the parameters chosen as τ = 0.5 , ζ = 0.2 , x 0 = 0.2 , σ 0 = 0.2 .
Figure 1. Modulus of the soliton q ( x , t ) in Equation (85) with the parameters chosen as τ = 0.5 , ζ = 0.2 , x 0 = 0.2 , σ 0 = 0.2 .
Symmetry 11 00826 g001
Figure 2. Modulus of the soliton q ( x , t ) in Equation (86) with the parameters chosen as ζ 1 = 1 , τ 1 = 0.75 , ζ 2 = 0.65 , τ 2 = 0.55 , c 1 = 3 2 , c 2 = 1 2 i .
Figure 2. Modulus of the soliton q ( x , t ) in Equation (86) with the parameters chosen as ζ 1 = 1 , τ 1 = 0.75 , ζ 2 = 0.65 , τ 2 = 0.55 , c 1 = 3 2 , c 2 = 1 2 i .
Symmetry 11 00826 g002
Figure 3. Modulus of the soliton q ( x , t ) in Equation (86) with the parameters chosen as ξ 1 = 0.8 , η 1 = 0.55 , ξ 2 = 0.8 , η 2 = 0.45 , c 1 = 1 , c 2 = 1 .
Figure 3. Modulus of the soliton q ( x , t ) in Equation (86) with the parameters chosen as ξ 1 = 0.8 , η 1 = 0.55 , ξ 2 = 0.8 , η 2 = 0.45 , c 1 = 1 , c 2 = 1 .
Symmetry 11 00826 g003

Share and Cite

MDPI and ACS Style

Lin, Y.; Dong, H.; Fang, Y. N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem. Symmetry 2019, 11, 826. https://doi.org/10.3390/sym11060826

AMA Style

Lin Y, Dong H, Fang Y. N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem. Symmetry. 2019; 11(6):826. https://doi.org/10.3390/sym11060826

Chicago/Turabian Style

Lin, Yuxin, Huanhe Dong, and Yong Fang. 2019. "N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem" Symmetry 11, no. 6: 826. https://doi.org/10.3390/sym11060826

APA Style

Lin, Y., Dong, H., & Fang, Y. (2019). N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem. Symmetry, 11(6), 826. https://doi.org/10.3390/sym11060826

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