Nonlocal symmetries,soliton-cnoidal wave solution and soliton molecules to a (2+1)-dimensional modified KdV system

Jianyong Wang and Bo Ren

1 Department of Mathematics and Physics,Quzhou University,Quzhou 324000,China

2 Department of Mathematics,Zhejiang University of Technology,Hangzhou 310014,China

Abstract A (2+1)-dimensional modified KdV (2DmKdV) system is considered from several perspectives.Firstly,residue symmetry,a type of nonlocal symmetry,and the Bäcklund transformation are obtained via the truncated Painlevé expansion method.Subsequently,the residue symmetry is localized to a Lie point symmetry of a prolonged system,from which the finite transformation group is derived.Secondly,the integrability of the 2DmKdV system is examined under the sense of consistent tanh expansion solvability.Simultaneously,explicit soliton-cnoidal wave solutions are provided.Finally,abundant patterns of soliton molecules are presented by imposing the velocity resonance condition on the multiple-soliton solution.

Keywords: soliton-cnoidal wave solution,nonlocal symmetry,soliton molecule,(2+1)-dimensional modified KdV system

1.Introduction

The study of interaction waves between solitons and other types of nonlinear waves is a topic of great interest in soliton theory.In recent years,there has been remarkable progress in finding the analytical solution of soliton-cnoidal waves due to the introduction of some new effective methods,such as the localization of nonlocal symmetry[1–6],the consistent Riccati expansion method [7,8],and the consistent tanh expansion(CTE)method[9–11].In fact,such types of interaction waves have been observed in real world and laboratory experiments.For example,it has been reported in an experimental observation that truly traveling waves,consisting of a strongly localized solitary wave residing on a small amplitude oscillatory tail,are generated from a woodpile lattice [12].

Recently,a bound state of multiple-soliton solution,called soliton molecule,has rapidly become a popular topic in both theoretical study and experimental observation,ranging from optical fibers [13,14] to Bose–Einstein condensation [15,16].To construct abundant patterns of soliton molecules,Lou proposed a velocity resonance mechanism for (1+1)-dimensional fluid systems [17–21].Based on this mechanism,soliton molecules have been discovered in several famous integrable systems [22–29],including the complex modified KdV equation [22],the (2+1)-dimensional bidirectional Sawada–Kotera system [23],and the Sharma–Tasso–Olver–Burgers system [24].In addition,dromion molecules described by the variable separation solution were reported for the first time by applying the multilinear variable separation approach to the(3+1)-dimensional Boiti–Leon–Manna–Pempinelli system [29].

In this work,we consider the following (2+1)-dimensional modified KdV (2DmKdV) system [30]

which reduces to the (2+1)-dimensional modified Bogoyavlenskii–Schiff equation [31,32] for α=β=0 and λ=–1.The existence of a three-soliton solution of the 2DmKdV system (1) is confirmed in [30].It is noteworthy that,as detailed in [33],Hu developed a bilinear Bäcklund transformation and a nonlinear superposition formula to the 2DmKdV system (1).

The paper is structured as follows.In section 2,the nonlocal residue symmetry and the Bäcklund transformation of the 2DmKdV system with λ=–1 are obtained from its truncated Painlevé expansion.Then,new auxiliary functions are introduced to localize the nonlocal symmetry to Lie point symmetry,and thus a finite transformation group is derived.In section 3,the CTE method is utilized to prove the integrablity of the 2DmKdV system with λ=–1 and to derive the interaction wave solution between one-soliton/kink and the surrounding periodic waves.In section 4,the 2DMKdV system with λ=1 is transformed to a bilinear form and the multiple-soliton solution is presented.By imposing the velocity resonance condition on the multiple-soliton solution,abundant patterns of soliton molecules are generated.Some concluding remarks are given in the last section.

2.The truncated painlevé expansion and its related nonlocal residue symmetry

According to the truncated Painlevé expansion method,the Laurent series expansions ofuandvin (1) read

whereu0,u1,v0,v1,v2are undetermined functions of(x,y,t),and φ=φ(x,y,t)is the singularity manifold.Substituting (2)into (1) and then eliminating the coefficients of all different powers of φ-1lead to the solutions of the unknown functions

where the variableCand the Schwarzian derivativeSare defined as

Consequently,the Bäcklund transformation of (1) is derived as

Putting equations (1) and (5) together,the 2DmKdV system is successfully transformed into the Schwarzian version

Obviously,the Schwarzian equation (6) is invariant under the Möbius transformation φ→(a+bφ)/(c+dφ),which implies φ possesses the following Lie point symmetry

In the truncated Painlevé analysis,it is readily verified that (u1,v2) is a solution of the 2DmKdV system.Additionally,the coefficients of φ-1,namely,(u0,v1),can be verified to be a nonlocal residue symmetry of the 2DmKdV system,given by

In order to find the finite transformation group associated with the nonlocal symmetry (8),one usually solves the corresponding initial value problem in terms of the Lie’s first principle.However,it is difficult to do so due to the existence of the partial derivatives of φ.Fortunately,the nonlocal symmetry can be localized to the following Lie point symmetry

by introducing some new auxiliary functions

to form a prolonged system constituted by equations(1),(5)and(10).Therefore,the corresponding initial value problem reads

Solving the above initial value problem,one can deduce the finite transformation group to the enlarged system as stated in the following theorem.

Theorem 1.If{u,v,φ,p,q,r}is a solution of the enlargedsystem of equations(1),(5),and(10),thenis given by

with∊being an arbitrary group parameter.

3.Soliton-cnoidal wave solution

Balancing the nonlinearity and dispersive terms in (1),it is clear that the solution of the 2DmKdV system has the following CTE

whereu0,u1,v0,v1,v2andware functions of (x,y,t) to be determined later.

Substituting the expansions (13) into (1) and vanishing the coefficients of different powers oftanh(w),we obtain twelve overdetermined equations for six unknown functions {u0,u1,v0,v1,v2,w}.Solving them leads to

where the variableCand the Schwarzian derivativeSare defined as

Consequently,the CTE (13) can be rearranged in terms ofw

followed by the associated compatibility condition ofw

upon the substitution of (16) into (1).

To find the soliton-cnoidal wave solution of the 2DmKdV system (17),we consider a special solution in the form of

The direct substitution of (18) into (17) yields an ordinary differential equation aboutW1(=Wξ),which can be identically satisfied by introducing the following elliptic function equation

Figure 1. Profiles of the kink-cnoidal wave solution u given by equation (22) with α=1,β=0.5,δ=1,k1=2, k2=-3,l1=1,l2=3,m=0.2,n=1,and ω1=2.(a)The plot of the wave at y=t=0;(b)Three-dimensional plot of the wave at t=0;(c)Space-time evolution of the wave in the x-t plane at y=0.

Figure 2. Profiles of the kink-cnoidal wave solution v in equation (22) with α=1,β=0.5,δ=–1,k1=1.25, k2=–0.75,l1=1,l2=2,m=0.5,n=1.25,and ω1=2.(a)The plot of the wave at y=t=0;(b)3D plot of the wave lies in x-y plane at t=0;(c)Spatiotemporal density plot of the wave interaction with y=0.

The solution of the elliptic function equation(19)can be expressed in different types of Jacobi elliptic functions.To show abundant interactions between one soliton and the surrounded cnoidal periodic waves,we take an ansatz ofwas

whereSn(nξ,m),Cn(nξ,m) andDn(nξ,m) are the Jacobi elliptic functions with an argumentnξ and a modulusm.

Substituting the ansatz (20) into the compatibility condition (17),and then eliminating the coefficients of different powers of the Jacobi elliptic functions,a set of overdetermined equations for ten wave parameters {k1,l1,ω1,k2,l2,ω2,a0,a2,m,n} is obtained.From these equations,a0,a1and ω2are determined

The combination of(16),(20),and(21)gives the explicit soliton-cnoidal wave solution

withSn≡Sn(nξ,m),Cn≡Cn(nξ,m),Dn≡Dn(nξ,m).

The wave structure of the solutionuin (22) is shown in figure 1,as can be observed that a kink resides on a cnoidal periodic wave background.While the wave solution ofvin(22),as exhibited in figure 2,displays a bell-shaped soliton core surrounded by a cnoidal periodic wave.It is clearly demonstrated in both figure 1(c) and figure 2(c) that the interactions between a kink/soliton and the surrounded cnoidal periodic wave are completely elastic except for a phase shift.

4.N-soliton solution and soliton molecules

Through a bi-logarithmic transformation

wheregandg*are complex functions of (x,y,t),the 2DmKdV system (1) with λ=1 can be converted into the bilinear form

where the Hirota derivative is defined as

To obtain the multiple soliton solution of the 2DmKdV system,we expandgandg*as power series with respect to a small parameter ∊

Substituting (25) into (24) and then eliminating the coefficients of different powers of ∊,we obtain the recursion relations(linear equation)ofgiandgi* .Finally,it is found thatNsoliton solutions can be obtained by taking

Case 1 forN=1.In this case,choosingg=1+arrives at a one-soliton solution

Case 2 forN=2.Ifg=1+,then a two-soliton solution can be produced as

Figure 3(a) presents a two-soliton structure of the solutionu.It is known that a soliton molecule can be obtained from the two-soliton solution by imposing a velocity resonance condition [17].Here,the following type of the resonance conditions is introduced

Figure 3.Profiles of the wave solution u given by (29) with α=0.5,β=0.25 for (a) The two-soliton solution with k1=0.25,k2=-0.2,l1=l2=1.25,ξ10=ξ20=0;(b) The soliton molecule with k1=1,k2=,ξ10=-ξ20=10;(c) An asymmetric soliton with k1=1,k2=,ξ10=ξ20=0.

Figure 4. (a) 3D plot of the three-soliton solution u with k1=1.5,k2=1,k3=0.5,l1=l2=l3=0.5,ξ10=ξ20=ξ30=0,α=0.2 and β=0.04.(b)3D plot of the soliton-asymmetric soliton solution with k1=0.25,k2=0.75,k3=–0.5,l3=0.5,ξ10=ξ20=0,ξ30=–3,α=2 and β=2.(c) Density plot of the space-time evolution of the soliton-asymmetric soliton interaction.

Figure 5.The soliton molecule patterns from the four soliton solution.(a) The density plot of the space-time evolution of two two-soliton molecules with k1=–1.5,k2=–0.4,k3=–0.9,k4=–0.3,ξ10=–ξ20=–15,ξ30=-ξ40=–10,α=0.5 and β=1.(b) 3D plot of the interaction between two asymmetric solitons with the same parameters as (a),except for ξ10=–ξ20=0.25,ξ30=0 and ξ40=4.5.(c) The density plot of the space-time evolution of a special soliton molecule with k1=1,k2=0.5,Kr=0.75,Ki=–0.25,ξ10=–ξ20=–10,α=0.5 and β=0.25

which requires

With the above parameter values (32),the two-soliton solution (29) becomes a soliton molecule,which is graphically displayed in figure 3(b)and(c).It is obviously observed from figure 3(b) that,despite the equal velocity,the height and width of the two peaks in the soliton molecule are visibly different.Additionally,the distance between the two peaks is determined by the phase constants ξi0.It is noted that,under suitable values of the phase constants,the two peaks of the molecule can get close enough to each other,and thus leads to the formation of an asymmetric soliton as exhibited in figure 3(c).

Case 3 forN=3.In this case,the explicit form ofgreads

By substituting(33)into(23),one can immediately obtain the three-soliton solution as shown in figure 4(a).Here,we omit the lengthy expressions of the wave solutions.Similarly,by imposing the velocity resonance condition (31) on the first two solitons,we findl1=-k1[α+] andl2=-k2[α+].Consequently,the interaction between one soliton and one two-soliton molecule is produced,as shown in figure 4(b),and the corresponding spacetime evolution is displayed in figure 4(c).

Case 4 forN=4.Now we have a four-soliton solution withggiven by

Clearly,the larger theN,the more abundant the patterns of soliton molecules provided by the solution.For instance,applying the velocity resonance condition(31)on two pairs of solitons,respectively,arrives at

Thereafter,the four-soliton solution degenerates into two two-soliton molecules,as shown in figure 5(a),where one group of parallel lines,representing a two-soliton molecule,intersects elastically with the other group of parallel lines representing the other two-soliton molecule.Similar to the previous situation,suitable choices of the phase constants ξi0will lead to an elastic interaction between two asymmetric solitons,as exhibited in figure 5(b).It is interesting to remark that two peaks embedded in an asymmetric soliton will get separated from each other after the interaction.In addition,one can also construct a special soliton molecule which is a bound state of two solitons and one breather,as shown in figure 5(c).To construct such a soliton molecule,we restrict

where Ωrand Ωisatisfy the following dispersion relation

The introduction of the velocity resonance conditions ω1/k1=ω2/k2=Ωr/Krand ω1/l1=ω2/l2=Ωr/Lrleads to

5.Conclusions and discussions

In summary,several aspects of the 2DmKdV system are studied.Based on the truncated Painlevé expansion method,we obtain the nonlocal residue symmetry and the Bäcklund transformation.The nonlocal symmetry is thus localized to the Lie point symmetry of the enlarged system via introducing several new auxiliary functions.By virtue of the CTE method,the soliton-cnoidal wave solutions are explicitly presented and graphically displayed.Finally,several types of soliton molecules are illustrated,such as the asymmetric soliton and the soliton-breather molecule,which are produced from theN-soliton solution through a velocity resonance mechanism.It is hoped that the explicit solutions obtained here may find their physical applications in the near future.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No.12375006).