On the soliton solutions of the generalized stochastic nonlinear Schrödinger equation with Kerr effect and higher order nonlinearity via two analytical methods

Ali Khalid Salih; Hajar Farhan Ismael

doi:10.2478/ijmce-2026-0007

Introduction

Nonlinear partial differential equations (NLPDEs) are extensively utilized across multiple fields, including mathematical biology, chemistry, engineering, plasma physics, quantum mechanics, and fluid dynamics [1]. In order to explore and explain nonlinear effects, researchers in mathematics and physics focus on obtaining exact solutions of the partial differential equations that govern such systems [2, 3]. The study of complex nonlinear partial differential equations has garnered significant academic interest, driven by the inherent challenges of solving them and their importance in understanding intricate natural phenomena [4]. Significant advancements have been observed in the development of effective methods for obtaining accurate solutions to NLPDEs in recent years, including the Hirota direct method [5, 6], the new extended auxiliary equation method [7], the new extended direct algebraic method [8], the sin-Gordon expansion method [9], the modified simple equation method [10], the sub-ODE method [11], the F-expansion method [12], the (m+G’/G) expansion method [13], the improved (φ(ξ)/2)-expansion method [14], the Bernoulli sub-equation function method [15], the Jacobi elliptic function method [16, 17], the unified Riccati equation expansion method [18, 19], the modified Sardar sub-equation approach [20, 21], and many more.

A stochastic process mathematically represents the evolution of a random phenomenon over time. Stochastic differential equations (SDEs) extend this idea by integrating random fluctuations into system dynamics. Similarly, stochastic partial differential equations (SPDEs) incorporate random variables or noise functions, providing suitable mathematical representations for complex systems influenced by uncertainty. In nonlinear optics, optical solitons, widely used in high-speed data transmission, are subject to stochastic disturbances from external influences. This requires the formulation of differential equations with stochastic elements for accurate modeling. The stochastic nonlinear Schrödinger equation (SNLSE) provides a mathematical framework for analyzing such systems, especially in quantum mechanics and nonlinear optics [22]. It describes the evolution of a quantum field that incorporates both nonlinear interactions and stochastic effects. Beyond optics, stochastic networked linear systems are also employed to capture random influences across diverse fields such as chemistry, physics, and electrical engineering. Increasingly, the SNLSE has been applied in areas including physics modeling, climate research, and information technology, where it serves as a powerful tool for constructing mathematical models of complex phenomena.

The governing model under consideration is the stochastic generalized nonlinear Schrödinger equation (SGNLSE), as described in [23]:

i q_{t} + q_{x x} + | q |^{2} q + i σ (q_{x x x} + γ_{1} | q |^{2} q_{x} + γ_{2} q^{2} q_{x}^{*}) + σ \frac{\partial W}{\partial t} q = 0 .

iq_{t}+q_{xx}+|q|^{2}q+i\sigma(q_{xxx}+\gamma_{1}|q|^{2}q_{x}+\gamma_{2}q^{2}q_{x}^{*})+\sigma\frac{\partial W}{\partial t}q=0.

Here, Eq.(1) describes the propagation of light in a medium by incorporating the Kerr effect, which causes intensity-dependent changes in the refractive index. This nonlinearity affects both the phase and intensity of the light, thereby altering the energy distribution of the wave. The generalized nonlinear Schrödinger equation (GNLSE) is often used to predict energy changes, such as gain or loss, that arise from these nonlinear interactions. As a wave travels through a medium, it may experience energy gain similar to amplification in electronic systems or energy loss due to mechanisms such as absorption, scattering, and transmission. A filter is a device that selectively permits certain frequency components of a wave to pass while blocking others, thereby manipulating the energy distribution of the wave and often causing a net change in its energy. In this framework, q = q(x,t) denotes the soliton pulse profile, with q_t representing temporal evolution, q_xxcorresponding to group velocity dispersion (GVD), and q_xxx accounting for third-order dispersion. The term ${| q |}^{2} q$ |q|^{2}q captures the Kerr effect, while q with an asterisk (q_*) denotes the complex conjugate of q. The expressions ${| q |}^{2} q_{x}$ |q|^{2}q_{x} and $q^{2} q_{x}^{*}$ q^{2}q_{x}^{*} represent higher-order nonlinear terms, with γ₁ and γ₂ as real-valued physical parameters. The term W(t) denotes a one-dimensional Wiener process, which introduces stochasticity into the system and is formally defined through the Itô integral as [24]:

W (t) = \int_{0}^{t} q (τ) d W (τ), where τ < t .

W(t)=\int_{0}^{t}q(\tau)dW(\tau),

The aim of this paper is to employ the $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2}) expansion and the (Gߣ/G,1/G)-expansion technique to Eq. (1) to find exact solutions and apply novel analytical methods to obtain exact soliton solutions. The $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2})-expansion technique, presented as an alternative and effective analytical method, has demonstrated its efficacy as a powerful instrument for deriving accurate solutions to NLPDEs. This method employs the hyperbolic tangent function to convert nonlinear equations into algebraic equations. For instance, the modified Zakharov-Kuznetsov equation [25] and Schrödinger type nonlinear evolution equations [26]. The two-variable (Gߣ/G,1/G)-expansion method is a recently developed analytical technique for constructing exact solutions of nonlinear evolution equations (NLEEs). This method involves introducing two independent auxiliary functions based on the ratio of derivatives and reciprocals of a function G, which satisfies a second-order linear ordinary differential equation. By converting complex NLPDEs into algebraic forms, the method facilitates the derivation of exact traveling wave solutions. It was first applied by Li et al. [27] to the Zakharov equations and has since been employed by several researchers to solve a wide range of NLEEs, including higher-order nonlinear Schrödinger and quantum Zakharov-Kuznetsov equations [28, 29].

In addition to classical expansion approaches, several modern analytical techniques have been developed for constructing exact solutions of NLPDEs. Kudryashov proposed effective methods for obtaining exact solitary solutions to the generalized Kuramoto-Sivashinsky equation, providing an important basis for further analytical investigations [30, 31]. Later, Kudryashov et al. [32] applied Painlevé analysis and the first integral method to the traveling wave reductions of a nonlinear equation, namely, the Radhakrishnan-Kundu-Lakshmanan equation. Demina and Kudryashov [33] presented a general method for constructing explicit meromorphic solutions of autonomous nonlinear ordinary differential equations based on the analysis of Laurent series and singularity structures. In a related but more specific study, they applied this approach to derive exact meromorphic solutions of the Kawahara equation, demonstrating its effectiveness for higher-order dispersive wave equations [34]. In another contribution, Kudryashov presented a robust approach for finding exact solutions to nonlinear differential equations, emphasizing its generality and efficiency [35]. Our study incorporates the $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2}) -expansion and ( Gߣ/G,1/G)-expansion methods while exploring their potential for generating physically meaningful solutions for the generalized stochastic nonlinear Schrödinger equation with Kerr effect and higher-order nonlinearity.

The rest of this paper is organized as follows. In Section 2, the core procedures of the $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2}) -expansion method and the (G’/G,1/G)-expansion method are outlined. Section 3 focuses on the application of both techniques to Eq.(1), followed by a detailed analysis of the resulting exact analytical solutions. Section 4 presents the stability analysis of the obtained soliton solution. Section 5 includes 2D and 3D visualizations and examines how the solution evolves over time. The results are discussed in Section 6. Section 7 concludes the paper with a summary of the main findings.

Description of the methods

2.1

The

t a n (\frac{ϖ (ξ)}{2})

tan(\frac{\varpi(\xi)}{2})

-expansion method

Suppose a NPDE follows the following form:

P (q, q_{x}, q_{t}, q_{t x}, \dots) = 0

P(q,q_{x},q_{t},q_{tx},\cdot\cdot\cdot)=0

where $q = q (x, t)$ q=q(x,t) and it can be transformed into an ordinary differential equation (ODE):

O (Q, Q^{'}, Q^{'}, Q^{″}, \dots) = 0 .

O(Q,Q^{\prime},Q^{\prime\prime},Q^{\prime\prime\prime},\cdot\cdot\cdot)=0.

By using the wave transformation $q = Q (ξ)$ q=Q(\xi) and $ξ = x - v t$ \xi=x-vt. The following expression is proposed as a traveling wave solution to Eq.(3):

Q (ξ) = W (ϖ) = \sum_{k = 0}^{n} a_{k} {(t a n (\frac{ϖ (ξ)}{2}))}^{k}, 0 \leq k \leq n,

Q(\xi)=W(\varpi)=\sum_{k=0}^{n}a_{k}(\frac{tan(\varpi(\xi))}{2})^{k},0\le k\le n.

with constants a_k≠0, which will be evaluated subsequently and must satisfy the governing ODE.

ϖ^{'} (ξ) = A s i n (ϖ (ξ)) + B c o s (ϖ (ξ)) + C .

\varpi^{\prime}(\xi)=A~sin(\varpi(\xi))+B~cos(\varpi(\xi))+C.

Here are a class of solutions to Eq.(5):

Family 1: When $A^{2} + B^{2} - C^{2} < 0$ A^{2}+B^{2}-C^{2}<0 and $B - C \neq 0$ B-C\ne0 the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{A}{B - C} - (\frac{\sqrt{- A^{2} - B^{2} + C^{2}}}{B - C}) t a n (\frac{\sqrt{- A^{2} - B^{2} + C^{2}}}{2} ξ)) .

\varpi(\xi)=2~tan^{-1}(\frac{A}{B-C}-(\frac{\sqrt{-A^{2}-B^{2}+C^{2}}}{B-C})tan(\frac{\sqrt{-A^{2}-B^{2}+C^{2}}}{2}\xi))

Family 2: In case $A^{2} + B^{2} - C^{2} > 0$ A^{2}+B^{2}-C^{2}>0 and $B - C \neq 0 .$ B-C\ne0, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{A}{B - C} + (\frac{\sqrt{A^{2} + B^{2} - C^{2}}}{B - C}) t a n h (\frac{\sqrt{A^{2} + B^{2} - C^{2}}}{2} ξ)) .

\varpi(\xi)=2~tan^{-1}(\frac{A}{B-C}+(\frac{\sqrt{A^{2}+B^{2}-C^{2}}}{B-C})tanh(\frac{\sqrt{A^{2}+B^{2}-C^{2}}}{2}\xi))

Family 3: In case A²+B²-C²>0, B= ≠ 0 and C = 0, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{A}{B - C} + \frac{\sqrt{B^{2} - A^{2}}}{B} t a n h (\frac{\sqrt{B^{2} - A^{2}}}{2} ξ)) .

\varpi(\xi)=2~tan^{-1}(\frac{A}{B-C}+\frac{\sqrt{B^{2}-A^{2}}}{B}tanh(\frac{\sqrt{B^{2}-A^{2}}}{2}\xi))

Family 4: In case A²+B²-C²<0, C ≠ 0 and B=0, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (- \frac{A}{C} + \frac{\sqrt{C^{2} - A^{2}}}{C} t a n (\frac{\sqrt{C^{2} - A^{2}}}{2} ξ)) .

\varpi(\xi)=2~tan^{-1}(-\frac{A}{C}+\frac{\sqrt{C^{2}-A^{2}}}{C}tan(\frac{\sqrt{C^{2}-A^{2}}}{2}\xi))

Family 5: In case A²+B²–C²>0, B –C ≠ 0 and A=0, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\sqrt{\frac{B + C}{B - C}} t a n h (\frac{\sqrt{B^{2} - C^{2}}}{2} ξ)) .

\varpi(\xi)=2~tan^{-1}(\sqrt{\frac{B+C}{B-C}}tanh(\frac{\sqrt{B^{2}-C^{2}}}{2}\xi))

Family 6: In case A = 0 and C = 0, the corresponding solution is given by

ϖ (ξ) = t a n^{- 1} (\frac{2 e^{B ξ}}{e^{2 B ξ} + 1}) .

\varpi(\xi)=tan^{-1}(\frac{2e^{B\xi}}{e^{2B\xi}+1})

Family 7: In case B = C = 0, the corresponding solution is given by

ϖ (ξ) = t a n^{- 1} (\frac{2 e^{A ξ}}{e^{2 A ξ} + 1}) .

\varpi(\xi)=tan^{-1}(\frac{2e^{A\xi}}{e^{2A\xi}+1})

Family 8: In case A²+B²=C², then the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{A ξ + 2}{(B - C) ξ}) . .

\varpi(\xi)=2~tan^{-1}(\frac{A\xi+2}{(B-C)\xi})

Family 9: In case A = B = C = kA, then the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (e^{k A ξ} - 1) .

\varpi(\xi)=2~tan^{-1}(e^{kA\xi}-1)

Family 10: In case A = C = kA and B = –kA, the corresponding solution is given by

ϖ (ξ) = - 2 t a n^{- 1} (\frac{e^{k A ξ}}{- 1 + e^{k A ξ}}) .

\varpi(\xi)=-2~tan^{-1}(\frac{e^{kA\xi}}{-1+e^{kA\xi}})

Family 11: In case C = A, the corresponding solution is given by

ϖ (ξ) = - 2 t a n^{- 1} (\frac{(A + B) e^{B ξ} - 1}{(A - B) e^{B ξ} - 1}) .

\varpi(\xi)=-2~tan^{-1}(\frac{(A+B)e^{B\xi}-1}{(A-B)e^{B\xi}-1})

Family 12: In case A = C, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{(B + A) e^{B ξ} + 1}{(B - A) e^{B ξ} - 1}) .

\varpi(\xi)=2~tan^{-1}(\frac{(B+A)e^{B\xi}+1}{(B-A)e^{B\xi}-1})

Family 13: In case C = –A, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{B - A + e^{B ξ}}{- B - A + e^{B ξ}}) .

\varpi(\xi)=2~tan^{-1}(\frac{B-A+e^{B\xi}}{-B-A+e^{B\xi}})

Family 14: In case B = –C, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{A e^{A ξ}}{1 - C e^{A ξ}}) .

\varpi(\xi)=2~tan^{-1}(\frac{Ae^{A\xi}}{1-Ce^{A\xi}})

Family 15: In case B = 0 and A = C, the corresponding solution is given by

ϖ (ξ) = - 2 t a n^{- 1} (\frac{C ξ + 2}{C ξ}) .

\varpi(\xi)=-2~tan^{-1}(\frac{C\xi+2}{C\xi})

Family 16: In case A = 0 and B = C, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (C ξ) .

\varpi(\xi)=2~tan^{-1}(C\xi)

Family 17: In case A = 0 while B = –C, the corresponding solution is given by

ϖ (ξ) = - 2 t a n^{- 1} (\frac{1}{C ξ})

\varpi(\xi)=-2~tan^{-1}(\frac{1}{C\xi})

Family 18: In case A = B = 0, the corresponding solution is given by

ϖ (ξ) = C ξ + ϵ .

\varpi(\xi)=C\xi+\epsilon

Family 19: In case B = C, the corresponding solution is given by

ϖ (ξ) = 2 t a n^{- 1} (\frac{e^{A ξ} - C}{A}) .

\varpi(\xi)=2~tan^{-1}(\frac{e^{A\xi}-C}{A})

The value of n is determined using the balance principle. By substituting Eq.(4) into Eq.(3) and collecting the coefficients of $(t a n (ϖ / 2))^{n}$ (tan(\varpi/2))^{n} for (n=1,2,3,...), a system of equations is obtained by setting each coefficient to zero. This results in an over-determined set of equations for A, B, C, and a_k (k = 1,2,...,n) Symbolic computations were carried out using the computational package program.

2.2

The (G’/G,1/G)- expansion method

This part focuses on elucidating fundamental ideas of mention technique, which is employed to determine precise solutions of non-linear equations. Let us examine second order of linear ODE.

G^{″} [ξ] + δ G [ξ] = ρ,

G^{\prime\prime}[\xi]+\delta G[\xi]=\rho.

set as,

Ψ = 1 / G, and ϕ = G^{'} / G .

\Psi=1/G,\text{ and }\phi=G^{\prime}/G.

It may be inferred from the aforementioned relationships,

ϕ^{'} = - ϕ^{2} + ρ Ψ - δ, and Ψ^{'} = - ϕ Ψ .

\phi^{\prime}=-\phi^{2}+\rho\psi-\delta,\text{ and }\psi^{\prime}=-\phi\psi.

Three types of solutions may be discovered for various values of δ:

Case 1: When δ < 0, the general solution corresponding to Eq.(6) is given by the following expression:

G (ξ) = A_{1} c o s h (\sqrt{- δ} ξ) + B_{2} s i n h (\sqrt{- δ} ξ) + \frac{ρ}{δ} .

G(\xi)=A_{1}cosh(\sqrt{-\delta}\xi)+B_{2}sinh(\sqrt{-\delta}\xi)+\frac{\rho}{\delta}

Therefore, we acquire

Ψ^{2} = \frac{- δ (ϕ^{2} - 2 ρ Ψ + δ)}{δ^{2} ({A_{1}}^{2} - {A_{2}}^{2}) + ρ^{2}} .

\psi^{2}=\frac{-\delta(\phi^{2}-2\rho~\psi+\delta)}{\delta^{2}({A_{1}}^{2}-{A_{2}}^{2})+\rho^{2}}.

Here A₁,A₂ are constants.

Case 2: When δ < 0, the general solution corresponding to Eq.(6) is given by the following expression:

G (ξ) = A_{1} c o s (\sqrt{δ} ξ) + A_{2} s i n (\sqrt{δ} ξ) + \frac{ρ}{δ} .

G(\xi)=A_{1}cos(\sqrt{\delta}\xi)+A_{2}sin(\sqrt{\delta}\xi)+\frac{\rho}{\delta}

and consequently, we obtain

Ψ^{2} = \frac{δ (ϕ^{2} - 2 ρ Ψ + δ)}{δ^{2} ({A_{1}}^{2} + {A_{2}}^{2}) - ρ^{2}} .

\psi^{2}=\frac{\delta(\phi^{2}-2\rho~\psi+\delta)}{\delta^{2}({A_{1}}^{2}+{A_{2}}^{2})-\rho^{2}}

Here A₁,A₂ are constants.

Case 3: When δ = 0, the general solution corresponding to Eq.(6) is given by the following expression:

G (ξ) = \frac{ρ}{2} ξ^{2} + A_{1} ξ + A_{2},

G(\xi)=\frac{\rho}{2}\xi^{2}+A_{1}\xi+A_{2},

and thus, we obtain

Ψ^{2} = \frac{(ϕ^{2} - 2 ρ Ψ)}{A_{1}^{2} - 2 ρ A_{2}} .

\psi^{2}=\frac{(\phi^{2}-2\rho~\psi)}{A_{1}^{2}-2\rho A_{2}}.

Here A₁,A₂ are constants.

The main steps of the (G’/G,1/G)-expansion method are as follows:

Step 1. Apply the coordinate transformation ξ = x – vt and assume q(x,t) = Q(ξ). Under this transformation, Eq.(2) is reduced to ODE in terms of Q(ξ), given by:

P (Q, - v Q^{'}, Q^{'}, v^{2} Q^{″}, - v Q^{″}, Q^{″}, \dots) = 0 .

P(Q,-vQ^{\prime},Q^{\prime},v^{2}Q^{\prime\prime},-vQ^{\prime\prime},Q^{\prime\prime},\cdot\cdot\cdot)=0.

Step 2. Assume that the solution to ODE Eq.(15) can be represented as a polynomial involving the functions ϕ and Ψ in the following form:

Q (ξ) = \sum_{i = 0}^{n} a_{i} ϕ^{i} + \sum_{i = 1}^{n} B_{i} ϕ^{i - 1} Ψ .

Q(\xi)=\sum_{i=0}^{n}a_{i}\phi^{i}+\sum_{i=1}^{n}B_{i}\phi^{i-1}\psi .

Here, the constants a_i(i=0,,...,n),B_i(i=1,,...,n), δ, and ρ are to be determined later. The positive integer n can be found by applying the homogeneous balance between the highest-order derivatives and the nonlinear terms present in ODE Eq. (15).

Step 3. By substituting Eq.(6) into Eq.(15) and applying Eq.(8) and Eq.(10), the left-hand side of Eq.(15) is transformed into a polynomial in ϕ and Ψ. Setting each coefficient of this polynomial to zero produces a system of algebraic equations for the constants a_i(i=0,,...,n),b_i(i=1,,...,n), δ (with δ < 0), ρ, A₁, and A₂.

Step 4. Using the computational program, the system of algebraic equations obtained in Step 3 is solved. By substituting the computed values of a_i(fori=0,,...,n),b_i(fori=1,,...,n) together with the parameters δ, ρ, A₁, and A₂ into Eq.(16), the travelling wave solutions are derived. These solutions are expressed in terms of hyperbolic functions, corresponding to the form presented in Eq.(15).

Step 5. By applying a similar procedure as in Steps 3 and 4, Eq.(16) is substituted into Eq.(15). Then, by making use of equations (8) and (12) (or alternatively, equations (8) and (14)), the travelling wave solutions of equation (15) are obtained. These solutions are expressed either in terms of trigonometric functions or as rational functions, depending on the specific case.

Applications

3.1

Application of

t a n (\frac{ϖ (ξ)}{2})

tan(\frac{\sigma(\xi)}{2})

-expansion method

This part of the study employs the subsequent transformation to get analytical solutions for Eq.(1)

q (x, t) = Q (ξ) e^{i [- κ x + ω t + σ W (t) - σ^{2} t]}, ξ = x - v t,

q(x,t)=Q(\xi)e^{i[-\kappa x+\omega t+\sigma W(t)-\sigma^{2}t]}, \xi=x-vt,

where κ,ω,σ and v are constants. By substituting Eq.(17) into Eq.(1), the real and imaginary parts are obtained as

(1 + k γ_{1} σ - k γ_{2} σ) Q^{3} - (k^{2} + k^{3} σ - σ^{2} + ω) Q + (1 + 3 k σ) Q^{″} = 0,

(1+k\gamma_{1}\sigma-k\gamma_{2}\sigma)Q^{3}-(k^{2}+k^{3}\sigma-\sigma^{2}+\omega)Q+(1+3k\sigma)Q^{\prime\prime}=0,

and

(- 2 k - v - 3 k^{2} σ) Q^{'} + (γ_{1} + γ_{2}) σ Q^{'} Q^{2} + σ Q^{‴} = 0 .

(-2k-v-3k^{2}\sigma)Q^{\prime}+(\gamma_{1}+\gamma_{2})\sigma Q^{\prime}Q^{2}+\sigma Q^{\prime\prime\prime}=0.

Integrating Eq.(19), and assuming the constant of integration C₀=0 for simplicity, we obtain

(- 2 k - v - 3 k^{2} σ) Q + \frac{1}{3} (γ_{1} + γ_{2}) σ Q^{3} + σ Q^{″} = 0 .

(-2k-v-3k^{2}\sigma)Q+\frac{1}{3}(\gamma_{1}+\gamma_{2})\sigma Q^{3}+\sigma Q^{\prime\prime}=0.

Balancing Q³ with Qⁿ obtaining n = 1 is feasible. Utilizing the derived value of n, the solution for Eq.(18) and Eq.(20) may be expressed as

Q (ξ) = a_{0} + a_{1} t a n (\frac{ϖ (ξ)}{2}),

Q(\xi)=a_{0}+a_{1}tan(\frac{\varpi(\xi)}{2})

a system of equations is generated by substituting Eq.(21) into equations (18) and (20). The solutions of this system yield the subsequent set of solutions as following

Set 1: When

\begin{array}{l} a_{0} = \frac{A \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}}, a_{1} = \frac{(- B + C) \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}}, γ_{2} = \frac{3 - γ_{1}}{1 + 6 k σ}, v = - 2 k - \frac{1}{2} (A^{2} + B^{2} - C^{2} + 6 k^{2}) σ, \\ ω = \frac{1}{2} (C^{2} - 2 k^{2} + σ (3 C^{2} k - 2 k^{3} + 2 σ) - A^{2} (1 + 3 k σ) - B^{2} (1 + 3 k σ)), \end{array}

a_{0}=\frac{A\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}},a_{1}=\frac{(-B+C)\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}},\gamma_{2}=\frac{3-\gamma_{1}}{1+6k\sigma},v=-2k-\frac{1}{2}(A^{2}+B^{2}-C^{2}+6k^{2})\sigma,\omega=\frac{1}{2}(C^{2}-2k^{2}+\sigma(3C^{2}k-2k^{3}+2\sigma)-A^{2}(1+3k\sigma)-B^{2}(1+3k\sigma)),

where $- 1 - 6 k σ > 0, 2 + 4 k γ_{1} σ > 0$ -1-6k\sigma>0,2+4k\gamma_{1}\sigma>0, the outcomes solutions are:

Solution 1: When –A²+B²–C²>0, the singular periodic solution corresponding to Family 1 is obtained

\begin{matrix} q (x, t) = & (\frac{\sqrt{(- 1 - 6 k σ) (- A^{2} - B^{2} + C^{2})}}{\sqrt{2 + 4 k γ_{1} σ}}) t a n (\frac{1}{4} \sqrt{- A^{2} - B^{2} + C^{2}} (4 k t + 2 x + (A^{2} + B^{2} - C^{2}) t σ + 6 k^{2} t σ)) \\ \times e^{- \frac{1}{2} i (B^{2} t - C^{2} t + 2 k^{2} t + 2 k x + 3 B^{2} k t σ - 3 C^{2} k t σ + 2 k^{3} t σ + A^{2} (t + 3 k t σ) - 2 σ W (t))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{\sqrt{(-1-6k\sigma)(-A^{2}-B^{2}+C^{2})}}{\sqrt{2+4k\gamma_{1}\sigma}})tan(\frac{1}{4}\sqrt{-A^{2}-B^{2}+C^{2}}(4kt+2x+(A^{2}+B^{2}-C^{2})t\sigma+6k^{2}t\sigma))\\ &\times e^{-\frac{1}{2}i(B^{2}t-C^{2}t+2k^{2}t+2kx+3B^{2}kt\sigma-3C^{2}kt\sigma+2k^{3}t\sigma+A^{2}(t+3kt\sigma)-2\sigma W(t))},\end{aligned}

Solution 2: When A²+B²–C²>0 the dark soliton solution corresponding to Family 2 is obtained

\begin{matrix} q (x, t) = & (- \frac{\sqrt{(- 1 - 6 k σ) (A^{2} + B^{2} - C^{2})}}{\sqrt{2 + 4 k γ_{1} σ}}) t a n h (\frac{1}{4} \sqrt{A^{2} + B^{2} - C^{2}} (4 k t + 2 x + (A^{2} + B^{2} - C^{2}) t σ + 6 k^{2} t σ)) \\ \times e^{- \frac{1}{2} i (- C^{2} t + 2 k^{2} t + 2 k x - 3 C^{2} k t σ + 2 k^{3} t σ - 2 W (t) σ + A^{2} (t + 3 k t σ) + B^{2} (t + 3 k t σ))} . \end{matrix}

\begin{aligned}q(x,t)=&(-\frac{\sqrt{(-1-6k\sigma)(A^{2}+B^{2}-C^{2})}}{\sqrt{2+4k\gamma_{1}\sigma}})tanh(\frac{1}{4}\sqrt{A^{2}+B^{2}-C^{2}}(4kt+2x+(A^{2}+B^{2}-C^{2})t\sigma+6k^{2}t\sigma))\\ &\times e^{-\frac{1}{2}i(-C^{2}t+2k^{2}t+2kx-3C^{2}kt\sigma+2k^{3}t\sigma-2W(t)\sigma+A^{2}(t+3kt\sigma)+B^{2}(t+3kt\sigma))}\end{aligned}

Solution 3: The dark soliton solution corresponding to Family 3 is obtained

\begin{matrix} q (x, t) = & (- \frac{\sqrt{A^{2} + B^{2}} \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}}) t a n h (\frac{1}{4} \sqrt{A^{2} + B^{2}} (4 k t + 2 x + (A^{2} + B^{2}) t σ + 6 k^{2} t σ)) \\ \times e^{- \frac{1}{2} i (A^{2} (t + 3 k t σ) + B^{2} (t + 3 k t σ) + 2 (k^{2} t + k x + k^{3} t σ - W (t) σ))} . \end{matrix}

\begin{aligned}q(x,t)=&(-\frac{\sqrt{A^{2}+B^{2}}\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}})tanh(\frac{1}{4}\sqrt{A^{2}+B^{2}}(4kt+2x+(A^{2}+B^{2})t\sigma+6k^{2}t\sigma))\\ &\times e^{-\frac{1}{2}i(A^{2}(t+3kt\sigma)+B^{2}(t+3kt\sigma)+2(k^{2}t+kx+k^{3}t\sigma-W(t)\sigma))}.\end{aligned}

Solution 4: The singular periodic solution corresponding to Family 4 is obtained

\begin{matrix} q (x, t) = & (\frac{\sqrt{(- 1 - 6 k σ) (- A^{2} + C^{2})}}{\sqrt{2 + 4 k γ_{1} σ}}) t a n (\frac{1}{4} \sqrt{- A^{2} + C^{2}} (4 k t + 2 x + (A^{2} - C^{2}) t σ + 6 k^{2} t σ)) \\ \times e^{- \frac{1}{2} i (A^{2} (t + 3 k t σ) - C^{2} (t + 3 k t σ) + 2 (k^{2} t + k x + k^{3} t σ - W (t) σ))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{\sqrt{(-1-6k\sigma)(-A^{2}+C^{2})}}{\sqrt{2+4k\gamma_{1}\sigma}})tan(\frac{1}{4}\sqrt{-A^{2}+C^{2}}(4kt+2x+(A^{2}-C^{2})t\sigma+6k^{2}t\sigma))\\ &\times e^{-\frac{1}{2}i(A^{2}(t+3kt\sigma)-C^{2}(t+3kt\sigma)+2(k^{2}t+kx+k^{3}t\sigma-W(t)\sigma))}.\end{aligned}

Solution 5: When B²–C²>0 the dark soliton solution corresponding to Family 5 is obtained

\begin{matrix} q (x, t) = & (\frac{\sqrt{(- 1 - 6 k σ) (B^{2} - C^{2})}}{\sqrt{2 + 4 k γ_{1} σ}}) t a n h (\frac{1}{4} \sqrt{B^{2} - C^{2}} (4 k t + 2 x + (B^{2} - C^{2}) t σ + 6 k^{2} t σ)) \\ \times e^{- \frac{1}{2} i (- C^{2} t + 2 k^{2} t + 2 k x - 3 C^{2} k t σ + 2 k^{3} t σ + B^{2} (t + 3 k t σ) - 2 σ W (t))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{\sqrt{(-1-6k\sigma)(B^{2}-C^{2})}}{\sqrt{2+4k\gamma_{1}\sigma}})tanh(\frac{1}{4}\sqrt{B^{2}-C^{2}}(4kt+2x+(B^{2}-C^{2})t\sigma+6k^{2}t\sigma))\\ &\times e^{-\frac{1}{2}i(-C^{2}t+2k^{2}t+2kx-3C^{2}kt\sigma+2k^{3}t\sigma+B^{2}(t+3kt\sigma)-2\sigma W(t))}.\end{aligned}

Solution 6: The singular soliton solution corresponding to Family 8 is obtained

q (x, t) = (- \frac{\sqrt{2} \sqrt{- 1 - 6 k σ}}{(2 k t + x + 3 k^{2} t σ) \sqrt{1 + 2 k γ_{1} σ}}) \times e^{- i k (k t + x + k^{2} t σ) + i σ W (t)} .

q(x,t)=(-\frac{\sqrt{2}\sqrt{-1-6k\sigma}}{(2kt+x+3k^{2}t\sigma)\sqrt{1+2k\gamma_{1}\sigma}})\times e^{-ik(kt+x+k^{2}t\sigma)+i\sigma W(t)}

Solution 7: The dark soliton solution corresponding to Family 11 is obtained

\begin{matrix} q (x, t) = & (\frac{A \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}} - \frac{(A - B) (- 1 + (A + B) e^{B (x - t (- 2 k - \frac{1}{2} (B^{2} + 6 k^{2}) σ))}) \sqrt{- 1 - 6 k σ}}{(- 1 + (A - B) e^{B (x - t (- 2 k - \frac{1}{2} (B^{2} + 6 k^{2}) σ))}) \sqrt{2 + 4 k γ_{1} σ}}) \\ \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2} - \frac{1}{2} B^{2} (1 + 3 k σ)))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}}-\frac{(A-B)(-1+(A+B)e^{B(x-t(-2k-\frac{1}{2}(B^{2}+6k^{2})\sigma))})\sqrt{-1-6k\sigma}}{(-1+(A-B)e^{B(x-t(-2k-\frac{1}{2}(B^{2}+6k^{2})\sigma))})\sqrt{2+4k\gamma_{1}\sigma}})\\ &\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}-\frac{1}{2}B^{2}(1+3k\sigma)))}.\end{aligned}

Solution 8: The dark soliton solution corresponding to Family 12 is obtained

\begin{matrix} q (x, t) = & (\frac{C \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}} + \frac{(- B + C) (1 + (B + C) e^{B (x - t (- 2 k - \frac{1}{2} (B^{2} + 6 k^{2}) σ))}) \sqrt{- 1 - 6 k σ}}{(- 1 + (B - C) e^{B (x - t (- 2 k - \frac{1}{2} (B^{2} + 6 k^{2}) σ))}) \sqrt{2 + 4 k γ_{1} σ}}) \\ \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2} - \frac{1}{2} B^{2} (1 + 3 k σ)))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{C\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}}+\frac{(-B+C)(1+(B+C)e^{B(x-t(-2k-\frac{1}{2}(B^{2}+6k^{2})\sigma))})\sqrt{-1-6k\sigma}}{(-1+(B-C)e^{B(x-t(-2k-\frac{1}{2}(B^{2}+6k^{2})\sigma))})\sqrt{2+4k\gamma_{1}\sigma}})\\ &\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}-\frac{1}{2}B^{2}(1+3k\sigma)))}.\end{aligned}

Solution 9: The singular soliton solution corresponding to Family 13 is obtained

\begin{matrix} q (x, t) = & (\frac{A \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}} - \frac{(A + B) (- A + B + e^{B (x - t (- 2 k - \frac{1}{2} (B^{2} + 6 k^{2}) σ))}) \sqrt{- 1 - 6 k σ}}{(- A - B + e^{B (x - t (- 2 k - \frac{1}{2} (B^{2} + 6 k^{2}) σ))}) \sqrt{2 + 4 k γ_{1} σ}}) \\ \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2} - \frac{1}{2} B^{2} (1 + 3 k σ)))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}}-\frac{(A+B)(-A+B+e^{B(x-t(-2k-\frac{1}{2}(B^{2}+6k^{2})\sigma))})\sqrt{-1-6k\sigma}}{(-A-B+e^{B(x-t(-2k-\frac{1}{2}(B^{2}+6k^{2})\sigma))})\sqrt{2+4k\gamma_{1}\sigma}})\\ &\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}-\frac{1}{2}B^{2}(1+3k\sigma)))}.\end{aligned}

Solution 10: The dark soliton solution corresponding to Family 14 is obtained

\begin{matrix} q (x, t) = & (\frac{A \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}} + \frac{2 A C e^{A (x - t (- 2 k - \frac{1}{2} (A^{2} + 6 k^{2}) σ))} \sqrt{- 1 - 6 k σ}}{(1 - C e^{A (x - t (- 2 k - \frac{1}{2} (A^{2} + 6 k^{2}) σ))}) \sqrt{2 + 4 k γ_{1} σ}}) \\ \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2} - \frac{1}{2} A^{2} (1 + 3 k σ)))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}}+\frac{2ACe^{A(x-t(-2k-\frac{1}{2}(A^{2}+6k^{2})\sigma))}\sqrt{-1-6k\sigma}}{(1-Ce^{A(x-t(-2k-\frac{1}{2}(A^{2}+6k^{2})\sigma))})\sqrt{2+4k\gamma_{1}\sigma}})\\ &\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}-\frac{1}{2}A^{2}(1+3k\sigma)))}.\end{aligned}

Solution 11: The singular soliton solution corresponding to Family 15 is obtained

\begin{matrix} q (x, t) = & (\frac{C \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}} - \frac{\sqrt{- 1 - 6 k σ} (2 + C (x + k t (2 + 3 k σ)))}{\sqrt{2 + 4 k γ_{1} σ} (x + k t (2 + 3 k σ))}) \\ \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2}))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{C\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}}-\frac{\sqrt{-1-6k\sigma}(2+C(x+kt(2+3k\sigma)))}{\sqrt{2+4k\gamma_{1}\sigma}(x+kt(2+3k\sigma))})\\ &\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}))}\end{aligned}

Solution 12: The singular soliton solution corresponding to Family 17 is obtained

q (x, t) = (- \frac{2 \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ} (x + k t (2 + 3 k σ))}) \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2}))} .

q(x,t)=(-\frac{2\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}(x+kt(2+3k\sigma))})\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}))}.

Solution 13: The singular periodic solution corresponding to Family 18 is obtained

\begin{matrix} q (x, t) = & (\frac{C \sqrt{- 1 - 6 k σ}}{\sqrt{2 + 4 k γ_{1} σ}}) t a n (\frac{1}{2} C (x - t (\frac{C^{2} σ}{2} - k (2 + 3 k σ)))) \\ \times e^{i (- k x + W (t) σ - t σ^{2} + t (- k^{2} - k^{3} σ + σ^{2} + \frac{1}{2} C^{2} (1 + 3 k σ)))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{C\sqrt{-1-6k\sigma}}{\sqrt{2+4k\gamma_{1}\sigma}})tan(\frac{1}{2}C(x-t(\frac{C^{2}\sigma}{2}-k(2+3k\sigma))))\\ &\times e^{i(-kx+W(t)\sigma-t\sigma^{2}+t(-k^{2}-k^{3}\sigma+\sigma^{2}+\frac{1}{2}C^{2}(1+3k\sigma)))}\end{aligned}

Set 2: In the case

\begin{matrix} a_{0} = & \frac{\sqrt{\frac{3}{2}} A}{\sqrt{- γ_{1} - γ_{2}}}, a_{1} = \frac{\sqrt{\frac{3}{2}} (- B + C)}{\sqrt{- γ_{1} - γ_{2}}}, k = - \frac{(γ_{1} + γ_{2} - 3)}{6 γ_{2} σ}, \\ w = & - \frac{(A^{2} + B^{2} - C^{2}) (3 - γ_{1} + γ_{2})}{4 γ_{2}} + \frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{3} σ^{2}} + σ^{2}, \\ v = & - \frac{9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) + 3 γ_{2} (2 + γ_{2} (- 1 + 2 (A^{2} + B^{2} - C^{2}) σ^{2}))}{12 γ_{2}^{2} σ}, \end{matrix}

\begin{aligned}a_{0}=&\frac{\sqrt{\frac{3}{2}}A}{\sqrt{-\gamma_{1}-\gamma_{2}}},a_{1}=\frac{\sqrt{\frac{3}{2}}(-B+C)}{\sqrt{-\gamma_{1}-\gamma_{2}}},k=-\frac{(\gamma_{1}+\gamma_{2}-3)}{6\gamma_{2}\sigma}.\\ w=&-\frac{(A^{2}+B^{2}-C^{2})(3-\gamma_{1}+\gamma_{2})}{4\gamma_{2}}+\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216{\gamma_{2}}^{3}\sigma^{2}}+\sigma^{2},\\ v=&-\frac{9+{\gamma_{1}}^{2}-2\gamma_{1}(3+\gamma_{2})+3\gamma_{2}(2+\gamma_{2}(-1+2(A^{2}+B^{2}-C^{2})\sigma^{2}))}{12{\gamma_{2}}^{2}\sigma}.\end{aligned}

where –γ₁–γ₂>0, γ₂≠0,σ≠0, the outcomes solutions are:

Solution 1: When −A²₁−B²₁+C²₁>0, the following singular periodic solution corresponding to Family 1 is obtained

\begin{matrix} q (x, t) = & (\sqrt{\frac{3}{2}} \frac{\sqrt{- A^{2} - B^{2} + C^{2}}}{\sqrt{- γ_{1} - γ_{2}}}) \times e^{\frac{1}{6} i (- \frac{3 (A^{2} + B^{2} - c^{2}) (3 - γ_{1} + γ_{2})}{2 γ_{2}} + \frac{t (- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{36 γ_{2}^{2} σ^{2}} + \frac{x (- 3 + γ_{1} + γ_{2})}{γ_{2} σ} + 6 W (t) σ)} \\ \times t a n (\frac{1}{2} \sqrt{- A^{2} - B^{2} + C^{2}} (x + \frac{t (9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) + 3 γ_{2} (2 + γ_{2} (- 1 + 2 (A^{2} + B^{2} - C^{2}) σ^{2})))}{12 γ_{2}^{2} σ})) . \end{matrix}

\begin{aligned}q(x,t)=&(\sqrt{\frac{3}{2}}\frac{\sqrt{-A^{2}-B^{2}+C^{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}})\times e^{\frac{1}{6}i(-\frac{3(A^{2}+B^{2}-c^{2})(3-\gamma_{1}+\gamma_{2})}{2\gamma_{2}}+\frac{t(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{36\gamma_{2}^{2}\sigma^{2}}+\frac{x(-3+\gamma_{1}+\gamma_{2})}{\gamma_{2}\sigma}+6W(t)\sigma)}\\ &\times tan(\frac{1}{2}\sqrt{-A^{2}-B^{2}+C^{2}}(x+\frac{t(9+\gamma_{1}^{2}-2\gamma_{1}(3+\gamma_{2})+3\gamma_{2}(2+\gamma_{2}(-1+2(A^{2}+B^{2}-C^{2})\sigma^{2})))}{12\gamma_{2}^{2}\sigma}))\end{aligned}

Solution 2: When A²+B²−C²>0, the dark soliton solution corresponding to Family 2 is obtained

\begin{matrix} q (x, t) = & (- \sqrt{\frac{3}{2}} \frac{\sqrt{A^{2} + B^{2} - C^{2}}}{\sqrt{- γ_{1} - γ_{2}}}) \times e^{\frac{1}{6} i (- \frac{3 (A^{2} + B^{2} - c^{2}) (3 - γ_{1} + γ_{2})}{2 γ_{2}} + \frac{t (- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{36 γ_{2}^{3} σ^{2}} + \frac{x (- 3 + γ_{1} + γ_{2})}{γ_{2} σ} + 6 W (t) σ)} \\ \times t a n h (\frac{1}{2} \sqrt{A^{2} + B^{2} - C^{2}} (x + \frac{t (9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) + 3 γ_{2} (2 + γ_{2} (- 1 + 2 (A^{2} + B^{2} - C^{2}) σ^{2})))}{12 γ_{2}^{2} σ})) . \end{matrix}

\begin{aligned}q(x,t)=&(-\sqrt{\frac{3}{2}}\frac{\sqrt{A^{2}+B^{2}-C^{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}})\times e^{\frac{1}{6}i(-\frac{3(A^{2}+B^{2}-c^{2})(3-\gamma_{1}+\gamma_{2})}{2\gamma_{2}}+\frac{t(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{36\gamma_{2}^{3}\sigma^{2}}+\frac{x(-3+\gamma_{1}+\gamma_{2})}{\gamma_{2}\sigma}+6W(t)\sigma)}\\ &\times tanh(\frac{1}{2}\sqrt{A^{2}+B^{2}-C^{2}}(x+\frac{t(9+{\gamma_{1}}^{2}-2\gamma_{1}(3+\gamma_{2})+3\gamma_{2}(2+\gamma_{2}(-1+2(A^{2}+B^{2}-C^{2})\sigma^{2})))}{12\gamma_{2}^{2}\sigma}))\end{aligned}

Solution 3: The dark soliton solution corresponding to Family 3 is obtained

\begin{matrix} q (x, t) = & (- \sqrt{\frac{3}{2}} \frac{\sqrt{A^{2} + B^{2}}}{\sqrt{- γ_{1} - γ_{2}}}) \times e^{\frac{1}{6} i (- \frac{3 (A^{2} + B^{2}) t (3 - γ_{1} + γ_{2})}{2 γ_{2}} + \frac{t (- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{36 γ_{2}^{3} σ^{2}} + \frac{x (- 3 + γ_{1} + γ_{2})}{γ_{2} σ} + 6 W (t) σ)} \\ \times t a n h (\frac{1}{2} \sqrt{A^{2} + B^{2}} (x + \frac{t (9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) + 3 γ_{2} (2 + γ_{2} (- 1 + 2 (A^{2} + B^{2}) σ^{2})))}{12 γ_{2}^{2} σ})) . \end{matrix}

\begin{aligned}q(x,t)=&(-\sqrt{\frac{3}{2}}\frac{\sqrt{A^{2}+B^{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}})\times e^{\frac{1}{6}i(-\frac{3(A^{2}+B^{2})t(3-\gamma_{1}+\gamma_{2})}{2\gamma_{2}}+\frac{t(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{36\gamma_{2}^{3}\sigma^{2}}+\frac{x(-3+\gamma_{1}+\gamma_{2})}{\gamma_{2}\sigma}+6W(t)\sigma)}\\ &\times tanh(\frac{1}{2}\sqrt{A^{2}+B^{2}}(x+\frac{t(9+\gamma_{1}^{2}-2\gamma_{1}(3+\gamma_{2})+3\gamma_{2}(2+\gamma_{2}(-1+2(A^{2}+B^{2})\sigma^{2})))}{12\gamma_{2}^{2}\sigma}))\end{aligned}

Solution 4: The singular periodic solution corresponding to Family 4 is obtained

\begin{matrix} q (x, t) = & (\sqrt{\frac{3}{2}} \frac{\sqrt{- A^{2} + C^{2}}}{\sqrt{- γ_{1} - γ_{2}}}) \times e^{\frac{1}{6} i (- \frac{3 (A^{2} + C^{2}) t (3 - γ_{1} + γ_{2})}{2 γ_{2}} + \frac{t (- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{36 γ_{2}^{3} σ^{2}} + \frac{x (- 3 + γ_{1} + γ_{2})}{γ_{2} σ} + 6 W (t) σ)} \\ \times t a n (\frac{1}{2} \sqrt{- A^{2} + C^{2}} (x + \frac{t (9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) + 3 γ_{2} (2 + γ_{2} (- 1 + 2 (A^{2} - C^{2}) σ^{2})))}{12 γ_{2}^{2} σ})) . \end{matrix}

\begin{aligned}q(x,t)=&(\sqrt{\frac{3}{2}}\frac{\sqrt{-A^{2}+C^{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}})\times e^{\frac{1}{6}i(-\frac{3(A^{2}+C^{2})t(3-\gamma_{1}+\gamma_{2})}{2\gamma_{2}}+\frac{t(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{36\gamma_{2}^{3}\sigma^{2}}+\frac{x(-3+\gamma_{1}+\gamma_{2})}{\gamma_{2}\sigma}+6W(t)\sigma)}\\ &\times tan(\frac{1}{2}\sqrt{-A^{2}+C^{2}}(x+\frac{t(9+\gamma_{1}^{2}-2\gamma_{1}(3+\gamma_{2})+3\gamma_{2}(2+\gamma_{2}(-1+2(A^{2}-C^{2})\sigma^{2})))}{12\gamma_{2}^{2}\sigma}))\end{aligned}

Solution 5: The dark soliton solution corresponding to Family 5 is obtained

\begin{matrix} q (x, t) = & (- \sqrt{\frac{3}{2}} \frac{\sqrt{B^{2} - C^{2}}}{\sqrt{- γ_{1} - γ_{2}}}) \times e^{\frac{1}{6} i (\frac{3 (A^{2} - C^{2}) t (3 - γ_{1} + γ_{2})}{2 γ_{2}} + \frac{t (- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{36 γ_{2}^{2} σ^{2}} + \frac{x (- 3 + γ_{1} + γ_{2})}{γ_{2} σ} + 6 W (t) σ)} \\ \times t a n h (\frac{1}{2} \sqrt{B^{2} - C^{2}} (x + \frac{t (9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) + 3 γ_{2} (2 + γ_{2} (- 1 + 2 (A^{2} - C^{2}) σ^{2})))}{12 γ_{2}^{2} σ})) . \end{matrix}

\begin{aligned}q(x,t)=&(-\sqrt{\frac{3}{2}}\frac{\sqrt{B^{2}-C^{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}})\times e^{\frac{1}{6}i(\frac{3(A^{2}-C^{2})t(3-\gamma_{1}+\gamma_{2})}{2\gamma_{2}}+\frac{t(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{36\gamma_{2}^{2}\sigma^{2}}+\frac{x(-3+\gamma_{1}+\gamma_{2})}{\gamma_{2}\sigma}+6W(t)\sigma)}\\ &\times tanh(\frac{1}{2}\sqrt{B^{2}-C^{2}}(x+\frac{t(9+{\gamma_{1}}^{2}-2\gamma_{1}(3+\gamma_{2})+3\gamma_{2}(2+\gamma_{2}(-1+2(A^{2}-C^{2})\sigma^{2})))}{12\gamma_{2}^{2}\sigma}))\end{aligned}

Solution 6: The singular soliton solution corresponding to Family 8 is obtained

\begin{matrix} q (x, t) = & (- \frac{12 \sqrt{6} γ_{2}^{2} σ}{\sqrt{- γ_{1} - γ_{2}} (t (9 + γ_{1}^{2} + 6 γ_{2} - 3 γ_{2}^{2} - 2 γ_{1} (3 + γ_{2})) + 12 x γ_{2}^{2} σ)}) \\ \times e^{\frac{i (- 3 + γ_{1} + γ_{2}) (t (9 + γ_{1}^{2} + 12 γ_{2} - 5 γ_{2}^{2} - 2 γ_{1} (3 + 2 γ_{2})) + 36 x γ_{2}^{2} σ)}{216 γ_{2}^{3} σ^{2}} + i σ W (t)} . \end{matrix}

\begin{aligned}q(x,t)=&(-\frac{12\sqrt{6}{\gamma_{2}}^{2}\sigma}{\sqrt{-\gamma_{1}-\gamma_{2}}(t(9+{\gamma_{1}}^{2}+6{\gamma_{2}}-3{\gamma_{2}}^{2}-2\gamma_{1}(3+\gamma_{2}))+12x{\gamma_{2}}^{2}\sigma)})\\ &\times e^{\frac{i(-3+\gamma_{1}+\gamma_{2})(t(9+{\gamma_{1}}^{2}+12\gamma_{2}-5\gamma_{2}^{2}-2\gamma_{1}(3+2\gamma_{2}))+36x{\gamma_{2}}^{2}\sigma)}{216{\gamma_{2}}^{3}\sigma^{2}}+i\sigma W(t)}\end{aligned}

Solution 7: The dark soliton solution corresponding to Family 11 is obtained

\begin{matrix} q (x, t) = & (\frac{A \sqrt{\frac{3}{2}}}{\sqrt{- γ_{1} - γ_{2}}} - \frac{\sqrt{\frac{3}{2}} (A - B) (- 1 + (A + B) e^{B (x - t v)})}{\sqrt{- γ_{1} - γ_{2}} (- 1 + (A - B) e^{B (x - t v)})}) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W (t) σ - t σ^{2} + t (- \frac{B^{2} (3 - γ_{1} + γ_{2})}{4 γ_{2}} + \frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{3} σ^{2}} + σ^{2}))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A\sqrt{\frac{3}{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}}-\frac{\sqrt{\frac{3}{2}}(A-B)(-1+(A+B)e^{B(x-tv)})}{\sqrt{-\gamma_{1}-\gamma_{2}}(-1+(A-B)e^{B(x-tv)})})\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W(t)\sigma-t\sigma^{2}+t(-\frac{B^{2}(3-\gamma_{1}+\gamma_{2})}{4\gamma_{2}}+\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{3}\sigma^{2}}+\sigma^{2}))}\end{aligned}

Solution 8: The dark soliton solution corresponding to Family 12 is obtained

\begin{matrix} q (x, t) = & (\frac{C \sqrt{\frac{3}{2}}}{\sqrt{- γ_{1} - γ_{2}}} + \frac{\sqrt{\frac{3}{2}} (- B + C) (1 + (B + C) e^{B (x - t v)})}{\sqrt{- γ_{1} - γ_{2}} (- 1 + (B - C) e^{B (x - t v)})}) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W (t) σ - t σ^{2} + t (- \frac{B^{2} (3 - γ_{1} + γ_{2})}{4 γ_{2}} + \frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{3} σ^{2}} + σ^{2}))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{C\sqrt{\frac{3}{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}}+\frac{\sqrt{\frac{3}{2}}(-B+C)(1+(B+C)e^{B(x-tv)})}{\sqrt{-\gamma_{1}-\gamma_{2}}(-1+(B-C)e^{B(x-tv)})})\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W(t)\sigma-t\sigma^{2}+t(-\frac{B^{2}(3-\gamma_{1}+\gamma_{2})}{4\gamma_{2}}+\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{2}\sigma^{2}}+\sigma^{2}))}\end{aligned}

Solution 9: The singular soliton solution corresponding to Family 13 is obtained

\begin{matrix} q (x, t) = & (\frac{A \sqrt{\frac{3}{2}}}{\sqrt{- γ_{1} - γ_{2}}} - \frac{\sqrt{\frac{3}{2}} (A + B) (- A + B + e^{B (x - t v)})}{\sqrt{- γ_{1} - γ_{2}} (- A - B + e^{B (x - t v)})}) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W (t) σ - t σ^{2} + t (- \frac{B^{2} (3 - γ_{1} + γ_{2})}{4 γ_{2}} + \frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{3} σ^{2}} + σ^{2}))} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A\sqrt{\frac{3}{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}}-\frac{\sqrt{\frac{3}{2}}(A+B)(-A+B+e^{B(x-tv)})}{\sqrt{-\gamma_{1}-\gamma_{2}}(-A-B+e^{B(x-tv)})})\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W(t)\sigma-t\sigma^{2}+t(-\frac{B^{2}(3-\gamma_{1}+\gamma_{2})}{4\gamma_{2}}+\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{2}\sigma^{2}}+\sigma^{2}))}\end{aligned}

Solution 10: The dark soliton solution corresponding to Family 14 is obtained

\begin{matrix} q (x, t) = & (\frac{A \sqrt{\frac{3}{2}}}{\sqrt{- γ_{1} - γ_{2}}} + \frac{\sqrt{\frac{3}{2}} 2 A C e^{A (x - t v)}}{(1 - C e^{A (x - t v)}) \sqrt{- γ_{1} - γ_{2}}}) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W (t) σ - t σ^{2} + t (- \frac{A^{2} (3 - γ_{1} + γ_{2})}{4 γ_{2}} + \frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{2} σ^{2}} + σ^{2}))} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A\sqrt{\frac{3}{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}}+\frac{\sqrt{\frac{3}{2}}2ACe^{A(x-tv)}}{(1-Ce^{A(x-tv)})\sqrt{-\gamma_{1}-\gamma_{2}}})\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W(t)\sigma-t\sigma^{2}+t(-\frac{A^{2}(3-\gamma_{1}+\gamma_{2})}{4\gamma_{2}}+\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{2}\sigma^{2}}+\sigma^{2}))}\end{aligned}

Solution 11: The singular soliton solution corresponding to Family 15 is obtained

\begin{matrix} q (x, t) = & (\frac{C \sqrt{\frac{3}{2}}}{\sqrt{- γ_{1} - γ_{2}}} - \frac{\sqrt{\frac{3}{2}} (2 + C (x - \frac{t (- 3 + γ_{1} + γ_{2}) (3 - γ_{1} + 3 γ_{2})}{12 γ_{2}^{2} σ}))}{\sqrt{- γ_{1} - γ_{2}} (x - \frac{t (- 3 + γ_{1} + γ_{2}) (3 - γ_{1} + 3 γ_{2})}{12 γ_{2}^{2} σ})}) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W σ - t σ^{2} + t (\frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{3} σ^{2}} + σ^{2}))} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{C\sqrt{\frac{3}{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}}-\frac{\sqrt{\frac{3}{2}}(2+C(x-\frac{t(-3+\gamma_{1}+\gamma_{2})(3-\gamma_{1}+3\gamma_{2})}{12\gamma_{2}^{2}\sigma}))}{\sqrt{-\gamma_{1}-\gamma_{2}}(x-\frac{t(-3+\gamma_{1}+\gamma_{2})(3-\gamma_{1}+3\gamma_{2})}{12\gamma_{2}^{2}\sigma})})\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W\sigma-t\sigma^{2}+t(\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{3}\sigma^{2}}+\sigma^{2}))}\end{aligned}

Solution 12: The singular soliton solution corresponding to Family 17 is obtained

\begin{matrix} q (x, t) = & (- \frac{\sqrt{6}}{\sqrt{- γ_{1} - γ_{2}} (x - \frac{t (- 3 + γ_{1} + γ_{2}) (3 - γ_{1} + 3 γ_{2})}{12 γ_{2}^{3} σ})}) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W σ - t σ^{2} + t (\frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{2} σ^{2}} + σ^{2}))} \end{matrix}

\begin{aligned}q(x,t)=&(-\frac{\sqrt{6}}{\sqrt{-\gamma_{1}-\gamma_{2}}(x-\frac{t(-3+\gamma_{1}+\gamma_{2})(3-\gamma_{1}+3\gamma_{2})}{12\gamma_{2}^{3}\sigma})})\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W\sigma-t\sigma^{2}+t(\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{2}\sigma^{2}}+\sigma^{2}))}\end{aligned}

Solution 13: The singular periodic solution corresponding to Family 18 is obtained

\begin{matrix} q (x, t) = & (\sqrt{\frac{3}{2}} \frac{C}{\sqrt{- γ_{1} - γ_{2}}}) t a n (\frac{1}{2} C (x + \frac{t (9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) - 3 γ_{2} (- 2 + γ_{2} + 2 C^{2} γ_{2} σ^{2}))}{12 γ_{2}^{2} σ})) \\ \times e^{i (\frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2} σ} + W (t) σ - t σ^{2} + t (\frac{C^{2} (3 - γ_{1} + γ_{2})}{4 γ_{2}} + \frac{(- 3 + γ_{1} - 5 γ_{2}) {(- 3 + γ_{1} + γ_{2})}^{2}}{216 γ_{2}^{2} σ^{2}} + σ^{2}))} \end{matrix}

\begin{aligned}q(x,t)=&(\sqrt{\frac{3}{2}}\frac{C}{\sqrt{-\gamma_{1}-\gamma_{2}}})tan(\frac{1}{2}C(x+\frac{t(9+\gamma_{1}^{2}-2\gamma_{1}(3+\gamma_{2})-3\gamma_{2}(-2+\gamma_{2}+2C^{2}\gamma_{2}\sigma^{2}))}{12\gamma_{2}^{2}\sigma}))\\ &\times e^{i(\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}\sigma}+W(t)\sigma-t\sigma^{2}+t(\frac{C^{2}(3-\gamma_{1}+\gamma_{2})}{4\gamma_{2}}+\frac{(-3+\gamma_{1}-5\gamma_{2})(-3+\gamma_{1}+\gamma_{2})^{2}}{216\gamma_{2}^{2}\sigma^{2}}+\sigma^{2}))}\end{aligned}

3.2

Application of (G’/G,1/G)-expansion method

In this part, we will apply the (G’/G,1/G)-expansion method to find the exact traveling wave solutions of Eq.(1).Balancing the terms of Q³ and Qⁿ in Eq.(18) and Eq.(20) gives n = 1 Putting n into Eq.(16), we derive

Q (ξ) = a_{0} + a_{1} ϕ (ξ) + B_{1} Ψ (ξ)

Q(\xi)=a_{0}+a_{1}\phi(\xi)+B_{1}\psi(\xi)

where a₀,a₁,B₁ are constants to be determined, such that $a_{1}^{2} + B_{1}^{2} \neq 0$ a_{1}^{2}+B_{1}^{2}\ne0. To analyze NLEE based on the sign of δ three distinct cases must be considered. Let us now elaborate on these three scenarios:

Case 1: When δ<0.

By putting Eq.(50) into Eq.(18) and Eq.(20) and applying Eq.(8) and Eq.(10), the lefthand side of equation Eq.(18) and Eq.(20) transforms into a polynomial in ϕ and Ψ. Collecting each coefficient of this polynomial to zero yields in a system of algebraic equations involving a₀,a₁,B₁,ρ,δ,A₁, and A₁. Utilizing package program, we obtain the desired solutions as follows

Set 1:

\begin{matrix} a_{0} = & 0, a_{1} = a_{1}, B_{1} = \frac{a_{1} δ \sqrt{ρ^{2} + (- A_{1}^{2} + A_{2}^{2}) δ^{2}}}{{(- δ)}^{3 / 2}}, γ_{1} = - \frac{6 (3 - 2 a_{1}^{2}) k^{2} + δ + 2 a_{1}^{2} δ + 12 k v}{8 a_{1}^{2} k (2 k + v)}, \\ ω = & \frac{- 24 k^{6} - 20 k^{4} δ + δ^{3} - 2 k^{2} (- 16 + δ^{2}) + 24 k^{5} v - 40 k^{3} δ v + 2 k (16 + 3 δ^{2}) v + 8 v^{2}}{2 {(- 6 k^{2} + δ)}^{2}}, \\ γ_{2} = & - \frac{(1 + 2 a_{1}^{2}) (6 k^{2} - δ)}{8 a_{1}^{2} k (2 k + v)}, σ = - \frac{2 (2 k + v)}{6 k^{2} - δ}, \end{matrix}

\begin{aligned}a_{0}=&0,a_{1}=a_{1},B_{1}=\frac{a_{1}\delta\sqrt{\rho^{2}+(-A_{1}^{2}+A_{2}^{2})\delta^{2}}}{(-\delta)^{3/2}},\gamma_{1}=-\frac{6(3-2a_{1}^{2})k^{2}+\delta+2a_{1}^{2}\delta+12kv}{8a_{1}^{2}k(2k+v)},\\ \omega=&\frac{-24k^{6}-20k^{4}\delta+\delta^{3}-2k^{2}(-16+\delta^{2})+24k^{5}v-40k^{3}\delta v+2k(16+3\delta^{2})v+8v^{2}}{2(-6k^{2}+\delta)^{2}},\\ \gamma_{2}=&-\frac{(1+2a_{1}^{2})(6k^{2}-\delta)}{8a_{1}^{2}k(2k+v)},\sigma=-\frac{2(2k+v)}{6k^{2}-\delta},\end{aligned}

where $ρ^{2} + (- A_{1}^{2} + A_{2}^{2}) δ^{2} > 0$ \rho^{2}+(-A_{1}^{2}+A_{2}^{2})\delta^{2}>0,a₁ ≠ 0, and k ≠ 0. Substituting the values of Eq.(51) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows

\begin{matrix} q (x, t) = & (\frac{a_{1} δ \sqrt{ρ^{2} + (- A_{1}^{2} + A_{2}^{2}) δ^{2}}}{{(- δ)}^{3 / 2} (\frac{ρ}{δ} + A_{1} c o s h (\sqrt{- δ} (x - t v)) + A_{2} s i n h (\sqrt{- δ} (x - t v)))} \\ + \frac{a_{1} (A_{2} \sqrt{- δ} c o s h (\sqrt{- δ} (x - t v)) + A_{1} \sqrt{- δ} s i n h (\sqrt{- δ} (x - t v)))}{\frac{ρ}{δ} + A_{1} c o s h (\sqrt{- δ} (x - t v)) + A_{2} s i n h (\sqrt{- δ} (x - t v))}) \\ \times e^{i (- k x - \frac{4 t {(2 k + v)}^{2}}{{(6 k^{2} - δ)}^{2}} + \frac{t (- 24 k^{6} - 20 k^{4} δ + δ^{3} - 2 k^{2} (- 16 + δ^{2}) + 24 k^{5} v - 40 k^{3} δ v + 2 k (16 + 3 δ^{2}) v + 8 v^{2})}{2 {(- 6 k^{2} + δ)}^{2}} - \frac{2 (2 k + v) W (t)}{6 k^{2} - δ})} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{a_{1}\delta\sqrt{\rho^{2}+(-A_{1}^{2}+A_{2}^{2})\delta^{2}}}{(-\delta)^{3/2}(\frac{\rho}{\delta}+A_{1}cosh(\sqrt{-\delta}(x-tv))+A_{2}sinh(\sqrt{-\delta}(x-tv)))}\\ &+\frac{a_{1}(A_{2}\sqrt{-\delta}cosh(\sqrt{-\delta}(x-tv))+A_{1}\sqrt{-\delta}sinh(\sqrt{-\delta}(x-tv)))}{\frac{\rho}{\delta}+A_{1}cosh(\sqrt{-\delta}(x-tv))+A_{2}sinh(\sqrt{-\delta}(x-tv))})\\ &\times e^{i(-kx-\frac{4t(2k+v)^{2}}{(6k^{2}-\delta)^{2}}+\frac{t(-24k^{6}-20k^{4}\delta+\delta^{3}-2k^{2}(-16+\delta^{2})+24k^{5}v-40k^{3}\delta v+2k(16+3\delta^{2})v+8v^{2})}{2(-6k^{2}+\delta)^{2}}-\frac{2(2k+v)W(t)}{6k^{2}-\delta})}\end{aligned}

Eq.(52) describes a singular soliton solution.

Set 2:

a_{0} = 0, a_{1} = \frac{1}{\sqrt{- 2 - \frac{8 γ_{2}}{3}}}, B_{1} = 0, γ_{1} = 3 (1 + γ_{2}), σ = \frac{2 v}{δ}, w = - \frac{δ}{2} - \frac{δ^{2}}{27 v^{2}} + \frac{4 v^{2}}{δ^{2}}, k = - \frac{δ}{3 v}, A_{1} = \frac{\sqrt{ρ^{2} + {A_{2}}^{2} δ^{2}}}{δ}

a_{0}=0,a_{1}=\frac{1}{\sqrt{-2-\frac{8\gamma_{2}}{3}}},B_{1}=0,\gamma_{1}=3(1+\gamma_{2}),\sigma=\frac{2v}{\delta},w=-\frac{\delta}{2}-\frac{\delta^{2}}{27v^{2}}+\frac{4v^{2}}{\delta^{2}},k=-\frac{\delta}{3v},A_{1}=\frac{\sqrt{\rho^{2}+{A_{2}}^{2}\delta^{2}}}{\delta}

where $- 2 - \frac{8 γ_{2}}{3} > 0$ -2-\frac{8\gamma_{2}}{3}>0, and v ≠ 0. Substituting the values of Eq.(53) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows:

\begin{matrix} q (x, t) = & (\frac{A_{2} \sqrt{- δ} c o s h (\sqrt{- δ} (x - t v)) + \frac{\sqrt{- δ} \sqrt{ρ^{2} + A_{2}^{2} δ^{2}} s i n h (\sqrt{- δ} (x - t v))}{δ}}{\sqrt{- 2 - \frac{8 γ_{2}}{3}} (\frac{ρ}{δ} + \frac{\sqrt{ρ^{2} + A_{2}^{2} δ^{2}} c o s h (\sqrt{- δ} (x - t v))}{δ} + A_{2} s i n h (\sqrt{- δ} (x - t v)))}) \\ \times e^{i (\frac{x δ}{3 v} - \frac{4 t v^{2}}{δ^{2}} + t (- \frac{δ}{2} - \frac{δ^{2}}{27 v^{2}} + \frac{4 v^{2}}{δ^{2}}) + \frac{2 v W (t)}{δ})} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A_{2}\sqrt{-\delta}cosh(\sqrt{-\delta}(x-tv))+\frac{\sqrt{-\delta}\sqrt{\rho^{2}+A_{2}^{2}\delta^{2}}sinh(\sqrt{-\delta}(x-tv))}{\delta}}{\sqrt{-2-\frac{8\gamma_{2}}{3}}(\frac{\rho}{\delta}+\frac{\sqrt{\rho^{2}+A_{2}^{2}\delta^{2}}cosh(\sqrt{-\delta}(x-tv))}{\delta}+A_{2}sinh(\sqrt{-\delta}(x-tv)))})\\ &\times e^{i(\frac{x\delta}{3v}-\frac{4tv^{2}}{\delta^{2}}+t(-\frac{\delta}{2}-\frac{\delta^{2}}{27v^{2}}+\frac{4v^{2}}{\delta^{2}})+\frac{2vW(t)}{\delta})}\end{aligned}

Eq.(54) describes a singular soliton solution.

Set 3:

\begin{matrix} a_{0} = & 0, a_{1} = \frac{\sqrt{- 6 k^{2} + δ}}{\sqrt{4 k^{2} (3 + 4 γ_{2}) - 2 δ + 8 k γ_{2} v}}, B_{1} = 0, γ_{1} = \frac{18 k^{2} (1 + γ_{2}) + (- 3 + γ_{2}) δ + 12 k γ_{2} v}{6 k^{2} - δ}, \\ ω = & \frac{- 24 k^{6} - 20 k^{4} δ + δ^{3} - 2 k^{2} (- 16 + δ^{2}) + 24 k^{5} v - 40 k^{3} δ v + 2 k (16 + 3 δ^{2}) v + 8 v^{2}}{2 {(- 6 k^{2} + δ)}^{2}}, \\ σ = & - \frac{2 (2 k + v)}{6 k^{2} - δ}, A_{1} = \frac{\sqrt{ρ^{2} + A_{2}^{2} δ^{2}}}{δ}, \end{matrix}

\begin{aligned}a_{0}=&0,a_{1}=\frac{\sqrt{-6k^{2}+\delta}}{\sqrt{4k^{2}(3+4\gamma_{2})-2\delta+8k\gamma_{2}v}},B_{1}=0,\gamma_{1}=\frac{18k^{2}(1+\gamma_{2})+(-3+\gamma_{2})\delta+12k\gamma_{2}v}{6k^{2}-\delta},\\ \omega=&\frac{-24k^{6}-20k^{4}\delta+\delta^{3}-2k^{2}(-16+\delta^{2})+24k^{5}v-40k^{3}\delta v+2k(16+3\delta^{2})v+8v^{2}}{2(-6k^{2}+\delta)^{2}},\\ \sigma=&-\frac{2(2k+v)}{6k^{2}-\delta},A_{1}=\frac{\sqrt{\rho^{2}+A_{2}^{2}\delta^{2}}}{\delta},\end{aligned}

where $- 6 k^{2} + δ > 0$ -6k^{2}+\delta>0, and $4 k^{2} (3 + 4 γ_{2}) - 2 δ + 8 k γ_{2} v > 0$ 4k^{2}(3+4\gamma_{2})-2\delta+8k\gamma_{2}v>0.

Substituting the values of Eq.(55) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows

\begin{matrix} q (x, t) = & (\frac{\sqrt{- δ (- 6 k^{2} + δ)} (A_{2} δ c o s h (\sqrt{- δ} (x - t v)) + \sqrt{ρ^{2} + A_{2}^{2} δ^{2}} s i n h (\sqrt{- δ} (x - t v)))}{\sqrt{4 k^{2} (3 + 4 γ_{2}) - 2 δ + 8 k γ_{2} v} (ρ + \sqrt{ρ^{2} + A_{2}^{2} δ^{2}} c o s h (\sqrt{- δ} (x - t v)) + A_{2} δ s i n h (\sqrt{- δ} (x - t v)))}) \\ \times e^{- \frac{i (4 k^{4} t + 4 k^{2} t δ + t δ^{2} - 2 k δ (x - 3 t v) + 4 k^{3} (3 x - t v) + 4 (2 k + v) W (t))}{12 k^{2} - 2 δ}} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{\sqrt{-\delta(-6k^{2}+\delta)}(A_{2}\delta~cosh(\sqrt{-\delta}(x-tv))+\sqrt{\rho^{2}+A_{2}^{2}\delta^{2}}sinh(\sqrt{-\delta}(x-tv)))}{\sqrt{4k^{2}(3+4\gamma_{2})-2\delta+8k\gamma_{2}v}(\rho+\sqrt{\rho^{2}+A_{2}^{2}\delta^{2}}cosh(\sqrt{-\delta}(x-tv))+A_{2}\delta~sinh(\sqrt{-\delta}(x-tv)))})\\ &\times e^{-\frac{i(4k^{4}t+4k^{2}t\delta+t\delta^{2}-2k\delta(x-3tv)+4k^{3}(3x-tv)+4(2k+v)W(t))}{12k^{2}-2\delta}}\end{aligned}

Eq.(56) describes a dark soliton solution.

Case 2: When δ>0.

By putting Eq.(50) into Eq.(18) and Eq.(20) and applying Eq.(8) and Eq.(12), the lefthand side of Eq.(18) and Eq.(20) transforms into a polynomial in Ψ and ϕ. Collecting each coefficient of this polynomial to zero yields in a system of algebraic equations involving a₀,a₁,B₁,ρ,δ,A₁, and A₂. Via package program, we obtain solutions as follows

Set 1:

\begin{matrix} a_{0} = & 0, a_{1} = - \frac{1}{\sqrt{- 2 + 4 k γ_{2}}}, B_{1} = - \frac{\sqrt{- ρ^{2} + A_{1}^{2} δ^{2} + A_{2}^{2} δ^{2}}}{\sqrt{- 2 δ + 4 k γ_{2} δ}}, γ_{1} = 3 - (1 + 6 k) γ_{2}, \\ σ = & 1, v = - k (2 + 3 k) + \frac{δ}{2}, w = \frac{1}{2} (2 + δ + k (- 2 k (1 + k) + 3 δ)), \end{matrix}

\begin{aligned}a_{0}=&0,a_{1}=-\frac{1}{\sqrt{-2+4k\gamma_{2}}},B_{1}=-\frac{\sqrt{-\rho^{2}+A_{1}^{2}\delta^{2}+A_{2}^{2}\delta^{2}}}{\sqrt{-2\delta+4k\gamma_{2}\delta}},\gamma_{1}=3-(1+6k)\gamma_{2},\\ \sigma=&1,v=-k(2+3k)+\frac{\delta}{2},w=\frac{1}{2}(2+\delta+k(-2k(1+k)+3\delta)),\end{aligned}

where $- 2 + 4 k γ_{2} > 0$ -2+4k\gamma_{2}>0 and $- ρ^{2} + A_{1}^{2} δ^{2} + A_{2}^{2} δ^{2} > 0$ -\rho^{2}+A_{1}^{2}\delta^{2}+A_{2}^{2}\delta^{2}>0.

Substituting the values of Eq.(57) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows:

\begin{matrix} q (x, t) = & (\frac{A_{1} \sqrt{δ} s i n (\sqrt{δ} (x - t (- k (2 + 3 k) + \frac{δ}{2}))) - A_{2} \sqrt{δ} c o s (\sqrt{δ} (x - t (- k (2 + 3 k) + \frac{δ}{2})))}{\sqrt{- 2 + 4 k γ_{2}} (\frac{ρ}{δ} + A_{1} c o s (\sqrt{δ} (x - t (- k (2 + 3 k) + \frac{δ}{2}))) + A_{2} s i n (\sqrt{δ} (x - t (- k (2 + 3 k) + \frac{δ}{2}))))} \\ - \frac{\sqrt{δ} \sqrt{- ρ^{2} + A_{1}^{2} δ^{2} + A_{2}^{2} δ^{2}}}{\sqrt{- 2 δ + 4 k γ_{2} δ} (\frac{ρ}{δ} + A_{1} c o s (\sqrt{δ} (x - t (- k (2 + 3 k) + \frac{δ}{2}))) + A_{2} s i n (\sqrt{δ} (x - t (- k (2 + 3 k) + \frac{δ}{2}))))}) \\ \times e^{i (- t - k x + \frac{1}{2} t (2 + δ + k (- 2 k (1 + k) + 3 δ)) + W (t))} . \end{matrix}

\begin{aligned}q(x,t)=&(\frac{A_{1}\sqrt{\delta}sin(\sqrt{\delta}(x-t(-k(2+3k)+\frac{\delta}{2})))-A_{2}\sqrt{\delta}cos(\sqrt{\delta}(x-t(-k(2+3k)+\frac{\delta}{2})))}{\sqrt{-2+4k\gamma_{2}}(\frac{\rho}{\delta}+A_{1}cos(\sqrt{\delta}(x-t(-k(2+3k)+\frac{\delta}{2})))+A_{2}sin(\sqrt{\delta}(x-t(-k(2+3k)+\frac{\delta}{2}))))}\\ &-\frac{\sqrt{\delta}\sqrt{-\rho^{2}+A_{1}^{2}\delta^{2}+A_{2}^{2}\delta^{2}}}{\sqrt{-2\delta+4k\gamma_{2}\delta}(\frac{\rho}{\delta}+A_{1}cos(\sqrt{\delta}(x-t(-k(2+3k)+\frac{\delta}{2})))+A_{2}sin(\sqrt{\delta}(x-t(-k(2+3k)+\frac{\delta}{2}))))})\\ &\times e^{i(-t-kx+\frac{1}{2}t(2+\delta+k(-2k(1+k)+3\delta))+W(t))}.\end{aligned}

Eq.(58) describes a singular periodic soliton solution.

Set 2:

\begin{matrix} a_{0} = & 0, a_{1} = \frac{\sqrt{\frac{3}{2}}}{\sqrt{- γ_{1} - γ_{2}}}, B_{1} = \frac{\sqrt{\frac{3}{2}} \sqrt{ρ^{2} - (A_{1}^{2} + A_{2}^{2}) δ^{2}}}{\sqrt{(γ_{1} + γ_{2}) δ}}, v = - \frac{9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) - 3 γ_{2} (- 2 + γ_{2} + 2 γ_{2} δ)}{12 γ_{2}^{2}}, \\ ω = & \frac{- 27 + γ_{1}^{3} - 3 γ_{1}^{2} (3 + γ_{2}) - 9 γ_{1} (- 3 + γ_{2} (- 2 + γ_{2} + 6 γ_{2} δ)) + γ_{2} (- 27 + γ_{2} (27 + 211 γ_{2} + 54 (3 + γ_{2}) δ))}{216 γ_{2}^{3}}, \\ σ = & 1, k = - \frac{- 3 + γ_{1} + γ_{2}}{6 γ_{2}}, \end{matrix}

\begin{aligned}a_{0}=&0,a_{1}=\frac{\sqrt{\frac{3}{2}}}{\sqrt{-\gamma_{1}-\gamma_{2}}},B_{1}=\frac{\sqrt{\frac{3}{2}}\sqrt{\rho^{2}-(A_{1}^{2}+A_{2}^{2})\delta^{2}}}{\sqrt{(\gamma_{1}+\gamma_{2})\delta}},v=-\frac{9+\gamma_{1}^{2}-2\gamma_{1}(3+\gamma_{2})-3\gamma_{2}(-2+\gamma_{2}+2\gamma_{2}\delta)}{12\gamma_{2}^{2}},\\ \omega=&\frac{-27+\gamma_{1}^{3}-3\gamma_{1}^{2}(3+\gamma_{2})-9\gamma_{1}(-3+\gamma_{2}(-2+\gamma_{2}+6\gamma_{2}\delta))+\gamma_{2}(-27+\gamma_{2}(27+211\gamma_{2}+54(3+\gamma_{2})\delta))}{216\gamma_{2}^{3}},\\ \sigma=&1,k=-\frac{-3+\gamma_{1}+\gamma_{2}}{6\gamma_{2}},\end{aligned}

where $- γ_{1} - γ_{2} > 0, ρ^{2} - (A_{1}^{2} + A_{2}^{2}) δ^{2} > 0$ \rho^{2}-(A_{1}^{2}+A_{2}^{2})\delta^{2}>0, and γ₂ ≠ 0. Substituting the values of Eq.(59) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows:

\begin{matrix} q (x, t) = & (\frac{\sqrt{\frac{3}{2}} \sqrt{- ρ^{2} + (A_{1}^{2} + A_{2}^{2}) δ^{2}}}{\sqrt{- (γ_{1} + γ_{2}) δ} (\frac{ρ}{δ} + A_{1} c o s (\sqrt{δ} (x + D_{1} t)) + A_{2} s i n (\sqrt{δ} (x + D_{1} t)))} \\ + \frac{\sqrt{\frac{3}{2}} (A_{2} \sqrt{δ} c o s (\sqrt{δ} (x + D_{1} t)) - A_{1} \sqrt{δ} s i n (\sqrt{δ} (x + D_{1} t)))}{\sqrt{- γ_{1} - γ_{2}} (\frac{ρ}{δ} + A_{1} c o s (\sqrt{δ} (x + D_{1} t)) + A_{2} s i n (\sqrt{δ} (x + D_{1} t)))}) \\ \times e^{i (- t + \frac{x (- 3 + γ_{1} + γ_{2})}{6 γ_{2}} + \frac{t (- 27 + γ_{1}^{3} - 3 γ_{1}^{2} (3 + γ_{2}) - 9 γ_{1} (- 3 + γ_{2} (- 2 + γ_{2} + 6 γ_{2} δ)) + γ_{2} (- 27 + γ_{2} (27 + 211 γ_{2} + 54 (3 + γ_{2}) δ)))}{216 γ_{2}^{3}} + W (t))} \end{matrix}

\begin{aligned}q(x,t)=&(\frac{\sqrt{\frac{3}{2}}\sqrt{-\rho^{2}+({A_{1}}^{2}+{A_{2}}^{2})\delta^{2}}}{\sqrt{-((\gamma_{1}+\gamma_{2})\delta)}(\frac{\rho}{\delta}+A_{1}cos(\sqrt{\delta}(x+D_{1}t))+A_{2}sin(\sqrt{\delta}(x+D_{1}t)))}\\ &+\frac{\sqrt{\frac{3}{2}}(A_{2}\sqrt{\delta}cos(\sqrt{\delta}(x+D_{1}t))-A_{1}\sqrt{\delta}sin(\sqrt{\delta}(x+D_{1}t)))}{\sqrt{-\gamma_{1}-\gamma_{2}}(\frac{\rho}{\delta}+A_{1}cos(\sqrt{\delta}(x+D_{1}t))+A_{2}sin(\sqrt{\delta}(x+D_{1}t)))})\\ &\times e^{i(-t+\frac{x(-3+\gamma_{1}+\gamma_{2})}{6\gamma_{2}}+\frac{t(-27+\gamma_{1}^{3}-3\gamma_{1}^{2}(3+\gamma_{2})-9\gamma_{1}(-3+\gamma_{2}(-2+\gamma_{2}+6\gamma_{2}\delta))+\gamma_{2}(-27+\gamma_{2}(27+211\gamma_{2}+54(3+\gamma_{2})\delta)))}{216\gamma_{2}^{3}}+W(t))}\end{aligned}

where $D_{1} = \frac{9 + γ_{1}^{2} - 2 γ_{1} (3 + γ_{2}) - 3 γ_{2} (γ_{2} - 2 + 2 γ_{2} δ)}{12 γ_{2}^{2}} .$ D_{1}=\frac{(9+{\gamma_{1}}^{2}-2\gamma_{1}(3+\gamma_{2})-3\gamma_{2}(\gamma_{2}-2+2\gamma_{2}\delta))}{12\gamma_{2}^{2}}.Eq.(60) describes a singular periodic soliton solution.

Case 3: When δ = 0. By putting Eq.(50) into Eq.(18) and Eq.(20) and applying Eq.(8) and Eq.(14), the lefthand side of equation Eq.(18) and Eq.(20) transforms into a polynomial in ϕ and Ψ. Collecting each coefficient of this polynomial to zero yields in a system of algebraic equations involving a₀,a₁,B₁,ρ,A₁ and A₂. By using package program, we obtain solutions as follows:

Set 1:

\begin{matrix} a_{0} = & 0, a_{1} = a_{1}, B_{1} = a_{1} \sqrt{A_{1}^{2} - 2 A_{2} ρ}, γ_{1} = \frac{3}{2} - \frac{3}{4 a_{1}^{2}}, γ_{2} = - \frac{3}{2} - \frac{3}{4 a_{1}^{2}}, \\ σ = & \frac{\sqrt{9 ω + \sqrt{24 + 81 ω^{2}}}}{3 \sqrt{2}}, v = \frac{\sqrt{2}}{\sqrt{9 ω + \sqrt{24 + 81 ω^{2}}}}, k = - \frac{\sqrt{2}}{\sqrt{9 ω + \sqrt{24 + 81 ω^{2}}}}, \end{matrix}

\begin{aligned}a_{0}=&0,a_{1}=a_{1},B_{1}=a_{1}\sqrt{A_{1}^{2}-2A_{2}\rho},\gamma_{1}=\frac{3}{2}-\frac{3}{4a_{1}^{2}},\gamma_{2}=-\frac{3}{2}-\frac{3}{4a_{1}^{2}},\\ \sigma=&\frac{\sqrt{9\omega+\sqrt{24+81\omega^{2}}}}{3\sqrt{2}},v=\frac{\sqrt{2}}{\sqrt{9\omega+\sqrt{24+81\omega^{2}}}},k=-\frac{\sqrt{2}}{\sqrt{9\omega+\sqrt{24+81\omega^{2}}}},\end{aligned}

where $A_{1}^{2} - 2 A_{2} ρ > 0$ A_{1}^{2}-2A_{2}\rho>0, and $9 ω + \sqrt{24 + 81 ω^{2}} > 0$ 9\omega+\sqrt{24+81\omega^{2}}>0. Substituting the values of Eq.(61) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows:

\begin{matrix} q (x, t) = & (- \frac{2 a_{1} \sqrt{9 ω + \sqrt{24 + 81 ω^{2}}} (- \sqrt{2} t ρ + (A_{1} + \sqrt{A_{1}^{2} - 2 A_{2} ρ} + ρ x) \sqrt{9 ω + \sqrt{24 + 81 ω^{2}}})}{- 2 t^{2} ρ + 2 \sqrt{2} t (A_{1} + ρ x) \sqrt{9 ω + \sqrt{24 + 81 ω^{2}}} - (2 A_{2} + x (2 A_{1} + ρ x)) (9 ω + \sqrt{24 + 81 ω^{2}})}) \\ \times e^{\frac{1}{18} i (9 t ω - t \sqrt{24 + 81 ω^{2}} + \frac{18 \sqrt{2} x}{\sqrt{9 ω + \sqrt{24 + 81 ω^{2}}}} + 3 \sqrt{2} \sqrt{9 ω + \sqrt{24 + 81 ω^{2}}} W (t))} \end{matrix}

\begin{aligned}q(x,t)=&(-\frac{2a_{1}\sqrt{9\omega+\sqrt{24+81\omega^{2}}}(-\sqrt{2}t\rho+(A_{1}+\sqrt{A_{1}^{2}-2A_{2}\rho}+\rho x)\sqrt{9\omega+\sqrt{24+81\omega^{2}}})}{-2t^{2}\rho+2\sqrt{2}t(A_{1}+\rho x)\sqrt{9\omega+\sqrt{24+81\omega^{2}}}-(2A_{2}+x(2A_{1}+\rho x))(9\omega+\sqrt{24+81\omega^{2}})})\\ &\times e^{\frac{1}{18}i(9t\omega-t\sqrt{24+81\omega^{2}}+\frac{18\sqrt{2}x}{\sqrt{9\omega+\sqrt{24+81\omega^{2}}}}+3\sqrt{2}\sqrt{9\omega+\sqrt{24+81\omega^{2}}}W(t))}\end{aligned}

Eq.(62) describes a singular soliton solution.

Set 2:

\begin{matrix} a_{0} = & 0, a_{1} = \frac{1}{\sqrt{- 2 - \frac{4 γ_{2}}{3}}}, B_{1} = \frac{\sqrt{- \frac{3 A_{1}^{2}}{2} + 3 A_{2} ρ}}{\sqrt{3 + 2 γ_{2}}}, γ_{1} = 3 + γ_{2}, \\ σ = & - \frac{1}{3 k}, v = - k, ω = \frac{1 - 6 k^{4}}{9 k^{2}} \end{matrix}

\begin{aligned}a_{0}=&0,a_{1}=\frac{1}{\sqrt{-2-\frac{4\gamma_{2}}{3}}},B_{1}=\frac{\sqrt{-\frac{3A_{1}^{2}}{2}+3A_{2}\rho}}{\sqrt{3+2\gamma_{2}}},\gamma_{1}=3+\gamma_{2},\\ \sigma=&-\frac{1}{3k},v=-k,\omega=\frac{1-6k^{4}}{9k^{2}}\end{aligned}

where $- 2 - \frac{4 γ_{2}}{3} > 0$ -2-\frac{4\gamma_{2}}{3}>0 and $- \frac{3 A_{1}^{2}}{2} + 3 A_{2} ρ > 0$ -\frac{3A_{1}^{2}}{2}+3A_{2}\rho>0. Substituting the values of Eq.(63) into Eq.(50) and applying Eq.(17), the solution for Eq.(1) is obtained as follows:

\begin{matrix} q (x, t) = & (\frac{A_{1} + ρ (k t + x)}{(A_{2} + A_{1} (k t + x) + \frac{1}{2} ρ {(k t + x)}^{2}) \sqrt{- 2 - \frac{4 γ_{2}}{3}}} + \frac{\sqrt{- \frac{3 A_{1}^{2}}{2} + 3 A_{2} ρ}}{(A_{2} + A_{1} (k t + x) + \frac{1}{2} ρ {(k t + x)}^{2}) \sqrt{3 + 2 γ_{2}}}) \\ \times e^{i (- \frac{t}{9 k^{2}} + \frac{(1 - 6 k^{4}) t}{9 k^{2}} - k x - \frac{W (t)}{3 k})} \end{matrix}

\begin{aligned}q(x,t)=&\left(\frac{A_{1}+\rho(kt+x)}{(A_{2}+A_{1}(kt+x)+\frac{1}{2}\rho(kt+x)^{2})\sqrt{-2-\frac{4\gamma_{2}}{3}}}+\frac{\sqrt{-\frac{3A_{1}^{2}}{2}+3A_{2}\rho}}{(A_{2}+A_{1}(kt+x)+\frac{1}{2}\rho(kt+x)^{2})\sqrt{3+2\gamma_{2}}}\right)\\ &\times e^{i\left(-\frac{t}{9k^{2}}+\frac{(1-6k^{4})t}{9k^{2}}-kx-\frac{W(t)}{3k}\right)}\end{aligned}

Eq.(64) describes a singular soliton solution.

Stability analysis of the solutions

The stability of the obtained traveling wave solutions was analyzed using the Hamiltonian technique [36,37].

U = \frac{1}{2} \int_{- \infty}^{\infty} {| q (x, t) |}^{2} d x,

U=\frac{1}{2}\int_{-\infty}^{\infty}|q(x,t)|^{2}dx,

where q(x,t) represents the momentum associated with the Hamiltonian system (HS). The corresponding stability criterion is given by

\frac{\partial U}{\partial v} > 0

\frac{\partial U}{\partial v}>0

The wave velocity is represented by v. By considering the parameter values γ₁=−1,γ₂=−1,γ=−1,σ=−1, A = 4, B=1,W(t)=sin(t),−10<x<10, and evaluating at t=1, we obtain:

\frac{\partial U}{\partial v} |_{v = 9} = 0.399894 > 0,

\frac{\partial U}{\partial v}|_{v=9}=0.399894>0,

which denotes that Eq.(39) is stable. And by considering the parameter values A₁=−1,A₂=0,δ=−1,ρ=3, k = 1, k=1,W(t)=sin(t),−10<x<10, and evaluating at t = 1, we obtain:

\frac{\partial U}{\partial v} |_{v = 4} = - 0.00336642 < 0,

\frac{\partial U}{\partial v}|_{v=4}=-0.00336642<0,

which denotes that Eq.(52) is unstable. However, other solutions' stability can be tested in a similar way.

Graphical representation of the results

This section presents 3D, 2D, and contour graphs for various obtained solutions to clarify the characteristics of the retrieved solutions. Eq.(23) possesses a singular periodic solution and illustrated in Figure 1 for σ = –1,k = 1, γ₁ = –1, A = 4, B = 1, C = 5, W(t) = sin(t). Eq.(24) possesses a dark soliton solution and illustrated in Figure 2, when σ = 1,k = –1, γ₁ = –1, A = 4, B = 1, C = 4, W(t) = sin(t). Eq.(28) possesses a singular soliton solution and illustrated in Figure 3, when σ = –1,k = 1, γ₁ = –1, W(t) = sin(t). A dark soliton solution for equation (56) with parameters k = 4, γ₂ = –4,A₂ = 1,δ = –1,ρ = 3, W(t) = sin(t) is illustrated in Figure 4. Eq.(60) possesses a singular periodic solution and illustrated in Figure 5, when γ₁ = –1 , γ₂ = –1,A₁ = 1,A₂ = 3,δ = 1,ρ = 1, W(t) = sin(t). Eq.(62) possesses a singular soliton solution, illustrated in Figure 6, when A₁ = –1,A₂ = 0,ρ = –10, a₁, ω = 1, W(t) = sin(t). These investigations demonstrate that the evolution of soliton solutions over time in optical fibers is influenced by the effect of dispersion. By selecting specific integer values for the system parameters, one can control the emergence and behavior of these solutions, which retain their shape and velocity during propagation.

Results and discussions

The $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2})-expansion method and the (G’/G,1/G)-expansion method are applied to the SGNLSE with Kerr effect and higher-order nonlinearity. The exact solutions obtained include dark solitons, singular periodic solutions, and singular solitons. Figures 1-6 show 3D and 2D plots of these solutions with appropriate parameter values. Eq.(23) gives a singular periodic solution shown in Figure (1a) as a 3D plot. Figure (1b) shows its time evolution at t=0,5,10, where the soliton shifts left over time. This solution is useful in signal processing and nonlinear optics. Eq.(24) describes dark soliton solutions shown in Figure (2a) as a 3D plot. Figure (2b) shows their evolution at t=0,30,60, with the soliton moving left. Dark solitons are localized dips in background intensity with a phase jump, maintaining shape and velocity during propagation. Eq.(28) represents singular soliton solutions shown in Figure (3a). Figure (3b) shows evolution at t=0,5,10, with soliton shifting right. Singular solitons have infinite peak intensity and help model extreme wave behavior. Eq.(56) also gives dark soliton solutions shown in Figure (4a). Figure (4b) shows evolution at t=0,5,10, soliton shifts right. Eq.(60) is a singular periodic solution shown in Figure (5a). Figure (5b) shows time evolution at t=0,5,10, soliton shifts left. Eq.(62) shows singular soliton solutions in Figure (6a), with evolution at t=0,20,40 in Figure (6b). Soliton shifts right over time.

In this study, the $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2})-expansion and the (G’/G,1/G)-expansion methods are employed to analyze the SGNLSE with Kerr effect and higher-order nonlinearity. A comparison with previously reported results [23] highlights the novelty of the present findings. The derived solutions consist of dark, singular, and singular periodic soliton solutions, which broaden the range of known solution types by introducing generalized trigonometric, hyperbolic, and rational forms. These results represent new contributions that have not been previously reported in literature.

Conclusions

In this work, we investigated the SGNLSE incorporating the Kerr effect and higher-order nonlinearity by employing the $t a n (\frac{ϖ (ξ)}{2})$ tan(\frac{\varpi(\xi)}{2})-expansion method and the (G’/G,1/G)-expansion method. A diverse set of exact solutions was successfully derived, including singular, singular periodic, and dark soliton solutions. To visualize the dynamics, we included two- and three-dimensional graphical representations that illustrate how nonlinear effects influence the behavior and evolution of these solutions over time. In particular, the figures clearly demonstrate the impact of the time variable on wave structure and stability. Overall, the findings offer new insights into the dynamics governed by NLPDEs and further establish the effectiveness of the proposed methods in mathematical physics.

On the soliton solutions of the generalized stochastic nonlinear Schrödinger equation with Kerr effect and higher order nonlinearity via two analytical methods

Full Article

Paradigm

My account