We use a generalized tanh function expansion method and a direct method to study the analytical solutions of the (1+2)-dimensional sine Gordon (2DsG) equation. We obtain some new interaction solutions among solitary waves and periodic waves, such as the kink-periodic wave interaction solution, two-periodic solitoff solution, and two-toothed-solitoff solution. We also investigate the propagation properties of these solutions.
1. Introduction
Sine-Gordon (sG) equation is one of the most famous partial differential equations that have been investigated by many physicists for decades years. The sG equation has played a central role in lots of different scientific fields, such as in differential geometry [1], plasma physics [2], nonlinear optics [3], condensed matter physics [4], quantum field theory [5, 6], and so forth. Researchers have been spending a great deal of effort to generalize (1+1)-dimensional soliton equations to (2+1)-dimensional equations. Remarkable of these equations, in the 1980s, the Nizhnik-Novikov-Veselov (NNV) equation [7–9] and the Davey-Stewartson (DS) equation [10–12] were found. The NNV equation and the DS equation are (2+1) dimensional generalizations of the Korteweg-de Vries (KdV) equation and nonlinear Schrödinger (NLS) equation, respectively. After that, in 1991, Konopelchenko and Rogers [13, 14] proposed a significant symmetry to generalize the (1+1)-dimensional sG equation to (2+1)-dimensional sG equation through a reinterpretation and generalization of a class of infinitesimal Bäcklund transformation. The well-known nonintegrable (2+1)-dimensional sine-Gordon (2DsG) equation is as follows:
(1)∂2Φ∂x2+∂2Φ∂y2-∂2Φ∂t2=msinΦ.
Various methods have been used to study this equation because of its rich symmetrical structure. The brief and effective methods for solving the 2DsG equation include the binary Darboux transformation [15, 16], the extensive symmetry group analysis [17, 18], Hirota’s method [19], Lamb’s method [20, 21], the Painlevé transcendents [22], and the Bäcklund transformation [23]. And researchers have found abundant types of solutions of 2DsG equation, such as the multisoliton solutions and vortex-like solution [24], line and ring solitons [25, 26], curve soliton, point instanton soliton and doubly periodic wave solutions [27–29], Solitoff structure solution, and snake-shape solitary wave solution [30].
Recently, some new useful and powerful methods have been proposed to search for the accurate solutions of nonlinear partial differential equations, such as the general algebra method for the coupled Schrödinger-Boussinesq equations [31], the general mapping deformation method for the generalized variable-coefficient Gardner equation with forcing term [32], the generalized tanh function expansion method for the Abowitz-Kaup-Nwell-Segur system [33], the bosonized supersymmetric KdV model [34], and the Broer-Kaup system [35]. Significantly, the generalized tanh function expansion method is an effective new technique for us to obtain some new interaction solutions of 2DsG equation. Also, we can solve the 2DsG equation by a direct method based on the mapping relations between 2DsG equation and the cubic nonlinear Klein-Gordon (CNKG) equation. This method can be also applied to solve the double sine-Gordon equation, the triple sine-Gordon equation, and the Ginzburg-Landau equation [36], and so forth. In this paper, we want to seek more interaction solutions of new types among solitary waves and periodic waves of the 2DsG equation by the generalized tanh function expansion method and the direct method.
This paper is organized as follows. In Section 2, a kink-periodic wave interaction solution of 2DsG equation is obtained by using of the generalized tanh function expansion method. In Section 3, two-periodic solitoff solution, periodic soliton-periodic travelling wave interaction solution, two-toothed-solitoff solution, and periodic solitoff-kink interaction solution of 2DsG equation are obtained by using the direct method. In Section 4, a short summary and discussions are given.
2. Kink-Periodic Wave Interaction Solutions
The 2DsG equation (1) cannot be solved directly by the generalized tanh function expansion method [33–35], and to find some soliton-periodic wave interaction solutions of 2DsG equation, we suppose
(2)Φ=-iln[W(X,T)],
and take the following coordinates transformation:
(3)X=x+α1y+β1t,T=α2x+y+β2t.
Then, we substitute (2) with (3) into (1) and arrive at
(4)WXWT-WWXT+W-W3=0,
with the constants α1, α2, β1, and β2 satisfying
(5)α22=β22-1,β12=1+α12,α1=β1β2-α2+m4.
It is worth noting that (4) can be solved by using the generalized tanh function expansion method. Firstly, we set
(6)W=u2tanh2(Ψ)+u1tanh(Ψ)+u0,
where u2, u1, u0, and Ψ are functions of variables (X,T). In order to obtain some soliton-periodic wave interaction solutions, let
(7)Ψ=ξ1+ψ1(ξ),
where ξ1=k1X+ω1T, ξ=kX+ωT, in which k, ω, k1, and ω1 are undetermined constants. Then, we substitute (6) and (7) into (4) and analyse the coefficients of function tanh(Ψ) order by order we get the expression of W(8)W=2[ω1+ψ2(ξ)ω][k1+ψ2(ξ)k]tanh2(Ψ)-2ψ3(ξ)kωtanh(Ψ)+ψ32(ξ)k2ω22[ω1+ψ2(ξ)ω][k1+ψ2(ξ)k],
where the functions ψ2(ξ) and ψ3(ξ) satisfy
(9)ψ2(ξ)=ψ1ξ(ξ),ψ3(ξ)=ψ2ξ(ξ).
Furthermore, ψ2(ξ) is a solution of the following Jacobi elliptic function equation:
(10)ψ2ξ2(ξ)=a0+a1ψ2(ξ)+a2ψ22(ξ)+a3ψ23(ξ)+4ψ24(ξ),
with these parameters a0, a1, a2, and a3 satisfying
(11)a0=2ω1k1(2γk1ω1g+γ2k1k-ω1ω)k2ω2γg,a1=2(4γk1ω1gh+γ2k1kq-ω1ωp)k2ω2γg,a2=2(4γk1kω1ωg+γ2k2p+2γgh2-ω2q)k2ω2γg,a3=2(γ2k2+4γgh-ω2)kωγg,
where g=kω1-k1ω, h=kω1+k1ω, p=(3h-g)/2, and q=(3h+g)/2, in which kω1≠k1ω and γ is a constant.
Now we choose the sine Jacobi elliptic function as a solution of (10),
(12)ψ2(ξ)=asn(bξ),
and the functions ψ1(ξ) and ψ3(ξ) are easy to be obtained
(13)ψ1(ξ)=aln[dn(bξ)-ncn(bξ)]bn,ψ3(ξ)=abcn(bξ)dn(bξ).
Then substituting (12), (13), and (11) into (10), relationships of these parameters are written as
(14)b=2an,ω1=-k1ω(γk2[γn2-4a2ω2(n2+1)]+ω2n2(8γk12+1))k(γk2[γn2+4a2ω2(n2+1)]-ω2n2(8γk12+3)),ω=nkγ-2(1+8γk12)[2a2k2γ(n2+1)-n2(4γk12+1)]+n2,γ=-n[-k1n+a2k2(n2+1)-k12n2]4k1[a2k2(n2+1)-2k12n2],k1=2ka2(n2+1)±aa2(n2+1)2-b2n42n,
where n is the modulus of the Jacobi elliptic function sn(z)=sn(z,n).
Finally, the accurate expression of Φ is gained:
(15)Φ=-iln[[[ξ1+aln(D-nC)bn]2(ω1+aωS)(k1+akS)mmmm×tanh[ξ1+aln(D-nC)bn]-abkωCD]2mmmm×(2(ω1+aωS)(k1+akS))-1[ξ1+aln(D-nC)bn]]],
where S=sn(bξ), C=cn(bξ), and D=dn(bξ). The solution of (15) denotes a kink-periodic wave interaction solution of 2DsG equation. Velocities of these two travelling waves are v1=(k1β1+ω1β2)/[(k1+ω1α2)2+(k1α1+ω1)2]1/2 and v2=(kβ1+ωβ2)/[(k+ωα2)2+(kα1+ω)2]1/2, respectively.
Figure 1 shows the density distribution of a kink-periodic wave interaction solution on the x-y plane given by [-exp(iΦ)] and (15) with these parameters
(16)a=45,n=910,k=72,b=169,k1=289,ω=-135108512,ω1=-27102128,m=14,α1=-2,α2=-33+51364,β1=5,β2=-335+6564.
at time t=1. This figure exhibits a special interaction structure of a kink and a periodic wave. Figure 2 shows the propagation of the kink-periodic wave solution at y=0 and t=1. In this figure, the soliton propagates along the negative direction of the x-axis, and its velocity is quicker than the one of the periodic wave, which also propagates along the negative x direction.
The density distribution of a kink-periodic wave interaction solution [-exp(iΦ)] and (15) with (16) on the x-y plane at time t=1.
The propagation of the kink-periodic interaction solution [-exp(iΦ)] and (15) with (16) at y=0 and t=1.
In this section, we use the direct method to study the 2DsG equation. Based on the Lamb substitution [20, 21], the solution of (1) can be set to the following form:
(17)Φ(x,y,t)=4arctan[M(x,y,t)],
in which the function M(x,y,t) is the solution of the CNKG equation [30, 36],
(18)Mxx+Myy-Mtt=λM+μM3,
under the constrained condition
(19)Mx2+My2-Mt2=λM2+μ2M4+μ2,
with m=λ-μ. Function M(x,y,t) can be various styles, such as exp, tanh, sn, and dn [36]. Here we take
(20)M=nsn(V~),
where function V~=(|m|V)/(n+1), in which V is a function of variables (x,y,t), and the constant n is the modulus of the Jacobi elliptic function. Then, we substitute (17) and (20) into 2DsG equation and get
(21)|m|((∇~V)2+m|m|)(n2+n)sn3(V~)+n|m|□Vdn(V~)cn(V~)sn2(V~)-|m|((∇~V)2+m|m|)(n+1)sn(V~)+n|m|□Vdn(V~)cn(V~)=0,
with the constrained conditions
(22)□V=0,(∇~V)2=±1.
Here we define
(23)□=∂x2+∂y2-∂t2,(∇~)2=(∂x)2+(∂y)2-(∂t)2;
then an arbitrary function v(ξ) can be included in the function V by solving (22), namely,
(24)V=v(ξ)+ξ0=v(k11x+k12y+ω1t)+k01x+k02y+ω0t,
and parameters k01, k02, k11, k12, ω0, and ω1 satisfy
(25)k012+k022-ω02=±1,k112+k122-ω12=0,k01k11+k02k12-ω0ω1=0,
where the sign “±” in (22) and (25) takes “−” when m>0 and takes “+” when m<0. Due to the existence of the arbitrary functions, abundant exact solutions of (1) will be obtained as long as the function v(ξ) is properly selected.
When we take
(26)v(ξ)=ξarctan(ξ),
a (2+1)-dimensional two-periodic solitoff solution of 2DsG equation can be obtained:
(27)Φ=4arctan(+k01x+k02y+ω0t((k11x+k12y+ω1t))|m|1+n))nsn(|m|1+n(..(k11x+k12y+ω1t)×arctan(k11x+k12y+ω1t)+k01x+k02y+ω0t..((k11x+k12y+ω1t))|m|1+n)).
We know that a solitoff is defined as a half line soliton. The solution of (27) indicates a solitoff type solution constructed by two travelling waves that propagate in different directions. Velocities of these two travelling waves are v1=ω0/k012+k022≠1 and v2=ω1/k112+k122=1, respectively.
Figure 3 shows a two-periodic solitoff solution (27) with these parameters
(28)k01=1,k02=3,k11=1.2,k12=1.6,ω0=3,ω1=2,m=-1,n=0.9
at time t=0. The angle of the two-periodic solitoff in this figure is actually an obtuse angle although it seems to be orthogonal. It is because k→1·k→2=k01k11+k02k12=ω0ω1≠0. Figure 4 shows more details of the two-periodic solitoff solution (27) with (28). The two-periodic solitoff solution with different wavelength has the same amplitude and keeps the peak unchanged during the propagation process. Their phase velocities are different, but their travelling directions are same; they propagate along the negative y-axis.
A two-periodic solitoff solution of 2DsG equation (27) with (28) at time t=0.
The propagation of two-periodic solitoff solution (27) with (28) along the y-axis when x=0 at (a) t=-5 and (b) t=25.
A periodic soliton-periodic travelling wave interaction solution of 2DsG equation can be obtained:
(29)Φ=4arctan(+3)|m|1+n))nsn(|m|1+n(..k01x+k02y+ω0t+7sech2(k11x+k12y+ω1t)+3..)|m|1+n)),
by choosing
(30)v(ξ)=7sech2(ξ)+3.
Figure 5(a) shows the periodic soliton-periodic travelling wave interaction solution (29) with these parameters
(31)k01=1,k02=3,k11=3,k12=4,ω0=3,ω1=5,m=-1,n=0.9
at time t=0. The solitoff-type structure solution does not appear, whereas these two travelling waves propagate in the different directions. The graph is similar to the soliton-periodic interaction wave in [33], but the soliton really has the periodicity and the peak of the soliton keeps periodically changing. Figure 5(b) shows the density distribution of Φ on the x-y plane.
(a) A periodic soliton-periodic travelling wave interaction solution (29) with (31) at time t=0. (b) The density of Φ on the x-y plane.
Furthermore, if we take
(32)v(ξ)=ξ2+1+54sin3(ξ),v(ξ)=3cos3(ξ)cn(ξ)+1,
then a two-sawtooth-solitoff solution and a periodic solitoff-kink interaction solution of 2DsG equation can be written as
(33)Φ=4arctan(+(k11x+k12y+ω1t)2+154)))nsn(|m|1+n(54sin3(k11x+k12y+ω1t)+k01x+k02y+ω0t+(k11x+k12y+ω1t)2+154)|m|1+n)),(34)Φ=4arctan(×cn(k11x+k12y+ω1t))|m|1+n))nsn(|m|1+n(..k01x+k02y+ω0t+1+3cos3(k11x+k12y+ω1t)×cn(k11x+k12y+ω1t)..)|m|1+n)),
respectively. Figure 6 shows a two-toothed-solitoff solution (33) with (31) in the limit case of the modulus n=1. The two-toothed-solitoff structure is constructed by a kink soliton and an antikink soliton. Their travelling velocities are different, but group velocities are the same. And travelling directions of these two solitoff waves construct a constant acute angle during the propagation process.
(a) A two-toothed-solitoff solution (33) with (31) at time t=1 except for the modulus n=1. (b) The density of Φ on the x-y plane.
Figure 7(a) displays a periodic solitoff-kink interaction solution constructed by a bright soliton and a kink soliton. Figures 7(b)–7(d) show that the bright soliton and the kink soliton have different travelling velocities, and they propagate along the negative x-axis. The peak of the bright soliton keeps increasing until it is arriving at the same amplitude of the kink soliton.
(a) A periodic solitoff-kink interaction solution (34) with (28) at time t=0 except for the modulus n=1. (b)–(d) show the propagation of the periodic solitoff-kink interaction solution along the x-axis when y=-1.3 at t=-3, t=-0.5, and t=1.
4. Summary and Discussion
First of all, we use the generalized tanh function expansion method to solve the 2DsG equation; a special new kink-periodic wave interaction solution is explicitly expressed both analytically and graphically. This interaction solution between tanh-type soliton and periodic wave of 2DsG equation is firstly obtained. Then, we use the direct method and obtain more new interaction solutions of the 2DsG equation, including the two-periodic solitoff solution (27), periodic soliton-periodic travelling wave interaction solution (29), two-toothed-solitoff solution (33), and periodic solitoff-kink interaction solution (34). The solution (34) is a generalization of a single straight-line kink soliton solution, while the solution (33) is an alternative generalization of periodic straight-line solitoff type of kink soliton solution. These types of interaction solutions are also firstly found for the 2DsG equation. All of these solutions indicate the interaction solution among solitary waves and periodic waves; their travelling velocities are different, but group velocities are same, and they propagate in different trajectories which contain linear shape, curve shape, and saw-tooth shape. In fact, the forms of (12) and (20) can be not only taken the sine Jacobi elliptic function (sn), more functions can be selected such as exp, cn, and cn/sn, and more explicit solutions can be gained. The abundant solutions solved by these two methods suggest that the rich structures of nonlinear systems do not only exist in the integrable systems but also in the nonintegrable systems. Furthermore, there are some types of localized solutions decaying in all directions, for instance, the dromions and ring solitons have not been found by these two methods; those will be left for us to do more research.
Conflict of Interests
The authors declare that they have no financial relationships with other people or organizations that can inappropriately influence this work or possible conflict of interests.
Acknowledgments
The work was supported by the National Natural Science Foundation of China no. 11175158 and by program for Innovative Research Team in Zhejiang Normal University.
EisenhartL. P.1960New York, NY, USADover PublicationsMR0115134WashimiH.TaniutiT.Propagation of ion-acoustic solitary waves of small amplitude196617199969982-s2.0-424384435710.1103/PhysRevLett.17.996LeblondH.MihalacheD.Ultrashort light bullets described by the two-dimensional sine-Gordon equation20108162-s2.0-7795370161910.1103/PhysRevA.81.063815063815LoutsenkoI.RoubtsovD.Critical velocities in exciton superfluidity19977815301130142-s2.0-000088648310.1103/PhysRevLett.78.3011SamuelS.Grand partition function in field theory with applications to sine-Gordon field theory1978186191619322-s2.0-3374473072710.1103/PhysRevD.18.1916BaroneA.EspositoF.MageeC. J.ScottA. C.Theory and applications of the sine-gordon equation1971122272672-s2.0-5124918880810.1007/BF02820622NovikovS. P.VeselovA. P.Two-dimensional Schrödinger operator: inverse scattering transform and evolutional equations1986181–326727310.1016/0167-2789(86)90187-9MR838332ZBL0609.35082AthorneC.NimmoJ. J. C.On the Moutard transformation for integrable partial differential equations199176, article 80910.1088/0266-5611/7/6/005MR1140316ZBL0737.35091RadhaR.LakshmananM.Singularity analysis and localized coherent structures in (2+1)-dimensional generalized Korteweg-de Vries equations19943594746475610.1063/1.530812MR1290899BoitiM.LeonJ. Jp.MartinaL.PempinelliF.Scattering of localized solitons in the plane19881328-943243910.1016/0375-9601(88)90508-7MR965704FokasA. S.SantiniP. M.Dromions and a boundary value problem for the Davey-Stewartson (2+1) equation1990441-29913010.1016/0167-2789(90)90050-YMR1069674ZBL0707.35144DaveyA.StewartsonK.On three-dimensional packets of surface waves1974338101110KonopelchenkoB. G.RogersC.On (2+1)-dimensional nonlinear systems of Loewner-type1991158839139710.1016/0375-9601(91)90680-7MR1130347KonopelchenkoB. G.RogersC.On generalized Loewner systems: novel integrable equations in 2+1 dimensions199334121424210.1063/1.530377MR1198634SchiefW. K.On the geometry of an integrable (2+1)-dimensional sine-Gordon system199745319631671168810.1098/rspa.1997.0090MR1469166SchiefW. K.On localized solitonic solutions of a (2+1)-dimensional sine-Gordon system19922524L1351L1354MR121087310.1088/0305-4470/25/24/002ClarksonP. A.MansfieldE. L.MilneA. E.Symmetries and exact solutions of a (2+1)-dimensional sine-Gordon system199635417131807183510.1098/rsta.1996.0079MR1400220LouS.-y.Symmetry analysis and exact solutions of the (2+1) dimensional sine-Gordon system20004196509652410.1063/1.1286770MR1779659ZBL0976.35067HirotaR.Exact three-soliton solution of the two-dimensional sine-Gordon equation197335515662-s2.0-5614912039710.1143/JPSJ.35.1566ZagrodzińskiJ.Particular solutions of the sine-Gordon equation in 2+1 dimensions1979724-528428610.1016/0375-9601(79)90469-9MR589817ChristiansenP. L.OlsenO. H.Return effect for rotationally symmetric solitary wave solutions to the sine-Gordon equation197868218518810.1016/0375-9601(78)90798-3MR602933KaliappanP.LakshmananM.Kadomtsev-Petviashvile and two-dimensional sine-Gordon equations: reduction to Painlevé transcendents19791210L249L252MR54538810.1088/0305-4470/12/10/002KonopelchenkoB. G.SchiefW.RogersC.A (2+1)-dimensional sine-Gordon system: its auto-Bäcklund transformation19921721-2394810.1016/0375-9601(92)90186-PMR1197747TakenoS.Multi-(resonant-soliton)-soliton solutions and vortex-like solutions to two- and three-dimensional sine-Gordon equations198268399299510.1143/PTP.68.992MR683543ZBL1194.35391DjidjeliK.PriceW. G.TwizellE. H.Numerical solutions of a damped sine-Gordon equation in two space variables199529434736910.1007/BF00042761MR1345523ZBL0841.65083DehghanM.MirzaeiD.The dual reciprocity boundary element method (DRBEM) for two-dimensional sine-Gordon equation20081976-847648610.1016/j.cma.2007.08.016MR2396038ZBL1169.76401LouS.-y.Symmetry analysis and exact solutions of the 2+1 dimensional sine-Gordon system20004196509652410.1063/1.1286770MR1779659ZBL0976.35067ChowK. W.A class of doubly periodic waves for nonlinear evolution equations2002351719010.1016/S0165-2125(01)00078-6MR1869308ZBL1163.74329TangX.-Y.LouS.-Y.ZhangY.Localized excitations in (2+1)-dimensional systems20026641704660110.1103/PhysRevE.66.046601MR1935202LouS. Y.HuH.-C.TangX.-Y.Interactions among periodic waves and solitary waves of the (N+1)-dimensional sine-Gordon field20057132-s2.0-4134909978210.1103/PhysRevE.71.036604036604HongB.LuD.New exact Jacobi elliptic function solutions for the coupled Schrödinger-Boussinesq equations20132013717083510.1155/2013/170835HongB.LuD.New exact solutions for the generalized variable-coefficient Gardner equation with forcing term201221952732273810.1016/j.amc.2012.08.104MR2988148LouS. Y.ChengX. P.TangX. Y.Interactions between solitons and other nonlinear Schrödinger waveshttp://arxiv.org/abs/1208.5314GaoX. N.LouS. Y.TangX. Y.Bosonization, singularity analysis, nonlocal symmetry reductions and exact solutions of supersymmetric KdV equation2013529MR3080912ChenC. L.LouS. Y.CTE solvability and exact solution to the Broer-Kaup system20133024110202LouS. Y.NiG. J.The relations among a special type of solutions in some (D+1)-dimensional nonlinear equations19893071614162010.1063/1.528294MR1002269