# Analytical Solutions of the Fractional Complex Ginzburg-Landau Model Using Generalized Exponential Rational Function Method with Two Different Nonlinearities

Authors: Ghazala Akram; Maasoomah Sadaf; M. Atta Ullah Khan; Hasan Hosseinzadeh
Journal: Advances in Mathematical Physics
 (2023)

---

## Abstract

The complex Ginzburg-Landau model appears in the mathematical description of wave propagation in nonlinear optics. In this paper, the fractional complex Ginzburg-Landau model is investigated using the generalized exponential rational function method. The Kerr law and parabolic law are considered to discuss the nonlinearity of the proposed model. The fractional effects are also included using a novel local fractional derivative of orderα. Many novel solutions containing trigonometric functions, hyperbolic functions, and exponential functions are acquired using the generalized exponential rational function method. The 3D-surface graphs, 2D-contour graphs, density graphs, and 2D-line graphs of some retrieved solutions are plotted using Maple software. A variety of exact traveling wave solutions are reported including dark, bright, and kink soliton solutions. The nature of the optical solitons is demonstrated through the graphical representations of the acquired solutions for variation in the fractional order of derivative. It is hoped that the acquired solutions will aid in understanding the dynamics of the various physical phenomena and dynamical processes governed by the considered model.


---

## Body

## 1. Introduction

Nonlinear partial differential equations (NPDEs) govern the majority of physical phenomena and dynamical processes. The study of nonlinear wave propagation problems provides motivations for developing NPDEs. The most significant task in nonlinear science is to obtain solutions for NPDEs, particularly solitary and soliton wave solutions. Nonlinear natural and physical phenomena are better understood with the aid of solutions. In chemistry, mathematical physics, mathematical biology, and many other basic sciences, NPDEs are frequently used to model different phenomena and dynamics. The NPDEs such as the Lakshmanan-Porsezian-Daniel model [1, 2], Navier-Stokes equations [3], shallow water-like equation [4], resonant nonlinear Schrödinger equation [5], and Korteweg-de Vries equation [6] are used to transform numerous natural processes in mathematical form. Numerical simulations, mathematical expressions, and formulations are increasingly used in various problems of science in recent years [7–10].In this paper, the complex Ginzburg-Landau (CGL) model can be expressed as(1)iUt+δ1Uxx+δ2FU2U−δ3U2U∗2U2U2xx−U2x2−δ4U=0,where U represents the profile of wave propagation and it depends upon the space variable x and time variable t. Moreover, δ1 and δ2 represent the group velocity dispersion and nonlinearity coefficient, respectively. The real parameters δ3 and δ4 are included to represent perturbation and detuning effects, respectively. Functional F is essential for complex function FU2U:C⟶C in two-dimensional linear space C. The k-times continuously differentiable function FU2U is of the form
(2)FU2U∈⋃m,n=1∞Ck−m,m×−n,n,R2.In scientific studies, fractional order derivatives [11–13] are very important and play an essential role in the investigation of NPDEs. Recently, the fractional CGL model is investigated by many researchers [14–16]. In this paper, the CGL model is investigated using the newly proposed definition of new local fractional derivative (NLFD) Nhyp,tα of order α defined in [17] on the basis of the study of local fractional derivatives introduced in [18].The purpose of this paper is to examine the fractional CGL model with newly proposed NLFDNhyp,tα of order α. A generalized exponential rational function (GERF) method is utilized to acquire the analytical solitary wave solutions of the fractional CGL equation. The proposed work involving the GERF method along with the local fractional derivative of order α is novel which has not been done previously in literature.The CGL model is examined using GERF along with the newly proposed local fractional differential operator for the first time in this work. A variety of novel interesting graphical observations are presented to exhibit the dynamical behavior of the CGL model. Bright, dark, and kink optical solitons are obtained which are graphically examined for the effect of variation in the value of the fractional parameter on the traveling wave solutions.In the following subsection, a description of the NLFD is illustrated.


### 1.1. The Description of NLFD Orderα

In this subsection, the basic properties of NLFDNhyp,tα of order α as defined in [18] are presented as follows.Definition 1.
Consider0<α<1, 0<t and define
(3)Nhyp,tαht=limε⟶0ht+εt1−α/2sech1−αt1+α/2−htε.
Forh:a,b⟶ℝ with a<t, if there exist limε⟶a+Nhypα, then
(4)Nhypαha=limε⟶a+Nhypαht.
In this case,h is known as α-differentiable at a with respect to Nhypαht. One significant outcome of the new fractional local derivative is of the form
(5)Nhypαht=t1−α/2sech1−αt1+α/2dftdt,where h:a,b⟶ℝ is differentiable with a<t. The following property which has been proved in [17] is given as
(6)Nhypα21−α2sinh1−αt1+α/2=1.
The chain rule is defined as(7)Nhypαg∘hc=g′hcNhypαhc.
Using Equations (6) and (7), it is obtained that
(8)NhypαHξ=c1H′ξ,ξ=21−α2sinh1−αt1+α/2,where c1 is constant.The fractional derivative defined in Equation (3) is a novel generalization of the classical integer order derivative which exhibits interesting mathematical properties. The application of this definition is easier due to the chain rule. Moreover, the properties of NLFD make it suitable to evaluate higher order derivatives. These results are essential to extract new solitary wave solutions and soliton solutions of nonlinear evolution equations.The rest of the research is structured as follows: the GERF method is explained in Section2. The implementation of the proposed method is given in Section 3. 3D-surface plots, 2D-contour plots, density plots, and comparison of 2D-line plots are given in Section 4. Section 5 includes the conclusion of the research paper.


## 1.1. The Description of NLFD Orderα

In this subsection, the basic properties of NLFDNhyp,tα of order α as defined in [18] are presented as follows.Definition 1.
Consider0<α<1, 0<t and define
(3)Nhyp,tαht=limε⟶0ht+εt1−α/2sech1−αt1+α/2−htε.
Forh:a,b⟶ℝ with a<t, if there exist limε⟶a+Nhypα, then
(4)Nhypαha=limε⟶a+Nhypαht.
In this case,h is known as α-differentiable at a with respect to Nhypαht. One significant outcome of the new fractional local derivative is of the form
(5)Nhypαht=t1−α/2sech1−αt1+α/2dftdt,where h:a,b⟶ℝ is differentiable with a<t. The following property which has been proved in [17] is given as
(6)Nhypα21−α2sinh1−αt1+α/2=1.
The chain rule is defined as(7)Nhypαg∘hc=g′hcNhypαhc.
Using Equations (6) and (7), it is obtained that
(8)NhypαHξ=c1H′ξ,ξ=21−α2sinh1−αt1+α/2,where c1 is constant.The fractional derivative defined in Equation (3) is a novel generalization of the classical integer order derivative which exhibits interesting mathematical properties. The application of this definition is easier due to the chain rule. Moreover, the properties of NLFD make it suitable to evaluate higher order derivatives. These results are essential to extract new solitary wave solutions and soliton solutions of nonlinear evolution equations.The rest of the research is structured as follows: the GERF method is explained in Section2. The implementation of the proposed method is given in Section 3. 3D-surface plots, 2D-contour plots, density plots, and comparison of 2D-line plots are given in Section 4. Section 5 includes the conclusion of the research paper.


## 2. The Description of GERF Method

The GERF method is a relatively new approach to NPDEs. Many researchers utilized this method to extract the solutions of NPDEs [19–23]. The NPDE is considered as
(9)HU,Ux,Ut,Uxx,Utt,⋯=0,where U represents the profile of wave propagation and it is dependent on space variable x and temporal variable t. Moreover, H represents a polynomial in U and partial derivatives of U. Using fractional derivatives as illustrated in Subsection 1.1, Equation (9) is transformed into fractional form as
(10)HNhyp,xαU,Nhyp,tαU,Nhyp,x2αU,Nhyp,t2αU,⋯=0,where 0<α<1. The fractional transformations are considered as
(11)Ux,t=Aξeiη,ξ=21−α2sinh1−αw1x1+α/2−k1t1+α/2,η=21−α2sinh1−αw2x1+α/2−k2t1+α/2.Using fractional transformations from Equation (8), Equation (10) is reduced to the following ODE:
(12)QA,A′,A″,⋯=0,where A′=dA/dξ, Q is the polynomial of A and its ordinary derivatives.The general solution for Equation (12) is considered as
(13)Aξ=a0+∑i=1Naiϕiξ+bjϕ−iξ,where a0, ai, and bi are unknown parameters to be evaluated. The function ϕξ is defined as
(14)ϕξ=τ1eκ1ξ+τ2eκ2ξτ3eκ3ξ+τ4eκ4ξ,where τj and κj are real or complex numbers j=1,2,3,4. For different values of τj and κj, Equation (13) together with Equation (14) provides the different exact soliton solutions of Equation (12). Among these solutions, the nontrivial solutions are of significance for the determination of the nontrivial solutions of Equation (10) by backward substitution.In formal solution Equation (13), N can be determined by the homogeneous balance principle (HBP). There are some instances in which N retains the negative or fractional values. The following substitutions are used in such cases:
(15)A=Vp/q,p,q=1,q≠0,A=V−N.By substituting the general solution Equation (13) together with Equation (14) into Equation (10) and coefficients of ejκξj=1,2,3,⋯, taken to be zero, in result, the system of linear equations is acquired. The unknown parameters can be determined by resolving the retrieved system of linear equations simultaneously with the aid of any Maple software. By substituting nontrivial values of solutions to Equation (13), the solution of Equation (10) can be obtained.


## 3. Mathematical Analysis

Using the NLFDNhyp,tα of order α, NPDE (1) is transformed to the following form:
(16)iNhyp,tαU+δ1Nhyp,x2αU+δ2FU2U−δ3U2U∗2U2Nhyp,x2αU2−Nhyp,xαU22−δ4U=0,where the time and space fractional derivatives are defined according to the definition of the novel local fractional derivative given by Equation (3) as
(17)Nhyp,tαUx,t=limε⟶0Ux,t+εt1−α/2sech1−αt1+α/2−Ux,tε,0<α<1,0<t,Nhyp,xαUx,t=limε⟶0Ux+εx1−α/2sech1−αx1+α/2,t−Ux,tε,0<α<1,0<x.Equation (16) is the fractional form of CGL Equation (1) with the newly proposed definition of local fractional derivatives. This is obtained by replacing integer order derivatives with fractional order derivatives. The fractional order derivatives are more generalized than the classical integer order derivatives. This definition of fractional derivative is proposed and utilized to report several interesting and useful results in recent research works [17, 24].The fractional transformations from Equation (8) are used to transform Equation (16) into ODE having the following real part:
(18)k2−δ1w22−δ4A+δ2FA2A+δ1w12−4δ3w12A″=0.The imaginary part is given as(19)−k1+2δ1w1w2A′=0.The following substitution is considered:(20)B=k2−δ1w22−δ4,C=δ1w12−4δ3w12.Hence, Equation (18) is simplified to the following form:
(21)BA+δ2FA2A+CA″=0.Equation (19) yields
(22)k1=−2δ1w1w2.The solution of Equation (21) is obtained subject to condition (22) by considering the Kerr law nonlinearity in the following subsection.


### 3.1. Kerr Law Nonlinearity

Due to the propagation of light waves, optical fibers exhibit many definite forms of nonlinearities. One of these nonlinearities is Kerr law nonlinearity which is mathematically described as(23)Fv=v.The Kerr law of nonlinearity arises due to the nonlinear responses from the nonharmonic motion of electrons bound in molecules faced by a light wave in an optical fiber, under the influence of an external electric field. The Kerr nonlinearity is observed in most optical fibers. The nonlinearity of Equation (21) with the Kerr law is considered as
(24)FAξ2=Aξ2.Using Equation (24), Equation (21) is transformed to the following form:
(25)BA+δ2A3+CA″=0.In Equation (13), formal solution N=1 is evaluated by HBL. Then the generalized solution of Equation (25) is transformed as
(26)Aξ=a0+a1ϕξ+b1ϕ−1ξ.The formal solution of Equation (1) is written as
(27)Ux,t=AξeiηorUx,t=a0+a1ϕξ+b1ϕ−1ξeiη,where a0, a1, and b1 are unknown constants to be determined. The solution of Equation (16) is investigated under the following cases.Case 1.
Ifτ1,τ2,τ3,τ4=1,−1,1,1, κ1,κ2,κ3,κ4=1,−1,1,−1, Equation (14) is simplified to the following form:
(28)ϕξ=tanhξ.
Equation (26) with Equation (28) is substituted in Equation (25). By taking the coefficients of ϕ’s zero, algebraic linear equations are acquired. By resolving these algebraic equations simultaneously, the following sets of solitary wave solutions are acquired as follows.
Set 1.(29)a0=0,a1=0,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=2δ1w12+δ1w22−8δ3w12+δ4.
Set 2.(30)a0=0,a1=±−2δ1−8δ3δ2w1,b1=0,w1=w1,w2=w2,k2=2δ1w12+δ1w22−8δ3w12+δ4.
Set 3.(31)a0=0,a1=±−2δ1−8δ3δ2w1,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=±32δ1−8δ3w12+2δ1w12+δ1w22−8δ3w12+δ4.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the solution of Equation (16) is determined as follows:
(32)U11x,t=±−2δ1−8δ3δ2w1cothξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the solution of Equation (16) is determined as follows:
(33)U21x,t=±−2δ1−8δ3δ2w1tanhξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 3, the solution of Equation (16) is determined as follows:
(34)U31x,t=±−2δ1−8δ3δ2w1tanhξ±−2δ1−8δ3δ2w1cothξeiη.Case 2.
Ifτ1,τ2,τ3,τ4=i,−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) is reduced to the following form:
(35)ϕξ=−tanξ.
Equation (26) with Equation (35) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(36)a0=0,a1=0,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
Set 2.(37)a0=0,a1=±−2δ1−8δ3δ2w1,b1=0,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
Set 3.(38)a0=0,a1=±−2δ1−8δ3δ2w1,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=±32δ1−8δ3w12−2δ1w12+δ1w22+8δ3w12+δ4.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the solution of Equation (16) is determined as
(39)U12x,t=±−2δ1−8δ3δ2w1cotξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the solution of Equation (16) is determined as
(40)U22x,t=±−2δ1−8δ3δ2w1tanξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 3, the following result for the solution of Equation (16) is obtained:
(41)U32x,t=±−2δ1−8δ3δ2w1tanξ±−2δ1−8δ3δ2w1cotξeiη.Case 3.
Ifτ1,τ2,τ3,τ4=1+i,1−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) is reduced to the form
(42)ϕξ=1−tanξ.
Equation (26) with Equation (42) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(43)a0=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,k2=−2δ1w12+δ1w22+8δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following solution of Equation (16) is obtained:
(44)U13x,t=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w11−tanξeiη.Case 4.
Ifτ1,τ2,τ3,τ4=1+i,1−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) is reduced to the form
(45)ϕξ=2+tanξ.
Equation (26) with Equation (45) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(46)a0=±20w1δ1−4δ3δ2−50δ1−200δ3/δ2,a1=0,b1=±−50δ1−200δ3δ2w1,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
Set 2.(47)a0=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the result for the solution of Equation (16) is obtained as
(48)U24x,t=±20w1δ1−4δ3δ2−50δ1−200δ3/δ2±−50δ1−200δ3δ2w112+tanξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the following result is determined for the solution of Equation (16):
(49)U14x,t=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w12+tanξeiη.Case 5.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of the considered equation. Therefore, with the values τ1,τ2,τ3,τ4=2,1,1,1 and κ1,κ2,κ3,κ4=1,0,1,0, Equation (14) is taken which reduces Equation (14) to the following form:
(50)ϕξ=2eξ+1eξ+1.
Equation (26) with Equation (50) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of eξ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(51)a0=±3δ1w22+δ4−k2δ2,a1=±−−4δ1w22−4δ4+4k2δ2,b0=0,w1=±−2δ1w22+2δ4−2k2δ1−4δ3.
Set 2.(52)a0=±3δ1w22+δ4−k2δ2,a1=0,b1=±−−16δ1w22−16δ4+16k2δ2,w1=±−2δ1w22+2δ4−2k2δ1−4δ3.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined.
(53)U15x,t=±3δ1w22+δ4−k2δ2±−−4δ1w22−4δ4+4k2δ22eξ+1eξ+1eiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(54)U25x,t=±3δ1w22+δ4−k2δ2±−−16δ1w22−16δ4+16k2δ2eξ+12eξ+1eiη.Case 6.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of the considered equation. Therefore, with the values τ1,τ2,τ3,τ4=2,0,1,1 and κ1,κ2,κ3,κ4=−1,0,1,−1, Equation (14) is taken which reduces Equation (14) to the following form:
(55)ϕξ=1−tanhξ.
Equation (26) with Equation (55) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(56)a0=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,k2=2δ1w12+δ1w22−8δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(57)U16x,t=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w11−tanhξeiη.Case 7.
Ifτ1,τ2,τ3,τ4=−3,−1,−1,1, κ1,κ2,κ3,κ4=−1,1,−1,1, Equation (14) is reduced to the following form:
(58)ϕξ=tanhξ−2.
Equation (26) with Equation (58) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(59)a0=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=−2δ1−8δ3δ2w1,b1=0,k2=2δ1w12+δ1w22−8δ3w12+δ4,w1=w1,w2=w2.
Set 2.(60)a0=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=0,b1=±3−2δ1−8δ3δ2w1,k2=2δ1w12+δ1w22−8δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(61)U17x,t=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w1tanhξ−2eiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(62)U27x,t=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2±3−2δ1−8δ3δ2w11tanhξ−2eiη.Case 8.
Ifτ1,τ2,τ3,τ4=1,0,1,1, κ1,κ2,κ3,κ4=0,0,1,0, Equation (14) is reduced to the following form:
(63)ϕξ=11+eξ.
Equation (26) with Equation (63) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(64)a0=±w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,k2=12δ1w12+δ1w22−2δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(65)U18x,t=±w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w111+eξeiη.


### 3.2. Parabolic Law Nonlinearity

The parabolic law nonlinearity is considered as(66)FAξ2=d1Aξ2+d2Aξ4.Using Equation (66), Equation (21) is reduced to the following form:
(67)BA+δ2d1A3+δ2d2A5+CA″=0.The value ofN=1/2 is extracted using the homogeneous balance principle. Since the value of N is not a positive integer, the following substitution is used:
(68)A=V1/2.Hence, Equation (67) is simplified to the following form:
(69)4BV2+4δ2d1V3+4δ2d2V4−CV′2+2CVV″=0.The positive integer value ofN=1 is obtained using the homogeneous balance principle. The generalized solution Equation (69) is written as
(70)Vξ=a0+a1ϕξ+b1ϕ−1ξ.The formal solution of Equation (16) is written as
(71)Ux,t=a0+a1ϕξ+b1ϕ−1ξ1/2eiη.In Equations (70) and (71), a0, a1, and b1 are unknown constants to be determined. The solution of Equation (13) is investigated in the following cases.Case 1.
Ifτ1,τ2,τ3,τ4=1,−1,1,1, κ1,κ2,κ3,κ4=1,−1,1,−1, Equation (14) is simplified to the following form:
(72)ϕξ=tanhξ.
Equation (70) with Equation (72) is substituted in Equation (69). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(73)a0=−38d1d2,a1=0,b1=±38d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(74)a0=−38d1d2,a1=±38d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 3.(75)a0=−38d1d2,a1=±316d1d2,b1=±316d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ264d2δ1−256d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(76)U19x,t=−38d1d2±38d1d2cothξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(77)U29x,t=−38d1d2±38d1d2tanhξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 3, the following result for the solution of Equation (16) is determined:
(78)U39x,t=−38d1d2±316d1d2tanhξ±316d1d2cothξ1/2eiη.Case 2.
Ifτ1,τ2,τ3,τ4=i,−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) reduces to the form
(79)ϕξ=−tanξ.
Equation (70) with Equation (79) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The system of algebraic linear equations is acquired by taking the coefficients of zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(80)a0=−38d1d2,a1=±316d1d2,b1=±316d1d2,k2=64δ1w22d2+15d12δ2+64δ4d264d2,w1=±−3δ264d2δ1−256d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(81)U110x,t=−38d1d2±316d1d2−tanξ±316d1d2−cotξ1/2eiη.Case 3.
Ifτ1,τ2,τ3,τ4=1+i,1−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) reduces to the following form:
(82)ϕξ=1−tanξ.
Equation (70) with Equation (82) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved system of equations simultaneously with the aid of Maple software.
Set 1.(83)a0=−34d1d2,a1=38d1d2,b1=34d1d2,k2=16δ1w22d2−3d12δ2+16δ4d2d2,w1=−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(84)a0=±3d144+25228d29+42,a1=±3+62d156d2,b1=3+62d128d2,k2=1784d2±401+2982233632δ1w22d22+33564d12δ22±314384δ1w22d2+233632δ42d22±2,43473d12±δ2,314384,w1=±±122δ2+27δ2−784d2δ1+3136d2δ3d1.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(85)U111x,t=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w11−tanξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(86)U211x,t=±3d144+25228d29+42±3+62d156d21−tanξ+3+62d128d21−tanξ1/2eiη.Case 4.
Ifτ1,τ2,τ3,τ4=2,1,1,1, κ1,κ2,κ3,κ4=1,0,1,0, Equation (14) is reduced to the following form:
(87)ϕξ=2eξ+1eξ+1.
Equation (70) with Equation (87) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(88)a0=−12d1d2,a1=14d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=δ220d2δ1−80d2δ3d1.
Set 2.(89)a0=−54d1d2,a1=54d1d2,b1=0,k2=−−16d2δ1w22+45d12δ2−16d2δ416d2,w1=5δ24d2δ1−16d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(90)U112x,t=−12d1d2+14d1d22eξ+1eξ+11/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(91)U212x,t=−54d1d2+54d1d22eξ+1eξ+11/2eiη.Case 5.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of considered equations. Therefore, τ1,τ2,τ3,τ4=2,0,1,1 and κ1,κ2,κ3,κ4=−1,0,1,−1 are taken which reduce Equation (14) to the following form:
(92)ϕξ=1−tanhξ.
Equation (70) with Equation (92) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(93)a0=0,a1=−38d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(94)a0=−34d1d2,a1=38d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=−3δ216d2δ1−64d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(95)U113x,t=−38d1d21−tanhξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(96)U213x,t=−34d1d2+38d1d21−tanhξ1/2eiη.Case 6.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of considered equations. Therefore, τ1,τ2,τ3,τ4=−3,−1,−1,1 and κ1,κ2,κ3,κ4=−1,1,−1,1 are taken which reduce Equation (14) to the following form:
(97)ϕξ=tanhξ−2.
Equation (70) with Equation (97) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(98)a0=38d1d2,a1=0,b1=98d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(99)a0=−98d1d2,a1=0,b1=−98d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 3.(100)a0=38d1d2,a1=38d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 4.(101)a0=−98d1d2,a1=−38d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 5.(102)a0=−34d1d2,a1=−316d1d2,b1=−916d1d2,k2=64δ1w22d2+3d12δ2+64δ4d264d2,w1=±−3δ264d2δ1−256d2δ3d1,w2=w2.
Set 6.(103)a0=3d1±9−5316d22+3,a1=−31±3d132d2,w1=−±3δ23+6δ2128d2δ1−512d2δ3d1,k2=1408δ1w22+1408δ4d2−219d12δ23±2432δ1w22−2432δ4d2+381d12δ214083d2±2432d2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(104)U114x,t=38d1d2+98d1d21tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(105)U214x,t=−98d1d2−98d1d21tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 3, the following result for the solution of Equation (16) is determined:
(106)U314x,t=38d1d2+38d1d2tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 4, the following result for the solution of Equation (16) is determined:
(107)U414x,t=−98d1d2−38d1d2tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 5, the following result for the solution of Equation (16) is determined:
(108)U514x,t=−34d1d2−316d1d2tanhξ−2−916d1d21tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 6, the following result for the solution of Equation (16) is determined:
(109)U614x,t=3d1±9−5316d22+3−31±3d132d2tanhξ−21/2eiη.Case 7.
Ifτ1,τ2,τ3,τ4=1,0,1,1, κ1,κ2,κ3,κ4=0,0,1,0, Equation (14) is reduced to the following form:
(110)ϕξ=11+eξ.
Equation (70) with Equation (110) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved system of equations simultaneously with the aid of Maple software.
Set 1.(111)a0=0,a1=−34d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=±−3δ24d2δ1−16d2δ3d1,w2=w2,
Set 2.(112)a0=−34d1d2,a1=34d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=−3δ24d2δ1−16d2δ3,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(113)U115x,t=−34d1d211+eξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(114)U215x,t=−34d1d2+34d1d211+eξ1/2eiη.


## 3.1. Kerr Law Nonlinearity

Due to the propagation of light waves, optical fibers exhibit many definite forms of nonlinearities. One of these nonlinearities is Kerr law nonlinearity which is mathematically described as(23)Fv=v.The Kerr law of nonlinearity arises due to the nonlinear responses from the nonharmonic motion of electrons bound in molecules faced by a light wave in an optical fiber, under the influence of an external electric field. The Kerr nonlinearity is observed in most optical fibers. The nonlinearity of Equation (21) with the Kerr law is considered as
(24)FAξ2=Aξ2.Using Equation (24), Equation (21) is transformed to the following form:
(25)BA+δ2A3+CA″=0.In Equation (13), formal solution N=1 is evaluated by HBL. Then the generalized solution of Equation (25) is transformed as
(26)Aξ=a0+a1ϕξ+b1ϕ−1ξ.The formal solution of Equation (1) is written as
(27)Ux,t=AξeiηorUx,t=a0+a1ϕξ+b1ϕ−1ξeiη,where a0, a1, and b1 are unknown constants to be determined. The solution of Equation (16) is investigated under the following cases.Case 1.
Ifτ1,τ2,τ3,τ4=1,−1,1,1, κ1,κ2,κ3,κ4=1,−1,1,−1, Equation (14) is simplified to the following form:
(28)ϕξ=tanhξ.
Equation (26) with Equation (28) is substituted in Equation (25). By taking the coefficients of ϕ’s zero, algebraic linear equations are acquired. By resolving these algebraic equations simultaneously, the following sets of solitary wave solutions are acquired as follows.
Set 1.(29)a0=0,a1=0,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=2δ1w12+δ1w22−8δ3w12+δ4.
Set 2.(30)a0=0,a1=±−2δ1−8δ3δ2w1,b1=0,w1=w1,w2=w2,k2=2δ1w12+δ1w22−8δ3w12+δ4.
Set 3.(31)a0=0,a1=±−2δ1−8δ3δ2w1,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=±32δ1−8δ3w12+2δ1w12+δ1w22−8δ3w12+δ4.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the solution of Equation (16) is determined as follows:
(32)U11x,t=±−2δ1−8δ3δ2w1cothξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the solution of Equation (16) is determined as follows:
(33)U21x,t=±−2δ1−8δ3δ2w1tanhξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 3, the solution of Equation (16) is determined as follows:
(34)U31x,t=±−2δ1−8δ3δ2w1tanhξ±−2δ1−8δ3δ2w1cothξeiη.Case 2.
Ifτ1,τ2,τ3,τ4=i,−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) is reduced to the following form:
(35)ϕξ=−tanξ.
Equation (26) with Equation (35) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(36)a0=0,a1=0,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
Set 2.(37)a0=0,a1=±−2δ1−8δ3δ2w1,b1=0,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
Set 3.(38)a0=0,a1=±−2δ1−8δ3δ2w1,b1=±−2δ1−8δ3δ2w1,w1=w1,w2=w2,k2=±32δ1−8δ3w12−2δ1w12+δ1w22+8δ3w12+δ4.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the solution of Equation (16) is determined as
(39)U12x,t=±−2δ1−8δ3δ2w1cotξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the solution of Equation (16) is determined as
(40)U22x,t=±−2δ1−8δ3δ2w1tanξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 3, the following result for the solution of Equation (16) is obtained:
(41)U32x,t=±−2δ1−8δ3δ2w1tanξ±−2δ1−8δ3δ2w1cotξeiη.Case 3.
Ifτ1,τ2,τ3,τ4=1+i,1−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) is reduced to the form
(42)ϕξ=1−tanξ.
Equation (26) with Equation (42) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(43)a0=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,k2=−2δ1w12+δ1w22+8δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following solution of Equation (16) is obtained:
(44)U13x,t=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w11−tanξeiη.Case 4.
Ifτ1,τ2,τ3,τ4=1+i,1−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) is reduced to the form
(45)ϕξ=2+tanξ.
Equation (26) with Equation (45) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(46)a0=±20w1δ1−4δ3δ2−50δ1−200δ3/δ2,a1=0,b1=±−50δ1−200δ3δ2w1,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
Set 2.(47)a0=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,w1=w1,w2=w2,k2=−2δ1w12+δ1w22+8δ3w12+δ4.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the result for the solution of Equation (16) is obtained as
(48)U24x,t=±20w1δ1−4δ3δ2−50δ1−200δ3/δ2±−50δ1−200δ3δ2w112+tanξeiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the following result is determined for the solution of Equation (16):
(49)U14x,t=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w12+tanξeiη.Case 5.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of the considered equation. Therefore, with the values τ1,τ2,τ3,τ4=2,1,1,1 and κ1,κ2,κ3,κ4=1,0,1,0, Equation (14) is taken which reduces Equation (14) to the following form:
(50)ϕξ=2eξ+1eξ+1.
Equation (26) with Equation (50) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of eξ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(51)a0=±3δ1w22+δ4−k2δ2,a1=±−−4δ1w22−4δ4+4k2δ2,b0=0,w1=±−2δ1w22+2δ4−2k2δ1−4δ3.
Set 2.(52)a0=±3δ1w22+δ4−k2δ2,a1=0,b1=±−−16δ1w22−16δ4+16k2δ2,w1=±−2δ1w22+2δ4−2k2δ1−4δ3.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined.
(53)U15x,t=±3δ1w22+δ4−k2δ2±−−4δ1w22−4δ4+4k2δ22eξ+1eξ+1eiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(54)U25x,t=±3δ1w22+δ4−k2δ2±−−16δ1w22−16δ4+16k2δ2eξ+12eξ+1eiη.Case 6.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of the considered equation. Therefore, with the values τ1,τ2,τ3,τ4=2,0,1,1 and κ1,κ2,κ3,κ4=−1,0,1,−1, Equation (14) is taken which reduces Equation (14) to the following form:
(55)ϕξ=1−tanhξ.
Equation (26) with Equation (55) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(56)a0=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,k2=2δ1w12+δ1w22−8δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(57)U16x,t=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w11−tanhξeiη.Case 7.
Ifτ1,τ2,τ3,τ4=−3,−1,−1,1, κ1,κ2,κ3,κ4=−1,1,−1,1, Equation (14) is reduced to the following form:
(58)ϕξ=tanhξ−2.
Equation (26) with Equation (58) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(59)a0=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=−2δ1−8δ3δ2w1,b1=0,k2=2δ1w12+δ1w22−8δ3w12+δ4,w1=w1,w2=w2.
Set 2.(60)a0=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=0,b1=±3−2δ1−8δ3δ2w1,k2=2δ1w12+δ1w22−8δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(61)U17x,t=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w1tanhξ−2eiη.
By putting the unknown constant values into Equation (27) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(62)U27x,t=±4w1δ1−4δ3δ2−2δ1−8δ3/δ2±3−2δ1−8δ3δ2w11tanhξ−2eiη.Case 8.
Ifτ1,τ2,τ3,τ4=1,0,1,1, κ1,κ2,κ3,κ4=0,0,1,0, Equation (14) is reduced to the following form:
(63)ϕξ=11+eξ.
Equation (26) with Equation (63) is substituted in Equation (25). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(64)a0=±w1δ1−4δ3δ2−2δ1−8δ3/δ2,a1=±−2δ1−8δ3δ2w1,b1=0,k2=12δ1w12+δ1w22−2δ3w12+δ4,w1=w1,w2=w2.
By putting the unknown constant values into Equation (27) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(65)U18x,t=±w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w111+eξeiη.


## 3.2. Parabolic Law Nonlinearity

The parabolic law nonlinearity is considered as(66)FAξ2=d1Aξ2+d2Aξ4.Using Equation (66), Equation (21) is reduced to the following form:
(67)BA+δ2d1A3+δ2d2A5+CA″=0.The value ofN=1/2 is extracted using the homogeneous balance principle. Since the value of N is not a positive integer, the following substitution is used:
(68)A=V1/2.Hence, Equation (67) is simplified to the following form:
(69)4BV2+4δ2d1V3+4δ2d2V4−CV′2+2CVV″=0.The positive integer value ofN=1 is obtained using the homogeneous balance principle. The generalized solution Equation (69) is written as
(70)Vξ=a0+a1ϕξ+b1ϕ−1ξ.The formal solution of Equation (16) is written as
(71)Ux,t=a0+a1ϕξ+b1ϕ−1ξ1/2eiη.In Equations (70) and (71), a0, a1, and b1 are unknown constants to be determined. The solution of Equation (13) is investigated in the following cases.Case 1.
Ifτ1,τ2,τ3,τ4=1,−1,1,1, κ1,κ2,κ3,κ4=1,−1,1,−1, Equation (14) is simplified to the following form:
(72)ϕξ=tanhξ.
Equation (70) with Equation (72) is substituted in Equation (69). The system of algebraic linear equations is acquired by taking the coefficients of ϕ’s zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(73)a0=−38d1d2,a1=0,b1=±38d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(74)a0=−38d1d2,a1=±38d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 3.(75)a0=−38d1d2,a1=±316d1d2,b1=±316d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ264d2δ1−256d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(76)U19x,t=−38d1d2±38d1d2cothξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(77)U29x,t=−38d1d2±38d1d2tanhξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 3, the following result for the solution of Equation (16) is determined:
(78)U39x,t=−38d1d2±316d1d2tanhξ±316d1d2cothξ1/2eiη.Case 2.
Ifτ1,τ2,τ3,τ4=i,−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) reduces to the form
(79)ϕξ=−tanξ.
Equation (70) with Equation (79) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The system of algebraic linear equations is acquired by taking the coefficients of zero. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(80)a0=−38d1d2,a1=±316d1d2,b1=±316d1d2,k2=64δ1w22d2+15d12δ2+64δ4d264d2,w1=±−3δ264d2δ1−256d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(81)U110x,t=−38d1d2±316d1d2−tanξ±316d1d2−cotξ1/2eiη.Case 3.
Ifτ1,τ2,τ3,τ4=1+i,1−i,1,1, κ1,κ2,κ3,κ4=i,−i,i,−i, Equation (14) reduces to the following form:
(82)ϕξ=1−tanξ.
Equation (70) with Equation (82) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved system of equations simultaneously with the aid of Maple software.
Set 1.(83)a0=−34d1d2,a1=38d1d2,b1=34d1d2,k2=16δ1w22d2−3d12δ2+16δ4d2d2,w1=−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(84)a0=±3d144+25228d29+42,a1=±3+62d156d2,b1=3+62d128d2,k2=1784d2±401+2982233632δ1w22d22+33564d12δ22±314384δ1w22d2+233632δ42d22±2,43473d12±δ2,314384,w1=±±122δ2+27δ2−784d2δ1+3136d2δ3d1.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(85)U111x,t=±2w1δ1−4δ3δ2−2δ1−8δ3/δ2±−2δ1−8δ3δ2w11−tanξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(86)U211x,t=±3d144+25228d29+42±3+62d156d21−tanξ+3+62d128d21−tanξ1/2eiη.Case 4.
Ifτ1,τ2,τ3,τ4=2,1,1,1, κ1,κ2,κ3,κ4=1,0,1,0, Equation (14) is reduced to the following form:
(87)ϕξ=2eξ+1eξ+1.
Equation (70) with Equation (87) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(88)a0=−12d1d2,a1=14d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=δ220d2δ1−80d2δ3d1.
Set 2.(89)a0=−54d1d2,a1=54d1d2,b1=0,k2=−−16d2δ1w22+45d12δ2−16d2δ416d2,w1=5δ24d2δ1−16d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(90)U112x,t=−12d1d2+14d1d22eξ+1eξ+11/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(91)U212x,t=−54d1d2+54d1d22eξ+1eξ+11/2eiη.Case 5.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of considered equations. Therefore, τ1,τ2,τ3,τ4=2,0,1,1 and κ1,κ2,κ3,κ4=−1,0,1,−1 are taken which reduce Equation (14) to the following form:
(92)ϕξ=1−tanhξ.
Equation (70) with Equation (92) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(93)a0=0,a1=−38d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(94)a0=−34d1d2,a1=38d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=−3δ216d2δ1−64d2δ3d1,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(95)U113x,t=−38d1d21−tanhξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(96)U213x,t=−34d1d2+38d1d21−tanhξ1/2eiη.Case 6.
According to this method, Equation (14) allows to select any complex or real values of τi,κi,i=1,2,3,4, which can provide the nontrivial wave solutions of considered equations. Therefore, τ1,τ2,τ3,τ4=−3,−1,−1,1 and κ1,κ2,κ3,κ4=−1,1,−1,1 are taken which reduce Equation (14) to the following form:
(97)ϕξ=tanhξ−2.
Equation (70) with Equation (97) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved homogeneous equations at once with the aid of Maple software.
Set 1.(98)a0=38d1d2,a1=0,b1=98d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 2.(99)a0=−98d1d2,a1=0,b1=−98d1d2,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 3.(100)a0=38d1d2,a1=38d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 4.(101)a0=−98d1d2,a1=−38d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=±−3δ216d2δ1−64d2δ3d1,w2=w2.
Set 5.(102)a0=−34d1d2,a1=−316d1d2,b1=−916d1d2,k2=64δ1w22d2+3d12δ2+64δ4d264d2,w1=±−3δ264d2δ1−256d2δ3d1,w2=w2.
Set 6.(103)a0=3d1±9−5316d22+3,a1=−31±3d132d2,w1=−±3δ23+6δ2128d2δ1−512d2δ3d1,k2=1408δ1w22+1408δ4d2−219d12δ23±2432δ1w22−2432δ4d2+381d12δ214083d2±2432d2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(104)U114x,t=38d1d2+98d1d21tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(105)U214x,t=−98d1d2−98d1d21tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 3, the following result for the solution of Equation (16) is determined:
(106)U314x,t=38d1d2+38d1d2tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 4, the following result for the solution of Equation (16) is determined:
(107)U414x,t=−98d1d2−38d1d2tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 5, the following result for the solution of Equation (16) is determined:
(108)U514x,t=−34d1d2−316d1d2tanhξ−2−916d1d21tanhξ−21/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 6, the following result for the solution of Equation (16) is determined:
(109)U614x,t=3d1±9−5316d22+3−31±3d132d2tanhξ−21/2eiη.Case 7.
Ifτ1,τ2,τ3,τ4=1,0,1,1, κ1,κ2,κ3,κ4=0,0,1,0, Equation (14) is reduced to the following form:
(110)ϕξ=11+eξ.
Equation (70) with Equation (110) is substituted in Equation (69). By putting the accumulated terms of different powers of ϕ’s to zero, homogeneous linear equations are acquired. The following set of solutions are obtained by resolving the retrieved system of equations simultaneously with the aid of Maple software.
Set 1.(111)a0=0,a1=−34d1d2,b1=0,k2=16d2δ1w22+3d12δ2+16d2δ416d2,w1=±−3δ24d2δ1−16d2δ3d1,w2=w2,
Set 2.(112)a0=−34d1d2,a1=34d1d2,b1=0,k2=16δ1w22d2+3d12δ2+16δ4d216d2,w1=−3δ24d2δ1−16d2δ3,w2=w2.
By putting the unknown constant values into Equation (70) from corresponding Set 1, the following result for the solution of Equation (16) is determined:
(113)U115x,t=−34d1d211+eξ1/2eiη.
By putting the unknown constant values into Equation (70) from corresponding Set 2, the following result for the solution of Equation (16) is determined:
(114)U215x,t=−34d1d2+34d1d211+eξ1/2eiη.


## 4. Graphical Explanation

This section presents graphical representations of some obtained solutions. The 3D-surface graphs, 2D-contour graphs, density plots, and comparison of 2D-line graphs of retrieved solutions are plotted using Maple software. In plotted graphs, (a), (b), (c), and (d) represent the3D-surface graphs, 2D-contour graphs, density plots, and comparison of 2D-line graphs, respectively. The graphs are plotted for absolute values of obtained solutions. Using Maple software, both ± values cannot be plotted simultaneously; therefore, plotting graphs only requires + or − instead of ±. The suitable numeric values are assigned to undetermined constants δii=1,2,3,4,d1,d2,w2, and fractional derivative order α, to generate the well-shaped graphs of obtained solutions. From the corresponding sets of obtained solutions, the values of a0, a1, b1, w1, k1, and κ can be taken for plotting the graphs of acquired solutions.The acquired solutions are graphically presented in Figures1–12. The graphical representations include dark, bright, kink, and solitary wave solutions. The graphs exhibit the effect of variation in the fractional parameter α on the obtained solutions. The evolution of a dark soliton is shown corresponding to the solution U21x,t through Figures 1–3 for increasing values of fractional parameter α. Figures 1–3 show the graphical simulations for α=0.5,α=0.75, and α=0.95, respectively. Graphical simulations for the solution U13x,t are presented in Figures 4–6. The evolution of a kink soliton is illustrated from Figures 7–9, which has been obtained corresponding to the solution U29x,t. It is evident that the amplitude of the wave varies from one asymptotic state to another in this case. The bright soliton solution is shown in Figures 9–12 for α=0.5,α=0.75,α=0.9. The bright soliton is characterized by a local amplitude increase in the wave profile.Figure 1
The graphical representations ofU21x,t for NLFD of order α and for parameters δ1=6,δ2=δ4=0.6,δ3=0.1,w1=1,w2=0.01.
(a)(b)(c)(d)Figure 2
The graphical representations ofU21x,t for NLFD of order α and for parameters δ1=6,δ2=δ4=0.6,δ3=0.1,w1=1,w2=0.01.
(a)(b)(c)(d)Figure 3
The graphical representations ofU21x,t for NLFD of order α and for parameters δ1=6,δ2=δ4=0.6,δ3=0.1,w1=1,w2=0.01.
(a)(b)(c)(d)Figure 4
The graphical simulation representation ofU13x,t for NLFD of order α and for parameters δ1=0.63,δ2=δ3=w1=w2=0.1,δ4=0.6.
(a)(b)(c)(d)Figure 5
The graphical simulation representation ofU13x,t for NLFD of order α and for parameters δ1=0.63,δ2=δ3=w1=w2=0.1,δ4=0.6.
(a)(b)(c)(d)Figure 6
The graphical simulation representation ofU13x,t for NLFD of order α and for parameters δ1=0.63,δ2=δ3=w1=w2=0.1,δ4=0.6.
(a)(b)(c)(d)Figure 7
The graphical simulation representation ofU29x,t for NLFD of order α and for parameters δ1=0.3,δ2=δ3=d1=d2=0.1,δ4=0.6,w2=0.01.
(a)(b)(c)(d)Figure 8
The graphical simulation representation ofU29x,t for NLFD of order α and for parameters δ1=0.3,δ2=δ3=d1=d2=0.1,δ4=0.6,w2=0.01.
(a)(b)(c)(d)Figure 9
The graphical simulation representation ofU29x,t for NLFD of order α and for parameters δ1=0.3,δ2=δ3=d1=d2=0.1,δ4=0.6,w2=0.01.
(a)(b)(c)(d)Figure 10
The graphical simulation representation ofU111x,t for NLFD of order α and for parameters δ1=0.3,δ2=δ3=d1=d2=0.1,δ4=0.6,w2=0.01.
(a)(b)(c)(d)Figure 11
The graphical simulation representation ofU111x,t for NLFD of order α and for parameters δ1=0.3,δ2=δ3=d1=d2=0.1,δ4=0.6,w2=0.01.
(a)(b)(c)(d)Figure 12
The graphical simulation representation ofU111x,t for NLFD of order α and for parameters δ1=0.3,δ2=δ3=d1=d2=0.1,δ4=0.6,w2=0.01.
(a)(b)(c)(d)


## 5. Conclusion

In this paper, the complex Ginzburg-Landau model is investigated using novel local fractional derivatives of order0<α<1. CGL model solutions are acquired using the GERF method. Two laws of nonlinearities, Kerr law and parabolic law, are used to examine the nonlinearity of the proposed model. Many solitary wave solutions are obtained using the proposed method for the CGL model. The trigonometric functions tan,cot, hyperbolic functions tanh,coth, and exponential functions are included in obtained soliton solutions of the CGL model. Many soliton solutions, such as dark, bright, kink, and many other solitary wave solutions, are extracted using the proposed method. The 3D-surface graphs, 2D-contour graphs, density plots, and comparison 2D-line graphs are included to represent the dynamical behavior of retrieved solutions. The results presented in this paper are novel as compared to the already reported results for the CGL model. The authors of [25] presented the solutions of a CGL model using an iterative technique which was successfully utilized to retrieve the approximate solutions corresponding to specific initial conditions. The dark-singular soliton solutions of the CGL model were constructed in [26], but no other traveling wave solutions were investigated. The CGL model with M-truncated, beta, and conformable derivatives was theoretically examined using different mathematical techniques in previous works [11, 14, 27]. However, this manuscript utilizes a newly proposed local fractional derivative along with the GERF method to report a variety of explicit exact closed form traveling wave solutions of the CGL model in terms of free parameters. Hence, the comparison of the obtained results for the already reported studies of the CGL model confirm the novelty and usefulness of the results presented in this paper [11, 14, 25–27]. The reported results are useful to understand the possible dynamical behaviors exhibited by the CGL model. The soliton solutions will be useful additions to the literature for understanding related physical systems. In the future, the fractional CGL model can also be investigated by other nonlinearity laws and exact methods.



---

*Source: 9720612-2023-03-13.xml*