The D'Alembert type waves and the soliton molecules in a (2+1)-dimensional Kadomtsev-Petviashvili with its hierarchy equation*

Hui-Ling Wu(吴慧伶),Sheng-Wan Fan(樊盛婉),Jin-Xi Fei(费金喜),Zheng-Yi Ma(马正义)

Communications in Theoretical Physics ›› 2021, Vol. 73 ›› Issue (10) : 105002.

PDF(621 KB)
Welcome to visit Communications in Theoretical Physics, May 6, 2025
PDF(621 KB)
Communications in Theoretical Physics ›› 2021, Vol. 73 ›› Issue (10) : 105002. DOI: 10.1088/1572-9494/ac11ef
Mathematical Physics

The D'Alembert type waves and the soliton molecules in a (2+1)-dimensional Kadomtsev-Petviashvili with its hierarchy equation*

Author information +
History +

Abstract

For a one (2+1)-dimensional combined Kadomtsev-Petviashvili with its hierarchy equation, the missing D'Alembert type solution is derived first through the traveling wave transformation which contains several special kink molecule structures. Further, after introducing the Bäcklund transformation and an auxiliary variable, the N-soliton solution which contains some soliton molecules for this equation, is presented through its Hirota bilinear form. The concrete molecules including line solitons, breathers and a lump as well as several interactions of their hybrid are shown with the aid of special conditions and parameters. All these dynamical features are demonstrated through the different figures.

Key words

Kadomtsev-Petviashvili equation / soliton molecule / breather/lump soliton / elastic interaction

Cite this article

Download Citations
Hui-Ling Wu(吴慧伶), Sheng-Wan Fan(樊盛婉), Jin-Xi Fei(费金喜), et al. The D'Alembert type waves and the soliton molecules in a (2+1)-dimensional Kadomtsev-Petviashvili with its hierarchy equation*[J]. Communications in Theoretical Physics, 2021, 73(10): 105002 https://doi.org/10.1088/1572-9494/ac11ef

1. Introduction

It is well known that the standard Kadomtsev-Petviashvili (KP) equation can be written as (vv(x, y, t))
vt3vvx14vxxx34vyydx=0,
(1)
which possesses the quadratic nonlinearity vvx and a weak dispersion vxxx [1]. Equation (1) admits the weakly dispersive waves in a paraxial wave approximation and can also be described as the evolution of long ion-acoustic waves of small amplitude propagating in plasma physics under the effect of long transverse perturbations [2, 3].
The KP hierarchy is a paradigm of the hierarchies of one integrable system and includes an infinite number of integrable nonlinear differential equations [46]. Generally speaking, integrable equations possess properties such as Lax pair, an infinite number of conservation laws and multiple soliton solutions. These nonlinear wave equations that model complex physical phenomena have attracted more and more attention from different fields [713]. The following member of the KP hierarchy [14]
vt12vxxy12vyyydxdx  2vxvydx4vvy=0,
(2)
is derived through the pseudo-differential formalism [15]. Its multiple soliton solutions are obtained mainly using a simplified version of Hirota's method [16].
In fact, an integrable nonlinear equation combined with another will often produce unexpected results [1720]. Inspired by this idea, we employ the following integrable system which is extended and combines the (2+1)-dimensional KP equation (1) with the KP hierarchy (2) (KP-h)
vt+a[12αvvx+αvxxx+12βvyydx] + b[2αvxxy+8βvyyydxdx + 8αvxvydx+16αvvy]=0,
(3)
to illustrate the related soliton structures and some nonlinear phenomena. Here, the parameters a, b, α and β are four real constants. When α=14,β=116 and a = 1, b = 0, equation (3) is just the KP equation (1), but a = 0, b = 1, equation (3) is the KP hierarchy (2).
This system (3) is integrable and its Lax pair can be derived as (ψψ(x, y, t), ΨΨ(x, y, t))
{ψy=12αβ(ψxx+2vψ),ψt=2ψ[α(3avx+2bvy)+2αβ(3aβvydx+2bβvyydxdx+2αbv2+αbvxx)]4αψx(3av+2bvydx+2bαβvx)2αbαβψxxxx4aαψxxx8αbαβvψxx,
(4)
and its dual
{Ψy=12αβ(Ψxx+2vΨ),Ψt=2Ψ[α(3avx+2bvy)2αβ(3aβvydx+2bβvyydxdx+2αbv2+αbvxx)]4αΨx(3av+2bvydx2bαβvx)+2αbαβΨxxxx4aαΨxxx+8αbαβvΨxx.
(5)
The outline of this paper is as follows: In section 2, after taking the traveling wave transformation, the missing D'Alembert type solution is derived which contains several special kink molecule structures for the (2+1)-dimensional Kadomtsev-Petviashvili with its hierarchy equation (3). In section 3, after introducing the Bäcklund transformation and an auxiliary variable τ, the Hirota bilinear form of this equation is conducted, the N-soliton solution which contains some soliton molecules is presented. These concrete molecules include the line solitons, the breathers and a lump as well as several interactions of their hybrid after studying some special conditions and parameters. These dynamic features are demonstrated using the different figures. In the last section, a brief summary is given for this paper.

2. The D'Alembert type solutions of the (2+1)-dimensional KP-h system

Recently, Lou studied the cKP3-4 equation and the Nizhnik-Novikov-Veselov (NNV) equation to investigate various types of solitons including the missing D'Alembert type solutions, the arbitrary traveling waves moving in one direction with a fixed model dependent velocity and soliton molecules [21, 22]. In fact, this (2+1)-dimensional KP-h system (3) may have the same properties. For this purpose, we rewrite equation (3) as (uu(x, y, t), ww(x, y, t))
{vt+a[12αvvx+αvxxx+12βwx]+2b[αuxxx+4βwy  + 4αuvx+8αuxv]=0,vy=ux,vyy=wxx.
(6)
After taking the traveling wave transformation
u=U(X), v=V(X), w=W(X),X=kx+py+qt,
(7)
Equation (6) becomes
{qVX+a[12αkVVX+αk3VXXX+12βkWX]  + 2b[αk3UXXX+4βpWX+4αkUVX+8αkUXV]=0,pVX=kUX,p2VXX=k2WXX,
(8)
which can be derived as
{(q+12aβp2k+8bβp3k2)VX+2α(ak+12bp)VVX  + αk2(ak+2bp)VXXX=0,U=pVk,W=p2Vk2.
(9)
This equation induces
p=ak2b, q=2a3βkb2,
(10)
and
u=aV(X)2b, v=V(X), w=a2V(X)4b2,X=xa2by2a3βb2t,
(11)
where the traveling wave V becomes an arbitrary D'Alembert type wave. This solution contains several special Kink molecule structures through equations (9)–(11). For example, when the related parameters are
a=1, b=α=1, β=2, k=k1=a12=1,k2=45, η1=0, η2=10,
(12)
while the solution v = V(X) of equation (3) is
V(X)=[ln(1+ek1X+η1+ek2X+η2+a12ek1X+k2X+η1+η2)]X,
(13)
V(X)=[ln(1+ek1X+η1+ek2X+η2+a12(1+15cos(X))ek1X+k2X+η1+η2)]X,
(14)
V(X)=[ln(1+ek1X+η1+(1+15sin(X))ek2X+η2+a12ek1X+k2X+η1+η2)]X,
(15)
V(X)=[ln(1+ek1X+η1+(1+15sin(X))ek2X+η2+a12(1+15cos(X))ek1X+k2X+η1+η2)]X,
(16)
we can present a kink molecule (figure 1(a)), a kink and a half periodic kink (HPK) molecule (figure 1(b)), a HPK-HPK molecule (figure 1(c)) and a HPK-PK molecule (figure 1(d)), respectively at time t = 0 [2022].
Figure 1. (a) Kink molecule expressed by equations (12) and (13); (b) Kink-HPK molecule given by equations (12) and (14); (c) HPK-HPK molecule described by equations (12) and (15); (d) HPK-PK molecule defined by equations (12) and (16) at time t = 0.

Full size|PPT slide

3. The soliton molecule solutions of the (2+1)-dimensional KP-h system

To derive the soliton molecule solutions of the (2+1)-dimensional KP-h system (3), we first introduce the following Bäcklund transformation (ff(x, y, t))
v=uxx, u=lnf.
(17)
Then, the trilinear form in terms of the auxiliary function f can be written out through equation (3)
a(2αffxxfxxx24βffyfxy12βffxfyy+24βfxfy2 + αf2fxxxxx5αffxfxxxx+8αfx2fxxx6αfxfxx2 +12βf2fxyy)+2b(4αfxfyfxxx4αfxfxxfxy+4αfx2fxxy  12βffyfyy2αfyfxx2+2αffxxfxxy+αf2fxxxxy  4αffxfxxxyαffyfxxxx+4βf2fyyy+8βfyyy3) + f2fxxt2ffxfxt+2ftfx2fftfxx=0.
(18)
Further, we introduce an auxiliary variable τ, such that this trilinear equation (18) can be bilinearized directly
[DxDτ+12βDy2+αDx4]ff=0.
(19)
After that, equation (18) becomes
(fx2fx)[4αbDx3Dy+3DxDt3aDxDτ   2bDyDτ]ff=0,
(20)
which can be solved through the bilinear equation
[4αbDx3Dy+3DxDt3aDxDτ2bDyDτ]ff=0,
(21)
with the D-operator defined by [23]
DxkDylDtmDτn(fg)=(xx)k(yy)l × (tt)m(ττ)n × f(x,y,t,τ)g(x,y,t,τ)|x=x,y=y,t=t,τ=τ.
(22)
Therefore, the N-soliton solutions described through equation (17) of the (2 + 1)-dimensional KP-h system (3) also have the following standard expression for the auxiliary function f
f=τ=0,1ei<jNτiτjAij+i=1Nτiξi,
(23)
from equation (21) [2428], while the relations of the parameters are
ξi=kix+piy+qit+ηi,qi=aαki5+12aβkipi2+8bβpi3+2bαpiki4ki2   (i=1,2,,N),
(24)
eAij=aij=αki2kj2(kikj)24β(kipjkjpi)2αki2kj2(ki+kj)24β(kipjkjpi)2(i<j, i,j=1,2,,N).
(25)
These multiple soliton solutions of equations (17) with (23) may contain many kinds of resonant excitations [29, 30]. In the following, we try our best to derive the soliton molecules, the breathers through the resonant condition, the rational lumps by further limit condition and their hybrid structures for the above conclusion. The concrete formula of the auxiliary function f of equation (23) for N = 4 can be rewritten as
f(x,y,t)=1+eξ1+eξ2+eξ3+eξ4+a12eξ1+ξ2 + a13eξ1+ξ3+a14eξ1+ξ4+a23eξ2+ξ3+a24eξ2+ξ4 + a34eξ3+ξ4+a12a13a23eξ1+ξ2+ξ3+a23a24a34eξ2+ξ3+ξ4 + a13a14a34eξ1+ξ3+ξ4+a12a14a24eξ1+ξ2+ξ4 + a12a13a14a23a24a34eξ1+ξ2+ξ3+ξ4,
(26)
where the phase shifts aij (i < j, i, j = 1, 2, 3, 4) and the parameters qi (i = 1, 2, 3, 4) satisfy equations (24) and (25).
In order to construct soliton molecules, we use the following velocity resonant conditions:
kikj=pipj=qiqj, ki|kj|, pi|pj|, qi|qj|, (ij).
(27)
(1) The molecules produced by the solitons

One two-soliton molecule

We first take the above phase shifts aij = 0(j = 3, 4) and ξi = 0(i = 3, 4) in equation (26), and require the velocity resonant condition (27) for i = 1, j = 2. The soliton molecule structure is satisfied with the condition
k1p1=k2p2=2ba,
(28)
from equations (24) and (25) for a12, q1 and q2.
Figure 2(a) exhibits the two-soliton molecule after taking the determining parameters
a=1, b=α=1, β=2, k1=32,k2=2, η1=0, η2=15,
(29)
which lead to the other parameters' values
p1=34, p2=1, q1=6, q2=8, a12=149,
(30)
In this case, the the auxiliary function f is given by
f=1+e32x+34y6t+e2x+y8t+15+149e72x+74y14t+15.
(31)
Figure 2. (a) The two-soliton molecule structure of the solution v from equations (17) and (31) at t = 0. (b) The three-soliton molecule structure of the solution v from equations (17) and (33) at t = 0. (c) The four-soliton molecule structure of the solution v from equations (17) and (35) at t = 0.

Full size|PPT slide

This soliton molecule structure of variable v of equation (17) for the (2+1)-dimensional KP-h system (3) is constructed by two line solitons in the molecule, which possess the same velocity, but their height and width are different (figure 2(a)).

One three-soliton molecule

We continue to take the phase shift ai4 = 0(i < 4) and ξ4 = 0 in equation (26), and require the velocity resonant condition (27) for i, j = 1, 2, 3 the three-soliton molecule structure may appear. For example, when we set the parameters in equation (29) and k3=52,η3=30, while the others are derived as
p3=54, q3=10,   a13=116, a23=181,
(32)
from equations (24), (25) and (27), then the transformation function f of equation (26) is given by
f=1+e32x+34y6t+e2x+y8t+15+e52x+54y10t+30+149e72x+74y14t+15+116e4x+2y16t+30+181e92x+94y18t+45+163504e6x+3y24t+45.
(33)
The three-soliton molecule structure of v for the (2+1)-dimensional KP-h system (3) is constructed directly from equations (17) and (33), where three line solitons in the molecule (figure 2(b)).

One four-soliton molecule

When taking the phase shifts aij ≠ 0(i < j, i, j = 1, 2, 3, 4) and the parameters qi ≠ 0(i = 1, 2, 3, 4) in equation (26), and expanding the resonant conditions (27) for i, j = 1,2,3,4, the four-soliton molecule can be produced.
For instance, we select the same parameters as equation (29) and
k3=52, k4=3, η3=30, η4=45,
(34)
then the function f of equation (26) is derived as
f=1+e32x+34y6t+e2x+y8t+15+e52x+54y10t+30+ e3x+32y12t+45+149e72x+74y14t+15+116e4x+2y16t+30+ 1081e92x+94y18t+45+125e5x+52y20t+60+ 1121e112x+114y22t+75+163504e6x+3y24t+45+ 1245025e152x+154y30t+90+117424e7x+72y28t+75+ 111025e132x+134y26t+60+11728896400e9x+92y36t+90.
(35)
Therefore, the molecule with four line solitons is produced through equations (17) and (35) (figure 2(c)).
(2) The molecules constructed by the solitons and a lump

The soliton molecule through a line soliton and a lump

For this purpose, after taking the long-wave limit
ki=λiεi, pi=λiliεi, ηi=λi0εi, εi0  (i=1,2),
(36)
the transformation function f of equation (26) is simplified
f=(θ1θ2+b12)(1+eξ3)+(C23θ1+C13θ2+C13C23)eξ3,
(37)
with the phase shift ai4 = 0(i < 4) and ξ4 = 0. Here,
b12=λ1λ2α(l1l2)2β,θi=λi(x+liy4βli2(3a+2bli)t)+λi0  (i=1,2),C13=4αλ1k32αk344β(l1k3p3)2,C23=4αλ2k32αk344β(l2k3p3)2.
(38)
If we take the conjugate constants
λ1=λ2¯=a1+a5i, l1=l2¯=a2+a6ia1+a5i,λ10=λ20¯=a4+a8i,
(39)
the supposing quadratic function
f0=(a1x+a2y+a3t+a4)2+(a5x+a6y+a7t+a8)2+a9,
(40)
is derived from the expression f0θ1θ2 + b12 of equation (37) with the constraint relations of the parameters ai(1 ≤ i ≤ 9) are
a3=4β[3a(a12+a52)(a1a22a1a62+2a2a5a6)+2b(a12a233a12a2a622a1a5a63+6a1a22a5a6+3a2a52a62a23a52)](a12+a52)2,a7=4β[3a(a12+a52)(a5a62a22a5+2a1a2a6)+2b(3a12a22a6a12a63+6a1a2a5a622a1a23a53a22a52a6+a52a63)](a12+a52)2,a9=α(a12+a52)34β(a1a6a2a5)2,
(41)
which the moving route of the lump is
x=(a2a7a3a6)t+a2a8a4a6a1a6a2a5,y=(a1a7a3a5)t+a1a8a4a5a1a6a2a5.
(42)
This solution may produce the soliton molecule through a line soliton and a lump with the condition
k3x+p3y+q3t=k3(xvxt)+p3(yvyt),
(43)
where vx2+vy2 is the moving velocity of this lump. Technically, after selecting the parameters a, b, α, β as equation (29) and
a1=12, a2=a5=1, a4=a8=0,a6=18, k3=2, p3=1, η3=20,
(44)
the transformation function f of equation (37) is
f=[(x2+y71t4)2+(x+y8+31t8)2+518]×(1+e2xy+8t20)20(2x+y8t16)e2xy+8t20.
(45)
Figure 3 is one soliton molecule constructed through a line soliton and a lump at different times with the parameters are equation (44), while the moving velocity of this lump is 483.25.
Figure 3. The soliton molecule consisting of a line soliton and a lump of the solution v from equations (17) and (45) at (a) t=12, (b) t = 0 and (c) t=12, respectively.

Full size|PPT slide

The soliton molecule through the line soliton molecule and a lump

After taking the parameters aij ≠ 0(i < j, i, j = 1, 2, 3, 4), qi ≠ 0(i = 1, 2, 3, 4) of equation (26), and the long-wave limit condition of equation (36), the transformation function f of equation (37) is expanded
f=(θ1θ2+b12)(1+eξ3+eξ4)+(C23θ1+C13θ2+C13C23)eξ3+(C24θ1+C14θ2+C14C24)eξ4+a34[(C23+C24)θ1+(C13+C14)θ2+(C13+C14)(C23+C24)+θ1θ2+b12]eξ3+ξ4,
(46)
with
C14=4αλ1k42αk444β(l1k4p4)2,C24=4αλ2k42αk444β(l2k4p4)2.
(47)
Through this auxiliary function, we can construct the soliton molecule structure of variable v of equation (17) for the (2+1)-dimensional KP-h system (3) by two line solitons and a lump in the molecule with the condition
kix+piy+qit=ki(xvxt)+pi(yvyt), (i=3,4),
(48)
and
k3k4=p3p4=q3q4, k3|k4|, p3|p4|, q3|q4|.
(49)
For example, after selecting the parameters a, b, α, β as equation (29) and
a1=12,a2=a5=1,a4=a8=15,a6=18,k3=32,k4=2,p3=34,p4=1,η3=10,η4=20,
(50)
the transformation function f of equation (46) becomes
f=15(1621189x274y+3334t+14xy2xt22132yt+14x2+1364y2+422564t2)×(1+e32x+34y6t10+e2x+y8t+20)15[275294+149(13x2512y+55112t+14xy2xt22132yt+14x2+1364y2+422564t2)]e72x+74y14t+1020(202xy+8t)e2x+y8t+20103(10032xy+8t)e32x+34y6t10.
(51)
Figure 4 shows the soliton molecule through the line soliton molecule and a lump at different times with the parameters are taken as equation (50).
Figure 4. The soliton molecule consisting of the line soliton molecule and a lump of the solution v from equations (17) and (51) at (a) t=35, (b) t = 0 and (c) t=35, respectively.

Full size|PPT slide

However, if we don't obey the rule of equation (48), but hold on equation (49), the above soliton molecule structure may be destroyed. For example, when taking the parameters a, b, α, β as equation (29), but
a1=12, a2=a5=a6=1, a4=a8=0, k3=32,k4=2, p3=34, p4=1, η3=0, η4=40,
(52)
the function f of equation (46) is
f=(125128+3xy76825xt83225yt+54x2+2y2 + 10245t2)(1+e32x+34y6t+e2x+y8t+40) [494924112544247661696+149(263675307054793x233886527054793y + 9948639232176369825t+3xy76825xt83225yt+54x2 + 2y2+10245t2)]e72x+74y14t+40+4012 × (250415x386y+15257625t)e2x+y8t+40  175934(45007335x6114y+285542425t)e32x+34y6t.
(53)
This conclusion induces the elastic interaction of the line soliton molecule and a lump of the solution v from equations (17) and (53) for the (2+1)-dimensional KP-h system (3). Figure 5 shows that the amplitude, the velocity and the shape of these structures not any change front and back their collision.
Figure 5. The elastic interaction between the line soliton molecule and a lump of the solution v from equations (17) and (53) at (a) t = −2, (b) t = −1.3, (c) t = 0 and (d) t = 1, respectively.

Full size|PPT slide

(3) The molecules induced by the solitons and the breathers

The soliton molecule through a line soliton and the breather

After taking the phase shift ai4 = 0(i < 4) and ξ4 = 0 of equation (26), we let the conjugate relation ξ1=ξ2¯ [31, 32], i.e. the variables ki, pi and ηi (i = 1, 2) satisfy
k1=k2¯=k+κi, p1=p2¯=ρ+σi, η1=η2¯=ζ+μi,
(54)
with i is an imaginary unit. Equation (26) is deduced
f=1+2cos(ϕ1)eϕ2+a12e2ϕ2+eξ3+2[M1cos(ϕ1)M2sin(ϕ1)]eϕ2+ξ3+a12(M12+M22)e2ϕ2+ξ3,
(55)
with
q1=q2¯=ω+Ωi, a13=a23¯=M1+M2i,ϕ1=κx+σy+Ωt+μ, ϕ2=kx+ρy+ωt+ζ,
(56)
and ξ3, a12 also obey equations (24), (25).
When taking the parameters a, b, α, β as equation (29), but
κ=2, σ=35, k=45, ρ=2, μ=ζ=0, η3=30,
(57)
and
k3k=p3ρ=q3ω, k3|k|, p3|ρ|, q3|ω|.
(58)
one can construct the soliton molecule expressed a line soliton and the breather of variable v of equation (17) for the (2+1)-dimensional KP-h system (3) through the function (55) (figure 6(a)).
Figure 6. (a) The soliton molecule expressed a line soliton and the breather for the variable v of equation (17) with the condition equations (55)–(58) at t = 0. (b) The soliton molecule constructed the line soliton molecule and the breather for the variable v of equation (17) with the condition equations (61) and (62) at t = 0. (c) The elastic interaction between the line soliton molecule and the breather with the condition equations (61) and (63) at t = 0. (d) The soliton molecule formed through two breathers with the condition equations (68) and (69) at t = 0.

Full size|PPT slide

The soliton molecule through the line soliton molecule and the breather

When equations (54) and (56) hold with
kik=piρ=qiω, ki|k|, pi|ρ|, qi|ω|, (i=3,4),
(59)
and
a12=M0, a14=a24¯=M3+M4i0, a34=M50,
(60)
equation (26) is induced
f=1+2cos(ϕ1)eϕ2+M0e2ϕ2+eξ3+eξ4+2[M1cos(ϕ1)M2sin(ϕ1)]eϕ2+ξ3+2[M3cos(ϕ1)M4sin(ϕ1)]eϕ2+ξ4+M0[(M12+M22)eξ3+(M32+M42)eξ4]e2ϕ2+2M5[(M1M3M2M4)cos(ϕ1)(M1M4+M2M3)sin(ϕ1)]eϕ2+ξ3+ξ4M5[1+M0(M12+M22)(M32+M42)e2ϕ2]eξ3+ξ4,
(61)
which is an expanded expression of equation (55).
Figure 6(b) shows the soliton molecule through the line soliton molecule and the breather at time t = 0 with the parameters are selected as
a=p4=1, b=α=1, β=3, κ=12,σ=14893400, k=45, ρ=25, μ=ζ=0,k3=32, k4=2, p3=34, p4=1, η3=15, η4=30.
(62)
However, if we don't obey the rule of equation (59), that is, when taking the parameters as equation (62) except for
β=κ=ρ=2, σ=35, η3=5, η4=10,
(63)
the above soliton molecule structure may be destroyed. The interaction of the line soliton molecule and the breather is elastic for the (2+1)-dimensional KP-h system (figure 6(c)).

The soliton molecule through two breathers

Based on equations (54) and (56), we further take
k3=k4¯=k0+κ0i, p3=p4¯=ρ0+σ0i,η3=η4¯=ζ0+μ0i,
(64)
and
q3=q4¯=ω0+Ω0i, ϕ3=κ0x+σ0y+Ω0t+μ0,ϕ4=k0x+ρ0y+ω0t+ζ0,
(65)
but
a12=M0, a13=M1+M2i, a14=M3+M4i,a23=M3M4i, a24=M1M2i, a34=M5,
(66)
with
kk0=ρρ0=ωω0, k|k0|, ρ|ρ0|,ω|ω0|,
(67)
equation (26) is reproduced as
f=1+2cos(ϕ1)eϕ2+M0e2ϕ2+2cos(ϕ3)eϕ4+M5e2ϕ4+2[cos(ϕ1+ϕ3)M1sin(ϕ1+ϕ3)M2+cos(ϕ1ϕ3)M3sin(ϕ1ϕ3)M4]eϕ2+ϕ4+2M0[(cos(ϕ3)M3+sin(ϕ3)M4)M1+(cos(ϕ3)M4sin(ϕ3)M3)M2]e2ϕ2+ϕ4+2M5[(cos(ϕ1)M3sin(ϕ1)M4)M1(cos(ϕ1)M4+sin(ϕ1)M3)M2]eϕ2+2ϕ4+M0M5(M12+M22)(M32+M42)e2(ϕ2+ϕ4).
(68)
For this time, the soliton molecule occurs after taking the following parameters
a=1, b=α=1, β=κ=ρ=2, σ=35, k=45,k0=35, κ0=2, ρ0=32, μ=μ0=ζ=0, ζ0=10,
(69)
which is constructed by two breathers (figure 6(d)).

4. Summary

We employ the integrable KP-h system (3) as the investigated subject which is the (2+1)-dimensional KP equation (1) combining with the KP hierarchy (2). After taking the traveling wave transformation (7), the missing D'Alembert type solution (11) for this equation (3) is first derived. This solution contains several special Kink molecule structures (13)–(16). Further, after introducing the Bäcklund transformation (17) and an auxiliary variable τ, the Hirota bilinear form (21) of equation (3) is conducted and the N-soliton solution (23) with equations (24) and (25) is presented. This N-soliton solution contains some typical soliton molecules including the line solitons, the breathers and a lump (figures 24, 6(a), (b) and (d)) as well as the elastic interaction between the line soliton molecule and a lump/the breather (figures 5 and 6(c)) after selecting some special conditions and parameters. Up to now, these dynamics features have not been reported for this (2+1)-dimensional KP-h system (3). We believe that these structures derived above would be worth underlying in the future research.

References

1
Kadomtsev B B Petviashvili V I 1970 Sov. Phys. Dokl. 15 539
2
Wazwaz A M 2009 Partial Differential Equations and Solitary Waves Theorem Beijing Higher Education Press
3
Wazwaz A M 2011 Appl. Math. Mech. 32 875
4
Chen M R Wang S K Wang X L Wu K Zhao W Z 2015 Nucl. Phys. B 891 655
5
Yamagishi K 1991 Phys. Lett. B 259 436
6
Yu F Wu Y S 1991 Phys. Lett. B 263 220
7
Liu X Zhou Q Biswas A Alzahrani A K Liu W 2020 J. Adv. Research 24 167
8
Han H B Li H J Dai C Q 2021 Appl. Math. Lett. 120 107302
9
Dai C Q Wang Y Y 2020 Nonlinear Dyn. 102 1733
10
Dai C Q Wang Y Y Zhang J F 2020 Nonlinear Dyn. 102 379
11
Li J Chen Y 2020 Commun. Theor. Phys. 72 105005
12
Li J Chen Y 2020 Commun. Theor. Phys. 72 115003
13
Li J Chen Y 2020 Commun. Theor. Phys. 73 015001
14
Wang X L Yu L Chen M R 2015 J. Nonlinear. Math. Phys. 22 194
15
Sato M Sato Y 1983 Nonlinear Partial Differential Equations in Applied Science Tokyo Kinokuniya/North-Holland Press
16
Wazwaz A M 2016 Appl. Math. Lett. 52 74
17
Hortoa R Ito M 1983 J. Phys. Soc. Jpn. 52 744
18
Wazwaz A M 2008 Appl. Math. Comput. 203 402
19
X Chen S T Ma W X 2016 Nonlinear Dyn. 86 523
20
Yan Z W Lou S Y 2020 Appl. Math. Lett. 104 106271
21
Lou S Y 2020 Chin. Phys. B 29 080502
22
Jia M Lou S Y 2020 Chaos Soliton Fract. 140 110135
23
Hirota R 2004 The Direct Method in Soliton Theory Cambridge Cambridge University Press
24
Zhang Z Yang S X Li B 2019 Chin. Phys. Lett. 36 120501
25
Yang S X Zhang Z Li B 2020 Adv. Math. Phys. 2020 2670710
26
Dong J J Li B Yuen M 2020 Commun. Theor. Phys. 72 025002
27
Yang X Y Fan R Li B 2020 Phys. Scr. 95 045213
28
Zhu J M Wang B Ma Z Y Fei J X 2020 Mod. Phys. Lett. B 34 2150115
29
Cao Y L He J S Mihalache D 2018 Nonlinear Dyn. 91 2593
30
Liu Y K Li B An H L 2018 Nonlinear Dyn. 92 2061
31
Dai Z D Liu Z J Li D L 2008 Chin. Phys. Lett. 25 1531
32
Wang C J Fang H Tang X X 2019 Nonlinear Dyn. 95 2943

Funding

* National Natural Science Foundation of China(11705077)
National Natural Science Foundation of China(11775104)
National Natural Science Foundation of China(11447017)

RIGHTS & PERMISSIONS

© 2021 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing
*" title="Share on Weibo" target="_blank">
PDF(621 KB)

258

Accesses

0

Citation

Detail

Sections
Recommended

/