Abstract

A Riemann-Hilbert approach is developed to the multicomponent Kaup-Newell equation. The formula is presented of -soliton solutions through an identity jump matrix related to the inverse scattering problems with reflectionless potential.

1. Introduction

Many nonlinear partial differential equations especially soliton equations have exact solutions [13]. There are a lot of methods to solve soliton equations such as the Hirota bilinear method [26], Wronskian technique (Casoratian technique) [59], and Darboux transformation [10, 11]. The inverse scattering transformation (IST) is one of the most powerful tools and closely connected with those methods mentioned above [1, 3]. It is also called nonlinear Fourier transform for its procedure to solve the nonlinear equations is similar to the linear Fourier transform. One advantage of the IST is that it can be applied to the whole soliton hierarchies [3]. Recently, researches show that the IST can solve not only classic soliton hierarchies but also soliton equations with self-consistent sources [12], nonisospectral soliton hierarchies [13], hierarchies mixed with isospectral and nonisospectral ones [14], and nonlocal soliton hierarchies [15]. Furthermore, it can generate both soliton and general matrix exponent solutions [16, 17].

The Riemann-Hilbert (RH) approach is another effective method to solve soliton equations. It actually shares a close relationship with the IST [1820]. Both of them start from same matrix spectral problems which possess bounded eigenfunctions analytically extendable to the upper or lower half-plane. To get scattering data, we must consider the asymptotic conditions at infinity on real axis by the IST to solve soliton equations. In fact, the considered conditions are used as the solutions to the corresponding RH problems. When the jump matrix is an identity matrix, the RH problem is equivalent to the IST with reflectionless potentials, and -soliton solutions can be generated [2123]. Recently, Ma has already used the method to solve multicomponent soliton equations such as multicomponent AKNS integrable hierarchies and a coupled mKdV equation [2426].

It is known to us all that the three famous derivative nonlinear Schrödinger equations, the Chen-Lee-Liu (CLL) equation [17, 27], Kaup-Newell (KN) equation [28], and Gerdjikov-Ivanov (GI) equation [29, 30], can be reduced from the Kundu equation by choosing different value of the arbitrary parameter [31, 32]. Many properties of them have been researched such as exact solutions [30, 33], conservation laws [34], multi-Hamilton structure [31], and -symmetry algebra [32, 34].

In this paper, we will present the multicomponent KN equation with its matrix Lax pairs. To formulate an RH problem of the equation, we consider a modify matrix Lax pairs. The formula of generating the -soliton solutions to multicomponent KN equation will be obtained through taking the identify jump matrix.

The paper is organized as follows. In Section 2, we will introduce the multicomponent isospectral KN equation and its Lax pairs. In Section 3, we will construct a multicomponent RH problem to the equation introduced in the previous section. In Section 4, the expression of -soliton solutions will be obtained. We conclude the paper in Section 5.

2. The Multicomponent KN Equation

In this section, we will present the isospectral multicomponent KN equation from a matrix spectral problem by the zero-curvature representation. To our knowledge, there is another powerful method to build soliton equation hierarchies through Kac-Moody algebra and principal gradation [35].

Suppose that and are smooth functions of variables and , denotes the transpose of matric, and is an identity matrix. Let us consider the following Lax pairs

where and are two real constants; , , and are potential functions; is a spectral parameter; and

Obviously, and are smooth component functions of variables and . Assume that , and their derivatives of any order with respect to vanish rapidly as .

The compatibility condition of (1), i.e., the zero curvature equation generates the multicomponent KN soliton equation

For example, when , the spectral problem (1a) becomes where

Its time evolution is with

The 4-component KN equation is

3. The RH Problem to the Multicomponent KN Equation

In this section, we will build the RH problem to the multicomponent KN equation (5). Here, we only focus on the positive flows. Constructing the RH problem from negative symmetry flows have already appeared in [36] for the homogeneous -hierarchy and its .

Setting

it is obvious that the trace of is zero, where

Thus, the equation (5) has the following Lax pairs

Next, we will present the scattering and inverse scattering methods for the multicomponent KN equation (5) by the RH approach. The resulting results will lay the groundwork for -soliton solutions in the next section. Suppose that all the potentials rapidly vanish when or and satisfy

In the RH approach, we treat in the spectral problem (1a) as a fundament matrix. From (12), we note, under (13), one has the asymptotic behavior: . This motivate us to introduce the variable transformation to have the canonical normalization for the associated RH problem:

where is the identity matrix. This way, the spectral problems in (12) equivalently lead to

where and . Noticing , we have by the Abel’s formula.

Let us now consider the formulation of an associated RH problem with the variable . In the scattering problem, we first introduce the matrix solutions of (16a) with the asymptotics conditions

respectively. The subscripts above refer to which end of the -axis the boundary conditions are required. Then, by (17), we have the determinant for all . Since are both solutions of (12), they must be linearly related, and so, we have where

is the scattering matrix. Note that since . Using the method of variation in parameters as well as the boundary condition (19), we can turn the -part of (12) into the following Volterra integral equation for :

Thus, allows analytical continuations off the real axis as long as the integrals on their right hand sides converge. Taking , it is direct to see that the integral equation for the first column of contains only the exponential factor . When is in the first or third quadrant, i.e., , let . Then, due to in the integral decays as , and the integral equation for the last columns of contains only the exponential factor , which due to in the integral, also decays when is in the upper half-plane . Thus, these columns can be analytically continued to the first or third quadrants. Similarly, we find that the last columns of and the first column of can be analytically continued to the second and fourth quadrants. Let us express

where stands for the th column of . Then, the matrix solution

is analytic in the first and third quadrants of , and the matrix solution

is analytic in the second and fourth quadrants of , where and . In addition, from the Volterra integral equation (21), we know that

Next, we construct the analytic counterpart of in the second and fourth quadrants of . Note that the adjoint equation of the -part of (12) and the adjoint equation of (16) read as

It is easy to see that the inverse matrices and solve these adjoint equations, respectively. If we express as follows:

where is the th row of . Then, by similar arguments, we can show that the adjoint matrix solution

is analytic when is in the second or fourth quadrants, and the other matrix solution

is analytic for in first and third quadrants. In the same way, we see that

Now, we have constructed two matrix functions and , which are analytic in the first or third quadrants and second or fourth quadrants, respectively. Defining

we can easily find that if on the real axis or imaginary axis, the two matrix functions and are related by where

Eq. (32) and Eq.(33) are exactly the associated matrix RH problem we wanted to present. The asymptotic conditions provide the canonical normalization condition for the established RH problem.

To finish the direct scattering transform, we take the derivative of (19) with time and use the vanishing conditions of the potentials; we can show that satisfies which gives the time evolution of the scattering coefficients:

and the other scattering data do not depend on time .

4. -Soliton Solutions

The RH problems with zeros can generate soliton solutions. The uniqueness of the associated RH problem (32) does not hold unless the zeros of det in the first or third quadrants and det in the second or fourth quadrants are specified and kernel structures of at these zeros are determined. Following the definitions of as well as the scattering relation between and , we find that

where we have used the fact . Similarly, we have and

Suppose that has zeros , and has zeros . For simplicity, we assume that all these zeros, and , are simple. Then, each of ker contains only a single column vector, denoted by , and each of ker contains a row vector, denoted by :

The RH problem (32) with the canonical normalization condition (35) and the zero structure (40) can be solved explicitly, and thus, one can readily reconstruct the potential as follows. Note that is a solution to the spectral problem (16). Therefore, as long as we expand at large as

inserting this expansion into (16) and comparing terms lead to which implies that where . Further, the potentials and , can be computed as

To obtain soliton solutions, we set in the RH problem (32). This can be achieved if we assume , which means that there is no reflection in the scattering problem. The solutions to this specific RH problem can be given as follows [24, 25] where is a square matrix whose entries read

Noting that the zeros and are constants, i.e., space and time independent, we can easily find the spatial and temporal evolutions for the vectors, and . For example, let us take the -derivative of both sides of the equation . By using (16) and , we get which implies

The time dependence of :

can be determined similarly through an associated RH problem with the variable . Summing up, we obtain where and , are arbitrary constant vectors. Finally, from (45), we get

and thus by (44), the -soliton solution to the system of multicomponent KN equations (4): where and , are arbitrary.

5. Conclusions

In general, we construct the RH problem for the multicomponent KN equation in this paper. To build the special RH problem with the identity jump matrix, we introduced a variable transformation to canonical normalization spectral problem. By recombining the solutions of the canonical spectral problem and its adjoint spectral problem, a general jump matrix to the special RH problem was constructed. Letting the general jump matrix to be identity jump matrix, the RH problem was solved. Finally, we obtained the expression of the -soliton solutions through power series expansion of the spectral parameter in the canonical normalization spectral problem.

In this method, the jump matrix is corresponding to the scattering matrix, and the identity jump matrix is equivalent to reflectionless coefficient of the IST. It is well known that there are not only soliton solutions to soliton equations but also rational solutions, Matveev solutions, complexiton solutions, and so on. Recently, there have been active studies on lumps and their interaction solutions with solitons [37, 38]. It would be very interesting to generalize this method to ()-dimensional equations and consider their lumps and interaction solutions. These will be our future projects.

Data Availability

The data that supports the findings of this study are available within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The work was supported in part by NSFC under the grants 11101350, 11671177, and 11771186; the Jiangsu Qing Lan Project (2014); and the Six Talent Peaks Project of Jiangsu Province (2016-JY-08).