New Optical Solutions of the Fractional Gerdjikov-Ivanov Equation With Conformable Derivative

Finding exact analytic solutions to the partial equations is one of the most challenging problems in mathematical physics. Generally speaking, the exact solution to many categories of such equations can not be found. In these cases, the use of numerical and approximate methods is inevitable. Nevertheless, the exact PDE solver methods are always preferred because they present the solution directly without any restrictions to use. This article aims to examine the perturbed Gerdjikov-Ivanov equation in an exact approach point of view. This equation plays a significant role in non-linear fiber optics. It also has many important applications in photonic crystal fibers. To this end, firstly, we obtain some novel optical solutions of the equation via a newly proposed analytical method called generalized exponential rational function method. In order to understand the dynamic behavior of these solutions, several graphs are plotted. To the best of our knowledge, these two techniques have never been tested for the equation in the literature. The findings of this article may have a high significance application while handling the other non-linear PDEs.


INTRODUCTION
Non-linear Schrödinger equations (NLSE) are often studied from different points of view. In recent years a great variety of analytical and numerical methods have been proposed for solving these equations [1][2][3][4]. The most studied NLSE equation is that which has a cubic non-linearity. In the present paper, we will explore an NLSE that has a quintic non-linearity, namely the perturbed Gerdjikov-Ivanov (pGI) equation.
The main achievement of this research is to utilize a new method to derive some novel solutions to a variant form of NLSE. In particular, we consider the pGI equation is given by [5][6][7][8][9][10][11] i ∂q ∂t + a ∂ 2 q ∂x 2 + b|q| 4 q = i cq 2 ∂q ∂x provided that q(x, t) indicates the macroscopic complex-valued wave profile of temporal and spatial independent variables of t and x, respectively. In this equation, ∂q/∂t is linear temporal evolution, ∂ 2 q/∂x 2 stands for the group velocity dispersion (GVD), and |q| 4 q is the present quintic non-linearity of the model. The parameters a, b are the coefficients of these quantities, respectively. Moreover c is the non-linear dispersion coefficient. Finally, the constants λ 1 , λ 2 , and θ are known parameters related to perturbative effects. More information on this model can be found in the references as mentioned earlier.
In recent years, because of its high importance, the model has attracted the attention of many researchers. For instance, Biswas and Alqahtani [6] have presented two varieties of bright soliton solutions by the use of the semi-inverse variational principle. The sine-Gordon equation approach has been used to extract the dark, bright, dark-bright, singular, and combined singular optical solitons of the equation in Yaşar et al. [7]. Biswas et al. [8] have retrieved some bright and singular optical soliton solutions to the pGI equation by the implementation of the extended trial equation method. The exp(φ(ξ ))-Expansion and the Kudryashov methods are two reliable techniques that have been used in Arshed [9] to investigate some solitary wave solutions of the equation. In Kaur and Wazwaz [10], several hyperbolic, trigonometric or rational function solutions have been proposed using two efficient techniques, namely exp(φ(ξ ))-Expansion and G ′ G 2 -expansion methods. Very recently, Hosseini et al. [11] have listed several Kink, bright, and dark optical solitons of the model by the aid of the exp a -function method and a new version of the Kudryashov method.
In light of previous work, we will apply the generalized exponential rational function method (GERFM) to retrieve some new analytical optical solutions of the fractional pGI equation with the conformable derivative [12]. This new definition of derivative is based on the basic limit definition of the derivative that has been successfully tackled in solving many different problems [13][14][15][16][17][18][19][20][21][22]. The main structure of the present article is as outlined. In the second section of this paper, some mathematical preliminaries have been reviewed. This section includes the necessary steps of applying GERFM, and the definition and basic properties of the conformable derivative will be presented. The main results of this article are achieved by following these steps in section 3 of this contribution. In section 4, we have performed some numerical simulations of the obtained results. These graphs can help us in better understanding of their dynamic properties. Finally, the article concludes with some conclusions.

MATHEMATICAL PRELIMINARIES AND BACKGROUNDS
This section first deals with the structure of the GERFM. Then in the next subsection, the basic concepts of the conformable derivative are expressed.

Analysis of GERFM
GERFM is a newly developed method introduced by Ghanbari and Inc [23] to solve the resonance non-linear Schrödinger equation [23]. Other successful applications of the technique in solving different types of PDEs have also been reported in references [24][25][26][27].
We will review how to use the method below.
1. Let us consider a typical non-linear PDE for q = q(x, t), giving by Under the wave transformations of q(x, t) = Q(ξ ) and ξ = σ x−lt, Equation (2) becomes an ordinary differential equation given by: 2. Now, we assume that Equation (3) admits the exact solution giving by where (ξ ) = m 1 e n 1 ξ + m 2 e n 2 ξ m 3 e n 3 ξ + m 4 e n 4 ξ .
and m i , n i 's and A 0 , A k , and B k 's are disposal parameters. Finally, N is a constant, which is evaluated by applying the homogeneous balance to Equation (3).

Inserting Equation (4) into (3) with Equation (5), and then gathering all possible powers of
Equating coefficients of P to zero, one derives a simultaneous system of equations regarding m i , . 4. Finally, solving the non-linear system and substituting the obtained solutions in Equations (4) and (5), the explicit form of the solutions of (2) will be extracted.

The Conformable Derivative
Definition: Let q : R + → R, then the conformable derivative of q of order α, is giving by [12] D α Theorem: For any α ∈ (0, 1], and two α-differentiable functions p, q, the following propositions hold dt holds. Theorem [14]: Let p :(0, 1] → R be a function such that p is classical, and α-conformable differentiable. Moreover, consider q as a differentiable function defined in the range of p. Thus, we have where prime stands for standard derivatives respect to t. Some of the benefits of the conformable derivative compared to other new definitions for the derivative are as follows: • According to this definition of the operator, the derivative of a constant function is zero. This feature is not available in many other definitions. • Unlike many existing definitions, this definition satisfies the known formula of the derivative of the product of two functions. • The conformable derivative does satisfy the known formula of the derivative of the quotient of two functions. • The conformable derivative does satisfy the well-known chain rule. • The conformable derivative satisfies the well-known semi group property.
The mentioned properties are very important and valuable features for any derivative definition that the conformable derivative has all of them.

MATHEMATICAL ANALYSIS
The main contribution in this paper is to consider the derivatives in the Equation (1) with the conformable derivative defined by (6), as follows The main assumption is to taking the stationary soliton solution form of where ν, k, and ω are the phase component, the frequency of solitons, and the wavenumber, respectively. Substituting the stationary soliton solution form (8) into Equation (7), we arrive at a complex equation whose real part is as follows So, we will have From the real part, the following formula is also extracted Thus, in the following, we focus our attention on deriving solutions of Equation (11). Now balancing between two terms of Q 5 and Q ′′ in Equation (11) suggests N = 1 2 . If we want to get a closed-form solution, we need to define a new variable of Q(ξ ) = R 2 (ξ ). This substitution leads us to Now, the homogeneous balance in Equation (12) suggests N = 1. Setting N = 1 along with Equation (4), one gets Inserting (13) into (12) and pursuing the steps outlined for the method, the analytical solutions for the Equation (7) will be determined consequently.
Case 1: Inserting these values in Equation (13), yields Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and Case 2: Inserting these values in Equation (13), yields Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and

Category 2:
It is attained m = [2, 0, 1, −1] and n = [1, 0, 1, −1], which offers Case 1: Inserting these values in Equation (13), yields Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and Case 2: Inserting these values in Equation (13), yields Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and Case 1: Accordingly, we derive a soliton solution of given PDE in (7) as Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and Case 2: Accordingly, we derive a soliton solution of given PDE in (7) as e ξ 6e 2ξ + 10e ξ + 4 . Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and Category 4: It is attained m = [−1, 0, 1, 0] and n = [0, 1, 0, 1], which offers Case 1: Inserting these values in Equation (13), yields Accordingly, we derive a soliton solution of given PDE in (7) as provided that ab < 0, and Comparing our acquired solutions with other existing reported in the literature shows that ours are different and new. All the acquired solutions are new and have not been reported in the previous papers. Particularly, since the form of the equation and derivative considered in this article are the same as those given in reference [7], we can check that the results and the requirements for their existence in the two papers are quite different. Furthermore, we have checked the correctness of all obtained solutions, and found they satisfy the original equation.

NUMERICAL SIMULATIONS
In this section, we have presented several numerical simulations using the algorithm proposed in subsection 2.1. To illustrate the dynamic behaviors of the analytical results obtained in section 3, Figures 1-7 have been depicted. Figure 1 shows the dynamic behavior of the solution q 1 (x, t) defined in (15) (18), for the given values a = 1, b = −3, c = 0.5, λ 1 = 1, λ 2 = 4, α = 0.9 is plotted in Figure 3. Moreover, the diagram of q 4 (x, t) is displayed in Figure 4 corresponding to the choices of a = 1.5, b = −2.5, c = 1, λ 1 = 7, λ 2 = 2, α = 0.9. Taking the parameters a = 0.2, b = −2, c = 0.5, λ 1 = 2, λ 2 = 2, α = 0.99 into consideration, the graph of the solution q 5 (x, t) presented in the Equation (21) is plotted in Figure 5. Moreover, Figure 6 illustrate the dynamic behaviors of the analytical solution q 6 (x, t) obtained in Equation (22) by taking a = 1.5, b = −5, c = 65, λ 1 = 0.3, λ 2 = 0.7, α = 0.9. And finally, the profiles of the exact solution q 7 (x, t) presented in Equation (24) is displayed in Figure 7, when a = −2, b = 3, c = 3.5, λ 1 = 2, λ 2 = 0.5, α = 0.95 are chosen as parameters in the main PDE of (7). The performed numerical simulations admit that the solutions are of kinky and anti-kinky, and the trigonometric classifications. Also, by carefully looking at the structure of the obtained solutions, it can be seen that the corresponding conformable derivative parameter of α appears in the formula of all the solutions.

CONCLUSIONS
Partial differential equation is a powerful and effective tool for modeling non-linear systems. Finding the exact solution to such equations is one of the most challenging problems in mathematics. There is also no specific way of solving many of these equations. In these cases, we must resort to the approximate analytical methods due to the limitations of exact solver methods. According to what stated above, new approaches to solving PDE equations are of great importance and application. The main objective of this paper is to employ a well-known technique called GERFM to solve the perturbed Gerdjikov-Ivanov equation with the comfortable derivative. One of the outstanding features of the model considered in this article is the use of the definition of the comfortable derivative in the structure of the model. This definition is one of the most interesting definitions for a derivative that has many ideal features for a derivative. Applying this definition to the model will provide us with many advantages compared to the standard derivative. One of the advantages of the method used in this article is the determination of various categories of solutions during the method. Several numerical simulations are presented to gain a better understanding of the properties of the acquired solutions. By comparing the obtained results with the results of the present papers, it can be seen that the obtained results are not reported in any of the previous literature. It is worth mentioning that GERFM is capable of reducing the volume of needed computational compared to some other analytical method. The straightforward application is another advantage of the technique compared to other known techniques. This method can also be utilized to solve many other similar problems.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.