Next Article in Journal
Two Forms of the Integral Representations of the Mittag-Leffler Function
Previous Article in Journal
Existence of Nonnegative Solutions of Linear Autonomous Functional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Solutions of an Extended Nonlocal Nonlinear Schrödinger Equation in Plasmas

1
School of Mathematics and Physics, North China Electric Power University, Beijing 100083, China
2
Department of Mathematics and Statistics, University of South Florida, Tampa, FL 33620, USA
3
Department of Applied Mathematics, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(7), 1099; https://doi.org/10.3390/math8071099
Submission received: 9 June 2020 / Revised: 29 June 2020 / Accepted: 1 July 2020 / Published: 5 July 2020
(This article belongs to the Section Mathematical Physics)

Abstract

:
The parity-time symmetric nonlocal nonlinear Schrödinger equation with self-consistent sources (PTNNLSESCS) is used to describe the interaction between an high-frequency electrostatic wave and an ion-acoustic wave in plasmas. In this paper, the soliton solutions, rational soliton solutions and rogue wave solutions are derived for the PTNNLSESCS via the generalized Darboux transformation. We find that the soliton solutions can exhibit the elastic interactions of different type of solutions such as antidark-antidark, dark-antidark, and dark-dark soliton pairs on a continuous wave background. Also, we discuss the degenerate case in which only one antidark or dark soliton remains. The rogue wave solution is derived in some specially chosen situations.

1. Introduction

In 1998, the parity-time (PT) symmetry firstly appeared in quantum mechanics since Bender and Boettcher pointed out that non-Hermitian Hamiltonians exhibit entirely real spectra, provided that they respect both the parity and time-reversal symmetries (usually called the parity-time symmetry) [1]. Since then, physical systems exhibiting non-Hermitian systems with PT-symmetry have been the subject of intense investigation. In general, a necessary condition for a Hamiltonian H = p ^ 2 + V ( x ) to be PT-symmetric is that the complex potential satisfies V ( x ) = V * ( x ) , where p ^ denotes the momentum operator and V(x) is the complex potential [1,2]. Also, the notion of PT-symmetry has been applied to other areas of theoretical physics. In optical, R. El-Ganainy et al. developed a formalism suitable for describing coupled optical parity-time symmetric systems [3]. K. G. Makris et al. investigated the possibility of PT-symmetric periodic potentials within the context of optics [4]. A. Guo et al. demonstrated experimentally passive PT-symmetry breaking within the realm of optics [5]. Y. J. He et al. reported the existence and stability of lattice solitons in PT-symmetric mixed linear-nonlinear optical lattices in Kerr media [6].
In 2013, Ablowitz and Musslimani proposed and studied the following PT-symmetric nonlocal nonlinear Schrödinger (NLS) equation [7]:
i u z ( x , z ) = u x x ( x , z ) + 2 σ u ( x , z ) u * ( x , z ) u ( x , z ) ( σ = ± 1 ) ,
where the asterisk * denotes the complex conjugate, σ = ± 1 is the sign of nonlinearity (with the plus sign being the focusing case and minus sign being the defocusing case). The nonlinear term in Equation (1) brings a self-induced potential of the form V ( x , z ) = u ( x , z ) u * ( x , z ) satisfying V ( x , z ) = V * ( x , z ) . So the Hamiltonian of the Equation (1) satisfies the PT-symmetric condition. The Equation (1) is nonlocal as it has a term u ( x , z ) , which means that the time evolution of u does not only depends on the u ( x , z ) , but u ( x , z ) as well. Following the studying of this nonlocal PT-symmetric NLS equation, M. Li et al. applied the Nth iterated Darboux transformation to derive a chain of nonsingular localized-wave solutions of the PT-symmetric nonlocal NLS equation that can describe the soliton interactions on the continuous-wave (cw) background [8]. M. Li et al. derived the rational solutions of the PT-symmetric nonlocal NLS equation by generalized Darboux transformation [9]. In addition, some other PT-symmetric nonlocal integrable models have also been proposed. C. Q. Dai and W. H. Huang studied the coupled nonlinear Schrödinger equation in PT-symmetric coupled waveguides by means of the modified Darboux transformation method [10]. A. K. Sarma et al. investigated the continuous and discrete Schrödinger systems with PT-symmetric nonlinearities [11]. X. Y. Wen et al. studied novel higher-order rational solitons of the integrable nonlocal nonlinear Schrödinger equation with the self-induced PT-symmetric potential by the generalized perturbation N-1 fold Darboux transformation [12]. Z. X. Zhou studied global explicit solutions for nonlocal Davey-Stewartson I equation by Darboux transformation [13]. T. Xu et al. studied the nonsingular localized wave solutions of the partially PT-symmetric nonlocal Davey Stewartson I equation with zero background via the elementary Darboux transformation [14]. M. Li et al. analyzed the generation mechanism of rogue waves for the discrete nonlinear Schrödinger equation from the viewpoint of structural discontinuities [15]. K. Chen et al. showed a reduction technique that enables us to obtain solutions for the reduced local and nonlocal equation from the Ablowitz-Kaup-Newell-Suger hierarchy [16]. M. Duanmu et al. studied dynamical systems of linear and nonlinear PT-symmetric oligomers [17].
There are many methods in solving the integrable equations, including the Darboux transformation and Hirota bilinear method and so on. Yujia Zhang et al. investigated the interactions of vector anti-dark solitons for variable coefficients coupled NLS equation, found the double-S structure interactions [18]. H. Q. Zhang studied a modified NLS equation in inhomogeneous fibers and obtain the dark and antidark soliton solutions with Hirota bilinear method [19]. X. Z. Zhang propose a generalized long-water wave system study its invariant solutions and conservation laws are studied to guarantee its integrability [20]. R. A. El-Nabulsi used matrix Lie algebra to derive some non-standard higher-order equations [21]. Many equations posses peakons and smooth periodic waves. The most important model is the Camassa-Holm equation, which is studied by R. Camassa and D. D. Holm in 1993 [22]. They could be seen as the negative flows of some integrable hierarchy [23]. C. Z. Qu considered the stability of the peakon solitons of a modified Camassa-Holm equation with cubic nonlinearity [24]. S. Y. Lou extended the Camassa-Holm type equations to nonlocal case and constructed the AB peakon equations including the AB Camassa-Holm equation and the AB Degasperis-Procesi equation [25]. The nonlocal case of the integrable system attract lots of attentions in recent years.
The soliton equations with self-consistent sources (SESCS) play an important role in many fields of physics. The nonlinear Schrödinger equation with self-consistent sources (NLSESCS) describes the soliton propagation in a medium with both resonant and nonresonant nonlinearities and it also describes the nonlinear interaction of high-frequency electrostatic waves with ion acoustic waves in plasma [26].
In Reference [26], Claude studied the following system of coupled equation
i q t + 1 2 q x x q | q | 2 = i Γ a 1 a ¯ 2 , a 1 , x = q a 2 , a 2 , x 2 i k 0 a 2 = q ¯ a 1 ,
with the initial-boundary value problem:
q ( x , 0 ) = q 0 ( x ) L 1 ( R ) , a 1 ( x , t ) 1 , a s x + , a 2 ( x , t ) 0 , a s x .
This system occurs in plasma physics as the small amplitude limit of the general wave equations in a fluid-type warm electros/cold ions plasma.
To study the solution of NLSESCS, Y. J. Shao constructed a generalized Darbous transformation with an arbitrary function of t to study the soliton and positon solution of NLSESCS [27]. It has abundant applications in nonlinear envelope pulsed in fibers, pressure pulses in artery vessels, nonlinear Rossby waves in atmosphere, and matter waves in dilute-gas Bose-Einstein condensates [28]. The SESCS were first studied by Mel’nikov [29]. Since then, the SESCS have attracted some attention. A systematic way to construct the SESCS and their zero-curvature representations was proposed. Y. B. Zeng et al. studied the mKdV hierarchy with self-consistent sources by integral-type Darboux transformation [30]. Y. H. Huang et al. studied Camassa-Holm equation with self-consistent sources and its solutions [31]. Y. Q. Yao et al. studied the Qiao-Liu equation with self-consistent sources and its solutions [32].
In this paper, we construct the generalized Darboux transformation for the parity-time symmetric nonlocal nonlinear Schrödinger equation with self-consistent sources (PTNNLSESCS) and further reveal the soliton solution and the rational soliton phenomena on the cw background. The soliton solutions obtained in this work can display the profiles of the dark and antidark solitons. Also, we dicuss the degenerate cases in which only one dark or antidark soliton survives. We get the rogue wave solution of the PTNNLSESCS in some specially chose situations.
This paper is organized as follows: In Section 2, we review elementary Darboux transformation of the PT-symmetric NLS equation and derive the generalized Darboux transformation of the PTNNLSESCS. In Section 3, we derive the soliton solutions and discuss the degenerate cases for the PTNNLESCS. In Section 4, we derive the rational soliton solutions, rogue wave solutions and analysis the degenerate case for the PTNNLESCS. In Section 5, the conclude is given.

2. Elementary Darboux Transformation for the PTNNLSESCS

In this section we recall the successively iterated Darboux transformation of Equation (1) and construct the first-order generalized Darboux transformation with an arbitrary function of c(z) for the PTNNLSESCS. The Darboux transformation, which comprises the eigenfunction and potential transformations, can be used to recursively generate solutions including the soliton solutions, rational solutions from a trivial solution which can further provide an algebraic basis to analyze the asymptotic behavior of the solutions.
The Lax pair of Equation (1) can be written in the form [8]
Ψ x = U Ψ = λ u ( x , z ) σ u * ( x , z ) λ Ψ ,
Ψ z = V Ψ = 2 i λ 2 i σ u ( x , z ) u * ( x , z ) 2 i λ u ( x , z ) i u x ( x , z ) 2 i σ λ u * ( x , z ) i σ u x * ( x , z ) 2 i λ 2 + i σ u ( x , z ) u * ( x , z ) Ψ ,
where Ψ = ( f , g ) T (the superscript T represents the vector transpose) is the vector eigenfunction, λ is the spectral parameter, and Equation (1) satisfies the compatibility condition U z V x + U V V U = 0 .
In Reference [8], the Nth iterated elementary Darboux transformation for Equation (1) can be constituted by the eigenfunction
Ψ [ N ] = T [ N ] Ψ , T [ N ] = λ N n = 1 N a n ( x , z ) λ n 1 n = 1 N b n ( x , z ) ( λ ) n 1 n = 1 N c n ( x , z ) λ n 1 λ N n = 1 N d n ( x , z ) ( λ ) n 1 ,
and the Darboux transformation
u [ N ] ( x , z ) = u ( x , z ) + 2 ( 1 ) N 1 b N ( x , z ) ,
where N represents the iterated time. The new eigenfunction Ψ [ N ] is required to satisfy the Lax pair in Equations (2a) and (2b) with u [ N ] ( x , z ) and u [ N ] * ( x , z ) instead of u ( x , z ) and u * ( x , z ) , respectively. The functions a n ( x , z ) , b n ( x , z ) , c n ( x , z ) and d n ( x , z ) ( 1 n N ) can be determined from
T [ N ] | λ = λ k Ψ k = 0 , T [ N ] | λ = λ k * Ψ ¯ k = 0 ( 1 n N ) ,
where Ψ k = [ f k ( x , z ) , g k ( x , z ) T ] and Ψ ¯ k = [ g k * ( x , z ) , σ f k * ( x , z ) T ] are the solutions of Equations (2a) and (2b) with λ = λ k and λ = λ k * , respectively. In particular, the functions b N ( x , z ) and c N ( x , z ) can be obtained in the determinant form
b N ( x , z ) = ( 1 ) N 1 τ N + 1 , N 1 τ N , N ,
with
τ N , N = F N × N G N × N σ G ¯ N × N F ¯ N × N ,
where the block matrices
F N × N = [ λ k m 1 f k ( x , z ) ] 1 k , m N , G N × N = [ ( λ k ) m 1 g k ( x , z ) ] 1 k , m N , F ¯ N × N = [ ( λ k * ) m 1 f k * ( x , z ) ] 1 k , m N , G ¯ N × N = [ ( λ k * ) m 1 g k * ( x , z ) ] 1 k , m N .
It is well known that the AKNS equation with self-consistent sources (AKNSESCS) is defined as [27]
q t = i ( q x x 2 q 2 r ) + j = 1 n ( φ j ( 1 ) ) 2 , r t = i ( r x x 2 q r 2 ) + j = 1 n ( φ j ( 2 ) ) 2 ,
φ j , x = φ 1 φ 2 x = λ j q r λ j φ j , j = 1 , , n ,
where λ j are n distinct complex constants, φ j = ( φ j ( 1 ) , φ j ( 2 ) ) T (hereafter, we use superscripts (1) and (2) to denote the first and second elements of a two-dimensional vector respectively). In this paper, we only consider the PT-symmetric nonlocal NLS equation with one source.
With the reduction q = u ( x , z ) , r = u * ( x , z ) and n = 1 , the PTNNLSESCS can be constructed from AKNSESCS
i u z ( x , z ) u x x ( x , z ) + 2 u ( x , z ) u * ( x , z ) u ( x , z ) = 2 i ( ( φ 2 * ( x , z ) ) 2 ( φ 1 ( x , z ) ) 2 ) ,
φ x ( x , z ) = φ 1 ( x , z ) φ 2 ( x , z ) x = λ 1 u ( x , z ) u * ( x , z ) λ 1 φ ( x , z ) .
The Lax pair of Equation (9) can be written in the form
Ψ x = λ u ( x , z ) u * ( x , z ) λ Ψ , Ψ z = 2 i λ 2 + i u ( x , z ) u * ( x , z ) 2 i λ u ( x , z ) i u x ( x , z ) 2 i λ u * ( x , z ) + i u x * ( x , z ) 2 i λ 2 i u ( x , z ) u * ( x , z ) Ψ + 1 λ λ 1 φ 1 ( x , z ) φ 2 ( x , z ) ( φ 1 ( x , z ) ) 2 ( φ 2 ( x , z ) ) 2 φ 1 ( x , z ) φ 2 ( x , z ) Ψ
1 λ λ 1 * φ 1 * ( x , z ) φ 2 * ( x , z ) ( φ 2 * ( x , z ) ) 2 ( φ 1 * ( x , z ) ) 2 φ 1 * ( x , z ) φ 2 * ( x , z ) Ψ .
Based on the Darboux transformation for the AKNS equation, we derive the generalized Darboux transformations of Equation (9).
Theorem 1.
Let f 1 , φ = 0 , and g 1 , φ = 0 be two solutions of Equation (10) with λ = λ 1 and F = c ( z ) f 1 + g 1 , then the Darboux transformation for Equation (9) is defined as
Ψ [ 1 ] = T [ 1 ] Ψ ,
u [ 1 ] = u + 2 ( λ 1 λ 1 * ) ( P [ 1 ] ) 12 ,
φ [ 1 ] = c ( z ) ( λ 1 λ 1 * ) d e t ( f 1 , g 1 ) T [ 1 ] ( λ 1 ) f 1 ,
where
T [ 1 ] = λ λ 1 + ( λ 1 λ 1 * ) P [ 1 ] , P [ 1 ] = F F ¯ F ¯ F , F = F 1 F 2 , F ¯ = F 1 * ( x , z ) F 2 * ( x , z ) ,
and ( P [ 1 ] ) 12 represents the entry of matrix P [ 1 ] of the first row and second column.
Proof. 
Let u [ 1 ] z = I 1 + I 2 c ( z ) , according to Equation (9) we know
i I 1 u [ 1 ] x x + 2 u [ 1 ] 2 u [ 1 ] * ( x , z ) = 0 ,
and
i I 2 c ( z ) = 2 i c ( z ) ( λ 1 λ 1 * ) F 1 F 2 * ( x , z ) F ¯ F = 2 i c ( z ) ( λ 1 λ 1 * ) ( F ¯ F ) 2 f 1 F 2 * ( x , z ) F ¯ F + F 1 f 2 * ( x , z ) F ¯ F F 1 F 2 * ( x , z ) f ¯ F F 1 F 2 * ( x , z ) F ¯ f .
After some calculations we have
i I 2 c ( z ) = 2 i c ( z ) ( λ 1 λ 1 * ) ( F ¯ F ) 2 ( F 2 * ( x , z ) ) 2 d e t ( f 1 , g 1 ) F 1 2 d e t ( f 1 * ( x , z ) , g 1 * ( x , z ) ) ,
from Equation (11c), we obtain
φ [ 1 ] = c ( z ) ( λ 1 λ 1 * ) d e t ( f 1 , g 1 ) λ 1 λ 1 * F ¯ F | F 2 | 2 F 1 F 2 * ( x , z ) F 1 * ( x , z ) F 2 | F 1 | 2 f 1 = c ( z ) ( λ 1 λ 1 * ) d e t ( f 1 , g 1 ) λ 1 λ 1 * F ¯ F d e t ( f 1 , g 1 ) F 2 * ( x , z ) F 1 * ( x , z ) ,
we find that
i I 2 c ( z ) = 2 i ( ( φ 2 [ 1 ] * ( x , z ) ) 2 φ 1 [ 1 ] 2 ) ,
so
i u [ 1 ] z u [ 1 ] x x + 2 u [ 1 ] u [ 1 ] * ( x , z ) u [ 1 ] = 2 i ( ( φ 2 [ 1 ] * ( x , z ) ) 2 φ 1 [ 1 ] 2 ) .
 □
From Theorem 1, we can obtain different types solutions of Equation (9). Next, we will use the generalized Darboux transformation to construct the soliton solutions and rational solutions on a cw background.

3. Soliton Solutions of the PTNNLSESCS

In this section we construct the soliton solutions on the cw background. It is not difficult to find that Equation (9) admits the plane-wave solution
u = ρ e 2 i ρ 2 z + i ϕ , φ = 0 ,
where ρ and ϕ are two real parameters. In the following, we take λ = λ 1 , using λ 1 R and λ 1 I to denote the real and imaginary parts of λ = λ 1 , respectively. We take λ 1 R = 0 and 0 < | λ 1 I | < ρ to derive the soliton solutions of the PTNNLSESCS. Inserting Equation (17) into the Lax pair (10), we obtain
f 1 = e 2 i ρ 2 z + i ϕ 2 α 1 e 2 i ρ 2 z + i ϕ 2 ( s 1 λ 1 ) α 1 ρ e s 1 χ 1 ,
g 1 = e 2 i ρ 2 z + i ϕ 2 β 1 e 2 i ρ 2 z + i ϕ 2 ( s 1 + λ 1 ) β 1 ρ e s 1 χ 1 ,
where s 1 = λ 1 2 + ρ 2 , χ 1 = x 2 i λ 1 z , α 1 and β 1 are arbitrary nonzero complex parameters. Then, with substitution of (18) into the Darboux transformations (11), the solution can be written as
u [ 1 ] = ρ e 2 i ρ 2 z + i ϕ 1 2 λ 1 I ( c ( z ) e 2 s 1 χ 1 + γ 1 ) ( κ 1 c ( z ) e 2 s 1 ω 1 + γ 1 * κ 1 * ) c 2 ( z ) ρ 2 e 2 s 1 ( χ 1 ω 1 ) + c ( z ) λ 1 I ( κ 1 γ 1 e 2 s 1 ω 1 + κ 1 * γ 1 * e 2 s 1 χ 1 ) + | γ 1 | 2 ρ 2 ,
φ 1 [ 1 ] = i ρ β 1 λ 1 I s 1 c ( z ) α 1 e 2 i ρ 2 z + i ϕ 2 + s 1 χ 1 β 1 s 1 ( κ 1 c ( z ) e 2 s 1 ω 1 + γ 1 * κ 1 * ) c 2 ( z ) ρ 2 e 2 s 1 ( χ 1 ω 1 ) + c ( z ) λ 1 I ( κ 1 γ 1 e 2 s 1 ω 1 + κ 1 * γ 1 * e 2 s 1 χ 1 ) + | γ 1 | 2 ρ 2 ,
φ 2 [ 1 ] = i ρ i ρ β 1 λ 1 I s 1 c ( z ) α 1 e 2 i ρ 2 z + i ϕ 2 + s 1 χ 1 β 1 s 1 ρ ( c ( z ) e 2 s 1 ω 1 + γ 1 * ) c 2 ( z ) ρ 2 e 2 s 1 ( χ 1 ω 1 ) + c ( z ) λ 1 I ( κ 1 γ 1 e 2 s 1 ω 1 + κ 1 * γ 1 * e 2 s 1 χ 1 ) + | γ 1 | 2 ρ 2 ,
where χ 1 = x + 2 λ 1 I z , ω 1 = x 2 λ 1 I z , s 1 = ρ 2 λ 1 I 2 , κ 1 = λ 1 I i s 1 and γ 1 = β 1 α 1 . In this section, we take c(z) as a real arbitrary function of z. It can be proved that the solution (19) has no singularity if and only if the following condition is satisfied:
I m ( γ 1 κ 1 ) 0 or c ( z ) λ 1 I R e ( κ 1 γ 1 ) > 0 .
We make an asymptotic analysis to better understand the solitonic behavior in the solution (19) under the condition (20) as follows.
(i) Along the line x + 2 λ 1 I z = 0 as z , we have
u [ 1 ] u 1 ± = ρ e 2 i ρ 2 z + i ϕ 1 2 λ 1 I ( c ( z ) e 2 s 1 χ 1 + γ 1 ) μ 1 ± c ( z ) e 2 s 1 χ 1 + ν 1 ± ,
| u [ 1 ] | 2 | u 1 ± | 2 = ρ 2 1 2 γ 1 I s 1 R e ( γ 1 κ 1 ) + s g n ( λ 1 I c ( z ) ) | γ 1 | ρ cosh ( 2 s 1 χ 1 + Δ 1 ± ) ,
with μ 1 = λ 1 I , ν 1 = γ 1 κ 1 , μ 1 + = κ 1 * , ν 1 + = γ 1 λ 1 I , Δ 1 = ln | c ( z ) λ 1 I | ρ | γ 1 | and Δ 1 + = ln ρ | c ( z ) | | γ 1 | | λ 1 I | , where the plus sign corresponds to λ 1 I > 0 as z or λ 1 I < 0 as z and the minus sign corresponds to λ 1 I < 0 as z or λ 1 I > 0 as z .
(ii) Along the line x 2 λ 1 I z = 0 as z , we have
u [ 1 ] u 2 ± = ρ e 2 i ρ 2 z + i ϕ 1 2 λ 1 I ( c ( z ) κ 1 e 2 s 1 ω 1 + γ 1 * κ 1 * ) c ( z ) ξ 1 ± e 2 s 1 ω 1 + η 1 ± ,
| u [ 1 ] | 2 | u 2 ± | 2 = ρ 2 + 2 s 1 I m ( λ 1 I κ 1 2 ) R e ( γ 1 κ 1 ) + s g n ( λ 1 I c ( z ) ) | γ 1 | ρ cosh ( 2 s 1 ω 1 Δ 1 ± ) ,
with ξ 1 = λ 1 I κ 1 , η 1 = ρ 2 γ 1 * , ξ 1 + = ρ 2 and η 1 + = λ 1 I κ 1 * γ 1 * , where the plus sign corresponds to λ 1 I < 0 as z or λ 1 I > 0 as z and the minus sign corresponds to λ 1 I > 0 as z or λ 1 I < 0 as z .
The above asymptotic analysis shows three different types of elastic interactions. The associated parametric conditions are given in Table 1. The asymptotic soliton u 1 ± represents the antidark soliton for c ( z ) λ 1 I γ 1 I < 0 or dark soliton for c ( z ) λ 1 I γ 1 I > 0 on a cw background and both the antidark soliton and dark solitons are localized along the line x + 2 z ρ = 0 , while the asymptotic soliton u 2 ± also represents the antidark soliton for c ( z ) I m ( γ 1 I κ 1 2 ) > 0 or dark soliton for c ( z ) I m ( γ 1 I κ 1 2 ) < 0 on the same cw background. In this case, both the antidark soliton and dark soliton solitons are localized along the line x 2 z ρ = 0 . In particular, for the degenerate case γ 1 I = 0 , the asymptotic soliton u 1 ± disappear as z ± . Similarly, for the degenerate case λ 1 I I m ( ( γ 1 κ 1 2 ) = 0 , the only surviving asymptotic soliton is u 1 ± .
The heights of the antidark solitons or the depth of the dark solitons from the cw background for | u 1 ± | 2 and | u 2 ± | 2 , that is
A 1 ± = 2 ρ 2 | s 1 γ 1 I | R e ( γ 1 κ 1 ) + s g n ( λ 1 I c ( z ) ) ρ | γ 1 I | ,
A 2 ± = 2 | s 1 I m ( γ 1 I κ 1 2 ) | R e ( γ 1 κ 1 ) + s g n ( λ 1 I c ( z ) ) ρ | γ 1 I | ,
respectively. The envelope velocity of u 1 ± and u 2 ± , that is v 1 = v 1 + = 2 λ 1 I and v 2 = v 2 + = 2 λ 1 I , respectively. There exists a phase shift between u i and u i + ( i = 1 , 2 ) , and its absolute value can be calculated by | Δ 1 + Δ 1 | = 2 | ln ρ | λ 1 I | | which means that the phase shift depends on both ρ and λ 1 I .
In any of these interactions, the interacting solitons can completely recover their individual shapes and velocities upon an interaction and have only the phase shifts for their envelopes. The Figure 1, Figure 2 and Figure 3 suggest that the solution (19a) under the condition (20) can describe the elastic soliton interactions on the cw background. The elastic interactions may occur between two antidark solitons, two dark solitons, or dark and antidark solitons. The existence of the antidark soliton or the dark soliton can be decided by the sign of c(z) [see Figure 1a–c]. Particularly taking γ 1 I = 0 , the asymptotic solitons ( u 1 , u 1 + ) disappears as z ± , but ( u 2 , u 2 + ) still exists and represent antidark solitons for c ( z ) γ 1 R < 0 or dark solitons for c ( z ) γ 1 R > 0 [see Figure 4]. However, the solution (19a) in this degenerate case cannot be regarded as the conventional single soliton because there exists a phase shift between u 2 and u 2 + . Similarly, for the degenerate case I m ( γ 1 κ 1 2 ) = 0 , one can find that only the pair of asymptotic solitons ( u 1 , u 1 + ) exists as | z | and there is also a phase shift between u 1 and u 1 + [see Figure 5]. Associated with c ( z ) λ 1 I γ 1 I < 0 and c ( z ) λ 1 I γ 1 I > 0 , ( u 1 , u 1 + ) can represent antidark solitons and dark solitons, respectively.

4. Rational Solutions of the PTNNLSESCS

In this section we construct the rational solutions of the PTNNLSESCS on the cw background. Equation (9) admits plane-wave solution u = ρ e 2 i ρ 2 z + i ϕ , φ = 0 . We take c ( z ) = c 1 ( z ) + i c 2 ( z ) , λ 1 = i ρ to derive the rational solitons and rogue wave solutions. Inserting them into the Lax pair (10), by eigenvalue method we obtain
f 1 = e 2 i ρ 2 z + i ϕ 2 i e 2 i ρ 2 z + i ϕ 2 ,
g 1 = 2 i ρ ( x + 2 z ρ + s 1 ) e 2 i ρ 2 z + i ϕ 2 2 ( i + ρ ( x + 2 z ρ + s 1 ) ) e 2 i ρ 2 z + i ϕ 2 ,
where s 1 is an arbitrary complex parameter.
Then, with substitution of (25) into the Darboux transformations (11), we obtain the first-order rational soliton solution as follows:
u [ 1 ] = ρ e 2 i ρ 2 z + i ϕ 1 ( 2 ξ + 2 K 1 2 i ) ( 2 η + 2 K 1 * 2 i ) 2 ξ η + 2 K 1 * ξ + 2 K 1 η + | K 1 | 2 + 1 ,
φ 1 [ 1 ] = i ρ c 2 ( z ) 2 ( 2 η + 2 K 1 * 2 i ) e 2 i ρ 2 z + i ϕ 2 2 ξ η + 2 K 1 * ξ + 2 K 1 η + | K 1 | 2 + 1 ,
φ 2 [ 1 ] = i i ρ c 2 ( z ) 2 ( 2 η + 2 K 1 * + 2 i ) e 2 i ρ 2 z + i ϕ 2 2 ξ η + 2 K 1 * ξ + 2 K 1 η + | K 1 | 2 + 1 ,
where ξ = i c 1 ( z ) + c 2 ( z ) + 2 ρ ( x + 2 z ρ ) , η = i c 1 ( z ) + c 2 ( z ) + 2 ρ ( x + 2 z ρ ) , s 1 = a + i b and K 1 = 2 ρ s 1 + i . The solution in Equation (26) has no singularity if and only if the following condition is satisfied:
c 1 ( z ) 2 ρ b + 1 2 .
Under this condition, we perform an asymptotic analysis of the solution in Equation (26) so as to clarify the dynamical behavior underlying the solution.
(i) when c 1 ( z ) = 0 . Firstly, we obtain the asymptotic expression of the solution in Equation (26) along the line c 2 ( z ) + 2 ρ ( x + 2 z ρ ) 0 as | z | as follows:
u [ 1 ] u I = ρ e 2 i ρ 2 z + i ϕ 1 2 ( 2 ξ + K 1 i ) 2 ξ + K 1 ,
| u [ 1 ] | 2 | u I | 2 = ρ 2 1 8 ρ I m ( s 1 ) | 2 ξ + K 1 | 2 ,
φ 1 [ 1 ] φ 1 I = 2 i ρ c 2 ( z ) e 2 i ρ 2 z + i ϕ 2 2 ξ + K 1 ,
| φ 1 [ 1 ] | 2 | φ 1 I | 2 = 2 ρ | c 2 ( z ) | | 2 ξ + K 1 | 2 ,
φ 2 [ 1 ] φ 2 I I = i 2 i ρ c 2 ( z ) e 2 i ρ 2 z + i ϕ 2 2 ξ + K 1 ,
| φ 2 [ 1 ] | 2 | φ 2 I I | 2 = 2 ρ | c 2 ( z ) | | 2 ξ + K 1 | 2 .
Secondly, we derive the asymptotic expression of the solution in Equation (26) along the line c 2 ( z ) + 2 ρ ( x + 2 z ρ ) 0 as | z | as follows:
u [ 1 ] u I I = ρ e 2 i ρ 2 z + i ϕ 1 2 ( 2 η + K 1 * i ) 2 η + K 1 * ,
| u [ 1 ] | 2 | u I I | 2 = ρ 2 1 + 8 ( 1 + ρ I m ( s 1 ) ) | 2 η + K 1 * | 2 ,
φ 1 [ 1 ] = i ρ c 2 ( z ) 2 ( 2 η + 2 K 1 * 2 i ) e 2 i ρ 2 z + i ϕ 2 2 ξ η + 2 K 1 * ξ + 2 K 1 η + | K 1 | 2 + 1 0 ,
φ 2 [ 1 ] = i i ρ c 2 ( z ) 2 ( 2 η + 2 K 1 * + 2 i ) e 2 i ρ 2 z + i ϕ 2 2 ξ η + 2 K 1 * ξ + 2 K 1 η + | K 1 | 2 + 1 0 .
The above asymptotic analysis shows three different types of elastic interactions. The associated parametric conditions are given in Table 2. For the cases I m ( s 1 ) > 0 and I m ( s 1 ) < 0 , the intensity | u I | 2 can respectively exhibit the rational dark (RD) soliton beneath the cw background u = ρ e 2 i ρ 2 z + i ϕ , φ = 0 and the rational antidark (RAD) soliton on top of the same background, and the valley and peak are both localized along the line 2 c 2 ( z ) + 2 ρ ( x + 2 z ρ ) + R e ( K 1 ) = 0 . The intensity | u 1 I I | 2 display the RD and RAD soliton profiles which are associated with 1 + ρ I m ( s 1 ) < 0 and 1 + ρ I m ( s 1 ) > 0 respectively. In this case, both the RD and RAD solitons are localized along the line 2 c 2 ( z ) + 2 ρ ( x + 2 z ρ ) + R e ( K 1 ) = 0 . In particular, for the degenerate case I m ( s 1 ) = 0 , the asymptotic soliton u I disappear as z ± . Similarly, for the degenerate case 1 + ρ I m ( s 1 ) = 0 , the only surviving asymptotic soliton is u I .
The heights of the RAD soliton or the depth of the RD soliton from the cw background for | u I | and | u I I | , that is A 1 = 8 ρ 3 I m ( s 1 ) 2 ρ I m ( s 1 ) + 1 and A 2 = 8 ρ 2 ( 1 + ρ I m ( s 1 ) ) 2 ρ I m ( s 1 ) + 1 , respectively. The velocities of the RD and RAD solitons from the cw background for | u I | and | u I I | , that is v 1 = 2 ρ c 2 ( z ) 2 ρ and v 2 = 2 ρ + c 2 ( z ) 2 ρ , respectively.
The asymptotic analysis implies that Equation (26) can describe the elastic interactions of rational solitons as two interacting solitons retain their individual shapes, intensities and velocities as z . However, different from the standard elastic interaction in the PTNNLSESCS, each soliton experiences no phase shift upon the interaction.
In general, Equation (26) exhibits three different types of elastic interactions between two rational solitons on a cw background, as shown in Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. More specifically, with I m ( s 1 ) < 0 and 1 + ρ I m ( s 1 ) > 0 , the asymptotic soliton u I and u I I both displays RAD solitons profile [see Figure 6a]. With I m ( s 1 ) < 0 and 1 + ρ I m ( s 1 ) < 0 , the asymptotic soliton u I and u I I displays RAD and RD solitons profile [see Figure 7a]. With I m ( s 1 ) > 0 and 1 + ρ I m ( s 1 ) > 0 , the asymptotic soliton u I and u I I displays RD and RAD solitons profile [see Figure 8a]. The function c(z) of the solutions (26) can change the shape of the RAD or RD [see Figure 9, Figure 10 and Figure 11]. The solution of φ 1 [ 1 ] and φ 2 [ 1 ] always diplay RAD soliton solution along the line c 2 ( z ) + 2 ρ ( x + 2 z ρ ) = 0 as | z | and has a constant height as c 2 ( z ) = z [see Figure 6b,c], but the height change as c 2 ( z ) = z 2 [see Figure 9b,c]. The solution of φ 1 [ 1 ] and φ 2 [ 1 ] disappear along the line c 2 ( z ) + 2 ρ ( x + 2 z ρ ) = 0 as | z | .
In particular, with I m ( s 1 ) = 0 the asymptotic soliton u I disappear as z ± while u I I displays a RAD soliton profile [see Figure 12]. Similarly, for the degenerate case 1 + ρ I m ( s 1 ) = 0 , the only surviving asymptotic soliton is u I and it takes the shape of the RAD type [see Figure 13]. In either of the two degenerate cases, one asymptotic soliton disappears in the far-field region, but it still affects the other one in the near-field region, that is, the surviving soliton is segmented into two pieces at some finite value of z. Therefore, such two degenerate cases of the solution in Equation (26a) cannot be simply regarded as the conventional single soliton.
(ii) when c 1 ( z ) 0 , Equation (26) can be written as follows:
u [ 1 ] = ρ e 2 i ρ 2 z + i ϕ ( 1 + G 1 H 1 + i H ) ,
φ 1 [ 1 ] = G 2 H 1 + i H ,
φ 2 [ 1 ] = G 3 H 1 + i H ,
where
G 1 = 2 ( 2 2 i ( c 2 ( z ) + 2 2 ρ 2 z ) + 4 i ρ a + 2 ) , G 2 = 2 i ρ c 2 ( z ) ( 2 i c 1 ( z ) + 2 c 2 ( z ) + 2 ρ ( x + 2 z ρ ) + 2 ( a i b ) 2 i ) e 2 i ρ 2 z + i ϕ 2 , G 3 = i 2 i ρ c 2 ( z ) ( 2 i c 1 ( z ) + 2 c 2 ( z ) + 2 ρ ( x + 2 z ρ ) + 2 ( a i b ) ) e 2 i ρ 2 z + i ϕ 2 , H 1 = 2 ( 2 ρ a + c 2 ( z ) + 2 2 ρ 2 z ) 2 + ( 2 ρ b c 1 ( z ) ) 2 + ( 2 ρ b + 2 c 1 ) 2 4 ρ 2 x 2 , H 2 = 4 ρ x ( 2 ρ b 2 c 1 ( z ) + 1 ) .
We can see that the solution of (34a) tends to background in the z direction as | z | because the denominator containing z 2 and ( c 2 ( z ) ) 2 has higher exponential than the molecule containing z and c 2 ( z ) . The solution of (34a) tends to background in the x direction as | x | because the denominator contains x 2 and x, but the molecule does not contain the x. So the solution of (34a) would tend to background in all direction, which means that we get the rogue wave solution as shown in Figure 14a. The solution of (34b) tends to background in the z direction as | z | because the denominator containing z 2 , ( c 1 ( z ) ) 2 and ( c 2 ( z ) ) 2 has higher exponential than the molecule containing z, c 1 ( z ) and c 2 ( z ) . The solution of (34b) tends to background in the x direction as | x | because the denominator contains x 2 and x, but the molecule only contains the x. So the solution of (34b) would tend to background in all direction, which means that we get the rogue wave solution as shown in Figure 14b. The solution of (34c) has the same reason as the solution of (34b) as shown in Figure 14c.

5. Conclusions

In this paper, we have studied the nonlinear localized wave on the cw background for the PTNNLSESCS. The generalized Darboux transformation for the PTNNLSESCS is derived. We have derived the nonsingular first-order soliton solution, the first-order rational soliton solution and rogue wave solution from a cw solution by the generalized Darboux transformation. We have revealed the elastic antidark-antidark soliton, antidark-dark soliton, and dark-dark soliton interactions, in which there is a phase shift for interacting solitons. we have also revealed the elastic RAD-RAD and RAD-RD soliton interactions, in which there is no phase shift for interacting solitons. Meanwhile, we have discussed the degenerate case in which only one soliton remains.

Author Contributions

Formal analysis, M.L. and Z.Y.; supervision, Y.Y.; writing—original draft, H.J.; writing—review & editing, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Fundamental Research Funds of the Central Universities with the [grant number 2020MS043]. Min Li thanks the financial support by the National Natural Science Foundation of China [grant number 261505054]. Yuqin Yao thanks the financial support by Beijing Natural Science Foundation [grant number 1182009].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bender, C.M.; Boettcher, S. Real Spectra in Non-Hermitian Hamiltonians Having PT-Symmetry. Phys. Rev. Lett. 1998, 80, 5243. [Google Scholar] [CrossRef] [Green Version]
  2. Bender, C.M.; Brody, D.C.; Jones, H.F.; Meister, B.K. Faster than Hermitian Quantum Mechanics. Phys. Rev. Lett. 2007, 98, 040403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. El-Ganainy, R.; Makris, K.G.; Christodoulides, D.N.; Musslimani, Z.H. Theory of coupled optical PT-symmetric structures. Opt. Lett. 2007, 32, 2632–2634. [Google Scholar] [CrossRef]
  4. Makris, K.G.; El-Ganainy, R.; Christodoulides, D.N.; Musslimani, Z.H. Beam Dynamics in PT Symmetric Optical Lattices. Phys. Rev. Lett. 2008, 100, 103904. [Google Scholar] [CrossRef]
  5. Guo, A.; Salamo, G.J.; Duchesne, D.; Morandotti, R.; Volatier-Ravat, M.; Aimez, V.; Siviloglou, G.A.; Christodoulides, D.N. Observation of PT-Symmetry Breaking in Complex Optical Potentials. Phys. Rev. Lett. 2009, 103, 093902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. He, Y.J.; Zhu, X.; Mihalache, D.; Liu, J.L.; Chen, Z.X. Lattice solitons in PT-symmetric mixed linear-nonlinear optical lattices. Phys. Rev. A 2012, 85, 013831. [Google Scholar] [CrossRef]
  7. Ablowitz, M.J.; Musslimani, Z.H. Integrable Nonlocal Nonlinear Schrödinger Equation. Phys. Rev. Lett. 2013, 110, 064105. [Google Scholar] [CrossRef] [Green Version]
  8. Li, M.; Xu, T. Dark and antidark soliton interactions in the nonlocal nonlinear Schrödinger equation with the self-induced parity-time-symmetric potential. Phys. Rev. E 2015, 91, 033202. [Google Scholar] [CrossRef]
  9. Li, M.; Xu, T.; Meng, D.X. Rational Solitons in the Parity-Time-Symmetric Nonlocal Nonlinear Schrödinger Model. J. Phys. Soc. Jpn. 2016, 85, 124001. [Google Scholar] [CrossRef] [Green Version]
  10. Dai, C.Q.; Huang, W.H. Multi-rogue wave and multi-breather solutions in PT-symmetric coupled waveguides. Appl. Math. Lett. 2014, 32, 35–40. [Google Scholar] [CrossRef]
  11. Sarma, A.K.; Miri, M.A.; Musslimani, Z.H.; Christodoulides, D.N. Continuous and discrete Schrödinger systems with parity-time-symmetric nonlinearities. Phys. Rev. E 2014, 89, 052918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wen, X.Y.; Yan, Z.Y.; Yang, Y.Q. Dynamics of higher-order rational solitons for the nonlocal nonlinear Schrödinger equation with the self-induced parity-time-symmetric potential. Chaos 2016, 6, 063123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhou, Z.X. Darboux transformations and global explicit solutions for nonlocal Davey-Stewartson I equation. arXiv 2018, arXiv:1612.05689. [Google Scholar] [CrossRef]
  14. Xu, T.; Li, M.; Huang, Y.H.; Chen, Y.; Yu, C. Nonsingular localized wave solutions for the nonlocal Davey-Stewartson I equation with zero background. Mod. Phys. Lett. B 2017, 35, 1750338. [Google Scholar] [CrossRef]
  15. Li, M.; Shui, J.J.; Xu, T. Generation mechanism of rogue waves for the discrete nonlinear Schrödinger equation. Appl. Math. Lett. 2018, 83, 110–115. [Google Scholar] [CrossRef]
  16. Chen, K.; Deng, X.; Lou, S.Y.; Zhang, D.J. Solutions of local and nonlocal equations reduced from the AKNS hierarchy. Stud. Appl. Math. 2018, 141, 113–141. [Google Scholar] [CrossRef]
  17. Duanmu, M.; Li, K.; Horne, R.L.; Kevrekidis, P.G.; Whitaker, N. Linear and nonlinear parity-time-symmetric oligomers: A dynamical systems analysis. Philos. Trans. R. Soc. A 2013, 1989, 20120171. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, Y.J.; Yang, C.Y.; Yu, W.T.; Mirzazadeh, M.; Zhou, Q.; Liu, W.J. Interactions of vector anti-dark solitons for the coupled nonlinear Schrödinger equation in inhomogeneous fibers. Nonlinear Dyn. 2018, 94, 1351–1360. [Google Scholar] [CrossRef]
  19. Zhang, H.Q.; Zhai, B.G.; Wang, X.L. Dark and antidark soliton solutions in the modified nonlinear Schrödinger equation with distributed coefficients in inhomogeneous fibers. Phys. Scr. 2012, 85, 015006. [Google Scholar] [CrossRef]
  20. Zhang, X.Z.; Zhang, Y.F. Some invariant solutions and conservation laws of a type of long-water wave system. Adv. Differ. Equ. 2019, 2019, 496. [Google Scholar] [CrossRef] [Green Version]
  21. El-Nabulsi, R.A. Non-standard higher-order G-strand partial differential equations on matrix Lie algebra. J. Niger. Math. Soc. 2017, 36, 101–112. [Google Scholar]
  22. Camassa, R.; Holm, D.D. An integrable shallow water equation with peaked solitons. Phys. Rev. Lett. 1993, 71, 1661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Novikov, V. Generalizations of the Camassa-Holm equation. J. Phys. A Math. Theor. 2009, 42, 342002. [Google Scholar] [CrossRef]
  24. Qu, C.Z.; Liu, X.C.; Liu, Y. Stability of peakons for an integrable modified Camassa-Holm equation with cubic nonliearity. Commun. Math. Phys. 2013, 322, 967–997. [Google Scholar] [CrossRef]
  25. Lou, S.Y.; Qiao, Z.J. Alice-Bob Peakon Systems. Chin. Phys. Lett. 2017, 34, 100201. [Google Scholar]
  26. Claude, C.; Latifi, A.; Leon, J. Nonlinear resonant scattering and plasma instability: An integrable model. J. Math. Phys. 1991, 32, 3321. [Google Scholar] [CrossRef]
  27. Shao, Y.J.; Zeng, Y.B. The solutions of the NLS equations with self-consistent sources. J. Phys. A Math. Gen. 2005, 38, 2441–2467. [Google Scholar] [CrossRef]
  28. Meng, G.Q.; Qin, J.L.; Yu, G.L. Breather and rogue wave solutions for a nonlinear Schrödinger-type systems in plasmas. Nonlinear Dyn. 2015, 81, 739. [Google Scholar] [CrossRef]
  29. Mel’nikov, V.K. On equations for wave intections. Lett. Math. Phys. 1983, 7, 129–136. [Google Scholar] [CrossRef]
  30. Zeng, Y.B.; Shao, Y.J.; Ma, W.X. Integral-Type Darboux Transformation for the mKdV Hierarchy with Self-Consistent Sources. Commun. Theor. Phys. 2002, 38, 641–648. [Google Scholar]
  31. Huang, Y.H.; Yao, Y.Q.; Zeng, Y.B. On Camassa-Holm Equation with Self-Consistent Sources and Its Solutions. Commun. Theor. Phys. 2010, 53, 403–412. [Google Scholar]
  32. Yao, Y.Q.; Huang, Y.H.; Zeng, Y.B. The Qiao-Liu Equation with Self-Consistent Sources and Its Solutions. Commun. Theor. Phys. 2012, 57, 909–913. [Google Scholar] [CrossRef]
Figure 1. Elastic interactions via the solution (19a): (a) dark and antidark solitons with ρ = 1 , λ 1 = 0.1 i , c ( z ) = z and γ 1 = 0.4 0.5 i ; (b) two dark solitons with ρ = 1 , λ 1 = 0.1 i , c ( z ) = z 2 and γ 1 = 0.4 0.5 i ; (c) two dark solitons with ρ = 1 , λ 1 = 0.1 i , c ( z ) = e z and γ 1 = 0.4 0.5 i .
Figure 1. Elastic interactions via the solution (19a): (a) dark and antidark solitons with ρ = 1 , λ 1 = 0.1 i , c ( z ) = z and γ 1 = 0.4 0.5 i ; (b) two dark solitons with ρ = 1 , λ 1 = 0.1 i , c ( z ) = z 2 and γ 1 = 0.4 0.5 i ; (c) two dark solitons with ρ = 1 , λ 1 = 0.1 i , c ( z ) = e z and γ 1 = 0.4 0.5 i .
Mathematics 08 01099 g001
Figure 2. Elastic interactions via the solution (19a): (a) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z and γ 1 = 2 0.6 i ; (b) two antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z 2 and γ 1 = 2 0.6 i ; (c) two antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = e z and γ 1 = 2 0.6 i .
Figure 2. Elastic interactions via the solution (19a): (a) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z and γ 1 = 2 0.6 i ; (b) two antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z 2 and γ 1 = 2 0.6 i ; (c) two antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = e z and γ 1 = 2 0.6 i .
Mathematics 08 01099 g002
Figure 3. Elastic interactions via the solution (19a): (a) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z and γ 1 = 2 0.5 i ; (b) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z 2 and γ 1 = 2 0.5 i ; (c) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = e z and γ 1 = 2 0.5 i .
Figure 3. Elastic interactions via the solution (19a): (a) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z and γ 1 = 2 0.5 i ; (b) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = z 2 and γ 1 = 2 0.5 i ; (c) dark and antidark solitons with ρ = 1 , λ 1 = 0.4 i , c ( z ) = e z and γ 1 = 2 0.5 i .
Mathematics 08 01099 g003
Figure 4. Degenerate two-solution interaction via the solution (19a): (a) dark and antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z and γ 1 = 1 ; (b) antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z 2 and γ 1 = 1 ; (c) antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = e z and γ 1 = 1 .
Figure 4. Degenerate two-solution interaction via the solution (19a): (a) dark and antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z and γ 1 = 1 ; (b) antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z 2 and γ 1 = 1 ; (c) antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = e z and γ 1 = 1 .
Mathematics 08 01099 g004
Figure 5. Degenerate two-solution interaction via the solution (19a): (a) dark and antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z and γ 1 = 3 3 + i ; (b) dark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z 2 and γ 1 = 3 3 + i ; (c) dark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = e z and γ 1 = 3 3 + i .
Figure 5. Degenerate two-solution interaction via the solution (19a): (a) dark and antidark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z and γ 1 = 3 3 + i ; (b) dark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = z 2 and γ 1 = 3 3 + i ; (c) dark soliton with a phase shift between two soliton segments as z ± , where ρ = 1 , λ 1 = 0.5 i , c ( z ) = e z and γ 1 = 3 3 + i .
Mathematics 08 01099 g005
Figure 6. (ac) are the first-order rational soliton solutions via Equations (26a)–(26c) with ρ = 0.5 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 0.1 i , respectively.
Figure 6. (ac) are the first-order rational soliton solutions via Equations (26a)–(26c) with ρ = 0.5 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 0.1 i , respectively.
Mathematics 08 01099 g006
Figure 7. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 3 i , respectively.
Figure 7. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 3 i , respectively.
Mathematics 08 01099 g007
Figure 8. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 + 3 i , respectively.
Figure 8. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 + 3 i , respectively.
Mathematics 08 01099 g008
Figure 9. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 0.1 i , respectively.
Figure 9. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 0.1 i , respectively.
Mathematics 08 01099 g009
Figure 10. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 3 i , respectively.
Figure 10. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 3 i , respectively.
Mathematics 08 01099 g010
Figure 11. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 + 3 i , respectively.
Figure 11. (ac) are the first-order rational soliton solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 + 3 i , respectively.
Mathematics 08 01099 g011
Figure 12. Degenerate two-solution interaction via the solution in Equation (26a): (a) antidark soliton as z ± , where ρ = 0.5 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 ; (b) antidark soliton as z ± , where ρ = 0.5 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 .
Figure 12. Degenerate two-solution interaction via the solution in Equation (26a): (a) antidark soliton as z ± , where ρ = 0.5 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z and γ 1 = 1 ; (b) antidark soliton as z ± , where ρ = 0.5 , ϕ = 0 , c 1 ( z ) = 0 , c 2 ( z ) = z 2 and γ 1 = 1 .
Mathematics 08 01099 g012
Figure 13. Degenerate two-solution interaction via the solution Equation (26a): (a) antidark soliton as z ± , where ρ = 0.5 , ϕ = 0 , c ( z ) = z and γ 1 = 1 2 i ; (b) antidark soliton as z ± , where ρ = 0.5 , c ( z ) = z 2 and γ 1 = 1 2 i .
Figure 13. Degenerate two-solution interaction via the solution Equation (26a): (a) antidark soliton as z ± , where ρ = 0.5 , ϕ = 0 , c ( z ) = z and γ 1 = 1 2 i ; (b) antidark soliton as z ± , where ρ = 0.5 , c ( z ) = z 2 and γ 1 = 1 2 i .
Mathematics 08 01099 g013
Figure 14. (ac) are the rogue wave solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = z 2 + 3 2 , c 2 ( z ) = z 2 and γ 1 = 1 + 2 i , respectively.
Figure 14. (ac) are the rogue wave solution via Equations (26a)–(26c) with ρ = 1 , ϕ = 0 , c 1 ( z ) = z 2 + 3 2 , c 2 ( z ) = z 2 and γ 1 = 1 + 2 i , respectively.
Mathematics 08 01099 g014
Table 1. Asymptotic patterns of the solution (19a) under different parametric conditions.
Table 1. Asymptotic patterns of the solution (19a) under different parametric conditions.
Parametric ConditionsAsymptotic Soliton u 1 ± Asymptotic Soliton u 2 ±
c ( z ) λ 1 I γ 1 I > 0 , c ( z ) λ 1 I I m ( γ 1 κ 1 2 ) < 0 dark solitondark soliton
c ( z ) λ 1 I γ 1 I < 0 , c ( z ) λ 1 I I m ( γ 1 κ 1 2 ) > 0 antidark solitonantidark soliton
c ( z ) λ 1 I γ 1 I < 0 , c ( z ) λ 1 I I m ( γ 1 κ 1 2 ) < 0 antidark solitondark soliton
c ( z ) λ 1 I γ 1 I > 0 , c ( z ) λ 1 I I m ( γ 1 κ 1 2 ) > 0 dark solitonantidark soliton
γ 1 I = 0 , c ( z ) γ 1 R > 0 disappeardark soliton
γ 1 I = 0 , c ( z ) γ 1 R < 0 disappearantidark soliton
λ 1 I I m ( γ 1 κ 1 2 ) = 0 , c ( z ) λ 1 I γ 1 I < 0 antidark solitondisappear
λ 1 I I m ( γ 1 κ 1 2 ) = 0 , c ( z ) λ 1 I γ 1 I > 0 dark solitondisappear
Table 2. Asymptotic patterns of the solution (26a) under different parametric conditions.
Table 2. Asymptotic patterns of the solution (26a) under different parametric conditions.
Parametric ConditionsAsymptotic Soliton u I Asymptotic Soliton u II
I m ( s 1 ) < 0 , 1 + ρ I m ( s 1 ) > 0 RAD solitonRAD soliton
I m ( s 1 ) < 0 , 1 + ρ I m ( s 1 ) < 0 RAD solitonRD soliton
I m ( s 1 ) > 0 , 1 + ρ I m ( s 1 ) > 0 RD solitonRAD soliton
I m ( s 1 ) = 0 disappearRAD soliton
1 + ρ I m ( s 1 ) = 0 RAD solitondisappear

Share and Cite

MDPI and ACS Style

Huang, Y.; Jing, H.; Li, M.; Ye, Z.; Yao, Y. On Solutions of an Extended Nonlocal Nonlinear Schrödinger Equation in Plasmas. Mathematics 2020, 8, 1099. https://doi.org/10.3390/math8071099

AMA Style

Huang Y, Jing H, Li M, Ye Z, Yao Y. On Solutions of an Extended Nonlocal Nonlinear Schrödinger Equation in Plasmas. Mathematics. 2020; 8(7):1099. https://doi.org/10.3390/math8071099

Chicago/Turabian Style

Huang, Yehui, Hongqing Jing, Min Li, Zhenjun Ye, and Yuqin Yao. 2020. "On Solutions of an Extended Nonlocal Nonlinear Schrödinger Equation in Plasmas" Mathematics 8, no. 7: 1099. https://doi.org/10.3390/math8071099

APA Style

Huang, Y., Jing, H., Li, M., Ye, Z., & Yao, Y. (2020). On Solutions of an Extended Nonlocal Nonlinear Schrödinger Equation in Plasmas. Mathematics, 8(7), 1099. https://doi.org/10.3390/math8071099

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