Solutions of the nonlocal (2+1)-D breaking solitons hierarchy and the negative order AKNS hierarchy

Jing Wang,Hua Wu,Da-jun Zhang

Communications in Theoretical Physics ›› 2020, Vol. 72 ›› Issue (4) : 45002.

PDF(1592 KB)
Welcome to visit Communications in Theoretical Physics, June 15, 2025
PDF(1592 KB)
Communications in Theoretical Physics ›› 2020, Vol. 72 ›› Issue (4) : 45002. DOI: 10.1088/1572-9494/ab7705
Mathematical Physics

Solutions of the nonlocal (2+1)-D breaking solitons hierarchy and the negative order AKNS hierarchy

Author information +
History +

Abstract

The (2+1)-dimensional nonlocal breaking solitons AKNS hierarchy and the nonlocal negative order AKNS hierarchy are presented. Solutions in double Wronskian form of these two hierarchies are derived by means of a reduction technique from those of the unreduced hierarchies. The advantage of our method is that we start from the known solutions of the unreduced bilinear equations, and obtain solitons and multiple-pole solutions for the variety of classical and nonlocal reductions. Dynamical behaviors of some obtained solutions are illustrated. It is remarkable that for some real nonlocal equations, amplitudes of solutions are related to the independent variables that are reversed in the real nonlocal reductions.

Key words

nonlocal / (2+1)-dimensional breaking solitons AKNS hierarchy / negative order AKNS hierarchy / double Wronskian solutions / reduction

Cite this article

Download Citations
Jing Wang, Hua Wu, Da-jun Zhang. Solutions of the nonlocal (2+1)-D breaking solitons hierarchy and the negative order AKNS hierarchy[J]. Communications in Theoretical Physics, 2020, 72(4): 45002 https://doi.org/10.1088/1572-9494/ab7705

1. Introduction

Recently, PT-symmetric integrable nonlocal nonlinear Schrödinger equation
iqt(x,t)+qxx(x,t)±q2(x,t)q(x,t)=0,
(1)
was proposed by Ablowitz and Musslimani [1], which is reduced from the second-order Ablowitz–Kaup–Newell–Segur (AKNS) system with reduction r(x,t)=q(x,t). It is PT-symmetric because of the potential V(x,t)=±q(x,t)q(x,t)=V(x,t) (see [2, 3]). Since then, more and more integrable nonlocal systems, such as nonlocal Davey–Stewartson model, nonlocal Korteweg–de Vries model and so forth, were proposed (e.g. [49]). Classical solving methods, such as the inverse scattering transform and Darboux transformation (e.g. [1, 5, 1014]), have been successfully used to find solutions of nonlocal integrable systems. However, the bilinear method cannot be directly applied to the nonlocal case as it is difficult to define a ‘nonlocal’ Hirota’s bilinear operator. Recently, Chen et al [1518] developed a reduction technique of double Wronskians on unreduced bilinear equations to obtain solutions to those reduced equations, including nonlocal ones. One advantage of this method is that one only needs to make use of the solutions (most of them are known) of the unreduced bilinear equations. The other advantage is one can obtain N-soliton solutions for the whole hierarchy, rather than implementing reductions solution by solution and equation by equation (see [19]).
The (2+1)-dimensional breaking solitons equations ware first systematically investigated by Bogoyavlenskii [2026]. Such models can be used to describe two-dimensional interaction of a Riemann wave with transverse long waves (see [20, 21]). Mathematically, these equations can be constructed based on (1+1)-D Lax pairs by imposing an evolution in y-direction in both spectral parameters and eigenfunctions (see [21, 24]). The (2+1)-D breaking AKNS(BAKNS in brief) equations were first proposed in [24], and has been investigated from many aspects (e.g. [2731]). And recently, bilinear form and double Wronskian solutions of (2+1)-D BAKNS hierarchy were worked out [32]. In this paper, we will take into account of bilinearization-reduction technique and derive double Wronskian solutions for the classical and nonlocal (2+1)-D BAKNS hierarchy and the negative order AKNS hierarchy. Note that the negative order AKNS equations also exhibits interesting dynamical behaviors (see [33]).
This paper is organized as follows. In section 2 we recall the unreduced (2+1)-D BAKNS hierarchy, its bilinear form and double Wronskian solutions. In section 3 we present classical and nonlocal reductions of the (2+1)-D BAKNS hierarchy. In section 4, solutions in double Wronskian form of the reduced (2+1)-D BAKNS hierarchies are derived by means of the reduction technique from those of the unreduced hierarchy. Section 5 is contributed to the negative order AKNS hierarchy. Finally, section 6 consists of concluding remarks.

2. Unreduced (2+1)-D BAKNS hierarchy and solutions

Let us recall the unreduced (2+1)-D BAKNS hierarchy, its bilinear form and double Wronskian solutions.
The unreduced (2+1)-D BAKNS hierarchy reads [27]
utn=Kn=(K1,nK2,n)=Lnuy,n=1,2,3,,
(2)
where u=(q,r)T, L is a recursion operator
L=(x+2qx1r2qx1q2rx1rx2rx1q),
(3)
in which x=x and xx1=x1x=1. This hierarchy is related to the AKNS spectral problem [24]
ϕx=(λqrλ)ϕ,
(4)
and
ϕtn=(2λ)nϕy+(AnBnCnAn)ϕ,
where λtn=(2λ)nλy and
An=1(r,q)(Bn,Cn)T,(Bn,Cn)T=j=1n2njLj1(qy,ry)Tλnj.
The hierarchy (2) can be rewritten as
qtn=qx,tn1+2qx1(qr)tn1,
(5a)
rtn=rx,tn12rx1(qr)tn1
(5b)
with n = 1, 2, 3, ⋯ and setting t0 = y. Introducing transformation
q=gf,r=hf,
(6)
the hierarchy is written as the following bilinear form [32]
(Dtn+DxDtn1)gf=0,
(7a)
(DtnDxDtn1)hf=0,
(7b)
Dx2ff=2gh,
(7c)
for n ≥ 1, where Hirota’s bilinear operator D is defined as [34]
DxmDynf(x,y)g(x,y)=(xx)m(yy)nf(x,y)g(x,y)|x=x,y=y.

System (7) admits double Wronskian solutions [32]

f=Wn+1,m+1(φ,ψ)=|φ^(n);ψ^(m)|,
(8a)
g=2Wn+2,m(φ,ψ)=2|φ^(n+1);ψ^(m1)|,
(8b)
h=2Wn,m+2(φ,ψ)=2|φ^(n1);ψ^(m+1)|.
(8c)
Here, φ and ψ are respectively (n+m+2)th order column vectors
φ=(φ1,φ2,,φn+m+2)T,ψ=(ψ1,ψ2,,ψn+m+2)T
(9)
defined by
φ=exp(Ax+A2y+j=12jAj+2tj)C+,ψ=exp(AxA2yj=12jAj+2tj)C,
(10)
where A is a (n+m+2)×(n+m+2) constant matrix in C(m+n+2)×(m+n+2) and
C±=(c1±,c2±,,cn+m+2±)T,ci±C,
(11)
φ^(n) and ψ^(m) respectively denote (n+m+2)×(n+1) and (n+m+2)×(m+1) Wronski matrices
φ^(n)=(φ,xφ,x2φ,,xnφ),ψ^(m)=(ψ,xψ,x2ψ,,xmψ).
(12)

It is easy to find that for the odd members in the unreduced (2+1)-D BAKNS hierarchy (2)
ut2l+1=K2l+1,l=0,1,2,,
(13)
their solutions are given through (6) and (8) where
φ=exp(Ax+A2y+j=022j+1A2j+3t2j+1)C+,ψ=exp(AxA2yj=022j+1A2j+3t2j+1)C,
(14)
and for the even members
ut2l=K2l,l=1,2,,
(15)
their solutions are given through (6) and (8) where
φ=exp(Ax+A2y+j=122jA2j+2t2j)C+,ψ=exp(AxA2yj=122jA2j+2t2j)C.
(16)

3. Reductions of the (2+1)-D BAKNS hierarchy

3.1. Classical reductions

The odd hierarchy (13) allows a classical complex reduction
r(x,y,t)=δq(x,y,t),δ=±1,yiy,
(17)
where i is the imaginary unit and ∗ stands for complex conjugate. And the representative one is
qt1=iqxy2iδqx1(qq)y,δ=±1.
(18)
The even hierarchy (15) allows classical complex reduction
r(x,y,t)=δq(x,y,t)δ=±1,yiy,t2lit2l,
(19)
when l = 1, the corresponding equation is
qt2=qxxy4δqqqy2δqxx1(qq)y2δqx1(qxqyqxqy),δ=±1,
(20)

3.2. Nonlocal reductions

For the odd hierarchy (13), first, it allows a real nonlocal reduction
r(x,y,t)=δq(σx,y,σt),δ,σ=±1,x,y,tR,
(21)
where the simplest one-component equation reads
qt1=qxy+2δqx1(qq(σx,y,σt)y,δ,σ=±1.
(22)
Equation (13) also allows a complex nonlocal reduction
r(x,y,t)=δq(σx,y,σt),δ,σ=±1,x,y,tR,
(23)
and the resulted simplest one-component equation is
qt1=qxy+2δqx1(qq(σx,y,σt)y,δ,σ=±1.
(24)
The even hierarchy (15) also allows two types of reductions. One is real
r(x,y,t)=δq(σx,y,t),δ,σ=±1,
(25)
and the resulted simplest one-component equation is
qt2=qxxy4δq~qqy2δqxx1(qq~)y2δqx1(qxq~yq~xqy),δ,σ=±1,
(26)
where
q~=q(σx,y,t).
(27)
The other reduction is complex
r(x,y,t)=δq(σx,y,t),δ,σ=±1,
(28)
and the resulted simplest one-component equation is
qt2=qxxy4δq~qqy2δqxx1(qq~)y2δqx1(qxq~yq~xqy),δ,σ=±1,
(29)
where q~ is given as (27).
It is remarkable that the integration operator x1 should specially take the form
x1=12(xx+)dx,
(30)
which has played an important role in nonlocal reduction involved with x (see [16, 18]). For example, for the integration term x1(qr)y in (5) with n = 1, one can handle it as the following. First, introduce
z(x,y,t)=x1(qr)y.
Then we have
z(σx,y,σt)=12(σxσx+)[q(x,y,σt)r(x,y,σt)]ydx=σ2(xx+)[q(σx,y,σt)r(σx,y,σt)]ydx=σ2(xx+)[r(x,y,t)q(x,y,t)]ydx=σz(x,y,t),
which provides a clear expression for the integration term with nonlocal reduction.

4. Reduction of solutions

4.1. Nonlocal cases

In the following we implement the reduction procedure on the double Wronskian solutions, by which solutions for the reduced nonlocal hierarchies can be obtained from those of the unreduced (2+1)-D BAKNS hierarchy (2) presented in theorem 1.

4.1.1. Reduction (21)

The nonlocal hierarchy

qt2l+1=K1,2l+1|(21),l=0,1,2,
(31)
allows the following solution
q(x,y,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(32)
where φ=(φ1,φ2,,φ2n+2)T, ψ=(ψ1,ψ2,,ψ2n+2)T(i.e. m = n in (9)), defined by (14) and satisfy
ψ(x,y,t)=Tφ(σx,y,σt),
(33a)
C=TC+,
(33b)
in which T is a constant matrix determined through
AT+σTA=0,
(34a)
T2=σδI,σ,δ=±1.
(34b)

It follows from (14), (33b) and (34a) that

ψ(x,y,t)=exp(AxA2yj=022j+1A2j+3t2j+1)C=exp(σ(TAT1)xA2yj=022j+1(σTA2j+3T1)t2j+1)TC+=Texp(A(σx)+A2(y)+j=022j+1A2j+3(σt2j+1)C+=Tφ(σx,y,σt),
which gives rise to the the constraint (33a). Thus, f, g, h in (8) (with m = n) are expressed as
f(x,y,t)=|φ^(n);ψ^(n)|=|φ^(n)(x,y,t)[x];Tφ^(n)(σx,y,σt)[x]|,
(35a)
g(x,y,t)=2|φ^(n+1);ψ^(n1)|=2|φ^(n+1)(x,y,t)[x];Tφ^(n1)(σx,y,σt)[x]|,
(35b)
h(x,y,t)=2|φ^(n1);ψ^(n+1)|=2|φ^(n1)(x,y,t)[x];Tφ^(n+1)(σx,y,σt)[x]|,
(35c)
where we employed a notation (see [15])
φ^(n)(ax)[bx]=(φ(ax),bxφ(ax),bx2φ(ax),,bxnφ(ax).
Next, with assumption (34b), we have
f(σx,y,σt)=|φ^(n)(σx,y,σt)[σx];Tφ^(n)(x,y,t)[σx]|=|T||T1φ^(n)(σx,y,σt)[σx];φ^(n)(x,y,t)[σx]|=(1)(n+1)2|T||φ^(n)(x,y,t)[x];T1φ^(n)(σx,y,σt)[x]|=(1)(n+1)2|T||φ^(n)(x,y,t)[x];σδTφ^(n)(σx,y,σt)[x]|=(1)(n+1)2(σδ)n+1|T||φ^(n)(x,y,t)[x];Tφ^(n)(σx,y,σt)[x]|=(1)(n+1)2(σδ)n+1|T|f(x,y,t),
and in a similar way
g(σx,y,σt)=(1)(n+2)nσn+1δn|T|h(x,y,t),
which then give rise to
δq(σx,y,σt)=δg(σx,y,σt)f(σx,y,σt)=h(x,y,t)f(x,y,t)=r(x,y,t).
This is nothing but the reduction (21). Thus we complete the proof.

Solutions of (34) can be written out by assuming T and A to be block matrices
T=(T1T2T3T4),A=(K100K4),
(36)
where Ti and Ki are (n + 1) × (n + 1) matrices and Ki is a complex matrix. We list them in table 1.
One can prove that A and any of its similar forms lead to same solution (32). We only need to consider canonical form of A. When
Kn+1=Diag(k1,k2,,kn+1),
(37)
we get
φ=(c1eζ(k1),c2eζ(k2),,cn+1eζ(kn+1),d1eζ(k1),d2eζ(k2),,dn+1eζ(kn+1))T.
(38)
When Kn+1 is a (n+1)×(n+1) Jordan matrix Jn+1(k)
Jn+1(k)=(k001k001k)(n+1)×(n+1),
(39)
we get
φ=(ceζ(k),k1!(ceζ(k)),,knn!(ceζ(k)),deζ(k),k1!(deζ(k)),,knn!(deζ(k))T,
(40)
where
ζ(ki)=kix+ki2y+j=022j+1ki2j+3t2j+1.
(41)
As examples, for equation (22) with different (σ, δ), its one-soliton solution (1SS) are
qσ=1,δ=1=4cdke2k2yc2e2kx+4k3t+d2e2kx4k3t,
(42a)
qσ=1,δ=1=4cdke2k2yc2e2kx+4k3td2e2kx4k3t,
(42b)
qσ=1,δ=1=4ke2k2ye2kx+4k3te2kx4k3t,
(42c)
qσ=1,δ=1=4ke2k2yie2kx+4k3tie2kx4k3t,
(42d)
where we have taken k=k1,t=t1,C+=(c,d)T.
Next, let us quickly investigate dynamics of q(σ,δ)=q(1,1) which is governed by equation
qt1=qxy+2δqx1(qq(x,y,t)y.
(43)
Its 1SS (42a) is rewritten as
q(x,y,t)=2kcde2k2y|cd|sech(2kx+4k3t+ln|c||d|),
(44)
which is a moving wave with an initial phase ln|c||d| and a y-dependent amplitude 2cdke2k2y|cd|. The top trajectory is x(t)=2k2t+ln|cd|2k, and travel velocity is 2k2. Fixing y, the trajectory of one-solton in coordinate frame {x, t} is depicted as figure 1(a).
Figure 1. (a). Shape and motion of 1SS (42a) for equation (43), in which k = 1, c = 1, d = 1, y = 0. (b). Shape and motion of 2SS (45) for equation (43), in which k1=0.8,k2=0.6,c1=1, c2=1,d1=1,d2=2,y=0. (c). Shape and motion of 2SS (45) for equation (43), in which k1=1.2,k2=0.8,c1=1,c2=1,d1=1,d2=2,y=0. 1(d). Breather provided by (45) for equation (43), in which k1=0.8+0.6i,k2=0.80.6i,c1=1, c2=1,d1=1,d2=1,y=0.

Full size|PPT slide

The 2SS of equation (43) can be written as
q(x,y,t)=AB,
(45)
where
A=4(k12k22)e2(k12+k22)y(c2c12d2k2e2k2(2k22t+k2y+x)+8k13t+c2d12d2k2e4k23t+2k2x+4k1x+2k22y)4c1d1k1(k12k22)(c22e8k23t+d22e4k2x)e(2k1(2k12t+k1y+x)+2(k12+k22)y),
B=c22(e2k22(4k2t+y)+2k12y)×(c12(k1k2)2e8k13t+d12(k1+k2)2e4k1x)+4c1c2d1d2k1k2(e4k12y+e4k22y)e4k13t+4k23t+2k1x+2k2xd22e4k2x+2k12y+2k22y(c12(k1+k2)2e8k13t+d12(k1k2)2e4k1x),
and we note that 2SS is non-singular if c1c2d1d2k1k2<0. As shown in figures 1(b)–(d), there are three types of two-soliton interactions, respectively are soliton-soliton interaction, soliton-anti-soliton interaction and breather soliton. To investigate asymptotic behavior of the first two types solution interactions, we make asymptotic analysis for t going to infinity in the case of k1>k2>0. To do that, we first rewrite the 2SS (45) in the following coordinate
(X1=x2k12t,t),
(46)
fix X1, let t±, and we get
q(x,y,t){2c1d1k1e2k12y|c1||d1|sech(2k1X1+ln|d1|(k1+k2)|c1|(k1k2)),t+,2c1d1k1e2k12y|c1||d1|sech(2k1X1+ln|d1|(k1k2)|c1|(k1+k2)),t.
(47)
It shows that along the line X1 = constant, there is a soliton; when t±, the soliton asymptotically follows
toptrajectory:x(t)=2k12tln|d1c1|2k1±lnk1+k2k1k22k1,
(48a)
amplitude:2c1d1k1e2k12y|c1||d1|,
(48b)
and phase shift after interaction is lnk1+k2k1k2k1.
We can also rewrite the 2SS (45) in the coordinate (X2=x2k22t,t). By fixing X2, we can obtain
q(x,y,t){2c2d2k2e2k22y|c2||d2|sech(2k2X2+ln|d2|(k1+k2)|c2|(k1k2)),t+,2c2d2k2e2k22y|c2||d2|sech(2k2X2+ln|d2|(k1k2)|c2|(k1+k2)),t.
(49)
It indicates that along the line X2= constant, there is a soliton; when t±, the soliton asymptotically follows
toptrajectory:x(t)=2k22tln|d2c2|2k2±lnk1+k2k1k22k2,
(50a)
amplitude:2c2d2k2e2k22y|c2||d2|,
(50b)
and the phase shift after interaction is lnk1+k2k1k2k2.
Supposing k1>k2>0,c1c2d1d2k1k2<0, from the asymptotic analysis, we can obtain that if sgn[c1d1]=sgn[c2d2], it will appear soliton-soliton interaction, such as figure 1(b); if sgn[c1d1]=sgn[c2d2], it will appear soliton-anti-soliton interaction, such as figure 1(c).
We list out solutions for other cases of nonlocal reductions without giving proof.

4.1.2. Reduction (23)

For the reduction (23), the reduced nonlocal hierarchy

qt2l+1=K1,2l+1|(23),l=0,1,2,
(51)
allow a solution
q(x,y,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(52)
where φ and ψ are (2n+2)th order column vectors (i.e. m=n in (9)), defined by (14) and satisfy
ψ(x,y,t)=Tφ(σx,y,σt),
(53a)
C=TC+,
(53b)
in which T is a constant matrix determined through
AT+σTA=0,
(54a)
TT=σδI,σ,δ=±1.
(54b)

When T and A are block matrices (36), solutions to equations (54) are given in table 2.
Table 2. T and A for (54).
σ,δ T A
(1, −1) T1=T4=0n+1,T3=T2=In+1 K1=K4=Kn+1
(1, 1) T1=T4=0n+1,T3=T2=In+1 K1=K4=Kn+1
(−1, −1) T1=T4=0n+1,T3=T2=In+1 K1=K4=Kn+1
(−1, 1) T1=T4=0n+1,T3=T2=In+1 K1=K4=Kn+1
When Kn+1 is a diagonal matrix (37), the vector φ in (52) can be given as
φ=(c1eθ(k1),c2eθ(k2),,cn+1eθ(kn+1),d1eθ(σk1),d2eθ(σk2),,dn+1eθ(σkn+1))T,
(55)
where
θ(ki)=kix+ki2y+j=022j+1ki2j+3t2j+1.
(56)
When Kn+1=Jn+1(k) (39), we have
φ=(ceθ(k),k1!(ceθ(k)),,knn!(ceθ(k)),deθ(σk),k1!(deθ(σk)),,knn!(deθ(σk))T.
(57)
For the equation (24) with different (σ, δ), its 1SS are
qσ=1,δ=1=2cd(k+k)|c|2e2kx2k2y+4k3t+|d|2e2kx2k2y4k3t,
(58a)
qσ=1,δ=1=2cd(k+k)|c|2e2kx2k2y+4k3t|d|2e2kx2k2y4k3t,
(58b)
qσ=1,δ=1=2cd(kk)|c|2e2kx2k2y+4k3t|d|2e2kx2k2y4k3t,
(58c)
qσ=1,δ=1=2cd(kk)|c|2e2kx2k2y+4k3t+|d|2e2kx2k2y4k3t,
(58d)
where k=k1,t=t1,C+=(c,d)T.

4.1.3. Reduction (25)

For the reduction (25), the reduced nonlocal hierarchy

qt2l=K1,2l|(25),l=1,2,
(59)
allow a solution
q(x,y,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(60)
where φ=(φ1,φ2,,φ2n+2)T, ψ=(ψ1,ψ2,,ψ2n+2)T, defined by (16) and satisfy
ψ(x,y,t)=Tφ(σx,y,t),
(61a)
C=TC+,
(61b)
in which T is a constant matrix satisfying (34).

If T and A are block matrices (36), solutions to (34) haven been given by table 1. In the case that Kn+1 is a diagonal matrix (37), the vector φ in (60) can be given as
φ=(c1eξ(k1),c2eξ(k2),,cn+1eξ(kn+1),d1eξ(k1),d2eξ(k2),,dn+1eξ(kn+1))T.
(62)
where
ξ(ki)=kix+ki2y+j=122jki2j+2t2j.
(63)
When Kn+1=Jn+1(k) (39), we have
φ=(ceξ(k),k1!(ceξ(k)),,knn!(ceξ(k)),deξ(k),k1!(deξ(k)),,knn!(de(k))T.
(64)
Table 1. T and A for (34).
(σ, δ) T A
(1, −1) T1=T4=0n+1,T3=T2=In+1 K1=K4=Kn+1
(1, 1) T1=T4=0n+1,T3=T2=In+1 K1=K4=Kn+1
(−1, −1) T1=T4=In+1,T3=T2=0n+1 K1=K4=Kn+1
(−1, 1) T1=T4=iIn+1,T3=T2=0n+1 K1=K4=Kn+1
For the equation (26) with different (σ, δ), its 1SS are
qσ=1,δ=1=4cdke2k2y+8k4tc2e2kx+d2e2kx,
(65a)
qσ=1,δ=1=4cdke2k2y+8k4tc2e2kxd2e2kx,
(65b)
qσ=1,δ=1=4ke2k2y+8k4te2kxe2kx,
(65c)
qσ=1,δ=1=4ike2k2y+8k4te2kxe2kx,
(65d)
where k=k1,t=t2,C+=(c,d)T.
As an example we consider the dynamics of q(σ,δ)=q(1,1), which is governed by equation
qt2=qxxy+4q~qqy+2qxx1(qq~)y+2qx1(qxq~yq~xqy),
(66)
where q~=q(x,y,t). Its 1SS (65a) can be rewritten as
q(x,y,t)=2cdke2k2y+8k4t|cd|sech(2kx+ln|cd|).
(67)
Depicted as figure 2(a), this is a stationary wave with an initial phase ln|cd| and an amplitude that exponentially increases with time t and y, and the top trajectory is x(t)=ln|cd|2k. It is interesting that the amplitude is y and t-dependent, which looks like in nonisospectral case (see [35]) where amplitudes are changed due to time-dependent eigenvalues (spectral parameters). However, here the amplitude changes might be caused by the reversed y and t in the nonlocal reduction, rather than nonspectral parameters as the eigenvalues (see A in table 1) are still constant.
Figure 2. (a). Shape and motion of 1SS (65a) for equation (66), in which k=0.67,c=1,d=1,y=0. (b). Shape and motion of 2SS (68) for equation (66), in which k1=0.8,k2=0.6,c1=2,c2=1, d1=1,d2=1,y=0. (c). Shape and motion of e8t(k1k2)2q1,1(x,t) where q1,1(x,t) is depicted in (b).

Full size|PPT slide

The 2SS of equation (66) can be written as
q(x,y,t)=GF,
(68)
where
G=4c1d1k1(k12k22)(c22+d22e4k2x)e8k14t+8k24t+2k12y+2k22y+2k1(4k13t+k1y+x)+4(k12k22)e8k14t+8k24t+2k12y+2k22y(c2c12d2k2e2k2(4k23t+k2y+x)+c2d12d2k2e8k24t+2k2x+4k1x+2k22y),
F=d12(c22(k1+k2)2+d22(k1k2)2e4k2x)e(8k14t+8k24t+4k1x+2k12y+2k22y)+4c1c2d1d2k1k2e2(k1+k2)x(e4k12(4k12t+y)+e4k22(4k22t+y))c12(e8k14t+8k24t+2k12y+2k22y)×(c22(k1k2)2+d22(k1+k2)2e4k2x).

4.1.4. Reduction (28)

For the reduction (28), the reduced nonlocal hierarchy

qt2l=K1,2l|(28),l=1,2,
(69)
allow a solution
q(x,y,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(70)
where φ and ψ are defined by (9) with m=n and (16), satisfying
ψ(x,y,t)=Tφ(σx,y,t),
(71a)
C=TC+,
(71b)
in which T is a constant matrix satisfied (54).

When T and A are block matrices (36), solutions to (54) are already given by table 2. In the case that Kn+1 is a diagonal matrix (37), the vector φ in (5) can be written as
φ=(c1eη(k1),c2eη(k2),,cn+1eη(kn+1),d1eη(σk1),d2eη(σk2),,dn+1eη(σkn+1))T,
(72)
where
η(ki)=kix+ki2y+j=122jki2j+2t2j.
(73)
When Kn+1=Jn+1(k) (39), we get
φ=(ceη(k),k1!(ceη(k)),,knn!(ceη(k)),deη(σk),k1!(deη(σk)),,knn!(deη(σk))T.
(74)
For equation (29) with different (σ, δ), its 1SS are
qσ=1,δ=1=2cd(k+k)|c|2e2kx2k2y8k4t+|d|2e2kx2k2y8k4t,
(75a)
qσ=1,δ=1=2cd(k+k)|c|2e2kx2k2y8k4t|d|2e2kx2k2y8k4t,
(75b)
qσ=1,δ=1=2cd(kk)|c|2e2kx2k2y8k4t|d|2e2kx2k2y8k4t,
(75c)
qσ=1,δ=1=2cd(kk)|c|2e2kx2k2y8k4t+|d|2e2kx2k2y8k4t,
(75d)
where k=k1,t=t2,C+=(c,d)T.

4.2. Classical cases

For the classical reduction (17), the reduced local hierarchy

qt2l+1=K1,2l+1|(17),l=0,1,2,
(76)
allow a solution
q(x,y,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(77)
where φ and ψ are defined by (9) with m=n and (14), satisfying
ψ(x,y,t)=Tφ(x,y,t),
(78a)
C=TC+,
(78b)
in which T is a constant matrix determined through
AT+TA=0,
(79a)
TT=δI,δ=±1.
(79b)

When T and A are block matrices (36), solutions to equations (79) are given in table 2 with σ = 1. When Kn+1 is a diagonal matrix (37), the vector φ in (77) can be given as
φ=(c1eθ(k1),c2eθ(k2),,cn+1eθ(kn+1),d1eθ(k1),d2eθ(k2),,dn+1eθ(kn+1))T,
(80)
where
θ(kj)=kjx+ikj2y+l=022l+1kj2l+3t2l+1.
(81)
When Kn+1=Jn+1(k) (39), we have
φ=(ceθ(k),k1!(ceθ(k)),,knn!(ceθ(k)),deθ(k),k1!(deθ(k)),,knn!(deθ(k))T.
(82)
In addition, we present the 1SS of the equation (18)
qδ=1=2cd(k+k)|c|2e2kx2ik2y+4k3t|d|2e2kx2ik2y4k3t,
(83a)
qδ=1=2cd(k+k)|c|2e2kx2ik2y+4k3t+|d|2e2kx2ik2y4k3t,
(83b)
where k=k1,t=t1,C+=(c,d)T.

For the reduction (19), the reduced local hierarchy

qt2l=iK1,2l|(19),l=1,2,
(84)
allow a solution
q(x,y,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(85)
where φ=(φ1,φ2,,φ2n+2)T, ψ=(ψ1,ψ2,,ψ2n+2)T expressed as (16) and satisfy the condition (78), in which T is a constant matrix satisfied (79).

When T and A are block matrices (36), solutions to (79) are already given by table 2. In the case that Kn+1 is a diagonal matrix (37), the vector φ in (85) can be written as
φ=(c1eη(k1),c2eη(k2),,cn+1eη(kn+1),d1eη(k1),d2eη(k2),,dn+1eη(kn+1))T,
(86)
where
η(kj)=kjx+ikj2y+l=122likj2l+2t2l.
(87)
When Kn+1=Jn+1(k) (39), we get
φ=(ceη(k),k1!(ceη(k)),,knn!(ceη(k)),deη(k),k1!(deη(k)),,knn!(deη(k))T.
(88)
For equation (20), its 1SS with different δ are given as
qδ=1=2cd(k+k)|c|2e2kx2ik2y8ik4t|d|2e2kx2ik2y8ik4t,
(89a)
qδ=1=2cd(k+k)|c|2e2kx2ik2y8ik4t+|d|2e2kx2ik2y8ik4t,
(89b)
where k=k1,t=t2,C+=(c,d)T.

5. Negative order AKNS hierarchy

In this section, we list double Wronskian solutions for the negative order AKNS hierarchy and consider several cases of reductions.

5.1. Solutions

For the negative order AKNS hierarchy [33]
(qtnrtn)=Kn=(K1,nK2,n)=Ln(qr),n=1,2,3,,
(90)
where L is defined by (3). The hierarchy can be rewritten as [36]
q=qx,t12qx1(qr)t1,
(91a)
r=rx,t12rx1(qr)t1,
(91b)
qtn1=qx,tn+2qx1(qr)tn,
(91c)
rtn1=rx,tn2rx1(qr)tn,(n=2,3,).
(91d)
By the transformation
q=gf,r=hf,
(92)
its bilinear form turns out to be [36]
Dx2ff=2gh,
(93a)
DxDt1gf=gf,
(93b)
DxDt1hf=hf,
(93c)
(Dtn1+DxDtn)gf=0,
(93d)
(Dtn1DxDtn)hf=0,(n=2,3),
(93e)
which allow double Wronskian solution [32]
f=|φ^(n);ψ^(m)|,g=2|φ^(n+1);ψ^(m1)|,h=2|φ^(n1);ψ^(m+1)|.
(94)
where φ and ψ are defined as
φ=exp(Axn=12(n+1)Antn)C+,ψ=exp(Ax+n=12(n+1)Antn)C,
(95)
AC(n+m+2)×(n+m+2) and C± are (n+m+2)th-order constant column vector.

5.2. Reductions and solutions

The odd members in the hierarchy (91) admit a real reduction
r(x,t)=δq(σx,σt),σ,δ=±1,
(96)
under which we have
qt2l+1=K1,2l+1|(96),l=0,1,2,
(97)
and the representative one
q=qx,t12δqx1(qq(σx,σt)t1.
(98)
In addition, the odd members in the hierarchy (91) also admit a complex reduction
r(x,t)=δq(σx,σt),σ,δ=±1,
(99)
under which we have
qt2l+1=K1,2l+1|(99),l=0,1,2,
(100)
and
q=qx,t12δqx1(qq(σx,σt)t1
(101)
as the first member.
For those even members in the hierarchy (91), first, they allow a real reduction
r(x,t)=δq(σx,t),σ,δ=±1,
(102)
which generate
qt2l=K1,2l|(102),l=1,2,,
(103)
the simplest one is
q=qxxt2+4δq^qqt2+2δqxx1(qq^)t2+2δqx1(qxq^t2q^xqt2),δ,σ=±1,
(104)
where q^=q(σx,t). And second, a complex reduction
r(x,t)=δq(σx,t),σ,δ=±1,t2jit2j,
(105)
leads to
 iqt2l=K1,2l|(105),l=1,2,,
(106)
in which the representative one is
 iq=qxxt2+4δq^qqt2+2δqxx1(qq^)t2+2δqx1(qxq^t2q^xqt2),δ,σ=±1,
(107)
where q^=q(σx,t).
In the following we list out double Wronskian solutions and we skip proof.

The system (97) admits the following solutions

q(x,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(108)
where φ=(φ1,φ2,,φ2n+2)T, ψ=(ψ1,ψ2,,ψ2n+2)T, are defined by
φ=exp(Axj=02(2j+2)A(2j+1)t2j+1)C+,ψ=exp(Ax+j=02(2j+2)A(2j+1)t2j+1)C
(109)
and satisfy
ψ(x,t)=Tφ(σx,σt),
(110a)
C=TC+,
(110b)
in which T is a constant matrix determined by (34).

For equation (98) with different (σ, δ), its 1SS are
qσ=1,δ=1=4cdkc2e2kxt2k+d2e2kx+t2k,
(111a)
qσ=1,δ=1=4cdkd2e2kx+t2kc2e2kxt2k,
(111b)
qσ=1,δ=1=4ke2kxt2k+e2kx+t2k,
(111c)
qσ=1,δ=1=4kie2kxt2k+e2kx+t2k,
(111d)
where k=k1,t=t1,C+=(c,d).

The system (100) admits solutions

q(x,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(112)
where φ=(φ1,φ2,,φ2n+2)T, ψ=(ψ1,ψ2,,ψ2n+2)T, are defined by (109) and satisfy
ψ(x,t)=Tφ(σx,σt),
(113a)
C=TC+,
(113b)
in which T is a constant matrix determined by (54).

1SS of the equation (101) with different (σ, δ) are
qσ=1,δ=1=2cd(k+k)|c|2e2kxt2k+|d|2e2kx+t2k,
(114a)
qσ=1,δ=1=2cd(k+k)|d|2e2kx+t2k|c|2e2kxt2k,
(114b)
qσ=1,δ=1=2cd(kk)|c|2e2kx+t2k|d|2e2kx+t2k,
(114c)
qσ=1,δ=1=2cd(kk)|c|2e2kx+t2k+|d|2e2kx+t2k,
(114d)
where k=k1,t=t1,C+=(c,d)T.

The hierarchy (103) admit

q(x,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(115)
where φ=(φ1,φ2,,φ2n+2)T, ψ=(ψ1,ψ2,,ψ2n+2)T are defined by
φ=exp(Axj=12(2j+1)A2jt2j)C+,ψ=exp(Ax+j=12(2j+1)A2jt2j)C
(116)
and satisfy
ψ(x,t)=Tφ(σx,t),
(117a)
C=TC+,
(117b)
in which T is a constant matrix determined by (34).

For equation (104), its 1SS with different (σ, δ) are
qσ=1,δ=1=4cdkc2e2kx+t4k2+d2e2kx+t4k2,
(118a)
qσ=1,δ=1=4cdkd2e2kx+t4k2c2e2kx+t4k2,
(118b)
qσ=1,δ=1=4ke2kx+t4k2+e2kx+t4k2,
(118c)
qσ=1,δ=1=4ike2kx+t4k2+e2kx+t4k2,
(118d)
where k=k1,t=t2,C+=(c,d)T.
As an example, let us illustrate dynamics of q(σ,δ)=q(1,1) which is governed by the equation
q=qxxt24q^qqt22qxx1(qq^)t22qx1(qxq^t2q^xqt2),
(119)
where q^=q(x,t), (118a) describes a single soliton
q(x,t)=2cdk|cd|et4k2sech(2kx+ln|cd|),
(120)
with an initial phase ln|cd|, top trajectory x(t)=ln|cd|2k, and an amplitude that decreases with time, see figure 3(a).
Figure 3. (a). Shape and motion of 1SS (118a) for equation (119), in which k=0.9,c=5,d=1. (b). Shape and motion of 2SS (121) for equation (119), in which k1=3,k2=1,c1=0.5, c2=1,d1=1,d2=1.

Full size|PPT slide

The 2SS of equation (119) is
q(x,t)=GF,
(121)
where
G=4(k12k22)(c1d1k1e2k1x(c22+d22et4k22+4k2x)+c2d12d2k2et4k12+4k1x+2k2x)4(k12k22)c2c12d2k2e2k2x,
F=c12(c22(k1k2)2+d22(k1+k2)2et4k22+4k2x)4c2c1d1d2k1k2(et4k12+et4k22)e2(k1+k2)x+d12et4k12+4k1x(c22(k1+k2)2+d22(k1k2)2et4k22+4k2x),
which we depict in figure 3(b).

The hierarchy (106) admit solutions

q(x,t)=2|φ^(n+1);ψ^(n1)||φ^(n);ψ^(n)|,
(122)
where φ and ψ are (2n+2)th order column vectors (i.e. m=n in (9)), are defined by (116) and satisfy
ψ(x,t)=Tφ(σx,t),
(123a)
C=TC+,
(123b)
in which T is a constant matrix determined by (54).

The equation (107) with different (σ, δ) has the following 1SS:
qσ=1,δ=1=2cd(k+k)|c|2e2kx+it4k2+|d|2e2kx+it4k2,
(124a)
qσ=1,δ=1=2cd(k+k)|d|2e2kx+it4k2|c|2e2kx+it4k2,
(124b)
qσ=1,δ=1=2cd(kk)|c|2e2kx+it4k2|d|2e2kx+it4k2,
(124c)
qσ=1,δ=1=2cd(kk)|c|2e2kx+it4k2+|d|2e2kx+it4k2,
(124d)
where k=k1,t=t2,C+=(c,d)T.

6. Conclusions

We have discussed possible nonlocal reductions of the (2+1)-D BAKNS hierarchy and the negative order AKNS hierarchy. By means of the reduction technique (see [1518]), we derived N-soliton solutions and multiple-pole solutions in double Wronskian form of the reduced hierarchies from those of the unreduced hierarchies. Note that these hierarchies contain integration w.r.t. x, the integration operator with form of (30) is necessary and helpful in considering nonlocal reductions. In addition, this reduction method is based on bilinear forms and double Wronskians, As we can see, it allows us to obtain N-soliton solutions and multiple-pole solutions for the whole reduced hierarchy. Compared with [14, 16, 19] where only some single equations in the (2+1)-D BAKNS hierarchy and the negative order AKNS hierarchy were considered, here we solved the two hierarchies and obtained more solutions. Finally, it is remarkable that for some real nonlocal equations, amplitudes of solutions are related to the independent variables that are reversed in the real nonlocal reductions, e.g. (44), (67) and (120).

References

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

Acknowledgments

This work was supported by the NSF of China [grant numbers 11 875 040, 11 631 007, 11 571 225].

RIGHTS & PERMISSIONS

© 2020 Chinese Physical Society and IOP Publishing Ltd
PDF(1592 KB)

388

Accesses

0

Citation

Detail

Sections
Recommended

/