Research Article | Open Access
Prolongation Structures and -Soliton Solutions for a New Nonlinear Schrödinger-Type Equation via Riemann-Hilbert Approach
In this paper, a new integrable nonlinear Schrödinger-type (NLST) equation is investigated by prolongation structures theory and Riemann-Hilbert (R-H) approach. Via prolongation structures theory, the Lax pair of the NLST equation, a matrix spectral problem, is derived. Depending on the analysis of red the spectral problem, a R-H problem of the NLST equation is formulated. Furthermore, through a specific R-H problem with the vanishing scattering coefficient, -soliton solutions of the NLST equation are expressed explicitly. Moreover, a few key differences are presented, which exist in the implementation of the inverse scattering transform for NLST equation and cubic nonlinear Schrödinger (NLS) equation. Finally, the dynamic behaviors of soliton solutions are shown by selecting appropriate spectral parameter , respectively.
In recent decades, nonlinear partial differential equations (PDEs) play a significant role in mathematics and theoretical physics, which have attracted much attentions in soliton theory and integrable system [1, 2]. Since exact solutions [3, 4] of the integrable equations [5, 6] can describe and explain many natural phenomena , the study of integrable equations has become a hot topic. It is well known that inverse scattering transform [8, 9] is a very important theory for solving exact solutions of integrable equations. There are two typical techniques for gaining inverse scattering: one is using Gelfand-Levitan-Marchenko equations and the other is formatting a Riemann-Hilbert(R-H) problem. The former has a complicated calculation procedure, while the latter can provide an equivalent and simpler method for solving integrable equations, especially soliton solutions [10–14]. R-H approach [15–19] is extensively applied in lots of nonlinear PDEs, for example, the coupled mKdV equation , the generalized Sasa-Satsuma equation , the general coupled nonlinear Schrödinger equation , which plays a vital role in dealing with initial boundary value problems [23–25], discussing long-time asymptotic behavior  and investigating the lump solutions .
We find that construction of the Lax pair is critical for solving the R-H problem. However, there is no unified way to build Lax pair so far. In 1975, Wahlquist and Estabrook  proposed prolongation structures theory [29–31]. This method is one of the effective skills to verify whether the equation is integrable or not. The theory has been widely applied to many integrable equations, such as coupled KdV equation [32, 33], coupled nonlinear Schröinger equation , and Heisenberg ferromagnet equation .
In this paper, we get Lax pair of NLST equation via prolongation structures theory, where a and b are real parameters describing the measure of derivative cubic nonlinearity, the “–” denotes complex conjugation. Moreover, we formulate a R-H problem by the obtained Lax pair and get the N-soliton solutions of (1). As we all know, the NLS equationis a vital model in physics, which plays a significant part in describing soliton propagation in water waves, nonlinear fiber optics plasma, etc. Similarly, system (1) can also describe different phenomena, such as collision and bound state. In particular, when setting , (1) is reduced to classical NLS equation (2). For convenience, we set in what follows.
This paper is structured as follows. In Section 2, Lax pair related to (1) is obtained via prolongation structures theory. In Section 3, the analytical properties for an equivalent spectral problem are analyzed; besides, a R-H problem related to a newly introduced spatial matrix spectral problem is formulated. Through the formalism of R-H problem, the construction of N-soliton solutions is discussed in Section 4. Finally, a few conclusions and some discussions are given.
2. Prolongation Structures of the NLST Equation
In this part, we study prolongation structures of equationwhich successfully gives the Lax pair of the NLST equation (3). The specific steps are as follows.
Firstly, we introduce an important proposition in the representation theory of Lie algebra.
Proposition 1 (see ). Let and be two elements of Lie , such that , and range ad . Then we may identity with and with , where are the nilpotent elements of and is the neutral element of .
Secondly, a new series of independent variables for (3) are defined by . Then (3) can be written as follows:And a set of differential 2-forms are defined on a differential manifold , whereIt is easy to testify the ideal is a closed ideal, that is , and limiting to solution manifold, we can get (3) and (4). Then, we introduce the differential 1-formswhere is linearly dependent on as well as , that is . For simplicity, is written as and is written as in what follows. According to the general theory of the exterior differential system, we obtain as a closed ideal, namelywhere are functions to be ensured and is differential 1-forms. Then, we obtain a suit of nonlinear PDEs of and :where . Then set in (8) as follows:where are matrices. Then, taking (9) into , we can get a form of Substituting (10) into (8), we gainby collecting the coefficients of . Setting , (11) can be written asThen we haveLetting , we finally obtainwhere constitutes a continuation algebra and they have the following relationships:Finally, we embed continuation algebra into Lie algebra . Through (15) and Proposition 1, we can get that and are nilpotent elements, and is neutral element, namely,Then, we take (16) into (15) yieldingMoreover, we obtain specific expressions for and :In summary, the Lax pair of (3) is
3. The Riemann-Hilbert Problem
Before we formulate a R-H problem, we provide a definition of Riemann-Hilbert problem.
Definition 2 (Riemann-Hilbert problem). Let the contour be the union of a finite number of smooth and oriented curves on the Riemann sphere , such that has only a finite number of connected components. Let be a matrix defined on the contour . The RHP is the problem of finding a matrix-valued function that satisfies the following:(i) is analytic for all and extends continuously to the contour ;(ii);(iii)
The Lax pair of the NLST equation (3) readswhere is a spectral parameter,and the “−” represents complex conjugation.
In this part, the inverse scattering transforms of the NLST equation (3) are given by formulating a R-H problem. The results will lay a foundation for the derivation of -soliton solutions in Section 4.
We suppose that the potential is smooth enough and decays to zero fast when In what follows, we treat as a matrix function. From (20a) and (20b), as , we notice that one possesses the asymptotic behavior: Thus it is convenient to introduce a new matrix spectral function which is defined as Inserting (22) into (20a) and (20b), the original form of Lax pair (20a) and (20b) becomeswithwhere , andwhere “” represents the Hermitian of a matrix.
Let us now consider formulating a correlated R-H problem with variable Firstly, in the scattering process, we introduce matrix Jost solutions for spectral problem (23a)which satisfy the asymptotic conditionsrespectively, where is a unit matrix, each (k=1,2) represents the -th column of the matrices , respectively. The subscript of indicates at which end of the -axis boundary conditions are set. Utilizing the parameter variation method and the boundary condition (28), the matrix spectral problem (23a) is transformed into the Volterra integral equations of Via a direct analysis of (29), due to the structure of the potential , it is easy to see the second column of involving the exponential factor , which, since , decays when . In addition, the first column of contains the exponential factor , which, since , also decays when . Therefore, we believe these two columns are analytic for and continuous for Similarly, the first column of and the second column of also can be analytic for and continuous for .
Remark 3. This implementation of the inverse scattering transform for the NLST is very similar to that for the cubic NLS  except for a few key differences. The difference here is that one instead needs to distinguish between the upper and lower half planes of And the other difference is the continuous of
Notice that , and we can use Abel’s formula to getThen, through the asymptotic condition (28), we gain for all Then, we introduce the notation Since is solution of (20a) and is solution as well, they must be linearly related, that is,where is the scattering matrix. Notice that because
Moreover, by analyzing the properties of , we find accepts analytic extensions to and extends to analytically.
The R-H problem requires (3) to contain two matrix functions: one is analytic in , and the other is analytic in . We construct the first matrix functionwhich is analytic in . Here, , . Furthermore, through the Volterra integral equations (29)-(30), we have
In order to formulate the R-H problem of (3), it is necessary to construct a matrix function , which is analytic in for . Actually, we consider the inverse matrices of where represent the -th row of , respectively. Moreover, satisfy the equationwhich is called adjoint scattering equation of (23a). The matrix function is introduced aswhich is analytic in . With a similar process as , we haveWe introducerelated towhereEquations (43)-(44) accurately give a R-H problem of a correlation matrix. From (35) and (40), we get the asymptotic properties of the above R-H problemand the canonical normalization condition
4. N-Soliton Solutions and Their Dynamics
In this part, we derive -soliton solutions for NLST equation (3). As is well known that R-H problem with zeros produces soliton solutions, the uniqueness of the solution of associated R-H problem defined in (43) and the zeros of are specified in their analytic domains, and the structure of at these zeros can be determined. From the definitions of equations (34) and (39) as well as the scattering relation (32), one getswhere
Let us assume that has zeros and has zeros . For convenience, we suppose that all zeros are simple zeros of . In what follows, each of , which only includes a single column vector ; each of , which only includes a single row vector , namely,
A key step for solving soliton solutions is to calculate the potential matrix through . Notice that is the solution of the spectral problem (23a); consequently, we consider the asymptotic expansion of ,submitting (49) into (23a) and comparing term, we obtainThus, a direct calculation shows that the potentials is represented asThe potential matrix has symmetry property, which produces in scattering matrix and Jost functions. In addition, the scattering equation (23a) has Hermitian property and we getNotice that meet the adjoint scattering equation (38); recalling also meets (38), and using the boundary conditions (28) we obtainThrough this involution property and definitions of (34) as well as (39) for , we find the involution property is also fit for the analytic solutions :
Furthermore, from the scattering relationship (32) between and , the involution property also applies to :In view of the involution property, we gain the relationfor the zeros of and . Without loss of generality, taking equations (48a), (54), (56), we haveThe derived special R-H problem is solved by referring to references and the results can be expressed as is an matrix given by
Notice that zeros and are constants; that is, spatial variable and time variable are independent. Thus, it is easy to determine the spatial and temporal evolution of vectors and We calculate the derivatives of on both sides of equation Taking advantage of equation (23a), one getsthus In a completely similar way, the time dependence of :can be determined by matrix spectral problem (20b). Summing up, we havewhere is an arbitrary constant column vector, and is an arbitrary constant row vector.
Finally, from (58a) and (58b), we obtainand from (51), the general -soliton solutions of NLST equation readwhere vectors are given by (63), , and matrix is given by (59). In order to get the soliton solutions of the NLST equation (3) conveniently, we let Furthermore, we introduce the notation Then the above solution is expressed explicitly aswhere
4.1. Single-Soliton Solutions
Solution (70) is a solitary wave in the NLST equation (1). It is easy to see that the shape of the amplitude function is hyperbolic secant, the peak amplitude is , and the velocity is . The phase of this solution depends linearly on time and space . In addition, the parameters and represent the initial position and the phase of the solitary wave, respectively. Notice that the soliton power depends on only; therefore, it does not change after soliton collisions. The solution is shown graphically in Figure 1.
Remark 4. One can see that the single-soliton solution (70) of the NLST equation and the single-soliton solution of the cubic NLS are different. From the difference in the Lax pair, one can follow through to see that eigenvalues of the cubic NLS (I will call them ) corresponding to in this NLST. Thus since in the cubic NLS corresponds to a pulse width of and velocity of , in the NLST we have corresponding to , so the pulse width will be and the velocity will be .
4.2. Two-Soliton Solutions
We set in formula (67), and the two-soliton solutions also can be expressed explicitly, while it is quite complicated. Using algebraic theory, two-soliton solutions are represented aswhereEquation (71) is read aswhereDenoting , is the real part of and is the imaginary part of Two cases are given below: one case is , which means that the two solitons have different velocities, and they form a collision, and the other case is , which means that the two solitons have same velocities, and they form a bound state.
Case I. We choose
In this case, when , this solution consists of two single-solitons that are far apart and move toward each other. They interact when they collide. When , these solitons reappear from the interaction without any change in shape and velocity. The two-soliton solutions are shown graphically in Figure 2.
Case II. We choose
In this case, we assume , the two solitons have same velocity; therefore, they will stay together and form a bound state. In addition, the amplitude function has periodic oscillations with time as shown in Figure 3.
In summary, we obtained the N-soliton solutions (67) by the R-H approach. In addition, as an example, the expressions of the single-soliton solutions (69) and the two-soliton solutions (71) are given explicitly. Moreover, we discuss the dynamic behaviors of two solitons. By giving two examples, two different phenomena are shown, one is , which describes a collision; the other is , which describes a bound state. Similarly, we can get three-soliton solutions, four-soliton solutions, etc. It is also interesting to study their dynamic behaviors in the future.
In this work, a new integrable NLST equation is investigated via prolongation structures theory and R-H approach. We apply the prolongation structures theory to the NLST equation; by discussing the continuation algebra of the NLST equation, a Lax pair is successfully derived. Then, we begin with the spectral analysis of the Lax pair and formulate a R-H problem. Through a specific R-H problem with vanishing scattering coefficient, which is under the condition of the reflection-less case, -soliton solutions are obtained explicitly. In addition, we present a few key differences between the NLST equation and the cubic NLS equation in the implementation of inverse scattering transform. Moreover, we analyze the dynamic behaviors between two solitons, which describe different phenomena, such as collision, bound state, and the dynamic behaviors of the one-soliton solution and two-soliton solutions are shown by graphs, respectively.
Recently, we notice that there are many other approaches to obtain exact solutions in the field of integrable systems, like Hirota’s bilinear method , Darboux transformation method , Wronskian technique , Frobenius integrable decompositions , etc. . Thus we consider whether the NLST equation can be solved with these approaches. These questions will be discussed in our future work.
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was supported by the Shandong Provincial Natural Science Foundation (Grant No. ZR2019QD018), National Natural Science Foundation of China (Grant No. 61602188), and Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (Grant Nos. 2017RCJJ068 and 2017RCJJ069).
- W. X. Ma, H. Wu, and J. He, “Partial differential equations possessing Frobenius integrable decompositions,” Physics Letters A, vol. 364, no. 1, pp. 29–32, 2007.
- X. Wang, Y. Fang, and H. Dong, “Component-trace identity for Hamiltonian structure of the integrable couplings of the Giachetti-Johnson (GJ) hierarchy and coupling integrable couplings,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 7, pp. 2680–2688, 2011.
- L. Zhang, Y. Wang, C. M. Khalique, and Y. Bai, “Peakon and cuspon solutions of a generalized Camassa-Holm-Novikov equation,” Journal of Applied Analysis and Computation, vol. 8, no. 6, pp. 1938–1958, 2018.
- Y. Liu, H. Dong, and Y. Zhang, “Solutions of a discrete integrable hierarchy by straightening out of its continuous and discrete constrained flows,” Analysis and Mathematical Physics, vol. 9, pp. 465–481, 2019.
- M. Guo, H. Dong, J. Liu, and H. Yang, “The time-fractional mZK equation for gravity solitary waves and solutions using sech-tanh and radial basis function method,” Nonlinear Analysis: Modelling and Control, vol. 24, no. 1, 19 pages, 2019.
- M. Tao and H. Dong, “Algebro-geometric solutions for a discrete integrable equation,” Discrete Dynamics in Nature and Society, vol. 2017, 9 pages, 2017.
- M. Tao, N. Zhang, D. Gao, and H. Yang, “Symmetry analysis for three-dimensional dissipation Rossby waves,” Advances in Difference Equations, vol. 300, 2018.
- M. J. Ablowitz and Z. H. Musslimani, “Inverse scattering transform for the integrable nonlocal nonlinear Schrödinger equation,” Nonlinearity, vol. 29, no. 3, pp. 915–946, 2016.
- H. C. Morris, “Prolongation structures and a generalized inverse scattering problem,” Journal of Mathematical Physics, vol. 17, no. 10, pp. 1867–1869, 1976.
- H. Yang, J. Sun, and C. Fu, “Time-fractional Benjamin-Ono equation for algebraic gravity solitary waves in baroclinic atmosphere and exact multi-soliton solution as well as interaction,” Communications in Nonlinear Science and Numerical Simulation, vol. 71, pp. 187–201, 2019.
- M. Guo, C. Fu, Y. Zhang, J. Liu, and H. Yang, “Study of ion-acoustic solitary waves in a magnetized plasma using the three-dimensional time-space fractional Schamel-KdV equation,” Complexity, vol. 2018, Article ID 6852548, 2018.
- M. Han, L. Zhang, Y. Wang, and C. M. Khalique, “The effects of the singular lines on the traveling wave solutions of modified dispersive water wave equations,” Nonlinear Analysis: Real World Applications, vol. 47, pp. 236–250, 2019.
- W. Q. Peng, S. F. Tian, T. T. Zhang, and Y. Fang, “Rational and semi-rational solutions of a nonlocal (2+1)-dimensional nonlinear Schrodinger equation,” Mathematical Methods in the Applied Sciences, 2019.
- H. W. Yang, M. Guo, and H. He, “Conservation laws of space-time fractional mZK equation for Rossby solitary waves with complete Coriolis force,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 20, pp. 17–32, 2018.
- V. S. Shchesnovich and J. Yang, “General soliton matrices in the Riemann-Hilbert problem for integrable nonlinear equations,” Journal of Mathematical Physics, vol. 44, no. 10, pp. 4604–4639, 2003.
- D.-S. Wang, D.-J. Zhang, and J. Yang, “Integrable properties of the general coupled nonlinear Schrödinger equation,” Journal of Mathematical Physics, vol. 51, no. 2, 17 pages, 2010.
- Z. Kang, T. Xia, and X. Ma, “Multi-soliton solutions for the coupled modified nonlinear Schrödinger equations via Riemann–Hilbert approach,” Chinese Physics B, vol. 27, no. 7, p. 070201, 2018.
- X. G. Geng and H. Nie, “N solutions for a derivative nonlinear Schrodinger-type equation via Riemann-Hilbert approach,” Mathematical Methods in the Applied Sciences, vol. 41, pp. 1653–1660, 2018.
- Y. Lin, H. Dong, and Y. Fang, “N-Soliton Solutions for the NLS-Like Equation and Perturbation Theory Based on the Riemann–Hilbert Problem,” Symmetry, vol. 11, no. 6, p. 826, 2019.
- W. X. Ma, “Riemann-Hilbert problems and N-soliton solutions for a coupled mKdV system,” Journal of Geometry and Physics, vol. 132, pp. 45–54, 2018.
- X. Geng and J. Wu, “Riemann-Hilbert approach and N-soliton solutions for a generalized Sasa-Satsuma equation,” Wave Motion. An International Journal Reporting Research on Wave Phenomena, vol. 60, pp. 62–72, 2016.
- D.-S. Wang, D.-J. Zhang, and J. Yang, “Integrable properties of the general coupled nonlinear Schrödinger equation,” Journal of Mathematical Physics, vol. 51, no. 2, Article ID 023510, 17 pages, 2010.
- J. Xu and E. Fan, “The initial-boundary value problem for the ostrovsky-vakhnenko equation on the half-line,” Mathematical Physics, Analysis and Geometry, vol. 19, no. 3, p. 17, 2016.
- B. Hu, T. Xia, and W. Ma, “Riemann–Hilbert approach for an initial-boundary value problem of the two-component modified Korteweg-de Vries equation on the half-line,” Applied Mathematics and Computation, vol. 332, pp. 148–159, 2018.
- T. Liu and H. Dong, “The prolongation structure of the modified nonlinear schrödinger equation and its initial-boundary value problem on the half line via the riemann-hilbert approach,” Mathematics, vol. 7, no. 2, 2019.
- J. Xu, E. G. Fan, and Y. Chen, “Long-time asymptotic for the derivative nonlinear Schroinger equation with step-like initial value,” Mathematical Physics, Analysis and Geometry, vol. 16, pp. 253–288, 2013.
- Y. Zhang, H. Dong, X. Zhang, and H. Yang, “Rational solutions and lump solutions to the generalized (3+1)-dimensional Shallow Water-like equation,” Journal of Computers & Mathematics with Applications, vol. 73, no. 2, pp. 246–252, 2017.
- H. D. Wahlquist and F. B. Estabrook, “Prolongation structures of nonlinear evolution equations,” Journal of Mathematical Physics, vol. 16, pp. 1–7, 1975.
- H. C. Morris, “A prolongation structure for the AKNS system and its generalization,” Journal of Mathematical Physics, vol. 18, no. 3, pp. 533–536, 1977.
- D. Wang, Y. Ma, and X. Li, “Prolongation structures and matter-wave solitons in F=1 spinor Bose-Einstein condensate with time-dependent atomic scattering lengths in an expulsive harmonic potential,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 10, pp. 3556–3569, 2014.
- Z. Yan and C. Li, “Fermionic covariant prolongation structure theory for a (2+1)-dimensional super nonlinear evolution equation,” International Journal of Geometric Methods in Modern Physics, vol. 15, no. 07, p. 1850114, 2018.
- Y. Jia and W. Duan, “The Prolongation Structure of a Coupled Kdv Equation,” Acta Mathematicae Applicatae Sinica, vol. 34, no. 2, pp. 430–437, 2018.
- Y. Cao and D. Wang, “Prolongation structures of a generalized coupled Korteweg-de Vries equation and Miura transformation,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2344–2349, 2010.
- D. Wang, S. Yin, Y. Tian, and Y. Liu, “Integrability and bright soliton solutions to the coupled nonlinear Schrödinger equation with higher-order effects,” Applied Mathematics and Computation, vol. 229, pp. 296–309, 2014.
- W.-Z. Zhao, Y.-Q. Bai, and K. Wu, “Generalized inhomogeneous Heisenberg ferromagnet model and generalized nonlinear Schrödinger equation,” Physics Letters A, vol. 352, no. 1-2, pp. 64–68, 2006.
- M. Wadati, K. Konno, and Y. Ichikawa, “Lax representation for the generalized nonlinear Schrodinger equation,” Journal of the Physical Society of Japan, vol. 24, p. 1965, 1979.
- J. E. Humphreys, Introduction to Lie Algebras and Representation Theory, Springer, New York, NY, USA, 1972.
- S. Cuccagna and D. E. Pelinovsky, “The asymptotic stability of solitons in the cubic NLS equation on the line,” Applicable Analysis: An International Journal, vol. 93, no. 4, pp. 791–822, 2012.
- W.-X. Ma and Y. Zhou, “Lump solutions to nonlinear partial differential equations via Hirota bilinear forms,” Journal of Differential Equations, vol. 264, no. 4, pp. 2633–2659, 2018.
- N. Zhang, T. Xia, and Q. Jin, “N-Fold Darboux transformation of the discrete Ragnisco–Tu system,” Advances in Difference Equations, vol. 2018, no. 1, 2018.
- S. Jia, Y. Gao, C. Ding, and G. Deng, “Solitons for a (2+1)-dimensional Sawada–Kotera equation via the Wronskian technique,” Applied Mathematics Letters, vol. 74, pp. 193–198, 2017.
- Y. Fang, H. Dong, Y. Hou, and Y. Kong, “Frobenius integrable decompositions of nonlinear evolution equations with modified term,” Applied Mathematics and Computation, vol. 226, no. 1, pp. 435–440, 2014.
- B. Song and Y. Kong, “Solution to PSPACE-complete problem using p systems with active membranes with time-freeness,” Mathematical Problems in Engineering, vol. 2019, 8 pages, 2019.
Copyright © 2019 Yuxin Lin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.