^{1}

^{*}

^{1}

^{1}

^{1}

In this paper, we investigate an integrable high order nonlocal coupled Ablowitz-Kaup-Newell-Segur (AKNS) system for the first time. With the aid of Lax pair of this nonlocal system, Darboux transformation (DT) and new soliton-like solutions are obtained. Different from local equations, Darboux transformation of nonlocal systems needs to meet certain conditions. In this article, under the condition of symmetry reduction, the components of Darboux transformation need to satisfy . In order to study the dynamic information of the solutions, the images of the solutions are given.

In the past decades, a number of researches involved in applied mathematics and theoretical physics are concerned with integrable systems. The KdV, the nonlinear Schrödinger, the mKdV, the sine-Gordon equations, etc. are well-known integrable equations, and they can all be derived from the AKNS system. Many scholars have studied these equations from different fields and got many excellent results. An important research topic is to find the exact solutions of these systems, especially their soliton solutions or soliton-like solutions. So far, many effective methods have been developed, such as, inverse scattering method [

It is well known that symmetry is one of the important characteristics of integrable systems. By using Lie group method, the system can be reduced and solved, the conservation law of the system can also be constructed at the same time. Lie symmetry is also called continuous symmetry because its variables are transformed in the following form, x * = x + ε ξ , t * = t + ε τ , u * = u + ε ψ , where ε is infinitesimal. However, many facts show that two or more events occur at different places and times but are closely related to each other. For example, the reduction of Arctic sea ice in the summer of 2007 is considered to be the reason for the heavy winter freezing rain in southern China in 2008. To explain these events, nonlocal models are proposed, such as nonlocal nonlinear Schrödinger equation [

i q t ( x , t ) + q x x ( x , t ) + 2 q 2 ( x , t ) q * ( − x , t ) = 0, (1.1)

was proposed by Ablowitz, where ∗ represents complex conjugation. One can see from Equation (1.1) that q is not only a function of x but also a function of −x, that is to say, the solution q ( x , t ) evaluated at x requires information from −x via the term q * ( − x , t ) . For system (1.1), the exact solution and the conservation law are studied by the inverse scattering method. Before long, the nonlocal modified KdV systems, the discrete nonlocal NLS system, nonlocal sine-Gordon equation and the nonlocal Davey-Stewartson systems were also proposed. In general, nonlocal equations fall into three categories: 1) full PT symmetry; 2) partial PT symmetry; 3) partial reverse space-time symmetry. For this new type of system, many scholars have studied the integrability, exact solution [

Constructing the exact solutions of these nonlinear nonlocal equations can help to verify the numerical solutions and help to analyze the stability of the solutions. However, it is very difficult to solve them because of their high nonlinearity. Over the past few decades, significant progress has been made in many methods. The Darboux transform based on Lax pair has been proved to be one of the most effective algorithms for solving explicit solutions of nonlinear evolution equations. By using Darboux transformation method, soliton solutions, rogue-wave solutions of those nonlocal integrable equations are given [

This paper is arranged as follows: In Section 2, the high order nonlocal coupled AKNS system is constructed by using the linear spectral problem. In Section 3, the process of Darboux transformations is introduced in detail. Exact 1-fold and N-fold solutions are obtained. In Section 4, using the trivial seed solutions, the exact solutions of the high order nonlocal coupled AKNS system are constructed. In order to study the properties of the solution, the appropriate parameters are selected and the corresponding images are made. Finally, some conclusions and discussions are given in Section 5.

In this section, we will construct the high order nonlocal coupled AKNS system. In fact, the starting point is the AKNS scattering problem [

ϕ x = U ϕ , ϕ t = V ϕ , (2.1)

ϕ = ( ϕ 1 ( x , t ) , ϕ 2 ( x , t ) , ϕ 3 ( x , t ) , ϕ 4 ( x , t ) ) T , and the matrices U , V have the following form matrixes,

U = ( − λ P − λ R Q λ S λ 0 0 − λ P 0 0 Q λ ) , V = ( A B E F C − A G − E 0 0 A B 0 0 C − A ) , (2.2)

which was first introduced by Ma et al., A , B , C , E , F , G have the following form,

A = α ( λ 3 − 1 4 p q x + 1 4 q p x − 1 2 λ p q ) , B = α ( − λ 2 p + 1 2 λ p x − 1 4 p x x + 1 2 p 2 q ) , C = α ( − λ 2 q − 1 2 λ q x − 1 4 q x x + 1 2 p q 2 ) , E = α ( λ ( p q − p s 2 − q r 2 ) + 3 4 p q x − 3 4 q p x − 1 4 p s x + 1 4 q r x − 1 4 r q x + 1 4 s p x ) + β ( λ 3 − 1 4 p q x + 1 4 q p x − 1 2 λ p q ) ,

F = α ( λ 2 ( p − r ) − λ ( p x − 1 2 r x + 3 4 p x x − 1 4 r x x − 3 2 p 2 q + 1 2 p 2 s + p q r ) ) + β ( − λ 2 p + 1 2 λ p x − 1 4 p x x + 1 2 p 2 q ) , G = α ( λ 2 ( q − s ) − λ ( q x − 1 2 s x + 3 4 q x x − 1 4 s x x − 3 2 p q 2 + 1 2 q 2 r + p q s ) ) + β ( − λ 2 q + 1 2 λ q x − 1 4 q x x + 1 2 p q 2 ) , (2.3)

where p , q , r , s are four potentials, α is arbitrary constant and λ is the spectral parameter.

The compatibility condition of the system (2.1) yields the zero curvature representation U t − V x + [ U , V ] = 0 , using this formula, one can obtain the following evolution equations,

p t = 3 2 α p q p x − 1 4 α p x x x , q t = 3 2 α p q q x − 1 4 α q x x x , r t = 3 4 α p x x x − 1 4 α r x x x − 1 4 β p x x x + 3 2 β p q p x − 9 2 α p q p x + 3 2 α p q r x + 3 2 α p s p x + 3 2 α q r p x ,

s t = 3 4 α q x x x − 1 4 α s x x x − 1 4 β q x x x + 3 2 β p q q x − 9 2 α p q q x + 3 2 α p q s x + 3 2 α q r q x + 3 2 α p s q x , (2.4)

system (2.4) is coupled mixed nonlinear Schrödinger equation, it possesses infinitely many symmetries, which constitute some infinite-dimensional Lie algebra.

In this section, we set a symmetry reduction to (2.4),

q ( x , t ) = p ( − x , − t ) , s ( x , t ) = r ( − x , − t ) , (2.5)

the high order nonlocal coupled AKNS systemn can be obtained,

p t = 3 2 α p p x p ( − x , − t ) − 1 4 α p x x x , r t = 1 4 ( 3 α p x x x − α r x x x − β p x x x ) + 3 2 p p x ( β p ( − x , − t ) − 3 α p ( − x , − t ) + α r ( − x , − t ) ) + 3 2 α p ( − x , − t ) ( p r ) x , (2.6)

The integrable high order nonlocal coupled AKNS system is constructed by using the symmetry reduction method. The following work is to construct the exact solution of the system (2.6). It is well known that the Darboux transform is an important method for solving the exact solution of integrable nonlinear systems. It is proved that this method is not only suitable for solving local equations, but also for solving nonlocal equations. In the next section, we will construct the Darboux transform and exact solution of the nonlocal integrable system (2.6).

In this section, we will construct the Darboux transformation (DT) by means of the lax pair of high order nonlocal coupled AKNS system. The idea of the DT is to construct a transformation such that the form of the lax pair remains unchanged. Similarly to classical DT, we take the following gauge transformation,

ϕ [ 1 ] = T [ 1 ] ϕ . (3.1)

Assume that the new lax pair has the following form under the gauge transformation (4.1),

ϕ x [ 1 ] = U [ 1 ] ϕ [ 1 ] , ϕ t [ 1 ] = V [ 1 ] ϕ [ 1 ] , (3.2)

where

U [ 1 ] = ( T x [ 1 ] + T [ 1 ] U ) ( T [ 1 ] ) − 1 , V [ 1 ] = ( T t [ 1 ] + T [ 1 ] V ) ( T [ 1 ] ) − 1 , (3.3)

the next work is to find a matrix T [ 1 ] such that U [ 1 ] , V [ 1 ] have the same forms as U , V , and p , q , r , s are mapped into new potentials p [ 1 ] , q [ 1 ] , r [ 1 ] , s [ 1 ] , i.e.

U [ 1 ] = ( − λ p [ 1 ] − λ r [ 1 ] q [ 1 ] λ s [ 1 ] λ 0 0 − λ p [ 1 ] 0 0 q [ 1 ] λ ) , V [ 1 ] = ( A [ 1 ] B [ 1 ] E [ 1 ] F [ 1 ] − A [ 1 ] G [ 1 ] − E [ 1 ] 0 0 A [ 1 ] B [ 1 ] 0 0 C [ 1 ] − A [ 1 ] ) , (3.4)

where A [ 1 ] , B [ 1 ] , C [ 1 ] , E [ 1 ] , F [ 1 ] , G [ 1 ] have the same forms as A , B , C , E , F , G (2.3), just like the forms of U [ 1 ] .

In order to construct the form of the gauge matrix T, we suppose,

T [ 1 ] = ( T 0 T 1 0 T 0 ) , (3.5)

where T 0 and T 1 have the following forms,

T 0 = ( λ + b 11 [ 1 ] ( x , t ) b 12 [ 1 ] ( x , t ) b 21 [ 1 ] ( x , t ) λ + b 22 [ 1 ] ( x , t ) ) , T 1 = ( λ + b 13 [ 1 ] ( x , t ) b 14 [ 1 ] ( x , t ) b 23 [ 1 ] ( x , t ) λ + b 24 [ 1 ] ( x , t ) ) ,

b i j [ 1 ] , i = 1 , 2 ; j = 1 , 2 , 3 , 4 are undetermined functions of x and t.

By substituting (3.5) into (3.3), we obtain the transformations between p , q , r , s and p [ 1 ] , q [ 1 ] , r [ 1 ] , s [ 1 ] by balancing the order of λ ,

p [ 1 ] ( x , t ) = p ( x , t ) + 2 b 12 [ 1 ] ( x , t ) , q [ 1 ] ( x , t ) = q ( x , t ) − 2 b 21 [ 1 ] ( x , t ) , r [ 1 ] ( x , t ) = r ( x , t ) + 2 b 14 [ 1 ] ( x , t ) , s [ 1 ] ( x , t ) = s ( x , t ) − 2 b 23 [ 1 ] ( x , t ) , (3.6)

and

∂ ∂ x b 11 = p [ 1 ] b 21 − b 12 q , ∂ ∂ x b 12 = p [ 1 ] b 22 − b 11 p , ∂ ∂ x b 13 = p [ 1 ] b 23 + r 1 b 21 − b 12 s − b 14 q , ∂ ∂ x b 14 = p [ 1 ] b 24 + r 1 b 22 − b 11 r − b 13 p , ∂ ∂ x b 21 = q [ 1 ] b 11 − b 22 q , ∂ ∂ x b 22 = q [ 1 ] b 12 − b 21 p , ∂ ∂ x b 23 = q [ 1 ] b 13 + s 1 b 11 − b 22 s − b 24 q , ∂ ∂ x b 24 = q [ 1 ] b 14 + s 1 b 12 − b 21 r − b 23 p , (3.7)

because q [ 1 ] ( x , t ) , p [ 1 ] ( x , t ) , r [ 1 ] ( x , t ) , s [ 1 ] ( x , t ) are still solutions to the Equation (2.4), so we assume that these four new potentials satisfy the symmetry reduction condition q [ 1 ] ( x , t ) = p [ 1 ] ( − x , − t ) , s [ 1 ] ( x , t ) = r [ 1 ] ( − x , − t ) . Using this condition, one can get,

b 21 [ 1 ] ( x , t ) = − b 12 [ 1 ] ( − x , − t ) , b 23 [ 1 ] ( x , t ) = − b 14 [ 1 ] ( − x , − t ) , (3.8)

so the N-fold Darboux transformation for nonlocal couplings of AKNS equations (2.6) can be constructed,

Φ [ n ] = T n ( λ ) Φ = T [ n ] ( λ ) T [ n − 1 ] ( λ ) ⋯ T [ k ] ( λ ) ⋯ T [ 1 ] ( λ ) Φ , (3.9)

the gauge transformation T [ k ] ( λ ) has a similar form to (2.4). Different from the Darboux transformation of local equations, the Darboux transformation of nonlocal equations needs to meet the following constraints,

b 21 [ k ] ( x , t ) = − b 12 [ k ] ( − x , − t ) , b 23 [ k ] ( x , t ) = − b 14 [ k ] ( − x , − t ) , k = 1 , 2 , ⋯ , N (3.10)

The relationship between new and old potentials is as follows,

p [ k ] = p + 2 b 12 [ k ] ( − x , − t ) , r [ k ] = r + 2 b 14 [ k ] ( − x , − t ) . (3.11)

Proposition 1. Under the action of transformation T [ 1 ] (3.5), the matrix U [ 1 ] determined by Equation (3.2) has the same form as U, that is,

U [ 1 ] = ( − λ p [ 1 ] − λ r [ 1 ] q [ 1 ] λ s [ 1 ] λ 0 0 − λ p [ 1 ] 0 0 q [ 1 ] λ ) ,

the transformations between p , q , r , s and p [ 1 ] , q [ 1 ] , r [ 1 ] , s [ 1 ] are given by Equation (3.6).

Proof: For convenience, we will write T [ 1 ] as T, according to formulas T − 1 = ( det T ) − 1 T * , and

( T x + T U ) T * = ( f 11 ( λ ) f 12 ( λ ) f 13 ( λ ) f 14 ( λ ) f 21 ( λ ) f 22 ( λ ) f 23 ( λ ) f 24 ( λ ) 0 0 f 33 ( λ ) f 34 ( λ ) 0 0 f 43 ( λ ) f 44 ( λ ) ) (3.12)

It can be known by calculation that f i j ( λ ) are four or five times polynomials about λ , and it can be proved that λ j , j = 1 , 2 are the roots of f i j ( λ ) .

Because λ j ( j = 1 , 2 ) are the roots of det T , we can rewrite (3.12) as

( T x + T U ) T * = ( det T ) P ( λ ) , (3.13)

with

P ( λ ) = ( ∑ i = 0 5 p 11 ( i ) λ i ∑ i = 0 4 p 12 ( i ) λ i ∑ i = 0 5 p 13 ( i ) λ i ∑ i = 0 4 p 14 ( i ) λ i ∑ i = 0 4 p 21 ( i ) λ i ∑ i = 0 5 p 22 ( i ) λ i ∑ i = 0 4 p 23 ( i ) λ i ∑ i = 0 5 p 24 ( i ) λ i 0 0 ∑ i = 0 5 p 33 ( i ) λ i ∑ i = 0 4 p 34 ( i ) λ i 0 0 ∑ i = 0 4 p 43 ( i ) λ i ∑ i = 0 5 p 44 ( i ) λ i ) , (3.14)

where undetermined P i j ( k ) are not functions of λ . Because T − 1 = T * ⋅ ( det T ) − 1 , so Equation (3.13) can be rewritten as,

T x + T U = P ( λ ) T . (3.15)

Substitute Equation (3.5) and Equation (3.14) into Equation (3.15), by comparing the coefficients of the same power of λ on both sides of Equation (3.15), the following system of equations can be obtained,

λ 6 : p 11 5 = 0 , p 13 5 = 0 , p 22 5 = 0 , p 24 5 = 0 , p 33 5 = 0 , p 44 5 = 0 ,

λ 5 : p 11 4 = 0 , p 12 4 = 0 , p 13 4 = 0 , p 21 4 = 0 , p 22 4 = 0 , p 34 4 = 0 , p 43 4 = 0 , p 13 4 = 0 , p 44 4 = 0 , p 14 4 = 0 , p 23 4 = 0 , p 24 4 = 0 ,

λ 4 : p 11 3 = 0 , p 12 3 = 0 , p 21 3 = 0 , p 22 3 = 0 , p 33 3 = 0 , p 34 3 = 0 , p 44 3 = 0 , p 43 3 = 0 , p 13 3 = 0 , p 14 3 = 0 , p 23 3 = 0 , p 24 3 = 0 ,

λ 3 : p 11 2 = 0 , p 12 2 = 0 , p 21 2 = 0 , p 22 2 = 0 , p 33 2 = 0 , p 34 2 = 0 , p 43 2 = 0 , p 44 2 = 0 , p 13 2 = 0 , p 14 2 = 0 , p 23 2 = 0 , p 24 2 = 0 ,

λ 2 : p 12 1 = 0 , p 21 1 = 0 , p 34 1 = 0 , p 43 1 = 0 , p 14 1 = 0 , p 23 1 = 0 , p 11 1 = − 1 , p 33 1 = − 1 , p 13 1 = − 1 , p 22 1 = 1 , p 44 1 = 1 , p 24 1 = 1 ,

λ 1 : p 12 0 = p + 2 b 12 = p [ 1 ] , p 34 0 = p + 2 b 12 = p [ 1 ] , p 21 0 = q − 2 b 21 = q [ 1 ] , p 43 0 = q − 2 b 21 = q [ 1 ] , p 14 0 = r + 2 b 14 = r [ 1 ] , p 23 0 = s − 2 b 23 = s [ 1 ] , p 22 0 = 0 , p 11 0 = 0 , p 33 0 = 0 , p 44 0 = 0 , p 13 0 = 0 , p 24 0 = 0 ,

and the coefficients of λ 0 satisfies Equation (3.7). By substituting the above results into Equation (3.14), we can get,

P ( λ ) = ( − λ p [ 1 ] − λ r [ 1 ] q [ 1 ] λ s [ 1 ] λ 0 0 − λ p [ 1 ] 0 0 q [ 1 ] λ ) ,

then from Equation (3.3), we can get U [ 1 ] = P ( λ ) .

Using the same process, we can prove that V [ 1 ] has the same form as V under the transformation.

Proposition 2. Under the action of transformation T [ 1 ] (3.5), the matrix V [ 1 ] determined by Equation (3.2) has the same form as V, that is,

V [ 1 ] = ( A [ 1 ] B [ 1 ] E [ 1 ] F [ 1 ] C [ 1 ] − A [ 1 ] G [ 1 ] − E [ 1 ] 0 0 A [ 1 ] B [ 1 ] 0 0 C [ 1 ] − A [ 1 ] )

where A [ 1 ] , B [ 1 ] , C [ 1 ] , E [ 1 ] , F [ 1 ] , G [ 1 ] have the form (3.4), and the transformations between p , q , r , s and p [ 1 ] , q [ 1 ] , r [ 1 ] , s [ 1 ] are given by Equation (3.6).

Proof: According to formulas T − 1 = ( det T ) − 1 T * , and

( T t + T V ) T * = ( g 11 ( λ ) g 12 ( λ ) g 13 ( λ ) g 14 ( λ ) g 21 ( λ ) g 22 ( λ ) g 23 ( λ ) g 24 ( λ ) 0 0 g 33 ( λ ) g 34 ( λ ) 0 0 g 43 ( λ ) g 44 ( λ ) ) (3.16)

where g i j ( λ ) are six or seven times polynomials about λ , and λ j , j = 1 , 2 are the roots of g i j ( λ ) . So we can rewrite (3.16) as

( T t + T U ) T * = ( det T ) K ( λ ) , (3.17)

with

K ( λ ) = ( ∑ i = 0 7 k 11 ( i ) λ i ∑ i = 0 6 k 12 ( i ) λ i ∑ i = 0 7 k 13 ( i ) λ i ∑ i = 0 6 k 14 ( i ) λ i ∑ i = 0 6 k 21 ( i ) λ i ∑ i = 0 7 k 22 ( i ) λ i ∑ i = 0 6 k 23 ( i ) λ i ∑ i = 0 7 k 24 ( i ) λ i 0 0 ∑ i = 0 7 k 33 ( i ) λ i ∑ i = 0 6 k 34 ( i ) λ i 0 0 ∑ i = 0 6 k 43 ( i ) λ i ∑ i = 0 7 k 44 ( i ) λ i ) , (3.18)

where undetermined k i j ( k ) are not functions of λ . The Equation (3.17) can be rewritten as,

T t + T U = K ( λ ) T . (3.19)

Substitute Equation (3.5) and Equation (3.18) into Equation (3.19). By comparing the coefficients of the same power of λ on both sides of Equation (3.19) we can get Proposition 2. Because the results are very complex, we do not give them here.

In this section, we construct the Darboux transform of the high order nonlocal coupled AKNS system and establish the relationship between the new and old potentials. In the next section we will construct the exact solution of (2.6) using the results of this Section.

According to the classical Darboux transformation, firstly, four basic solutions of lax pairs (2.1) are given,

Δ = ( ϕ ( λ j ) ψ ( λ j ) Φ ( λ j ) Ψ ( λ j ) ) = ( ϕ 1 ψ 1 Φ 1 Ψ 1 ϕ 2 ψ 2 Φ 2 Ψ 2 ϕ 3 ψ 3 Φ 3 Ψ 3 ϕ 4 ψ 4 Φ 4 Ψ 4 ) . (4.1)

We know that by the Darboux transformation,

φ [ 1 ] = T [ 1 ] φ = ( T 11 T 12 T 13 T 14 T 21 T 22 T 23 T 24 T 31 T 23 T 33 T 34 T 41 T 42 T 43 T 44 ) ( ϕ 1 ψ 1 Φ 1 Ψ 1 ϕ 2 ψ 2 Φ 2 Ψ 2 ϕ 3 ψ 3 Φ 3 Ψ 3 ϕ 4 ψ 4 Φ 4 Ψ 4 ) , (4.2)

when λ = λ j , the 4 columns of φ [ 1 ] are linear correlation. In other words, there are different constants γ j ( 1 ) , γ j ( 2 ) , γ j ( 3 ) that satisfy the following relationship,

{ ∑ j = 1 4 T 1 j ϕ j + γ j ( 1 ) ∑ j = 1 4 T 1 j ψ j + γ j ( 2 ) ∑ j = 1 4 T 1 j Φ j + γ j ( 3 ) ∑ j = 1 4 T 1 j Ψ j = 0 ∑ j = 1 4 T 2 j ϕ j + γ j ( 1 ) ∑ j = 1 4 T 2 j ψ j + γ j ( 2 ) ∑ j = 1 4 T 2 j Φ j + γ j ( 3 ) ∑ j = 1 4 T 2 j Ψ j = 0 ∑ j = 1 4 T 3 j ϕ j + γ j ( 1 ) ∑ j = 1 4 T 3 j ψ j + γ j ( 2 ) ∑ j = 1 4 T 3 j Φ j + γ j ( 3 ) ∑ j = 1 4 T 3 j Ψ j = 0 ∑ j = 1 4 T 4 j ϕ j + γ j ( 1 ) ∑ j = 1 4 T 4 j ψ j + γ j ( 2 ) ∑ j = 1 4 T 4 j Φ j + γ j ( 3 ) ∑ j = 1 4 T 4 j Ψ j = 0 (4.3)

Substitute (3.5) into (4.3), we get,

b 11 [ 1 ] + α j ( 1 ) b 12 [ 1 ] + α j ( 2 ) b 13 [ 1 ] + α j ( 3 ) b 14 [ 1 ] = − λ j ( 1 + α j ( 2 ) ) b 21 [ 1 ] + α j ( 1 ) b 22 [ 1 ] + α j ( 2 ) b 23 [ 1 ] + α j ( 3 ) b 24 [ 1 ] = − λ j ( α j ( 1 ) + α j ( 3 ) ) α j ( 2 ) b 11 [ 1 ] + α j ( 3 ) b 12 [ 1 ] = − λ j α j ( 2 ) α j ( 2 ) b 21 [ 1 ] + α j ( 3 ) b 22 [ 1 ] = − λ j α j ( 3 ) (4.4)

where

α j ( 1 ) = γ j ( 1 ) ψ 2 + γ j ( 2 ) Φ 2 + γ j ( 3 ) Ψ 2 + ϕ 2 γ j ( 1 ) ψ 1 + γ j ( 2 ) Φ 1 + γ j ( 3 ) Ψ 1 + ϕ 1 , α j ( 2 ) = γ j ( 1 ) ψ 3 + γ j ( 2 ) Φ 3 + γ j ( 3 ) Ψ 3 + ϕ 3 γ j ( 1 ) ψ 1 + γ j ( 2 ) Φ 1 + γ j ( 3 ) Ψ 1 + ϕ 1 , α j ( 3 ) = γ j ( 1 ) ψ 4 + γ j ( 2 ) Φ 4 + γ j ( 3 ) Ψ 4 + ϕ 4 γ j ( 1 ) ψ 1 + γ j ( 2 ) Φ 1 + γ j ( 3 ) Ψ 1 + ϕ 1 , j = 1 , 2. (4.5)

Expression (4.4) can be written as a product of determinants,

A X T = b T , (4.6)

where

A = ( 1 α 1 ( 1 ) α 1 ( 2 ) α 1 ( 3 ) 0 0 0 0 0 0 0 0 1 α 1 ( 1 ) α 1 ( 2 ) α 1 ( 3 ) α 1 ( 2 ) α 1 ( 3 ) 0 0 0 0 0 0 0 0 0 0 α 1 ( 2 ) α 1 ( 3 ) 0 0 1 α 2 ( 1 ) α 2 ( 2 ) α 3 ( 3 ) 0 0 0 0 0 0 0 0 1 α 2 ( 1 ) α 2 ( 2 ) α 2 ( 3 ) α 2 ( 2 ) α 2 ( 3 ) 0 0 0 0 0 0 0 0 0 0 α 2 ( 2 ) α 3 ( 3 ) 0 0 ) ,

and

X = ( b 11 [ 1 ] b 12 [ 1 ] b 13 [ 1 ] b 14 [ 1 ] b 21 [ 1 ] b 22 [ 1 ] b 23 [ 1 ] b 24 [ 1 ] ) , b = ( − λ 1 ( 1 + α 1 ( 2 ) ) − λ 1 ( α 1 ( 1 ) + α 1 ( 3 ) ) − λ 1 α 1 ( 2 ) − λ 1 α 1 ( 3 ) − λ 2 ( 1 + α 2 ( 2 ) ) − λ 2 ( α 2 ( 1 ) + α 2 ( 3 ) ) − λ 2 α 2 ( 2 ) − λ 2 α 2 ( 3 ) ) .

To get b 1 i [ 1 ] , b 2 i [ 1 ] , i = 1 , 2 , 3 , 4 , we need to solve the system (4.6), but, we need to give the basic solutions for Lax pair firstly,

ϕ ( λ j ) = ( c j ( 1 ) λ j ( β λ j 2 t − x ) e α λ j 3 t − λ j x , 0 , c j ( 1 ) e α λ j 3 t − λ j x , 0 ) T ,

ψ ( λ j ) = ( 0 , − c j ( 2 ) λ j ( β λ j 2 t − x ) e λ j x − α λ j 3 t , 0 , c j ( 2 ) e λ j x − α λ j 3 t ) T , Φ ( λ j ) = ( c j ( 3 ) e α λ j 3 t − λ j x , 0 , 0 , 0 ) T , Ψ ( λ j ) = ( 0 , c j ( 4 ) e λ j x − α λ j 3 t , 0 , 0 , 0 ) T , (4.7)

with c j ( k ) ( j = 1 , 2 ; k = 1 , 2 , 3 , 4 ) are arbitrary constants.

By solving system of linear Equation (4.6) and using Equations (4.5), (4.7), we can get the expression of b 12 [ 1 ] ( x , t ) , b 21 [ 1 ] ( x , t ) and b 14 [ 1 ] ( x , t ) , b 23 [ 1 ] ( x , t ) since only these four variables are included in the Darboux transformation.

where

Δ 1 = ( c 1 ( 1 ) ) 2 ( λ 2 − λ 1 ) [ 2 γ 1 ( 2 ) c 1 ( 2 ) c 2 ( 2 ) ( β λ 1 3 t − λ 1 x ) + γ 1 ( 2 ) c 2 ( 2 ) ( γ 2 ( 2 ) c 3 ( 2 ) + c 1 ( 2 ) ) − γ 3 ( 2 ) c 1 ( 2 ) c 4 ( 2 ) ] , Δ 2 = − ( c 1 ( 2 ) ) 2 ( λ 2 − λ 1 ) [ 2 γ 1 ( 1 ) c 1 ( 1 ) c 2 ( 1 ) ( β λ 2 3 t − λ 2 x ) + γ 1 ( 1 ) c 2 ( 1 ) ( γ 2 ( 1 ) c 3 ( 1 ) + c 1 ( 1 ) ) − γ 3 ( 1 ) c 1 ( 1 ) c 4 ( 1 ) ] ,

Δ 3 = ( c 2 ( 1 ) γ 1 ( 1 ) ) 2 ( λ 2 − λ 1 ) [ 2 γ 1 ( 2 ) c 1 ( 2 ) c 2 ( 2 ) ( β λ 1 3 t − λ 1 x ) + γ 1 ( 2 ) c 2 ( 2 ) ( γ 2 ( 2 ) c 3 ( 2 ) − c 1 ( 2 ) ) − γ 3 ( 2 ) c 1 ( 2 ) c 4 ( 2 ) ] , Δ 4 = − ( c 1 ( 2 ) γ 1 ( 2 ) ) 2 ( λ 2 − λ 1 ) [ 2 γ 1 ( 1 ) c 1 ( 1 ) c 2 ( 1 ) ( β λ 2 3 t − λ 2 x ) + γ 1 ( 1 ) c 2 ( 1 ) ( γ 2 ( 1 ) c 3 ( 1 ) − c 1 ( 1 ) ) − γ 3 ( 1 ) c 1 ( 1 ) c 4 ( 1 ) ] .

Using the symmetry reduction condition (3.8), we can get

γ 1 ( 2 ) = i c 1 ( 2 ) c 2 ( 2 ) , γ 3 ( 2 ) = i c 3 ( 2 ) γ 2 ( 2 ) c 4 ( 2 ) , γ 1 ( 1 ) = − i c 1 ( 1 ) c 2 ( 1 ) , γ 3 ( 1 ) = − i c 3 ( 1 ) γ 2 ( 1 ) c 4 ( 1 ) , (4.9)

where i 2 = − 1 . In this case, we obtain the new solution,

p [ 1 ] = 2 b 12 [ 1 ] ( x , t ) , r [ 1 ] = 2 b 14 [ 1 ] ( x , t ) , (4.10)

with λ 1 , λ 2 , c 1 ( 1 ) , c 1 ( 2 ) , c 2 ( 1 ) , c 3 ( 1 ) , c 2 ( 2 ) , c 3 ( 2 ) , c 4 ( 1 ) , c 4 ( 2 ) are arbitrary constants, γ 1 ( 2 ) , γ 3 ( 2 ) , γ 1 ( 1 ) , γ 3 ( 1 ) conditions satisfied (4.9).

Through verification, we can know that solutions (4.10) are exact solutions of high order nonlocal coupled AKNS system (2.6). It can be seen from the expression (4.8) that the solutions of the system are in the form of exponential function, and the soliton solutions of the nonlocal system can be constructed by selecting appropriate parameters. In order to study the characteristics of the solutions, the graphs of module p and r are illustrated in

We know that if this solution is taken further as a new solution, we can perform the Darboux transformation again to produce another new solution. This process can be carried out continuously and a series of soliton solutions for high order nonlocal coupled AKNS system (2.6) can be theoretically obtained.

In this paper, considering the practical importance of nonlocal equations, a high-order nonlocal coupled AKNS system is constructed and studied by using the symmetric reduction method for the first time. With the help of lax pairs, the gauge transformation of the system is constructed, and the Darboux transformation of the system is obtained. In this paper, a new exact solution of the system is constructed by using the trivial seed solution and Darboux transformation, and the soliton solution of the system is obtained by selecting appropriate parameters. In order to study the properties of soliton solution, the corresponding images are given.

It is very meaningful to study the coupled nonlocal integrable model. However, much remains to be done. For example, the construction of Lie symmetry, symmetry reduction and similar solutions of nonlocal equations should be studied. In addition, the rogue wave solution of nonlocal equations is one of the research hotspots in recent years. Therefore, it is very meaningful to construct the rogue wave solution of nonlocal equations with variable coefficients by using the Darboux transform method. The above topics will be discussed in a series of future research works.

This work is supported by the National Natural Science Foundation of China (Nos. 11505090, 1171041), the Natural Science Foundation of Shaanxi Province under Grant No. 2018JQ-1045, and the Science and Technology Innovation Foundation of Xian under Grant No. 2017CGWL06. We also acknowledge Research Award Foundation for Outstanding Young Scientists of Shandong Province (No. BS2015SF009) and the doctoral foundation of Liaocheng University under Grant No. 318051413.

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

Xin, X.P., Guo, Z.X., Hu, Y.X. and Zhang, L.L. (2021) Darboux Transformation and Exact Solutions for High Order Nonlocal Coupled AKNS System. Journal of Applied Mathematics and Physics, 9, 2592-2608. https://doi.org/10.4236/jamp.2021.911167