Open AccessArticle

Exploring Families of Solitary Wave Solutions for the Fractional Coupled Higgs System Using Modified Extended Direct Algebraic Method

Muhammad Bilal

¹,

Javed Iqbal

^1,*,

Rashid Ali

^2,*

Fuad A. Awwad

and

Emad A. A. Ismail

Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan

School of Mathematics and Statistics, Central South University, Changsha 410083, China

Department of Quantitative Analysis, College of Business Administration, King Saud University, P.O. Box 71115, Riyadh 11587, Saudi Arabia

Authors to whom correspondence should be addressed.

Fractal Fract. 2023, 7(9), 653; https://doi.org/10.3390/fractalfract7090653

Submission received: 6 August 2023 / Revised: 24 August 2023 / Accepted: 28 August 2023 / Published: 30 August 2023

(This article belongs to the Special Issue Fractional Differential Equations: Computation and Modelling with Applications)

Download

Browse Figures

Versions Notes

Abstract

In this paper, we suggest the modified Extended Direct Algebraic Method (mEDAM) to examine the existence and dynamics of solitary wave solutions in the context of the fractional coupled Higgs system, with Caputo’s fractional derivatives. The method begins with the formulation of nonlinear differential equations using a fractional complex transformation, followed by the derivation of solitary wave solutions. Two-dimensional, Three-dimensional and contour graphs are used to investigate the behavior of traveling wave solutions. The research reveals many families of solitary wave solutions as well as their deep interrelationships and dynamics. These discoveries add to a better understanding of the dynamics of the fractionally coupled Higgs system and have potential applications in areas that use nonlinear Fractional Partial Differential Equations (FPDEs).

Keywords:

FPDEs; mEDAM; variable transformation; NODE; kink waves; lump waves

1. Introduction

FPDEs are essential mathematical tools for correctly simulating systems with non-local behavior [1], anomalous diffusion [2], and long-range interactions [3]. They are also widely used in the fields of data analysis [4], image and signal processing [5] and materials research [6]. FPDEs are therefore receiving a lot of attention as a result of their capacity to give more accurate modelling of systems.

Because of their intrinsic non-local and non-linear properties, both Partial Differential Equations (PDEs) and FPDEs can pose major obstacles in their solutions. Though they are often employed, numerical approaches, such as finite difference method [7], finite element method [8] and many more [9,10,11,12,13,14,15,16,17], can be computationally costly for complicated FPDEs. When feasible, analytical approaches are chosen since they produce precise results. These techniques employ mathematical approaches to reduce the complexity of the FPDE such that it may be resolved analytically. Analytical techniques are useful because they shed light on the system’s behavior and connections between its variables. They also make it possible to examine unusual circumstances and spot singularities that numerical approaches could overlook. In general, analytical approaches are suggested for solving FPDEs due to their precision, effectiveness, and capacity to offer insightful information about the system under study.

Therefore, the Laplace transform method [18], Fourier transform method [19], Adomian decomposition method [20], direct algebraic method [21], homotopy perturbation method [22], variational iteration method [23] and other analytical techniques have all been employed by researchers to solve FPDEs. Each approach has benefits and drawbacks, and the selection of approach relies on the particular topic being investigated. However, where possible, analytical approaches are preferable over numerical methods because they offer precise solutions and insightful information about how the system behaves.

To construct solitary wave solutions for PDEs and FPDEs, one of the novel analytical methods is mEDAM. The technique entails transforming the FPDE into a nonlinear ODE using the appropriate transformation, which is then resolved using a series form solution. The nonlinear ODE that results from this solution is then utilized to create a set of algebraic equations, which are then solved to get the solitary wave solutions for the FPDE. A soliton, commonly referred to as a solitary wave, is a self-sustaining wave-like solution to certain nonlinear PDEs and FPDEs that retains its shape and speed while traveling with no alteration to its form. Among the different analytical approaches for developing soliton solutions, the mEDAM stands out as the most successful, producing a wider range of soliton solution families. For instance, Sayed et al., have employed mEDAM to construct a number of families of soliton solutions for the three different types of Tzitzeica type PDEs arising in nonlinear optics [24]. Similarly, Yasmin et al., have investigated 32 different families of symmetric soliton solutions for the fractional coupled Konno–Onno system in [25] and 33 families of 131 optical soliton solutions for a fractional perturbed Radhakrishnan-Kundu–Lakshmanan model in [26], respectively, using mEDAM.

The mEDAM is employed in this work to examine solitary wave solutions for the fractional coupled Higgs system. In a relativistic quantum field theory, the fractional coupled Higgs system is a mathematical model that describes the dynamics of two scalar fields. In this system, the interaction of the two scalar fields results in intricate behavior. We considered this system in fractional form to characterize anomalous diffusion and non-local behavior, which are known to occur in many physical systems. By accounting for non-local influences on the system, fractional derivatives allow for more detailed models of solitons’ behavior. The mathematical form of this system is given as follows [27]:

\{\begin{matrix} D_{t}^{γ} (D_{t}^{γ} u) - D_{x}^{δ} (D_{x}^{δ} u) + {| u |}^{2} u - 2 u v = 0, \\ D_{t}^{γ} (D_{t}^{γ} v) + D_{x}^{δ} (D_{x}^{δ} v) - D_{x}^{δ} (D_{x}^{δ} {| u |}^{2}) = 0 . \end{matrix}

(1)

where

0 < γ, δ \leq 1

and

u = u (x, t)

, and

v = v (x, t)

. It should be noted that the proposed model in [27] is time-fractional; however, we have generalized the model by replacing the space derivative by Caputo’s fractional derivative. The system is composed of two equations: the first represents how the amplitude and phase of the soliton have changed over time, and the second defines how the background field has changed over time. The system shows a wide range of events, including solitons production, scattering, and formation. The fractional Higgs system is especially important in particle physics since it controls particle interactions and behavior, including how particles are assigned mass. As a result of the system’s symmetry breakdown, topological solitons, and nontrivial scattering amplitudes, much theoretical and experimental study has been conducted.

Several researchers have already used various analytical methods to investigate the coupled Higgs system in integer and fractional orders. The Sin–Gorden method was utilized by Rezazadeh et al. in [27] to solve the coupled Higgs system in time-fractional form. Jabbari et al., obtained precise solutions to the coupled Higgs system (in integer order) and the Maccari system using He’s semi-inverse technique and a simple

(\frac{G^{'}}{G})

-expansion method [28]. Similar to this, Atas et al. [29] identified invariant optical soltons as coupled-Higgs equation solutions. The exp-expansion method was utilized by Seadway et al. [30] to get traveling wave solutions for the given system. Mu and Qin provided Rogue Waves for the Coupled Schrödinger–Boussinesq Equation and the Coupled Higgs Equation using the Hirota technique [31].

The fractional derivatives present in (1) are defined using Caputo’s derivative operator. This operator with some of its properties is defined as follows [32]:

^{C} D_{t}^{α} f (t) = \frac{1}{Γ (n - α)} \int_{0}^{t} \frac{\frac{\partial^{n} f (s)}{\partial s^{n}}}{{(t - s)}^{α - n + 1}} d s,

(2)

D_{y}^{α} y^{p} = \frac{Γ (1 + p)}{Γ (1 + p - α)} y^{p - α},

(3)

D_{y}^{α} g [h (y)] = g_{h}^{'} (h (y)) D_{y}^{α} h (y) = D_{h}^{α} g (h (y)) {[h^{'} (y)]}^{α} 0 < α \leq 1,

(4)

where

α

is the order of the derivative,

g (t)

and

h (t)

are the functions to be differentiated, n is the smallest integer greater than or equal to

α

, and

Γ (\cdot)

is the Gamma function.

2. Methodology of mEDAM

This section outlines the modified EDAM technique. Consider the FPDE of the form [25,26]:

P (w, \partial_{t}^{α} w, \partial_{y_{1}}^{β} w, \partial_{y_{2}}^{γ} w, w^{2}, \dots) = 0, 0 < α, β, γ \leq 1,

(5)

where w is a function of

y_{1}, y_{2}, y_{3}, \dots, y_{r}

and t.

The below steps are used to solve Equation (5):

1.: We start by transforming the variables $w (y_{1}, y_{2}, y_{3}, l d o t s, y_{r})$ into $W (x i)$ , where $ξ$ is specified in a variety of ways. Equation (5) is changed by this transformation into a nonlinear ODE of the following form:

$Q (W, W^{'} W, W^{'}, \dots) = 0,$

(6)

where W in (6) has derivatives with respect to $ξ$ . Equation (6) may be integrated once or more times to obtain the constant(s) of integration.
2.: Next, we assume that (6) has the following solution:

$W (ξ) = \sum_{l = - m}^{m} d_{l} {(Z (ξ))}^{l},$

(7)

where $d_{l} (l = - m, \dots, 0, 1, 2, \dots, m)$ are constants to be determined, and $Z (ξ)$ is the general solution of the following ODE:

$Z^{'} (ξ) = l n (A) (a + b Z (ξ) + c {(Z (ξ))}^{2}),$

(8)

where $a, b$ and c are constants and $A \neq 0, 1$ .
3.: Establishing the homogeneous balance between the greatest nonlinear term and the highest-order derivative in (6) yields the positive integer m given in (7). More precisely the balance number can be estimated by using the given two formulas [33]:

$D (\frac{d^{k} W}{d ξ^{k}}) = m + k and D (W^{j} {(\frac{d^{k} W}{d ξ^{k}})}^{l}) = m j + l (k + m),$

where D stands for degree of $W (ξ)$ as $D [W (ξ)] = m$ and $j, k$ and l are positive integers.
4.: Then, we put (7) into (6) or the equation that comes from integrating (6), and arrange all of the terms of $Z (ξ)$ in the same order. The subsequent polynomial’s coefficients are then all set to zero, resulting in a system of algebraic equations for $d_{l} (l = - m, \dots, 0, 1, 2, \dots, m)$ and additional parameters.
5.: To solve this set of algebraic equations, we use MAPLE.
6.: The analytical solutions to (5) are then obtained by determining the unknown values and inserting them into (7) together with the $Z (ξ)$ (solution of Equation (8)). Using the general solution of (8), we may produce the following families of solutions.

Family 1.

For

M < 0

and

c \neq 0

Z_{1} (ξ) = - \frac{b}{2 c} + \frac{\sqrt{- M} \tan_{A} (1 / 2 \sqrt{- M} ξ)}{2 c},

Z_{2} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ)}{2 c},

Z_{3} (ξ) = - \frac{b}{2 c} + \frac{\sqrt{- M} (\tan_{A} (\sqrt{- M} ξ) \pm (\sqrt{p q} \sec_{A} (\sqrt{- M} ξ)))}{2 c},

Z_{4} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{- M} (\cot_{A} (\sqrt{- M} ξ) \pm (\sqrt{p q} \csc_{A} (\sqrt{- M} ξ)))}{2 c},

and

Z_{5} (ξ) = - \frac{b}{2 c} + \frac{\sqrt{- M} (\tan_{A} (1 / 4 \sqrt{- M} ξ) - \cot_{A} (1 / 4 \sqrt{- M} ξ))}{4 c} .

Family 2.

For

M > 0

and

c \neq 0

Z_{6} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{M} \tanh_{A} (1 / 2 \sqrt{M} ξ)}{2 c},

Z_{7} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{M} \coth_{A} (1 / 2 \sqrt{M} ξ)}{2 c},

Z_{8} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{Z} (\tanh_{A} (\sqrt{M} ξ) \pm (\sqrt{p q} {sech}_{A} (\sqrt{M} ξ)))}{2 c},

Z_{9} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{M} (\coth_{A} (\sqrt{M} ξ) \pm (\sqrt{p q} {csch}_{A} (\sqrt{M} ξ)))}{2 c},

and

Z_{10} (ξ) = - \frac{b}{2 c} - \frac{\sqrt{M} (\tanh_{A} (1 / 4 \sqrt{M} ξ) - \coth_{A} (1 / 4 \sqrt{M} ξ))}{4 c} .

Family 3.

For

a c > 0

and

b = 0

Z_{11} (ξ) = \sqrt{\frac{a}{c}} \tan_{A} (\sqrt{a c} ξ),

Z_{12} (ξ) = - \sqrt{\frac{a}{c}} \cot_{A} (\sqrt{a c} ξ),

Z_{13} (ξ) = \sqrt{\frac{a}{c}} (\tan_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \sec_{A} (2 \sqrt{a c} ξ))),

Z_{14} (ξ) = - \sqrt{\frac{a}{c}} (\cot_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{a c} ξ))),

and

Z_{15} (ξ) = \frac{1}{2} \sqrt{\frac{a}{c}} (\tan_{A} (1 / 2 \sqrt{a c} ξ) - \cot_{A} (1 / 2 \sqrt{a c} ξ)) .

Family 4.

For

a c < 0

and

b = 0

, then

Z_{16} (ξ) = - \sqrt{- \frac{a}{c}} \tanh_{A} (\sqrt{- a c} ξ),

Z_{17} (ξ) = - \sqrt{- \frac{a}{c}} \coth_{A} (\sqrt{- a c} ξ)

Z_{18} (ξ) = - \sqrt{- \frac{a}{c}} (\tanh_{A} (2 \sqrt{- a c} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- a c} ξ))),

Z_{19} (ξ) = - \sqrt{- \frac{a}{c}} (\coth_{A} (2 \sqrt{- a c} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- a c} ξ))),

and

Z_{20} (ξ) = - \frac{1}{2} \sqrt{- \frac{a}{c}} (\tanh_{A} (1 / 2 \sqrt{- a c} ξ) + \coth_{A} (1 / 2 \sqrt{- a c} ξ)) .

Family 5.

For

c = a

and

b = 0

Z_{21} (ξ) = \tan_{A} (a ξ),

Z_{22} (ξ) = - \cot_{A} (a ξ),

Z_{23} (ξ) = \tan_{A} (2 a ξ) \pm (\sqrt{p q} \sec_{A} (2 a ξ)),

Z_{24} (ξ) = - \cot_{A} (2 a ξ) \pm (\sqrt{p q} \csc_{A} (2 a ξ)),

and

Z_{25} (ξ) = \frac{1}{2} \tan_{A} (1 / 2 a ξ) - 1 / 2 \cot_{A} (1 / 2 a ξ) .

Family 6.

For

c = - a

and

b = 0

Z_{26} (ξ) = - \tanh_{A} (a ξ),

Z_{27} (ξ) = - \coth_{A} (a ξ),

Z_{28} (ξ) = - \tanh_{A} (2 a ξ) \pm (i \sqrt{p q} {sech}_{A} (2 a ξ)),

Z_{29} (ξ) = - \coth_{A} (2 a ξ) \pm (\sqrt{p q} {csch}_{A} (2 a ξ)),

and

Z_{30} (ξ) = - \frac{1}{2} \tanh_{A} (1 / 2 a ξ) - 1 / 2 \coth_{A} (1 / 2 a ξ) .

Family 7.

For

M = 0

Z_{31} (ξ) = - 2 \frac{a (b ξ \ln (A) + 2)}{b^{2} ξ \ln (A)} .

Family 8.

For

b = λ

a = n λ (n \neq 0)

and

c = 0

Z_{32} (ξ) = A^{λ ξ} - n .

Family 9.

For

b = c = 0

Z_{33} (ξ) = a ξ \ln (A) .

Family 10.

For

b = a = 0

Z_{34} (ξ) = - \frac{1}{c ξ \ln (A)} .

Family 11.

For

a = 0

b \neq 0

and

c \neq 0

Z_{35} (ξ) = - \frac{p b}{c (\cosh_{A} (b ξ) - \sinh_{A} (b ξ) + p)},

and

Z_{36} (ξ) = - \frac{b (\cosh_{A} (b ξ) + \sinh_{A} (b ξ))}{c (\cosh_{A} (b ξ) + \sinh_{A} (b ξ) + q)} .

Family 12.

For

b = λ

c = n λ (n \neq 0)

and

a = 0

Z_{37} (ξ) = \frac{p A^{λ ξ}}{p - n q A^{λ ξ}} .

where

M = b^{2} - 4 a c

and

p, q > 0

and are referred to as deformation parameters. The generalized trigonometric and hyperbolic functions are described as follows:

\begin{matrix} \cos_{A} (ξ) = \frac{p A^{i ξ} + q A^{- i ξ}}{2}, \sin_{A} (ξ) = \frac{p A^{i ξ} - q A^{- i ξ}}{2 i}, \\ \cot_{A} (ξ) = \frac{\cos_{A} (ξ)}{\sin_{A} (ξ)}, \tan_{A} (ξ) = \frac{\sin_{A} (ξ)}{\cos_{A} (ξ)}, \\ \csc_{A} (ξ) = \frac{1}{\sin_{A} (ξ)}, \sec_{A} (ξ) = \frac{1}{\cos_{A} (ξ)} . \end{matrix}

Similarly,

\begin{matrix} \cosh_{A} (ξ) = \frac{p A^{ξ} + q A^{- ξ}}{2}, \sinh_{A} (ξ) = \frac{p A^{ξ} - q A^{- ξ}}{2}, \\ \coth_{A} (ξ) = \frac{\cosh_{A} (ξ)}{\sinh_{A} (ξ)}, \tanh_{A} (ξ) = \frac{\sinh_{A} (ξ)}{\cosh_{A} (ξ)}, \\ {csch}_{A} (ξ) = \frac{1}{\sinh_{A} (ξ)}, {sech}_{A} (ξ) = \frac{1}{\cosh_{A} (ξ)} . \end{matrix}

3. Execution of the Approach

In the present section, we use the suggested method for solving the fractional coupled Higgs system (1). We offer a transformation so that this approach may be expanded to solve (1):

\begin{matrix} u (x, t) = U (ξ) e^{i θ}, v (x, t) = V (ξ), \\ w h e r e ξ = \frac{x^{δ}}{Γ (1 + δ)} + \frac{k t^{γ}}{Γ (1 + γ)}, θ = \frac{a_{1} x^{δ}}{Γ (1 + δ)} + \frac{b_{1} t^{γ}}{Γ (1 + γ)}, \end{matrix}

(9)

in (9), the soliton’s phase component is denoted by

θ (x, t)

, its wavenumber is denoted by

b_{1}

, and its velocity is denoted by

a_{1}

. When (9) is substituted into (1), we have

\begin{matrix} (k^{2} - 1) U^{''} + (a_{1}^{2} - b_{1}^{2}) U - 2 U V + U^{3} = 0, \\ (k^{2} + 1) V^{''} - 2 {(U^{'})}^{2} - 2 U U^{''} = 0 . \end{matrix}

(10)

Ignoring the integration constant and integrating the second equation in (10), we get

V (1 + k^{2}) = U^{2} .

(11)

Putting (11) into the first equation of (10) and integrating gives us

(k^{4} - 1) U^{''} + (a_{1}^{2} - b_{1}^{2}) (k^{2} + 1) U + (k^{2} - 1) U^{3} = 0,

(12)

where prime represents differentiation with respect to

ξ

m = 1

is obtained by balancing the terms of

U^{''}

and

U^{3}

in (12).

Implementation of mEDAM to the Problem

To solve the NODE in (12) generated from the fractional coupled Higgs system using mEDAM, put

m = 1

in (7). Thus, we have

U (ξ) = \sum_{i = - 1}^{1} d_{i} {(Z (ξ))}^{i} = d_{- 1} {(Z (ξ))}^{- 1} + d_{0} + d_{1} {(Z (ξ))}^{1},

(13)

where

d_{- 1}

d_{0}

and

d_{1}

are coefficients that will be determined.

With the help of (8), we insert (13) into (12), which creates a polynomial in

Z (ξ)

by collecting the terms with the same power of

Z (ξ)

. By equating the coefficients of the polynomial to zero, we get a system of nonlinear algebraic equations. We use Maple to solve the system, leading to the identification of the following two distinct cases of solutions:

Case 1:

\begin{matrix} d_{1} = 2 \sqrt{\frac{- M {(\ln (A))}^{2} - Ω}{M}} c, d_{- 1} = 0, d_{0} = \sqrt{\frac{- M {(\ln (A))}^{2} - Ω}{M}} b, k = \sqrt{\frac{M {(\ln (A))}^{2} + 2 Ω}{M {(\ln (A))}^{2}}} . \end{matrix}

(14)

Case 2:

\begin{matrix} d_{1} = 0, d_{- 1} = 2 \sqrt{\frac{- M {(\ln (A))}^{2} - Ω}{M}} a, d_{0} = \sqrt{\frac{- M {(\ln (A))}^{2} - Ω}{M}} b, k = \sqrt{\frac{M {(\ln (A))}^{2} + 2 Ω}{M {(\ln (A))}^{2}}}, \end{matrix}

(15)

where

Ω = a_{1}^{2} - b_{1}^{2}

and

M = b^{2} - 4 a c

Assuming Case 1, we get the following families of solutions:

Family 13.

When

M < 0

and

a, b, c \neq 0

, then the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{1} (x, t) & = \sqrt{- \frac{Z {(\ln (A))}^{2} + Ω}{Z}} \sqrt{- Z} \tan_{A} (1 / 2 \sqrt{- Z} ξ) . \end{matrix}

(16)

Using (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{1} (x, t) & = e^{i θ} (\sqrt{- \frac{Z {(\ln (A))}^{2} + Ω}{Z}} \sqrt{- Z} \tan_{A} (1 / 2 \sqrt{- Z} ξ)), \end{matrix}

(17)

and

\begin{matrix} v_{1} (x, t) & = \frac{1}{1 + k^{2}} {(\sqrt{- \frac{Z {(\ln (A))}^{2} + Ω}{Z}} \sqrt{- Z} \tan_{A} (1 / 2 \sqrt{- Z} ξ))}^{2}, \end{matrix}

(18)

\begin{matrix} U_{2} (x, t) & = - \sqrt{- \frac{M {(\ln (A))}^{2} + Ω}{M}} \sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ) . \end{matrix}

(19)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{2} (x, t) & = e^{i θ} (- \sqrt{- \frac{M {(\ln (A))}^{2} + Ω}{M}} \sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ)), \end{matrix}

(20)

and

\begin{matrix} v_{2} (x, t) & = \frac{1}{1 + k^{2}} {(- \sqrt{- \frac{M {(\ln (A))}^{2} + Ω}{M}} \sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ))}^{2}, \end{matrix}

(21)

\begin{matrix} U_{3} (x, t) & = \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sin_{A} (\sqrt{- M} ξ) + \sqrt{p q})}{\cos_{A} (\sqrt{- M} ξ)} . \end{matrix}

(22)

Using (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{3} (x, t) & = e^{i θ} (\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sin_{A} (\sqrt{- M} ξ) + \sqrt{p q})}{\cos_{A} (\sqrt{- M} ξ)}), \end{matrix}

(23)

and

\begin{matrix} v_{3} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sin_{A} (\sqrt{- M} ξ) + \sqrt{p q})}{\cos_{A} (\sqrt{- M} ξ)})}^{2}, \end{matrix}

(24)

\begin{matrix} U_{4} (x, t) & = - \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\cos_{A} (\sqrt{- M} ξ) + \sqrt{p q})}{\sin_{A} (\sqrt{- M} ξ)} . \end{matrix}

(25)

Using (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{4} (x, t) & = e^{i θ} (- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\cos_{A} (\sqrt{- M} ξ) + \sqrt{p q})}{\sin_{A} (\sqrt{- M} ξ)}), \end{matrix}

(26)

and

\begin{matrix} v_{4} (x, t) & = \frac{1}{1 + k^{2}} {(- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\cos_{A} (\sqrt{- M} ξ) + \sqrt{p q})}{\sin_{A} (\sqrt{- M} ξ)})}^{2}, \end{matrix}

(27)

and

\begin{matrix} U_{5} (x, t) & = - 1 / 2 \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (- 1 + 2 {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2})}{\cos_{A} (1 / 4 \sqrt{- M} ξ) \sin_{A} (1 / 4 \sqrt{- M} ξ)} . \end{matrix}

(28)

Using (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{5} (x, t) & = e^{i θ} (- 1 / 2 \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (- 1 + 2 {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2})}{\cos_{A} (1 / 4 \sqrt{- M} ξ) \sin_{A} (1 / 4 \sqrt{- M} ξ)}), \end{matrix}

(29)

and

\begin{matrix} v_{5} (x, t) & = \frac{1}{1 + k^{2}} {(- 1 / 2 \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (- 1 + 2 {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2})}{\cos_{A} (1 / 4 \sqrt{- M} ξ) \sin_{A} (1 / 4 \sqrt{- M} ξ)})}^{2} . \end{matrix}

(30)

Family 14.

For

M > 0

and

a, b, c \neq 0

, then the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{6} (x, t) & = - \sqrt{M {(\ln (A))}^{2} + Ω} \tanh_{A} (1 / 2 \sqrt{M} ξ) . \end{matrix}

(31)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{6} (x, t) & = e^{i θ} (- \sqrt{M {(\ln (A))}^{2} + Ω} \tanh_{A} (1 / 2 \sqrt{M} ξ)), \end{matrix}

(32)

and

\begin{matrix} v_{6} (x, t) & = \frac{1}{1 + k^{2}} {(- \sqrt{M {(\ln (A))}^{2} + Ω} \tanh_{A} (1 / 2 \sqrt{M} ξ))}^{2}, \end{matrix}

(33)

\begin{matrix} U_{7} (x, t) & = - \sqrt{M {(\ln (A))}^{2} + Ω} \coth_{A} (1 / 2 \sqrt{M} ξ) . \end{matrix}

(34)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{7} (x, t) & = e^{i θ} (- \sqrt{M {(\ln (A))}^{2} + Ω} \coth_{A} (1 / 2 \sqrt{M} ξ)), \end{matrix}

(35)

and

\begin{matrix} v_{7} (x, t) & = \frac{1}{1 + k^{2}} {(- \sqrt{M {(\ln (A))}^{2} + Ω} \coth_{A} (1 / 2 \sqrt{M} ξ))}^{2}, \end{matrix}

(36)

\begin{matrix} U_{8} (x, t) & = - \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sinh_{A} (\sqrt{M} ξ) + \sqrt{p q})}{\cosh_{A} (\sqrt{M} ξ)} . \end{matrix}

(37)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{8} (x, t) & = e^{i θ} (- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sinh_{A} (\sqrt{M} ξ) + \sqrt{p q})}{\cosh_{A} (\sqrt{M} ξ)}), \end{matrix}

(38)

and

\begin{matrix} v_{8} (x, t) & = \frac{1}{1 + k^{2}} {(- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sinh_{A} (\sqrt{M} ξ) + \sqrt{p q})}{\cosh_{A} (\sqrt{M} ξ)})}^{2}, \end{matrix}

(39)

\begin{matrix} U_{9} (x, t) & = - \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\cosh_{A} (\sqrt{M} ξ) + \sqrt{p q})}{\sinh_{A} (\sqrt{M} ξ)} . \end{matrix}

(40)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{9} (x, t) & = e^{i θ} (- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\cosh_{A} (\sqrt{M} ξ) + \sqrt{p q})}{\sinh_{A} (\sqrt{M} ξ)}), \end{matrix}

(41)

and

\begin{matrix} v_{9} (x, t) & = \frac{1}{1 + k^{2}} {(- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\cosh_{A} (\sqrt{M} ξ) + \sqrt{p q})}{\sinh_{A} (\sqrt{M} ξ)})}^{2}, \end{matrix}

(42)

and

\begin{matrix} U_{10} (x, t) & = 1 / 2 \frac{\sqrt{M {(\ln (A))}^{2} + Ω}}{\cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ)} . \end{matrix}

(43)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{10} (x, t) & = e^{i θ} (1 / 2 \frac{\sqrt{M {(\ln (A))}^{2} + Ω}}{\cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ)}), \end{matrix}

(44)

and

\begin{matrix} v_{10} (x, t) & = \frac{1}{1 + k^{2}} {(1 / 2 \frac{\sqrt{M {(\ln (A))}^{2} + Ω}}{\cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ)})}^{2} . \end{matrix}

(45)

Family 15.

For

c a > 0

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{11} (x, t) = \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \tan_{A} (\sqrt{c a} ξ) . \end{matrix}

(46)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{11} (x, t) = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \tan_{A} (\sqrt{c a} ξ)), \end{matrix}

(47)

and

\begin{matrix} v_{11} (x, t) = \frac{1}{1 + k^{2}} {(\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \tan_{A} (\sqrt{c a} ξ))}^{2}, \end{matrix}

(48)

\begin{matrix} U_{12} (x, t) & = - \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \cot_{A} (\sqrt{c a} ξ) . \end{matrix}

(49)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{12} (x, t) = e^{i θ} (- \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \cot_{A} (\sqrt{c a} ξ)), \end{matrix}

(50)

and

\begin{matrix} v_{12} (x, t) & = \frac{1}{1 + k^{2}} {(- \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \cot_{A} (\sqrt{c a} ξ))}^{2}, \end{matrix}

(51)

\begin{matrix} U_{13} (x, t) & = \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \times \\ (\tan_{A} (2 \sqrt{c a} ξ) \pm (\sqrt{p q} \sec_{A} (2 \sqrt{c a} ξ))) . \end{matrix}

(52)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{13} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \times \\ (\tan_{A} (2 \sqrt{c a} ξ) \pm (\sqrt{p q} \sec_{A} (2 \sqrt{c a} ξ)))), \end{matrix}

(53)

and

\begin{matrix} v_{13} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \times \\ (\tan_{A} (2 \sqrt{c a} ξ) \pm (\sqrt{p q} \sec_{A} (2 \sqrt{c a} ξ))))^{2}, \end{matrix}

(54)

\begin{matrix} U_{14} (x, t) & = - \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \times \\ (\cot_{A} (2 \sqrt{c a} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{c a} ξ))) . \end{matrix}

(55)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{14} (x, t) & = e^{i θ} (- \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \times \\ (\cot_{A} (2 \sqrt{c a} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{c a} ξ)))), \end{matrix}

(56)

and

\begin{matrix} v_{14} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{- (4 {(\ln (A))}^{2} c a - Ω)} \times \\ (\cot_{A} (2 \sqrt{c a} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{c a} ξ))))^{2}, \end{matrix}

(57)

and

\begin{matrix} U_{15} (x, t) & = \frac{\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)}}{2} \times \\ (\tan_{A} (\frac{\sqrt{c a} ξ}{2}) - \cot_{A} (\frac{\sqrt{c a} ξ}{2})) . \end{matrix}

(58)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{15} (x, t) & = e^{i θ} (\frac{\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)}}{2} \times \\ (\tan_{A} (\frac{\sqrt{c a} ξ}{2}) - \cot_{A} (\frac{\sqrt{c a} ξ}{2}))), \end{matrix}

(59)

and

\begin{matrix} v_{15} (x, t) & = \frac{1}{1 + k^{2}} (\frac{\sqrt{- (4 {(\ln (A))}^{2} c a - Ω)}}{2} \times \\ (\tan_{A} (\frac{\sqrt{c a} ξ}{2}) - \cot_{A} (\frac{\sqrt{c a} ξ}{2})))^{2} . \end{matrix}

(60)

Family 16.

For

c a < 0

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{16} (x, t) = - \sqrt{4 {(\ln (A))}^{2} c a - Ω} \tanh_{A} (\sqrt{- c a} ξ) . \end{matrix}

(61)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{16} (x, t) = e^{i θ} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \tanh_{A} (\sqrt{- c a} ξ)), \end{matrix}

(62)

and

\begin{matrix} v_{16} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \tanh_{A} (\sqrt{- c a} ξ))}^{2}, \end{matrix}

(63)

\begin{matrix} U_{17} (x, t) = - \sqrt{4 {(\ln (A))}^{2} c a - Ω} \coth_{A} (\sqrt{- c a} ξ) \end{matrix} .

(64)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{17} (x, t) = e^{i θ} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \coth_{A} (\sqrt{- c a} ξ)), \end{matrix}

(65)

and

\begin{matrix} v_{17} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \coth_{A} (\sqrt{- c a} ξ))}^{2}, \end{matrix}

(66)

\begin{matrix} U_{18} (x, t) & = - \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\tanh_{A} (2 \sqrt{- c a} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- c a} ξ))) . \end{matrix}

(67)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{18} (x, t) & = e^{i θ} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\tanh_{A} (2 \sqrt{- c a} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- c a} ξ)))), \end{matrix}

(68)

and

\begin{matrix} v_{18} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\tanh_{A} (2 \sqrt{- c a} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- c a} ξ))))^{2}, \end{matrix}

(69)

\begin{matrix} U_{19} (x, t) & = - \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\coth_{A} (2 \sqrt{- c a} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- c a} ξ))) . \end{matrix}

(70)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{19} (x, t) & = e^{i θ} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\coth_{A} (2 \sqrt{- c a} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- c a} ξ)))), \end{matrix}

(71)

and

\begin{matrix} v_{19} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\coth_{A} (2 \sqrt{- c a} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- c a} ξ))))^{2}, \end{matrix}

(72)

and

\begin{matrix} U_{20} (x, t) = - \sqrt{4 {(\ln (A))}^{2} c a - Ω} (\tanh_{A} (\sqrt{- c a} ξ) + \coth_{A} (\sqrt{- c a} ξ)) . \end{matrix}

(73)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{20} (x, t) & = e^{i θ} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\tanh_{A} (\sqrt{- c a} ξ) + \coth_{A} (\sqrt{- c a} ξ))), \end{matrix}

(74)

and

\begin{matrix} v_{20} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{4 {(\ln (A))}^{2} c a - Ω} \times \\ (\tanh_{A} (\sqrt{- c a} ξ) + \coth_{A} (\sqrt{- c a} ξ)))^{2} . \end{matrix}

(75)

Family 17.

For

c = a

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{21} (x, t) = \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \tan_{A} (c ξ) . \end{matrix}

(76)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{21} (x, t) = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \tan_{A} (c ξ)), \end{matrix}

(77)

and

\begin{matrix} v_{21} (x, t) = \frac{1}{1 + k^{2}} {(\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \tan_{A} (c ξ))}^{2}, \end{matrix}

(78)

\begin{matrix} U_{22} (x, t) = - \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \cot_{A} (c ξ) . \end{matrix}

(79)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{22} (x, t) = e^{i θ} (- \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \cot_{A} (c ξ)), \end{matrix}

(80)

and

\begin{matrix} v_{22} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \cot_{A} (c ξ))}^{2}, \end{matrix}

(81)

\begin{matrix} U_{23} (x, t) & = \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} (\tan_{A} (2 c ξ) \pm (\sqrt{p q} \sec_{A} (2 c ξ))) . \end{matrix}

(82)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{23} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (\tan_{A} (2 c ξ) \pm (\sqrt{p q} \sec_{A} (2 c ξ)))), \end{matrix}

(83)

and

\begin{matrix} v_{23} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (\tan_{A} (2 c ξ) \pm (\sqrt{p q} \sec_{A} (2 c ξ))))^{2}, \end{matrix}

(84)

\begin{matrix} U_{24} (x, t) & = \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (- \cot_{A} (2 c ξ) \mp (\sqrt{p q} \csc_{A} (2 c ξ))) . \end{matrix}

(85)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{24} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (- \cot_{A} (2 c ξ) \mp (\sqrt{p q} \csc_{A} (2 c ξ)))), \end{matrix}

(86)

and

\begin{matrix} v_{24} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (- \cot_{A} (2 c ξ) \mp (\sqrt{p q} \csc_{A} (2 c ξ))))^{2}, \end{matrix}

(87)

and

\begin{matrix} U_{25} (x, t) = \sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} (\frac{1}{2} (\tan_{A} (\frac{c ξ}{2}) - \cot_{A} (\frac{c ξ}{2}))) . \end{matrix}

(88)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{25} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (\frac{1}{2} (\tan_{A} (\frac{c ξ}{2}) - \cot_{A} (\frac{c ξ}{2})))), \end{matrix}

(89)

and

\begin{matrix} v_{25} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} c^{2} - Ω)} \times \\ (\frac{1}{2} (\tan_{A} (\frac{c ξ}{2}) - \cot_{A} (\frac{c ξ}{2}))))^{2} . \end{matrix}

(90)

Family 18.

For

c = - a

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{26} (x, t) = - \sqrt{- 4 {(\ln (A))}^{2} c^{2} - Ω} \tanh_{A} (c ξ) . \end{matrix}

(91)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{26} (x, t) = e^{i θ} (- \sqrt{- 4 {(\ln (A))}^{2} c^{2} - Ω} \tanh_{A} (c ξ)), \end{matrix}

(92)

and

\begin{matrix} v_{26} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- 4 {(\ln (A))}^{2} c^{2} - Ω} \tanh_{A} (c ξ))}^{2}, \end{matrix}

(93)

\begin{matrix} U_{27} (x, t) = - \sqrt{- 4 {(\ln (A))}^{2} c^{2} - Ω} \coth_{A} (c ξ) . \end{matrix}

(94)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{27} (x, t) = e^{i θ} (- \sqrt{- 4 {(\ln (A))}^{2} c^{2} - Ω} \coth_{A} (c ξ)), \end{matrix}

(95)

and

\begin{matrix} v_{27} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- 4 {(\ln (A))}^{2} c^{2} - Ω} \coth_{A} (c ξ))}^{2}, \end{matrix}

(96)

\begin{matrix} U_{28} (x, t) = \sqrt{- {(2 c \ln (A))}^{2} - Ω} (- \tanh_{A} (2 c ξ) \mp (i \sqrt{p q} {sech}_{A} (2 c ξ))) . \end{matrix}

(97)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{28} (x, t) & = e^{i θ} (\sqrt{- {(2 c \ln (A))}^{2} - Ω} \\ (- \tanh_{A} (2 c ξ) \mp (i \sqrt{p q} {sech}_{A} (2 c ξ)))), \end{matrix}

(98)

and

\begin{matrix} v_{28} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- {(2 c \ln (A))}^{2} - Ω} \times \\ (- \tanh_{A} (2 c ξ) \mp (i \sqrt{p q} {sech}_{A} (2 c ξ))))^{2}, \end{matrix}

(99)

\begin{matrix} U_{29} (x, t) = \sqrt{- {(2 c \ln (A))}^{2} - Ω} (- \coth_{A} (2 c ξ) \mp (\sqrt{p q} {csch}_{A} (2 c ξ))) . \end{matrix}

(100)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{29} (x, t) & = e^{i θ} (\sqrt{- {(2 c \ln (A))}^{2} - Ω} \times \\ (- \coth_{A} (2 c ξ) \mp (\sqrt{p q} {csch}_{A} (2 c ξ)))), \end{matrix}

(101)

and

\begin{matrix} v_{29} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- {(2 c \ln (A))}^{2} - Ω} \times \\ (- \coth_{A} (2 c ξ) \mp (\sqrt{p q} {csch}_{A} (2 c ξ))))^{2}, \end{matrix}

(102)

and

\begin{matrix} U_{30} (x, t) = \sqrt{- {(2 c \ln (A))}^{2} - Ω} (- \frac{1}{2} (\tanh_{A} (\frac{c ξ}{2}) + \coth_{A} (\frac{c ξ}{2}))) . \end{matrix}

(103)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{30} (x, t) & = e^{i θ} (\sqrt{- {(2 c \ln (A))}^{2} - Ω} \times \\ (- \frac{1}{2} (\tanh_{A} (\frac{c ξ}{2}) + \coth_{A} (\frac{c ξ}{2})))), \end{matrix}

(104)

and

\begin{matrix} v_{30} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- {(2 c \ln (A))}^{2} - Ω} \times \\ (- \frac{1}{2} (\tanh_{A} (\frac{c ξ}{2}) + \coth_{A} (\frac{c ξ}{2}))))^{2} . \end{matrix}

(105)

Family 19.

For

a = 0

b \neq 0

and

c \neq 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{31} (x, t) & = \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \\ - 2 \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} p {(\cosh_{A} (b ξ) - \sinh_{A} (b ξ) + p)}^{- 1} . \end{matrix}

(106)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{31} (x, t) & = e^{i θ} (\sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \\ - 2 \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} p {(\cosh_{A} (b ξ) - \sinh_{A} (b ξ) + p)}^{- 1}), \end{matrix}

(107)

and

\begin{matrix} v_{31} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \\ - 2 \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} p {(\cosh_{A} (b ξ) - \sinh_{A} (b ξ) + p)}^{- 1})^{2}, \end{matrix}

(108)

and

\begin{matrix} U_{32} (x, t) & = \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \\ - 2 \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \frac{(\cosh_{A} (b ξ) + \sinh_{A} (b ξ))}{(\cosh_{A} (b ξ) + \sinh_{A} (b ξ) + q)} . \end{matrix}

(109)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{32} (x, t) & = e^{i θ} (\sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \\ - 2 \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \frac{(\cosh_{A} (b ξ) + \sinh_{A} (b ξ))}{(\cosh_{A} (b ξ) + \sinh_{A} (b ξ) + q)}), \end{matrix}

(110)

and

\begin{matrix} v_{32} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \\ - 2 \sqrt{- {(\ln (A))}^{2} b^{2} - Ω} \frac{(\cosh_{A} (b ξ) + \sinh_{A} (b ξ))}{(\cosh_{A} (b ξ) + \sinh_{A} (b ξ) + q)})^{2} . \end{matrix}

(111)

Family 20.

For

b = λ

c = n λ (n \neq 0)

and

a = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{33} (x, t) & = \sqrt{- {(\ln (A))}^{2} λ^{2} - Ω} \\ + 2 \sqrt{- {(\ln (A))}^{2} λ^{2} - Ω} n p A^{λ ξ} {(p - n q A^{λ ξ})}^{- 1} . \end{matrix}

(112)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{33} (x, t) & = e^{i θ} (\sqrt{- {(\ln (A))}^{2} λ^{2} - Ω} \\ + 2 \sqrt{- {(\ln (A))}^{2} λ^{2} - Ω} n p A^{λ ξ} {(p - n q A^{λ ξ})}^{- 1}), \end{matrix}

(113)

and

\begin{matrix} v_{33} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- {(\ln (A))}^{2} λ^{2} - Ω} \\ + 2 \sqrt{- {(\ln (A))}^{2} λ^{2} - Ω} n p A^{λ ξ} {(p - n q A^{λ ξ})}^{- 1})^{2}, \end{matrix}

(114)

where

ξ = \frac{x^{δ}}{Γ (1 + δ)} + \sqrt{\frac{b^{2} {(\ln (A))}^{2} + 2 {a_{1}}^{2} - 4 c {(\ln (A))}^{2} a - 2 {b_{1}}^{2}}{(b^{2} - 4 a c) {(\ln (A))}^{2}}} \frac{t^{γ}}{Γ (1 + γ)} a n d θ = \frac{a_{1} x^{δ}}{Γ (1 + δ)} + \frac{b_{1} t^{γ}}{Γ (1 + γ)} .

Now, assuming Case 2, we get the following families of solitary wave solutions.

Family 21.

When

M < 0

and

a, b, c \neq 0

, then the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{34} (x, t) & = \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} + b \tan_{A} (1 / 2 \sqrt{- M} ξ))}{- b + \sqrt{- M} \tan_{A} (1 / 2 \sqrt{- M} ξ)} \end{matrix}

(115)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{34} (x, t) & = e^{i θ} (\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} + b \tan_{A} (1 / 2 \sqrt{- M} ξ))}{- b + \sqrt{- M} \tan_{A} (1 / 2 \sqrt{- M} ξ)}), \end{matrix}

(116)

and

\begin{matrix} v_{34} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} + b \tan_{A} (1 / 2 \sqrt{- M} ξ))}{- b + \sqrt{- M} \tan_{A} (1 / 2 \sqrt{- M} ξ)})}^{2}, \end{matrix}

(117)

\begin{matrix} U_{35} (x, t) & = \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} + b \cot_{A} (1 / 2 \sqrt{- M} ξ))}{- b + \sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ)}, \end{matrix}

(118)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{35} (x, t) & = e^{i θ} (\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} + b \cot_{A} (1 / 2 \sqrt{- M} ξ))}{- b + \sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ)}), \end{matrix}

(119)

and

\begin{matrix} v_{35} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} + b \cot_{A} (1 / 2 \sqrt{- M} ξ))}{- b + \sqrt{- M} \cot_{A} (1 / 2 \sqrt{- M} ξ)})}^{2}, \end{matrix}

(120)

\begin{matrix} U_{36} (x, t) & = \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} \cos_{A} (\sqrt{- M} ξ) + b \sin_{A} (\sqrt{- M} ξ) + b \sqrt{p q})}{- b \cos_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sin_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sqrt{p q}}, \end{matrix}

(121)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{36} (x, t) & = e^{i θ} (\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} \cos_{A} (\sqrt{- M} ξ) + b \sin_{A} (\sqrt{- M} ξ) + b \sqrt{p q})}{- b \cos_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sin_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sqrt{p q}}), \end{matrix}

(122)

and

\begin{matrix} v_{36} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} \cos_{A} (\sqrt{- M} ξ) + b \sin_{A} (\sqrt{- M} ξ) + b \sqrt{p q})}{- b \cos_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sin_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sqrt{p q}})}^{2}, \end{matrix}

(123)

\begin{matrix} U_{37} (x, t) & = \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} \sin_{A} (\sqrt{- M} ξ) + b \cos_{A} (\sqrt{- M} ξ) + b \sqrt{p q})}{b \sin_{A} (\sqrt{- M} ξ) + \sqrt{- M} \cos_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sqrt{p q}}, \end{matrix}

(124)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{37} (x, t) & = e^{i θ} (\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} \sin_{A} (\sqrt{- M} ξ) + b \cos_{A} (\sqrt{- M} ξ) + b \sqrt{p q})}{b \sin_{A} (\sqrt{- M} ξ) + \sqrt{- M} \cos_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sqrt{p q}}), \end{matrix}

(125)

and

\begin{matrix} v_{37} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (\sqrt{- M} \sin_{A} (\sqrt{- M} ξ) + b \cos_{A} (\sqrt{- M} ξ) + b \sqrt{p q})}{b \sin_{A} (\sqrt{- M} ξ) + \sqrt{- M} \cos_{A} (\sqrt{- M} ξ) + \sqrt{- M} \sqrt{p q}})}^{2}, \end{matrix}

(126)

and

\begin{matrix} U_{38} (x, t) & = - \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (2 R + b - 2 b {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2})}{2 R - \sqrt{- M} + 2 \sqrt{- M} {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2}}, \end{matrix}

(127)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{38} (x, t) & = e^{i θ} (- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (2 R + b - 2 b {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2})}{2 b R - \sqrt{- M} + 2 \sqrt{- M} {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2}}), \end{matrix}

(128)

and

\begin{matrix} v_{38} (x, t) & = \frac{1}{1 + k^{2}} {(- \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (2 R + b - 2 b {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2})}{2 R - \sqrt{- M} + 2 \sqrt{- M} {(\cos_{A} (1 / 4 \sqrt{- M} ξ))}^{2}})}^{2}, \end{matrix}

(129)

where

R = \sqrt{- M} \cos_{A} (1 / 4 \sqrt{- M} ξ) \sin_{A} (1 / 4 \sqrt{- M} ξ)

Family 22.

For

M > 0

and

a, b, c \neq 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{39} (x, t) & = \frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} + b \tanh_{A} (1 / 2 \sqrt{M} ξ))}{b + \sqrt{M} \tanh_{A} (1 / 2 \sqrt{M} ξ)}, \end{matrix}

(130)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{39} (x, t) & = e^{i θ} (\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} + b \tanh_{A} (1 / 2 \sqrt{M} ξ))}{b + \sqrt{M} \tanh_{A} (1 / 2 \sqrt{M} ξ)}), \end{matrix}

(131)

and

\begin{matrix} v_{39} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} + b \tanh_{A} (1 / 2 \sqrt{M} ξ))}{b + \sqrt{M} \tanh_{A} (1 / 2 \sqrt{M} ξ)})}^{2}, \end{matrix}

(132)

\begin{matrix} U_{40} (x, t) & = \frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} + b \coth_{A} (1 / 2 \sqrt{M} ξ))}{b + \sqrt{M} \coth_{A} (1 / 2 \sqrt{M} ξ)}, \end{matrix}

(133)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{40} (x, t) & = e^{i θ} (\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} + b \coth_{A} (1 / 2 \sqrt{M} ξ))}{b + \sqrt{M} \coth_{A} (1 / 2 \sqrt{M} ξ)}), \end{matrix}

(134)

and

\begin{matrix} v_{40} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} + b \coth_{A} (1 / 2 \sqrt{M} ξ))}{b + \sqrt{M} \coth_{A} (1 / 2 \sqrt{M} ξ)})}^{2}, \end{matrix}

(135)

\begin{matrix} U_{41} (x, t) & = \frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} \cosh_{A} (\sqrt{M} ξ) + b \sinh_{A} (\sqrt{M} ξ) + b \sqrt{p q})}{b \cosh_{A} (\sqrt{M} ξ) + \sqrt{M} \sinh_{A} (\sqrt{M} ξ) + \sqrt{M} \sqrt{p q}}, \end{matrix}

(136)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{41} (x, t) & = e^{i θ} (\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} \cosh_{A} (\sqrt{M} ξ) + b \sinh_{A} (\sqrt{M} ξ) + b \sqrt{p q})}{b \cosh_{A} (\sqrt{M} ξ) + \sqrt{M} \sinh_{A} (\sqrt{M} ξ) + \sqrt{M} \sqrt{p q}}), \end{matrix}

(137)

and

\begin{matrix} v_{41} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} \cosh_{A} (\sqrt{M} ξ) + b \sinh_{A} (\sqrt{M} ξ) + b \sqrt{p q})}{b \cosh_{A} (\sqrt{M} ξ) + \sqrt{M} \sinh_{A} (\sqrt{M} ξ) + \sqrt{M} \sqrt{p q}})}^{2}, \end{matrix}

(138)

\begin{matrix} U_{42} (x, t) & = \frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} \sinh_{A} (\sqrt{M} ξ) + b \cosh_{A} (\sqrt{M} ξ) + b \sqrt{p q})}{b \sinh_{A} (\sqrt{M} ξ) + \sqrt{M} \cosh_{A} (\sqrt{M} ξ) + \sqrt{M} \sqrt{p q}} . \end{matrix}

(139)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{42} (x, t) & = e^{i θ} (\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} \sinh_{A} (\sqrt{M} ξ) + b \cosh_{A} (\sqrt{M} ξ) + b \sqrt{p q})}{b \sinh_{A} (\sqrt{M} ξ) + \sqrt{M} \cosh_{A} (\sqrt{M} ξ) + \sqrt{M} \sqrt{p q}}), \end{matrix}

(140)

and

\begin{matrix} v_{42} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{- M {(\ln (A))}^{2} - Ω} (\sqrt{M} \sinh_{A} (\sqrt{M} ξ) + b \cosh_{A} (\sqrt{M} ξ) + b \sqrt{p q})}{b \sinh_{A} (\sqrt{M} ξ) + \sqrt{M} \cosh_{A} (\sqrt{M} ξ) + \sqrt{M} \sqrt{p q}})}^{2}, \end{matrix}

(141)

and

\begin{matrix} U_{43} (x, t) & = \frac{\sqrt{M {(\ln (A))}^{2} + Ω} (2 \sqrt{- M} \cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ) + b)}{- 2 b \cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ) + \sqrt{M}} . \end{matrix}

(142)

Utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1) as

\begin{matrix} u_{43} (x, t) & = e^{i θ} (\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (2 \sqrt{- M} \cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ) + b)}{- 2 b \cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ) + \sqrt{M}}), \end{matrix}

(143)

and

\begin{matrix} v_{43} (x, t) & = \frac{1}{1 + k^{2}} {(\frac{\sqrt{M {(\ln (A))}^{2} + Ω} (2 \sqrt{- M} \cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ) + b)}{- 2 b \cosh_{A} (1 / 4 \sqrt{M} ξ) \sinh_{A} (1 / 4 \sqrt{M} ξ) + \sqrt{M}})}^{2} . \end{matrix}

(144)

Family 23.

For

a c > 0

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{44} (x, t) = \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} {(\tan_{A} (\sqrt{a c} ξ))}^{- 1} . \end{matrix}

(145)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{44} (x, t) = e^{i θ} (\sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} {(\tan_{A} (\sqrt{a c} ξ))}^{- 1}), \end{matrix}

(146)

and

\begin{matrix} v_{44} (x, t) = \frac{1}{1 + k^{2}} {(\sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} {(\tan_{A} (\sqrt{a c} ξ))}^{- 1})}^{2}, \end{matrix}

(147)

\begin{matrix} U_{45} (x, t) & = - \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} {(\cot_{A} (\sqrt{a c} ξ))}^{- 1} . \end{matrix}

(148)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{45} (x, t) = e^{i θ} (- \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} {(\cot_{A} (\sqrt{a c} ξ))}^{- 1}), \end{matrix}

(149)

and

\begin{matrix} v_{45} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} {(\cot_{A} (\sqrt{a c} ξ))}^{- 1})}^{2}, \end{matrix}

(150)

\begin{matrix} U_{46} (x, t) & = \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\tan_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \sec (2 \sqrt{a c} ξ)))}^{- 1} . \end{matrix}

(151)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{46} (x, t) & = e^{i θ} (\sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\tan_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \sec_{A} (2 \sqrt{a c} ξ)))}^{- 1}), \end{matrix}

(152)

and

\begin{matrix} v_{46} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\tan_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \sec_{A} (2 \sqrt{a c} ξ)))}^{- 1})^{2}, \end{matrix}

(153)

\begin{matrix} U_{47} (x, t) & = - \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\cot_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{a c} ξ)))}^{- 1} . \end{matrix}

(154)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{47} (x, t) & = e^{i θ} (- \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\cot_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{a c} ξ)))}^{- 1}), \end{matrix}

(155)

and

\begin{matrix} v_{47} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\cot_{A} (2 \sqrt{a c} ξ) \pm (\sqrt{p q} \csc_{A} (2 \sqrt{a c} ξ)))}^{- 1})^{2}, \end{matrix}

(156)

and

\begin{matrix} U_{48} (x, t) & = 2 \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\tan_{A} (\frac{\sqrt{a c} ξ}{2}) - \cot_{A} (\frac{\sqrt{a c} ξ}{2}))}^{- 1} . \end{matrix}

(157)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{48} (x, t) & = e^{i θ} (2 \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\tan_{A} (\frac{\sqrt{a c} ξ}{2}) - \cot_{A} (\frac{\sqrt{a c} ξ}{2}))}^{- 1}), \end{matrix}

(158)

and

\begin{matrix} v_{48} (x, t) & = \frac{1}{1 + k^{2}} (2 \sqrt{- (4 c {(\ln (A))}^{2} a - Ω)} \times \\ {(\tan_{A} (\frac{\sqrt{a c} ξ}{2}) - \cot_{A} (\frac{\sqrt{a c} ξ}{2}))}^{- 1})^{2} . \end{matrix}

(159)

Family 24.

For

a c < 0

and

b = 0

, then the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{49} (x, t) = - \sqrt{4 c {(\ln (A))}^{2} a - Ω} {(\tanh_{A} (\sqrt{- a c} ξ))}^{- 1} . \end{matrix}

(160)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{49} (x, t) = e^{i θ} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} {(\tanh_{A} (\sqrt{- a c} ξ))}^{- 1}), \end{matrix}

(161)

and

\begin{matrix} v_{49} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{4 c {(\ln (A))}^{2} a - Ω} {(\tanh_{A} (\sqrt{- a c} ξ))}^{- 1})}^{2}, \end{matrix}

(162)

\begin{matrix} U_{50} (x, t) = - \sqrt{4 c {(\ln (A))}^{2} a - Ω} {(\coth_{A} (\sqrt{- a c} ξ))}^{- 1} . \end{matrix}

(163)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{50} (x, t) = e^{i θ} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} {(\coth_{A} (\sqrt{- a c} ξ))}^{- 1}), \end{matrix}

(164)

and

\begin{matrix} v_{50} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{4 c {(\ln (A))}^{2} a - Ω} {(\coth_{A} (\sqrt{- a c} ξ))}^{- 1})}^{2}, \end{matrix}

(165)

\begin{matrix} U_{51} (x, t) & = - \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\tanh_{A} (2 \sqrt{- a c} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- a c} ξ)))}^{- 1} . \end{matrix}

(166)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{51} (x, t) & = e^{i θ} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\tanh_{A} (2 \sqrt{- a c} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- a c} ξ)))}^{- 1}), \end{matrix}

(167)

and

\begin{matrix} v_{51} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\tanh_{A} (2 \sqrt{- a c} ξ) \pm (i \sqrt{p q} {sech}_{A} (2 \sqrt{- a c} ξ)))}^{- 1})^{2}, \end{matrix}

(168)

\begin{matrix} U_{52} (x, t) & = - \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\coth_{A} (2 \sqrt{- a c} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- a c} ξ)))}^{- 1} . \end{matrix}

(169)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{52} (x, t) & = e^{i θ} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\coth_{A} (2 \sqrt{- a c} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- a c} ξ)))}^{- 1}), \end{matrix}

(170)

and

\begin{matrix} v_{52} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\coth_{A} (2 \sqrt{- a c} ξ) \pm (\sqrt{p q} {csch}_{A} (2 \sqrt{- a c} ξ)))}^{- 1})^{2}, \end{matrix}

(171)

and

\begin{matrix} U_{53} (x, t) & = - \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\tanh_{A} (\frac{\sqrt{- a c} ξ}{2}) + \coth_{A} (\frac{\sqrt{- a c} ξ}{2}))}^{- 1} . \end{matrix}

(172)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{53} (x, t) & = e^{i θ} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\tanh_{A} (\frac{\sqrt{- a c} ξ}{2}) + \coth_{A} (\frac{\sqrt{- a c} ξ}{2}))}^{- 1}), \end{matrix}

(173)

and

\begin{matrix} v_{53} (x, t) & = \frac{1}{1 + k^{2}} (- \sqrt{4 c {(\ln (A))}^{2} a - Ω} \times \\ {(\tanh_{A} (\frac{\sqrt{- a c} ξ}{2}) + \coth_{A} (\frac{\sqrt{- a c} ξ}{2}))}^{- 1})^{2} . \end{matrix}

(174)

Family 25.

For

c = a

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{54} (x, t) = \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\tan_{A} (a ξ))}^{- 1}, \end{matrix}

(175)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{54} (x, t) = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\tan_{A} (a ξ))}^{- 1}), \end{matrix}

(176)

and

\begin{matrix} v_{54} (x, t) = \frac{1}{1 + k^{2}} {(\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\tan_{A} (a ξ))}^{- 1})}^{2}, \end{matrix}

(177)

\begin{matrix} U_{55} (x, t) = - \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\cot_{A} (a ξ))}^{- 1} . \end{matrix}

(178)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{55} (x, t) = e^{i θ} (- \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\cot_{A} (a ξ))}^{- 1}), \end{matrix}

(179)

and

\begin{matrix} v_{55} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\cot_{A} (a ξ))}^{- 1})}^{2}, \end{matrix}

(180)

\begin{matrix} U_{56} (x, t) = \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} {(\tan_{A} (2 a ξ) \pm (\sqrt{p q} \sec_{A} (2 a ξ)))}^{- 1} . \end{matrix}

(181)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{56} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(\tan_{A} (2 a ξ) \pm (\sqrt{p q} \sec_{A} (2 a ξ)))}^{- 1}), \end{matrix}

(182)

and

\begin{matrix} v_{56} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(\tan_{A} (2 a ξ) \pm (\sqrt{p q} \sec_{A} (2 a ξ)))}^{- 1})^{2}, \end{matrix}

(183)

\begin{matrix} U_{57} (x, t) & = \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(- \cot_{A} (2 a ξ) \mp (\sqrt{p q} \csc_{A} (2 a ξ)))}^{- 1} . \end{matrix}

(184)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{57} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(- \cot_{A} (2 a ξ) \mp (\sqrt{p q} \csc_{A} (2 a ξ)))}^{- 1}), \end{matrix}

(185)

and

\begin{matrix} v_{57} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(- \cot_{A} (2 a ξ) \mp (\sqrt{p q} \csc_{A} (2 a ξ)))}^{- 1})^{2}, \end{matrix}

(186)

and

\begin{matrix} U_{58} (x, t) & = \sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(\frac{1}{2} (\tan_{A} (\frac{a ξ}{2}) - \cot_{A} (\frac{a ξ}{2})))}^{- 1} . \end{matrix}

(187)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{58} (x, t) & = e^{i θ} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(\frac{1}{2} (\tan_{A} (\frac{a ξ}{2}) - \cot_{A} (\frac{a ξ}{2})))}^{- 1}), \end{matrix}

(188)

and

\begin{matrix} v_{58} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- (4 {(\ln (A))}^{2} a^{2} - Ω)} \times \\ {(\frac{1}{2} (\tan_{A} (\frac{a ξ}{2}) - \cot_{A} (\frac{a ξ}{2})))}^{- 1})^{2} . \end{matrix}

(189)

Family 26.

For

c = - a

and

b = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{59} (x, t) = - \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} {(\tanh_{A} (a ξ))}^{- 1} . \end{matrix}

(190)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{59} (x, t) = e^{i θ} (- \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} {(\tanh_{A} (a ξ))}^{- 1}), \end{matrix}

(191)

and

\begin{matrix} v_{59} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} {(\tanh_{A} (a ξ))}^{- 1})}^{2}, \end{matrix}

(192)

\begin{matrix} U_{60} (x, t) = - \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} {(\coth_{A} (a ξ))}^{- 1} . \end{matrix}

(193)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{60} (x, t) = e^{i θ} (- \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} {(\coth_{A} (a ξ))}^{- 1}), \end{matrix}

(194)

and

\begin{matrix} v_{60} (x, t) = \frac{1}{1 + k^{2}} {(- \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} {(\coth_{A} (a ξ))}^{- 1})}^{2}, \end{matrix}

(195)

\begin{matrix} U_{61} (x, t) & = \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(- \tanh_{A} (2 a ξ) \mp (i \sqrt{p q} {sech}_{A} (2 a ξ)))}^{- 1} . \end{matrix}

(196)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{61} (x, t) & = e^{i θ} (\sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(- \tanh_{A} (2 a ξ) \mp (i \sqrt{p q} {sech}_{A} (2 a ξ)))}^{- 1}), \end{matrix}

(197)

and

\begin{matrix} v_{61} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(- \tanh_{A} (2 a ξ) \mp (i \sqrt{p q} {sech}_{A} (2 a ξ)))}^{- 1})^{2}, \end{matrix}

(198)

\begin{matrix} U_{62} (x, t) & = \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(- \coth_{A} (2 a ξ) \mp (\sqrt{p q} {csch}_{A} (2 a ξ)))}^{- 1} . \end{matrix}

(199)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{62} (x, t) & = e^{i θ} (\sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(- \coth_{A} (2 a ξ) \mp (\sqrt{p q} {csch}_{A} (2 a ξ)))}^{- 1}), \end{matrix}

(200)

and

\begin{matrix} v_{62} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(- \coth_{A} (2 a ξ) \mp (\sqrt{p q} {csch}_{A} (2 a ξ)))}^{- 1})^{2}, \end{matrix}

(201)

and

\begin{matrix} U_{63} (x, t) & = \sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(\frac{- 1}{2} (\tanh_{A} (\frac{a ξ}{2}) + \coth_{A} (\frac{a ξ}{2})))}^{- 1} . \end{matrix}

(202)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{63} (x, t) & = e^{i θ} (\sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(\frac{- 1}{2} (\tanh_{A} (\frac{a ξ}{2}) + \coth_{A} (\frac{a ξ}{2})))}^{- 1}), \end{matrix}

(203)

and

\begin{matrix} v_{63} (x, t) & = \frac{1}{1 + k^{2}} (\sqrt{- 4 {(\ln (A))}^{2} a^{2} - Ω} \times \\ {(\frac{- 1}{2} (\tanh_{A} (\frac{a ξ}{2}) + \coth_{A} (\frac{a ξ}{2})))}^{- 1})^{2} . \end{matrix}

(204)

Family 25.

For

b = λ

a = n λ (n \neq 0)

and

c = 0

, the following set of equations presents the corresponding family of solitary wave solutions:

\begin{matrix} U_{64} (x, t) & = 2 \sqrt{- b^{2} {(\ln (A))}^{2} - Ω} n {(A^{λ ξ} - n)}^{- 1} + \\ \sqrt{- b^{2} {(\ln (A))}^{2} - Ω} . \end{matrix}

(205)

As a result of utilizing (9) and (11), the subsequent solitary wave solutions are obtained for (1):

\begin{matrix} u_{64} (x, t) & = e^{i θ} (2 \sqrt{- b^{2} {(\ln (A))}^{2} - Ω} n {(A^{λ ξ} - n)}^{- 1} + \\ \sqrt{- b^{2} {(\ln (A))}^{2} - Ω}), \end{matrix}

(206)

and

\begin{matrix} v_{64} (x, t) & = \frac{1}{1 + k^{2}} (2 \sqrt{- b^{2} {(\ln (A))}^{2} - Ω} n {(A^{λ ξ} - n)}^{- 1} + \\ \sqrt{- b^{2} {(\ln (A))}^{2} - Ω})^{2}, \end{matrix}

(207)

where

ξ = \frac{x^{δ}}{Γ (1 + δ)} + \sqrt{\frac{b^{2} {(\ln (A))}^{2} + 2 {a_{1}}^{2} - 4 c {(\ln (A))}^{2} a - 2 {b_{1}}^{2}}{(b^{2} - 4 a c) {(\ln (A))}^{2}}} \frac{t^{γ}}{Γ (1 + γ)}

and

θ = \frac{a_{1} x^{δ}}{Γ (1 + δ)} + \frac{b_{1} t^{γ}}{Γ (1 + γ)} .

4. Discussion and Graphs

The findings of this work provide significant light on how the fractional coupled Higgs system behaves. By analyzing the graphs and traveling wave solutions, we were able to identify a number of crucial system characteristics. In the solutions, we discovered solitary waves, kink waves, rogue waves, periodic and hyperbolic solitary waves, all of which are crucial to the physical fields that the fractional coupled Higgs system represents.

Confined disturbances known as solitary waves maintain their shape and amplitude as they pass through a medium. These waves are typically seen in nonlinear systems, where the interaction between dispersive and nonlinear effects balances out to determine how they behave. Solitary waves suggest that the fractional coupled Higgs system may behave in a nonlinear and dispersive manner.

On the other hand, kink waves are disturbances that move over the boundary between two different media with different physical properties. These waves typically occur in systems with phase transitions and are characterized by their sudden change from one medium to another. Kink waves in the solutions demonstrate that the fractional Higgs system may undergo phase transitions under certain circumstances.

Strong and sporadic phenomena known as rogue waves appear unpredictably in otherwise well-behaved systems. These waves have the power to completely destroy physical systems, including fiber-optic communications and sea waves. The existence of rogue waves in the fractional Higgs system demonstrates that the system is vulnerable to catastrophic and unexpected phenomena, which might have important ramifications for comprehending the system’s behavior.

In the fractional coupled Higgs system, we mostly found periodic and hyperbolic solitary wave solutions. These responses are crucial because they provide a more thorough understanding of the system’s conduct in diverse circumstances.

This study demonstrates that mEDAM offers a reliable and flexible method for solving challenging mathematical models. The idea has been used to study a number of physical systems, and it seems to be a workable method for figuring out how nonlinear and dispersive systems behave. Investigating the solutions obtained with this method can provide us with information on the behavior of the fractional Higgs system and related phenomena, which may have significant ramifications for our knowledge of fundamental physics. Furthermore, our analytical methodology has the benefit of allowing us to derive solutions produced by other analytical methods such as the tan-method, (G

^{'}

/G)-expansion method, subequation method, and many others. The three families of solitary solutions produced by the (G

^{'}

/G)-expansion approach, for example, may be derived from our results. Similarly, by substituting generalized trigonometric and hyperbolic functions with ordinary trigonometric and hyperbolic functions, we may recover all families of solutions found using the tan-method. This adaptability in acquiring multiple answers broadens our understanding of the problem from numerous angles and broadens the usefulness of our discoveries in varied contexts.

Remark 1.

Overall, Figure 1 shows the solitary kink wave profile, which is a special and distinctive wave form distinguished by a steep bend or abrupt change in direction.

Remark 2.

Overall, Figure 2 shows a lump wave, which is a little bulge or disturbance in a waveform that exhibits a quick rise in amplitude followed by a steady fall to its initial level.

Remark 3.

Figure 3 shows the periodic kink wave profile, which is a special and distinctive wave form distinguished by a steep bend or abrupt change in direction.

Remark 4.

Overall, lone rogue wave, a rare and extreme oceanic occurrence in physics, is depicted in Figure 4. With a noticeably larger amplitude and steeper profile than the other waves, it stands out for its size and unpredictability.

Remark 5.

Overall, a singular kink wave is depicted in Figure 5. Singular kink waves are nonlinear wave system amplitude or phase variations that occur in a variety of scientific environments, including particle physics. They are especially important in theories including scalar fields, such as the Higgs field in the Standard Model, which influences particle mass. To solve unexplained events, extensions such as the “Fractional Coupled Higgs System” include more scalar fields. The mEDAM’s analysis of these waves provides insights into the interaction of new scalar fields with the Higgs field. This research extends particle physics beyond the Standard Model by improving scalar field comprehension in the Higgs region.

Remark 6.

Overall, a lump wave is depicted in Figure 6. A lump wave, also known as a compacton, is a specialized localized solution in nonlinear systems that, unlike other waves, retains a fixed spatial extent while propagating. Nonlinear and dispersive effects cause this behavior. The study of lump waves is relevant in the context of the “Fractional Coupled Higgs System”, an extension of the Standard Model of particle physics. The study of lump waves within this system, which includes various scalar fields in addition to the Higgs field, provides insights on scalar field behavior and interactions. Exploring lump waves contributes to a better understanding of scalar field dynamics and their influence on system behavior, which contributes to the larger objective of grasping extended particle physics models and hypothetical phenomena beyond the standard model.

5. Conclusions

We used mEDAM’s power to construct insightful solitary wave solutions for the fractional coupled Higgs system, adding critical generalized functions to describe its intricate physical processes. We investigated the dynamic behavior of adjustable free parameters and built practical linkages by incorporating them. This study expands our understanding of the system’s features and dynamic patterns, proving mEDAM’s resilience in decoding complicated mathematical models. The relevance of these results extends beyond engineering and scientific areas, promising useful paths for future research into the complicated dynamics of the fractional coupled Higgs system and its associated phenomena.

Author Contributions

The ideas of the present paper come from J.I. and M.B.; J.I., F.A.A. and R.A. wrote and completed the calculations; M.B., E.A.A.I. and R.A. checked all the results. All authors have read and agreed to the published version of the manuscript.

Funding

This project is funded by King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgments

Researchers supporting project number (RSPD2023R576), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

Abro, K.A.; Atangana, A. Dual fractional modeling of rate type fluid through non-local differentiation. Numer. Methods Partial Differ. Equ. 2020, 38, 390–405. [Google Scholar] [CrossRef]
Metzler, R.; Glöckle, W.G.; Nonnenmacher, T.F. Fractional model equation for anomalous diffusion. Phys. Stat. Mech. Its Appl. 1994, 211, 13–24. [Google Scholar] [CrossRef]
Tarasov, V.E.; Zaslavsky, G.M. Fractional dynamics of systems with long-range interaction. Commun. Nonlinear Sci. Numer. Simul. 2006, 11, 885–898. [Google Scholar] [CrossRef]
Su, N. Random fractional partial differential equations and solutions for water movement in soils: Theory and applications. Hydrol. Process. 2023, 37, e14844. [Google Scholar] [CrossRef]
Che, J.; Guan, Q.; Wang, X. Image denoising based on adaptive fractional partial differential equations. In Proceedings of the 2013 6th International Congress on Image and Signal Processing (CISP), Hangzhou, China, 16–18 December 2013; Volume 1, pp. 288–292. [Google Scholar]
Kachhia, K.B.; Prajapati, J.C. Solution of fractional partial differential equation aries in study of heat transfer through diathermanous materials. J. Interdiscip. Math. 2015, 18, 125–132. [Google Scholar] [CrossRef]
Zhang, Y. A finite difference method for fractional partial differential equation. Appl. Math. Comput. 2009, 215, 524–529. [Google Scholar] [CrossRef]
Ford, N.J.; Xiao, J.; Yan, Y. A finite element method for time fractional partial differential equations. Fract. Calc. Appl. Anal. 2011, 14, 454–474. [Google Scholar] [CrossRef]
Lyu, W.; Wang, Z.A. Global classical solutions for a class of reaction-diffusion system with density-suppressed motility. arXiv 2021, arXiv:2102.08042. [Google Scholar] [CrossRef]
Li, C.; Chen, A. Numerical methods for fractional partial differential equations. Int. J. Comput. Math. 2018, 95, 1048–1099. [Google Scholar] [CrossRef]
Yang, Q.; Liu, F.; Turner, I. Numerical methods for fractional partial differential equations with Riesz space fractional derivatives. Appl. Math. Model. 2010, 34, 200–218. [Google Scholar] [CrossRef]
Lin, X.; Liu, Y.; Yu, J.; Yu, R.; Zhang, J.; Wen, H. Stability analysis of three-phase grid-connected inverter under the weak grids with asymmetrical grid impedance by LTP theory in time domain. Int. J. Electr. Power Energy Syst. 2022, 142, 108244. [Google Scholar] [CrossRef]
Naeem, M.; Azhar, O.F.; Zidan, A.M.; Nonlaopon, K.; Shah, R. Numerical analysis of fractional-order parabolic equations via Elzaki transform. J. Funct. Spaces 2021, 2021, 3484482. [Google Scholar] [CrossRef]
Liu, X.; Shi, T.; Zhou, G.; Liu, M.; Yin, Z.; Yin, L.; Zheng, W. Emotion classification for short texts: An improved multi-label method. Humanit. Soc. Sci. Commun. 2023, 10, 306. [Google Scholar] [CrossRef]
Al-Sawalha, M.M.; Ababneh, O.Y.; Shah, R.; Nonlaopon, K. Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators. Aims Math. 2023, 8, 2308–2336. [Google Scholar] [CrossRef]
Liu, X.; Zhou, G.; Kong, M.; Yin, Z.; Li, X.; Yin, L.; Zheng, W. Developing Multi-Labelled Corpus of Twitter Short Texts: A Semi-Automatic Method. Systems 2023, 11, 390. [Google Scholar] [CrossRef]
Chen, D.; Wang, Q.; Li, Y.; Li, Y.; Zhou, H.; Fan, Y. A general linear free energy relationship for predicting partition coefficients of neutral organic compounds. Chemosphere 2020, 247, 125869. [Google Scholar] [CrossRef]
Elagan, S.K.; Sayed, M.; Higazy, M. An analytical study on fractional partial differential equations by Laplace transform operator method. Int. J. Appl. Eng. Res. 2018, 13, 545–549. [Google Scholar]
Mahor, T.C.; Mishra, R.; Jain, R. Analytical solutions of linear fractional partial differential equations using fractional Fourier transform. J. Comput. Appl. Math. 2021, 385, 113202. [Google Scholar] [CrossRef]
Thabet, H.; Kendre, S. New modification of Adomian decomposition method for solving a system of nonlinear fractional partial differential equations. Int. J. Adv. Appl. Math. Mech. 2019, 6, 1–13. [Google Scholar]
Duran, S. Exact solutions for time-fractional Ramani and Jimbo—Miwa equations by direct algebraic method. Adv. Sci. Eng. Med. 2020, 12, 982–988. [Google Scholar] [CrossRef]
Momani, S.; Odibat, Z. Homotopy perturbation method for nonlinear partial differential equations of fractional order. Phys. Lett. A 2007, 365, 345–350. [Google Scholar] [CrossRef]
Odibat, Z.; Momani, S. The variational iteration method: An efficient scheme for handling fractional partial differential equations in fluid mechanics. Comput. Math. Appl. 2009, 58, 2199–2208. [Google Scholar] [CrossRef]
Mirhosseini-Alizamini, S.M.; Rezazadeh, H.; Eslami, M.; Mirzazadeh, M.; Korkmaz, A. New extended direct algebraic method for the Tzitzica type evolution equations arising in nonlinear optics. Comput. Methods Differ. Equ. 2020, 8, 28–53. [Google Scholar]
Yasmin, H.; Aljahdaly, N.H.; Saeed, A.M.; Shah, R. Investigating Symmetric Soliton Solutions for the Fractional Coupled Konno–Onno System Using Improved Versions of a Novel Analytical Technique. Mathematics 2023, 11, 2686. [Google Scholar] [CrossRef]
Yasmin, H.; Aljahdaly, N.H.; Saeed, A.M.; Shah, R. Probing Families of Optical Soliton Solutions in Fractional Perturbed Radhakrishnan–Kundu–Lakshmanan Model with Improved Versions of Extended Direct Algebraic Method. Fractal Fract. 2023, 7, 512. [Google Scholar] [CrossRef]
Rezazadeh, H.; Mirhosseini-Alizamini, S.M.; Neirameh, A.; Souleymanou, A.; Korkmaz, A.; Bekir, A. Fractional Sine–Gordon equation approach to the coupled higgs system defined in time-fractional form. Iran. J. Sci. Technol. Trans. A Sci. 2019, 43, 2965–2973. [Google Scholar] [CrossRef]
Jabbari, A.; Kheiri, H.; Bekir, A. Exact solutions of the coupled Higgs equation and the Maccari system using He’s semi-inverse method and (G’/G)-expansion method. Comput. Math. Appl. 2011, 62, 2177–2186. [Google Scholar] [CrossRef]
Atas, S.S.; Ali, K.K.; Sulaiman, T.A.; Bulut, H. Invariant optical soliton solutions to the Coupled-Higgs equation. Opt. Quantum Electron. 2022, 54, 754. [Google Scholar] [CrossRef]
Seadawy, A.R.; Lu, D.; Khater, M.M. Bifurcations of traveling wave solutions for Dodd–Bullough–Mikhailov equation and coupled Higgs equation and their applications. Chin. J. Phys. 2017, 55, 1310–1318. [Google Scholar] [CrossRef]
Mu, G.; Qin, Z. Rogue waves for the coupled Schrödinger–Boussinesq equation and the coupled Higgs equation. J. Phys. Soc. Jpn. 2012, 81, 084001. [Google Scholar] [CrossRef]
Sikora, B. Remarks on the Caputo fractional derivative. MINUT 2023, 5, 76–84. [Google Scholar]
Khan, H.; Barak, S.; Kumam, P.; Arif, M. Analytical Solutions of Fractional Klein-Gordon and Gas Dynamics Equations, via the (G^′/G)-Expansion Method. Symmetry 2019, 11, 566. [Google Scholar] [CrossRef]

$Fractalfract 07 00653 g001$

Figure 1. Graph of (28) in which (a,b) present 3D and contour plots of real part respectively while (c,d) present 3D and contour plots of imaginary part respectively for

a = 3, b = 1, c = 5, A = e

a_{1} = 1000, b_{1} = 0, γ = δ = 1

Figure 1. Graph of (28) in which (a,b) present 3D and contour plots of real part respectively while (c,d) present 3D and contour plots of imaginary part respectively for

a = 3, b = 1, c = 5, A = e

a_{1} = 1000, b_{1} = 0, γ = δ = 1

$Fractalfract 07 00653 g001$

$Fractalfract 07 00653 g002$

Figure 2. Graph of (46) in which (a,b) present 3D and contour plots of real part respectively while (c,d) and present 3D and contour plots of imaginary part respectively for

a = 3, b = 0, c = 5, A = e

a_{1} = 1000, b_{1} = 0, γ = δ = 1

Figure 2. Graph of (46) in which (a,b) present 3D and contour plots of real part respectively while (c,d) and present 3D and contour plots of imaginary part respectively for

a = 3, b = 0, c = 5, A = e

a_{1} = 1000, b_{1} = 0, γ = δ = 1

$Fractalfract 07 00653 g002$

$Fractalfract 07 00653 g003$

Figure 3. Graph of the real part of (136) in which (a,b) present 3D and contour plots of real part respectively for

a = 3, b = 10, c = 5, A = 2, a_{1} = 1000, b_{1} = 0, γ = δ = 1

Figure 3. Graph of the real part of (136) in which (a,b) present 3D and contour plots of real part respectively for

a = 3, b = 10, c = 5, A = 2, a_{1} = 1000, b_{1} = 0, γ = δ = 1

$Fractalfract 07 00653 g003$

$Fractalfract 07 00653 g004$

Figure 4. Graph of (204) in which (a,b) present 3D and contour plots of real part respectively while (c,d) present 3D and contour plots of imaginary part respectively for

a = 3, b = 0, c = - 3, A = e

a_{1} = 1000, b_{1} = 10, 000, γ = δ = 1

Figure 4. Graph of (204) in which (a,b) present 3D and contour plots of real part respectively while (c,d) present 3D and contour plots of imaginary part respectively for

a = 3, b = 0, c = - 3, A = e

a_{1} = 1000, b_{1} = 10, 000, γ = δ = 1

$Fractalfract 07 00653 g004$

$Fractalfract 07 00653 g005$

Figure 5. Depicts graph of the squared norm of (44). The 3D graph in (a) is plotted for

a = 3

b = 10, c = 1, A = e, a_{1} = b_{1} = 25, γ = δ = 1

while the 2D graph in (b) is drawn for

t = 0

and the similar values of involved parameters.

Figure 5. Depicts graph of the squared norm of (44). The 3D graph in (a) is plotted for

a = 3

b = 10, c = 1, A = e, a_{1} = b_{1} = 25, γ = δ = 1

while the 2D graph in (b) is drawn for

t = 0

and the similar values of involved parameters.

$Fractalfract 07 00653 g005$

$Fractalfract 07 00653 g006$

Figure 6. Depicts graph of the squared norm of (62). The 3D graph in (a) is plotted for

a = 2

b = 0, c = - 1, A = e, a_{1} = 0, b_{1} = 0, γ = δ = 1

while the 2D graph in (b) is drawn for

t = 0

and the similar values of involved parameters.

Figure 6. Depicts graph of the squared norm of (62). The 3D graph in (a) is plotted for

a = 2

b = 0, c = - 1, A = e, a_{1} = 0, b_{1} = 0, γ = δ = 1

while the 2D graph in (b) is drawn for

t = 0

and the similar values of involved parameters.

$Fractalfract 07 00653 g006$

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Bilal, M.; Iqbal, J.; Ali, R.; Awwad, F.A.; Ismail, E.A.A. Exploring Families of Solitary Wave Solutions for the Fractional Coupled Higgs System Using Modified Extended Direct Algebraic Method. Fractal Fract. 2023, 7, 653. https://doi.org/10.3390/fractalfract7090653

AMA Style

Bilal M, Iqbal J, Ali R, Awwad FA, Ismail EAA. Exploring Families of Solitary Wave Solutions for the Fractional Coupled Higgs System Using Modified Extended Direct Algebraic Method. Fractal and Fractional. 2023; 7(9):653. https://doi.org/10.3390/fractalfract7090653

Chicago/Turabian Style

Bilal, Muhammad, Javed Iqbal, Rashid Ali, Fuad A. Awwad, and Emad A. A. Ismail. 2023. "Exploring Families of Solitary Wave Solutions for the Fractional Coupled Higgs System Using Modified Extended Direct Algebraic Method" Fractal and Fractional 7, no. 9: 653. https://doi.org/10.3390/fractalfract7090653

APA Style

Bilal, M., Iqbal, J., Ali, R., Awwad, F. A., & Ismail, E. A. A. (2023). Exploring Families of Solitary Wave Solutions for the Fractional Coupled Higgs System Using Modified Extended Direct Algebraic Method. Fractal and Fractional, 7(9), 653. https://doi.org/10.3390/fractalfract7090653

Article Menu

Exploring Families of Solitary Wave Solutions for the Fractional Coupled Higgs System Using Modified Extended Direct Algebraic Method

Abstract

1. Introduction

2. Methodology of mEDAM

3. Execution of the Approach

Implementation of mEDAM to the Problem

4. Discussion and Graphs

5. Conclusions

Author Contributions

Funding

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI