# Soliton Solutions of the Coupled Schrödinger-Boussinesq Equations for Kerr Law Nonlinearity

Authors: Anwar Ja’afar Mohamad Jawad; Mahmood Jawad Abu-AlShaeer
Journal: Abstract and Applied Analysis
 (2018)

---

## Abstract

In this paper, the coupled Schrödinger-Boussinesq equations (SBE) will be solved by the sech, tanh, csch, and the modified simplest equation method (MSEM). We obtain exact solutions of the nonlinear for bright, dark, and singular 1-soliton solution. Kerr law nonlinearity media are studied. Results have proven that modified simple equation method does not produce the soliton solution in general case. Solutions may find practical applications and will be important for the conservation laws for dispersive optical solitons.


---

## Body

## 1. Introduction

All optical communications are being used for transcontinental and transoceanic data transfer, through long-haul optical fibers, at the present time. There are various aspects of soliton communication that still need to be addressed. One of the features is the dispersive optical solitons. In presence of higher order dispersion terms, soliton communications are sometimes a hindrance as these dispersion terms produce soliton radiation. Nonlinear evolution equations have a major role in various scientific and engineering fields, such as optical fibers. Nonlinear wave phenomena of dispersion, dissipation, diffusion, reaction, and convection are very important in nonlinear wave equations. In recent years, quite a few methods for obtaining explicit traveling and solitary wave solutions of nonlinear evolution equations have been proposed. In recent years, exact homoclinic and heteroclinic solutions were proposed for some NEEs like nonlinear Schrödinger equation, Sine-Gordon equation, Davey-Stewartson equation, Zakharov equation, and Boussinesq equation [1–7].In particular, the study of the coupled Schrödinger-Boussinesq equations has attracted much attention of mathematicians and physicists [8–10]. The existence of the global solution of the initial boundary problem for the equations was investigated in [8]. The existence of a periodic solution for the equations was considered in [9]. Kilicman and Abazari [10] used the (G′/G)-expansion method to construct periodic and soliton solutions for the Schrödinger-Boussinesq. The investigation of nonlinear partial differential equations plays an important role in the study of nonlinear physical phenomena [9–12].The nonlinear coupled Schrödinger-Boussinesq equation (SBE) governs the propagation of optical solitons in a dispersive optical fiber and is a very important equation in the area of theoretical and mathematical physics. This paper is going to take a look at the bright, dark, and singular soliton solutions for Kerr law nonlinearity media.


## 2. Governing Equations

Consider the coupled Schrödinger-Boussinesq equations (SBE). They appeared in [13] as a special case of general systems governing the stationary propagation of coupled nonlinear upper hybrid and magneto sonic waves in magnetized plasma. These equations were in the form [14] (1)i∂E∂t+∂2E∂x2+α1E-NE=0,3∂2N∂t2-∂4N∂x4+3∂2∂x2N2+α2∂2N∂x2-∂2∂x2E2=0,where α1,α2 are real constants, E(x,t) is a complex function, and N(x,t) is a real function. The complete integrability of (1) was studied by Chowdhury et al. [15], andN-soliton solution, homoclinic orbit solution, and rogue solution were obtained by Hu et al. [16], Dai et al. [17–19], and Mu and Qin [20].


## 3. The Traveling Solution

Consider the nonlinear partial differential equation in the form(2)Fu,ut,ux,utt,uxx,…=0,where u(x,t) is a traveling wave solution of nonlinear partial differential equation (2). We use the transformations(3)ux,t=fξ,where ξ=x-ρt+χ. This enables us to use the following changes: (4)∂∂t·=-ρddξ·,∂∂x·=ddξ·.Using (4) to transfer the nonlinear partial differential equation (2) to nonlinear ordinary differential equation, (5)Qf,f′,f′′,f′′′,…=0.The ordinary differential equation (5) is then integrated as long as all terms contain derivatives, where we neglect the integration constants.


## 4. Hyperbolic Function Methods

The solutions of many nonlinear equations can be expressed in the following form.


### 4.1. Sech Function Method (Bright Soliton) [21]

(6)
f
ξ
=
A
sech
β
⁡
μ
ξ
,
f
′
ξ
=
-
A
β
μ
sech
β
⁡
μ
ξ
·
tanh
⁡
μ
ξ
,
f
′
′
ξ
=
-
A
β
μ
2
β
+
1
sech
β
+
2
⁡
μ
ξ
-
β
sech
β
⁡
μ
ξ
,
f
′
′
′
ξ
=
A
β
μ
3
β
+
1
β
+
2
sech
β
+
2
⁡
μ
ξ
-
β
2
sech
β
⁡
μ
ξ
tanh
⁡
μ
ξ
.


### 4.2. Tanh Function Method (Dark Soliton) [22]

(7)
f
ξ
=
A
tanh
β
⁡
μ
ξ
,
ξ
≤
π
2
μ
,
f
′
ξ
=
A
β
μ
tanh
β
-
1
⁡
μ
ξ
-
tanh
β
+
1
⁡
μ
ξ
,
f
′
′
ξ
=
A
β
μ
2
β
-
1
tanh
β
-
2
⁡
μ
ξ
-
2
β
tanh
β
⁡
μ
ξ
+
β
+
1
tanh
β
+
2
⁡
μ
ξ
,
f
′
′
′
ξ
=
A
β
μ
3
β
-
1
β
-
2
tanh
β
-
3
⁡
μ
ξ
-
β
-
1
β
-
2
+
2
β
tanh
β
-
1
⁡
μ
ξ
+
β
+
1
β
+
2
+
2
β
tanh
β
+
1
⁡
μ
ξ
-
β
+
1
β
+
2
tanh
β
+
3
⁡
μ
ξ
.


### 4.3. Csch Function Method (Singular Soliton) [21, 22]

(8)
f
ξ
=
A
csch
β
⁡
μ
ξ
,
f
′
ξ
=
-
A
β
μ
csch
β
⁡
μ
ξ
·
coth
⁡
μ
ξ
,
f
′
′
ξ
=
A
β
μ
2
β
+
1
csch
β
+
2
⁡
μ
ξ
+
β
csch
β
⁡
μ
ξ
,
f
′
′
′
ξ
=
-
A
β
μ
3
β
+
1
β
+
2
csch
β
+
2
⁡
μ
ξ
+
β
2
csch
β
⁡
μ
ξ
coth
⁡
μ
ξ
,where A represent the amplitudes of the solitons and μ represents the solitons width.We substitute (6), (7), or (8) into the reduced equation (5), balance the terms of the sech, tanh, and csch functions, and solve the resulting system of algebraic equations by using computerized symbolic packages. We next collect all terms with the same power in sechk⁡μξ, tanhk⁡μξ, or cschk⁡μξ, set to zero their coefficients to get a system of algebraic equations among the unknowns A, μ, and β, and solve the subsequent system.


## 4.1. Sech Function Method (Bright Soliton) [21]

(6)
f
ξ
=
A
sech
β
⁡
μ
ξ
,
f
′
ξ
=
-
A
β
μ
sech
β
⁡
μ
ξ
·
tanh
⁡
μ
ξ
,
f
′
′
ξ
=
-
A
β
μ
2
β
+
1
sech
β
+
2
⁡
μ
ξ
-
β
sech
β
⁡
μ
ξ
,
f
′
′
′
ξ
=
A
β
μ
3
β
+
1
β
+
2
sech
β
+
2
⁡
μ
ξ
-
β
2
sech
β
⁡
μ
ξ
tanh
⁡
μ
ξ
.


## 4.2. Tanh Function Method (Dark Soliton) [22]

(7)
f
ξ
=
A
tanh
β
⁡
μ
ξ
,
ξ
≤
π
2
μ
,
f
′
ξ
=
A
β
μ
tanh
β
-
1
⁡
μ
ξ
-
tanh
β
+
1
⁡
μ
ξ
,
f
′
′
ξ
=
A
β
μ
2
β
-
1
tanh
β
-
2
⁡
μ
ξ
-
2
β
tanh
β
⁡
μ
ξ
+
β
+
1
tanh
β
+
2
⁡
μ
ξ
,
f
′
′
′
ξ
=
A
β
μ
3
β
-
1
β
-
2
tanh
β
-
3
⁡
μ
ξ
-
β
-
1
β
-
2
+
2
β
tanh
β
-
1
⁡
μ
ξ
+
β
+
1
β
+
2
+
2
β
tanh
β
+
1
⁡
μ
ξ
-
β
+
1
β
+
2
tanh
β
+
3
⁡
μ
ξ
.


## 4.3. Csch Function Method (Singular Soliton) [21, 22]

(8)
f
ξ
=
A
csch
β
⁡
μ
ξ
,
f
′
ξ
=
-
A
β
μ
csch
β
⁡
μ
ξ
·
coth
⁡
μ
ξ
,
f
′
′
ξ
=
A
β
μ
2
β
+
1
csch
β
+
2
⁡
μ
ξ
+
β
csch
β
⁡
μ
ξ
,
f
′
′
′
ξ
=
-
A
β
μ
3
β
+
1
β
+
2
csch
β
+
2
⁡
μ
ξ
+
β
2
csch
β
⁡
μ
ξ
coth
⁡
μ
ξ
,where A represent the amplitudes of the solitons and μ represents the solitons width.We substitute (6), (7), or (8) into the reduced equation (5), balance the terms of the sech, tanh, and csch functions, and solve the resulting system of algebraic equations by using computerized symbolic packages. We next collect all terms with the same power in sechk⁡μξ, tanhk⁡μξ, or cschk⁡μξ, set to zero their coefficients to get a system of algebraic equations among the unknowns A, μ, and β, and solve the subsequent system.


## 5. The Application

The starting hypothesis for solving (1) by the aid of traveling waves solution is as follows: introduce the transformations(9)Ex,t=eiθx,tuξ,Nx,t=vξ,where(10)θ=kx+ωt+ϵ0,ξ=x-ρt+χ,where k,ω,ρ,ϵ0, and χ are real constants. The parameter ρ represents the soliton velocity.Substituting (9) and (10) into (1) and decomposing into real and imaginary parts leads to(11)u′′-ω+k2-α1u-uv=0,(12)2k-ρ=0,(13)12k2+α2v′′-v′′′′+3v2′′-u2′′=0.Substitute (12) in (10), then (14)ξ=x-2kt+χ.Integrating (13) twice with zero constant, (13) can be written as(15)12k2+α2v-v′′+3v2-u2=0.


### 5.1. Bright Soliton

Seeking the solution by sech function method as in (6)(16)uξ=A1sechβ1⁡μξ,vξ=A2sechβ2⁡μξ,the system of equations in (11) and (15) becomes, respectively,(17)β1μ2β1+1sechβ1+2⁡μξ-β1sechβ1⁡μξ+ω+k2-α1sechβ1⁡μξ-A2sechβ1+β2⁡μξ=0,(18)12k2+α2A2sechβ2⁡μξ+A2β2μ2β2+1sechβ2+2μξ-β2sechβ2⁡μξ+3A22sech2β2⁡μξ-A12sech2β1⁡μξ=0.Equating the exponents and the coefficients of each pair of the sech functions, we find (19)2β1=β2+2,β1+β2=β1+2,thenβ1=β2=2.Thus setting coefficients of (17)-(18) to zero yields set system of equations:(20)4μ2-ω+k2-α1=0,12k2+α2-4μ2=0,6μ2-A2=0,6A2μ2+3A22-A12=0.Solving the system of equations in (20), we get(21)A1=∓312k2+α2,A2=3212k2+α2,(22)μ=∓12k2+α24,ω=11k2+α1+α2,(23)E1x,t=∓3eikx+11k2+α1+α2t+ϵ012k2+α2sech2⁡12k2+α24x-2kt+χ,(24)N1x,t=3212k2+α2sech2⁡12k2+α24x-2kt+χ.For k=α1=1,α2=4,ϵ0=χ=0, the real part of E1x,t=48cos⁡x+16tsech2⁡{2x-2t}, and N1x,t=24sech2⁡{2x-2t}.Figures1 and 2 represent the solitary wave of the real part of E1(x,t) in (23) and N1(x,t) in (24), respectively, for -5≤x≤5, 0≤t≤1.Figure 1
The solitary wave of the real part ofE1(x,t) in (23) for -5≤x≤5, 0≤t≤1.Figure 2
The solitary waveN1(x,t) in (24) for -5≤x≤5, 0≤t≤1.


### 5.2. Dark Soliton

Seeking the solution by tanh function method as in (7)(25)uξ=A1tanhβ1⁡μξ,vξ=A2tanhβ2⁡μξ,the system of equations in (11) and (15) becomes, respectively,(26)β1μ2β1-1tanhβ1-2⁡μξ-2β1tanhβ1⁡μξ+β1+1tanhβ1+2⁡μξ-ω+k2-α1tanhβ1⁡μξ-A2tanhβ1+β2⁡μξ=0,(27)12k2+α2A2tanhβ2⁡μξ-A2β2μ2β2-1tanhβ2-2⁡μξ-2β2tanhβ2⁡μξ+β2+1tanhβ2+2⁡μξ+3A22tanh2β2⁡μξ-A12tanh2β1⁡μξ=0.Equating the exponents and the coefficients of each pair of the sech functions, we find (28)2β1=β2+2,β1+β2=β1+2,thenβ1=β2=2.Thus setting coefficients of (26)-(27) to zero yields set system of equations:(29)8μ2+ω+k2-α1=0,12k2+α2+8μ2=0,6μ2-A2=0,-6A2μ2+3A22-A12=0.Solving the system of equations in (29), we get(30)ω=11k2+α1+α2,μ=∓i12k2+α28,(31)A1=±32412k2+α2,A2=-3412k2+α2,(32)E2x,t=±eikx+11k2+α1+α2t+ϵ03212k2+α24sec212k2+α28x-2kt+χ,(33)N2x,t=-312k2+α24sec212k2+α28x-2kt+χ.For k=α1=1,α2=4,ϵ0=χ=0, the real part of E2x,t=122cos⁡x+16tsec2{2x-2t}, and N2x,t=-12sec2{2x-2t}.Figures3 and 4 represent the solitary wave of the real part of E2(x,t) in (32) and N2(x,t) in (33) for -5≤x≤5, 0≤t≤1.Figure 3
The solitary wave of the real part ofE2(x,t) in (32) for -5≤x≤5, 0≤t≤1.Figure 4
The solitary waveN2(x,t) in (33) for -5≤x≤5, 0≤t≤1.


### 5.3. Singular Soliton

Seeking the solution by sech function method as in (8)(34)uξ=A1cothβ1⁡μξ,vξ=A2cothβ2⁡μξ,the system of equations in (11) and (15) becomes, respectively,(35)β1μ2β1+1cothβ1+2⁡μξ+β1cothβ1⁡μξ+ω+k2-α1cothβ1⁡μξ-A2cothβ1+β2⁡μξ=0,(36)12k2+α2A2cothβ2⁡μξ+A2β2μ2β2+1cothβ2+2⁡μξ+β2cothβ2⁡μξ+3A22coth2β2⁡μξ-A12coth2β1⁡μξ=0.Equating the exponents and the coefficients of each pair of the sech functions, we find (37)2β1=β2+2,β1+β2=β1+2,thenβ1=β2=2.Thus setting coefficients of (35)-(36) to zero yields set system of equations:(38)4μ2+ω+k2-α1=0,12k2+α2+4μ2=0,6μ2-A2=0,6A2μ2+3A22-A12=0.Solving the system of equations in (38), we get(39)A1=±312k2+α2,A2=-312k2+α22,μ=∓i12k2+α24,ω=11k2+α1+α2,(40)E3x,t=∓eikx+11k2+α1+α2t+ϵ066k2+1cot26k2+12x-2kt+χ,(41)N3x,t=36k2+1cot26k2+12x-2kt+χ.For k=α1=1,α2=4,ϵ0=χ=0, the real part of E3x,t=42cos⁡x+16tcot2{7/2x-2t}, and N3x,t=21cot2{7/2x-2t}.Figures5 and 6 represent the solitary wave of the real part of E3(x,t) in (40) and N3(x,t) in (41) for -5≤x≤5, 0≤t≤1.Figure 5
The solitary wave of the real part ofE3(x,t) in (40) for -5≤x≤5, 0≤t≤1.Figure 6
The solitary wave ofN3(x,t) in (41) for -5≤x≤5, 0≤t≤1.


### 5.4. Modified Simple Equation Method

This section will analyze (11) and (15) by the modified simple equation method; assume that solutions are of the form [23](42)uξ=∑i=0mAiψξψi,vξ=∑i=0nBiψξψi,where the parameters m, n can be found by balancing the highest-order linear term with the nonlinear terms in (11) and (15), respectively.In (11), we balance u′′ with uv, to obtain m+2=m+n, and then n=2. While in (15), We balance v′′ with u2, to obtain n+2=2m, and then m=2.Then(43)uξ=A0+A1ψξψ+A2ψξ2ψ2,uξ=A1ψξξψ-A1ψξ2ψ2+2A2ψξψξξψ2-2A2ψξ3ψ3,uξξ=A1ψξξξψ+2A2ψξξ2ψ2+2A2ψξξξψξψ2-3A1ψξψξξψ2-10A2ψξξψξ2ψ3-2A1ψξ3ψ3-2A2ψξ4ψ4,vξ=B0+B1ψξψ+B2ψξ2ψ2,vξ=B1ψξξψ-B1ψξ2ψ2+2B2ψξψξξψ2-2B2ψξ3ψ3,vξξ=B1ψξξξψ+2B2ψξξ2ψ2+2B2ψξξξψξψ2-3B1ψξψξξψ2-10B2ψξξψξ2ψ3-2B1ψξ3ψ3-2B2ψξ4ψ4,where A0,A1,A2,B0,B1, and B2 are constants to be calculated.Substitute (43) in (11) and (15), respectively, to get(44)A1ψξξξψ+2A2ψξξ2ψ2+2A2ψξξξψξψ2-3A1ψξψξξψ2-10A2ψξξψξ2ψ3-2A1ψξ3ψ3-2A2ψξ4ψ4-ω+k2-α1A0+A1ψξψ+A2ψξ2ψ2-B0A0+A1ψξψ+A2ψξ2ψ2-B1A0ψξψ+A1ψξ2ψ2+A2ψξ3ψ3-B2A0ψξ2ψ2+A1ψξ3ψ3+A2ψξ4ψ4=0,12k2+α2B0+B1ψξψ+B2ψξ2ψ2-B1ψξξξψ+2B2ψξξ2ψ2+2B2ψξξξψξψ2-3B1ψξψξξψ2-10B2ψξξψξ2ψ3-2B1ψξ3ψ3-2B2ψξ4ψ4+3B0B0+B1ψξψ+B2ψξ2ψ2+3B1B0ψξψ+B1ψξ2ψ2+B2ψξ3ψ3+3B2B0ψξ2ψ2+B1ψξ3ψ3+B2ψξ4ψ4-A0A0+A1ψξψ+A2ψξ2ψ2-A1A0ψξψ+A1ψξ2ψ2+A2ψξ3ψ3-A2A0ψξ2ψ2+A1ψξ3ψ3+A2ψξ4ψ4=0.In (44) equating expressions at ψjj=0,-1,-2,-3,-4 to zero, we get the following system of equations:(45)ψξξ+2A1+A2B1+A1B210A2ψξ=0,2A2ψξξ2+2A2ψξξξψξ-3A1ψξψξξ-ω+k2-α1A2+A2B0+A1B1+A0B2ψξ2=0,ψξξξ-ω+k2-α1A1-B0A1-B1A0A1ψξ=0,ω+k2-α1+B0=0,12k2+α2B0+3B02-A02=0,2B2+3B22-A22=0,B2=-2,ψξξ+2B1+6B1B2-2A1A210B2ψξ=0,2B2ψξξ2+2B2ψξξξψξ-3B1ψξψξξ-12k2+α2B2+6B0B2+3B12-2A0A2-A12ψξ2=0,B1ψξξξ-12k2+α2B1+6B1B0-2A0A1ψξ=0.Obviously when solving the system of (45), we conclude that equations can be satisfied simultaneously for the following constraints. Hence, the modified simple equation method does not produce the soliton solution in general case:(46)A0=0,B0=0,B1=A1=0,A2=22,B2=-2,ω=α1-k2,α2=-12k2.Then we will solve the following ordinary differential equation:(47)ψξξ=0,and therefore(48)ψ=a0+a1ξ,where a0,a1 are arbitrary constants.And(49)ux,t=22a12a0+a1x-2kt+χ2,vx,t=-2a12a0+a1x-2kt+χ2.Finally solutions become(50)E4x,t=22eikx+α1-k2t+ϵ0a12a0+a1x-2kt+χ2,N4x,t=-2a12a0+a1x-2kt+χ2.


## 5.1. Bright Soliton

Seeking the solution by sech function method as in (6)(16)uξ=A1sechβ1⁡μξ,vξ=A2sechβ2⁡μξ,the system of equations in (11) and (15) becomes, respectively,(17)β1μ2β1+1sechβ1+2⁡μξ-β1sechβ1⁡μξ+ω+k2-α1sechβ1⁡μξ-A2sechβ1+β2⁡μξ=0,(18)12k2+α2A2sechβ2⁡μξ+A2β2μ2β2+1sechβ2+2μξ-β2sechβ2⁡μξ+3A22sech2β2⁡μξ-A12sech2β1⁡μξ=0.Equating the exponents and the coefficients of each pair of the sech functions, we find (19)2β1=β2+2,β1+β2=β1+2,thenβ1=β2=2.Thus setting coefficients of (17)-(18) to zero yields set system of equations:(20)4μ2-ω+k2-α1=0,12k2+α2-4μ2=0,6μ2-A2=0,6A2μ2+3A22-A12=0.Solving the system of equations in (20), we get(21)A1=∓312k2+α2,A2=3212k2+α2,(22)μ=∓12k2+α24,ω=11k2+α1+α2,(23)E1x,t=∓3eikx+11k2+α1+α2t+ϵ012k2+α2sech2⁡12k2+α24x-2kt+χ,(24)N1x,t=3212k2+α2sech2⁡12k2+α24x-2kt+χ.For k=α1=1,α2=4,ϵ0=χ=0, the real part of E1x,t=48cos⁡x+16tsech2⁡{2x-2t}, and N1x,t=24sech2⁡{2x-2t}.Figures1 and 2 represent the solitary wave of the real part of E1(x,t) in (23) and N1(x,t) in (24), respectively, for -5≤x≤5, 0≤t≤1.Figure 1
The solitary wave of the real part ofE1(x,t) in (23) for -5≤x≤5, 0≤t≤1.Figure 2
The solitary waveN1(x,t) in (24) for -5≤x≤5, 0≤t≤1.


## 5.2. Dark Soliton

Seeking the solution by tanh function method as in (7)(25)uξ=A1tanhβ1⁡μξ,vξ=A2tanhβ2⁡μξ,the system of equations in (11) and (15) becomes, respectively,(26)β1μ2β1-1tanhβ1-2⁡μξ-2β1tanhβ1⁡μξ+β1+1tanhβ1+2⁡μξ-ω+k2-α1tanhβ1⁡μξ-A2tanhβ1+β2⁡μξ=0,(27)12k2+α2A2tanhβ2⁡μξ-A2β2μ2β2-1tanhβ2-2⁡μξ-2β2tanhβ2⁡μξ+β2+1tanhβ2+2⁡μξ+3A22tanh2β2⁡μξ-A12tanh2β1⁡μξ=0.Equating the exponents and the coefficients of each pair of the sech functions, we find (28)2β1=β2+2,β1+β2=β1+2,thenβ1=β2=2.Thus setting coefficients of (26)-(27) to zero yields set system of equations:(29)8μ2+ω+k2-α1=0,12k2+α2+8μ2=0,6μ2-A2=0,-6A2μ2+3A22-A12=0.Solving the system of equations in (29), we get(30)ω=11k2+α1+α2,μ=∓i12k2+α28,(31)A1=±32412k2+α2,A2=-3412k2+α2,(32)E2x,t=±eikx+11k2+α1+α2t+ϵ03212k2+α24sec212k2+α28x-2kt+χ,(33)N2x,t=-312k2+α24sec212k2+α28x-2kt+χ.For k=α1=1,α2=4,ϵ0=χ=0, the real part of E2x,t=122cos⁡x+16tsec2{2x-2t}, and N2x,t=-12sec2{2x-2t}.Figures3 and 4 represent the solitary wave of the real part of E2(x,t) in (32) and N2(x,t) in (33) for -5≤x≤5, 0≤t≤1.Figure 3
The solitary wave of the real part ofE2(x,t) in (32) for -5≤x≤5, 0≤t≤1.Figure 4
The solitary waveN2(x,t) in (33) for -5≤x≤5, 0≤t≤1.


## 5.3. Singular Soliton

Seeking the solution by sech function method as in (8)(34)uξ=A1cothβ1⁡μξ,vξ=A2cothβ2⁡μξ,the system of equations in (11) and (15) becomes, respectively,(35)β1μ2β1+1cothβ1+2⁡μξ+β1cothβ1⁡μξ+ω+k2-α1cothβ1⁡μξ-A2cothβ1+β2⁡μξ=0,(36)12k2+α2A2cothβ2⁡μξ+A2β2μ2β2+1cothβ2+2⁡μξ+β2cothβ2⁡μξ+3A22coth2β2⁡μξ-A12coth2β1⁡μξ=0.Equating the exponents and the coefficients of each pair of the sech functions, we find (37)2β1=β2+2,β1+β2=β1+2,thenβ1=β2=2.Thus setting coefficients of (35)-(36) to zero yields set system of equations:(38)4μ2+ω+k2-α1=0,12k2+α2+4μ2=0,6μ2-A2=0,6A2μ2+3A22-A12=0.Solving the system of equations in (38), we get(39)A1=±312k2+α2,A2=-312k2+α22,μ=∓i12k2+α24,ω=11k2+α1+α2,(40)E3x,t=∓eikx+11k2+α1+α2t+ϵ066k2+1cot26k2+12x-2kt+χ,(41)N3x,t=36k2+1cot26k2+12x-2kt+χ.For k=α1=1,α2=4,ϵ0=χ=0, the real part of E3x,t=42cos⁡x+16tcot2{7/2x-2t}, and N3x,t=21cot2{7/2x-2t}.Figures5 and 6 represent the solitary wave of the real part of E3(x,t) in (40) and N3(x,t) in (41) for -5≤x≤5, 0≤t≤1.Figure 5
The solitary wave of the real part ofE3(x,t) in (40) for -5≤x≤5, 0≤t≤1.Figure 6
The solitary wave ofN3(x,t) in (41) for -5≤x≤5, 0≤t≤1.


## 5.4. Modified Simple Equation Method

This section will analyze (11) and (15) by the modified simple equation method; assume that solutions are of the form [23](42)uξ=∑i=0mAiψξψi,vξ=∑i=0nBiψξψi,where the parameters m, n can be found by balancing the highest-order linear term with the nonlinear terms in (11) and (15), respectively.In (11), we balance u′′ with uv, to obtain m+2=m+n, and then n=2. While in (15), We balance v′′ with u2, to obtain n+2=2m, and then m=2.Then(43)uξ=A0+A1ψξψ+A2ψξ2ψ2,uξ=A1ψξξψ-A1ψξ2ψ2+2A2ψξψξξψ2-2A2ψξ3ψ3,uξξ=A1ψξξξψ+2A2ψξξ2ψ2+2A2ψξξξψξψ2-3A1ψξψξξψ2-10A2ψξξψξ2ψ3-2A1ψξ3ψ3-2A2ψξ4ψ4,vξ=B0+B1ψξψ+B2ψξ2ψ2,vξ=B1ψξξψ-B1ψξ2ψ2+2B2ψξψξξψ2-2B2ψξ3ψ3,vξξ=B1ψξξξψ+2B2ψξξ2ψ2+2B2ψξξξψξψ2-3B1ψξψξξψ2-10B2ψξξψξ2ψ3-2B1ψξ3ψ3-2B2ψξ4ψ4,where A0,A1,A2,B0,B1, and B2 are constants to be calculated.Substitute (43) in (11) and (15), respectively, to get(44)A1ψξξξψ+2A2ψξξ2ψ2+2A2ψξξξψξψ2-3A1ψξψξξψ2-10A2ψξξψξ2ψ3-2A1ψξ3ψ3-2A2ψξ4ψ4-ω+k2-α1A0+A1ψξψ+A2ψξ2ψ2-B0A0+A1ψξψ+A2ψξ2ψ2-B1A0ψξψ+A1ψξ2ψ2+A2ψξ3ψ3-B2A0ψξ2ψ2+A1ψξ3ψ3+A2ψξ4ψ4=0,12k2+α2B0+B1ψξψ+B2ψξ2ψ2-B1ψξξξψ+2B2ψξξ2ψ2+2B2ψξξξψξψ2-3B1ψξψξξψ2-10B2ψξξψξ2ψ3-2B1ψξ3ψ3-2B2ψξ4ψ4+3B0B0+B1ψξψ+B2ψξ2ψ2+3B1B0ψξψ+B1ψξ2ψ2+B2ψξ3ψ3+3B2B0ψξ2ψ2+B1ψξ3ψ3+B2ψξ4ψ4-A0A0+A1ψξψ+A2ψξ2ψ2-A1A0ψξψ+A1ψξ2ψ2+A2ψξ3ψ3-A2A0ψξ2ψ2+A1ψξ3ψ3+A2ψξ4ψ4=0.In (44) equating expressions at ψjj=0,-1,-2,-3,-4 to zero, we get the following system of equations:(45)ψξξ+2A1+A2B1+A1B210A2ψξ=0,2A2ψξξ2+2A2ψξξξψξ-3A1ψξψξξ-ω+k2-α1A2+A2B0+A1B1+A0B2ψξ2=0,ψξξξ-ω+k2-α1A1-B0A1-B1A0A1ψξ=0,ω+k2-α1+B0=0,12k2+α2B0+3B02-A02=0,2B2+3B22-A22=0,B2=-2,ψξξ+2B1+6B1B2-2A1A210B2ψξ=0,2B2ψξξ2+2B2ψξξξψξ-3B1ψξψξξ-12k2+α2B2+6B0B2+3B12-2A0A2-A12ψξ2=0,B1ψξξξ-12k2+α2B1+6B1B0-2A0A1ψξ=0.Obviously when solving the system of (45), we conclude that equations can be satisfied simultaneously for the following constraints. Hence, the modified simple equation method does not produce the soliton solution in general case:(46)A0=0,B0=0,B1=A1=0,A2=22,B2=-2,ω=α1-k2,α2=-12k2.Then we will solve the following ordinary differential equation:(47)ψξξ=0,and therefore(48)ψ=a0+a1ξ,where a0,a1 are arbitrary constants.And(49)ux,t=22a12a0+a1x-2kt+χ2,vx,t=-2a12a0+a1x-2kt+χ2.Finally solutions become(50)E4x,t=22eikx+α1-k2t+ϵ0a12a0+a1x-2kt+χ2,N4x,t=-2a12a0+a1x-2kt+χ2.


## 6. Conclusion

In this paper the dispersive bright, dark, and singular soliton solutions to SBE with Kerr law of nonlinearity were studied. The sech, tanh, csch, and the modified simplest equation method have been successfully applied to find solitons solutions for the coupled Schrödinger-Boussinesq equations. Several constraint conditions were assuring the existence of such solitons with Kerr law nonlinearity. The modified simple equation method does not produce the soliton solution in general case. Solutions by three methods are plotted in figures for the real and imaginary parts forE(x,t) and N(x,t). Compatibility in figures shape between the solutions of E(x,t) and N(x,t) by the same method sometimes appeared. Solutions may be important for the conservation laws for dispersive optical solitons. Those research outcomes will be soon disseminated.



---

*Source: 8325919-2018-01-01.xml*