Exact Soliton Solutions to the Cubic-Quartic Non-linear Schrödinger Equation With Conformable Derivative

The research paper aims to investigate the space-time fractional cubic-quartic non-linear Schrödinger equation in the appearance of the third, and fourth-order dispersion impacts without both group velocity dispersion, and disturbance with parabolic law media by utilizing the extended sinh-Gordon expansion method. This method is one of the strongest methods to find the exact solutions to the non-linear partial differential equations. In order to confirm the existing solutions, the constraint conditions are used. We successfully construct various exact solitary wave solutions to the governing equation, for example, singular, and dark-bright solutions. Moreover, the 2D, 3D, and contour surfaces of all obtained solutions are also plotted. The finding solutions have justified the efficiency of the proposed method.


INTRODUCTION AND MOTIVATION
Non-linear partial differential equations have different types of equations, one of them is the non-linear Schrödinger equation (NLSE) that relevant to the classical and quantum mechanics. The non-linear Schrödinger equation is a generalized (1+1)-dimensional version of the Ginzburg-Landau equation presented in 1950 in their study on supraconductivity and has been specifically reported by Chiao et al. [1] in their research of optical beams. In the past several years, various methods have been proposed to obtain the exact optical soliton solutions of the non-linear Schrödinger equation [2][3][4][5][6][7][8][9][10][11][12]. Dispersion and non-linearity are two of the essential components for the distribution of solitons across inter-continental regions. Usually, group velocity dispersion (GVD) level with self-phase modulation in a sensible manner allows these solitons to sustain tall range travel. In fact it might happen that the GVD is tiny and thus totally ignored, in this case the dispersion effect is determined by third and fourth order dispersion effects. Subsequently, this equation has been studied in a variety of ways, such as the Lie symmetry [13], both the m + G ′ G -improved expansion, and the exp (−ϕ (ξ )) −expansion methods [14], and the semiinverse variation principle method [4]. In this study, the extended sinh-Gordon expansion method (ShGEM) is applied to the non-linear cubic-quartic Schrödinger equations with the Parabolic law of fractional order, which is given by where u (x, t) is the complex valued wave function. The operator D α of order α, where α ∈ (0, 1] is the fractional derivative, the parameters γ and β are real constants, a real-valued algebraic function F |u| 2 is p-times continuously differentiable, then By using the relation of on Equation (1), we obtain the fractional non-linear Schrödinger equations with Parabolic law as follows: The extended sinh-Gordon expansion method is intended to a generalization of the sine-Gordon expansion equation because it is based on an auxiliary equation namely the sine-Gordon equation (see previous studies [15,16] for details). Moreover, different computational and numerical methods have been utilized to constructed new solutions to the non-linear partial differential equations, such as the variable separated method [17], the auxiliary parameters and residual power series method [18], the Bernoulli sub-equation method [19,20], the modified auxiliary expansion method [21], the homotopy analysis transform method [22][23][24][25][26], the homotopy perturbation sumudu transform method [27], the shooting method with the explicit Runge-Kutta scheme [28,29], and the Adomian decomposition method [30]. Recently, several fractional operators have been applied to the mathematical models in order to seek their exact solutions, such as the Laplace transform [31,32], the Nabla operator [33][34][35].
The outline of paper are organize the paper as follows: A short review of the conformable derivative is presented in section 2. Section 3 deals with the analysis of the ShGEM. In section 4, the method is applied to solve the non-linear Schrödinger equation involving the fractional derivatives with the parabolic law. Eventually, in section 5, we presented our conclusion of this paper.
Definition 2.2. Assume that c ≥ 0 and t ≥ c, let h be a function defined on (c, t] as well as α ∈ R. Then, the α-fractional integral of h is given by  (26) where if the Riemann improper integral exists.
here "prime" is the classical derivatives with respect to t.

THE EXTENDED ShGEM
In the current section, we presented the main steps of the e ShGEM (see previous study [42,43]). Consider the following fractional non-linear PDE:  where p = p (x, t).
Consider the wave transformation by substitute relation Equation (8) into Equation (7), we obtain the following non-linear ODE: Consider the trial solution of Equation (9) of the form The parameters A j , B j , for (j = 1, 2, . . . , k) and A 0 are real constants, and θ is a function of η that hold the following ODE: The homogeneous balance principle is applied on Equation (9) to find the value of k. From the space-time fractional the sinh-Gordon equation, we have (see previous study [15,16]).

IMPLEMENT OF THE EXTENDED ShGEM
The implementation of the extended ShGEM to the cubic-quartic non-linear Schrödinger equation with conformable derivative is provided in this section. Consider the wave transformation In Equation (18), θ (x, t) represents the phase component of the soliton. The ω, κ, ν are the wave number, the soliton frequency, and the soliton velocity, respectively. Substituting wave transformation into Equation (2) and splitting the outcomes equation into real and imaginary parts, we gain Multiply both sides of Equation (19) by U ′ and integrate it, we obtain From Equation (20), we get constraint conditions ν = 4γ κ 3 − 3βκ 2 and β = 4γ κ. Balancing the terms U ′′′ U ′ and U 6 yields κ = 1. With κ = 1, Equations (10), (16), and (17) change to and ψ (ζ ) = ±m 1 csch (ζ ) ± n 1 coth (ζ ) + n 0 , respectively. Inserting Equation (22) along with Equation (12) into Equation (21), and using constraint conditions provides a non-linear algebraic system. Equaling each coefficient of sinh i (θ ) cosh j ( θ ) with the same power to zero, and finding the obtained system of algebraic equations, we gain the values of the parameters. Putting the obtained values of the parameters into Equations (23) and (24), give the solutions of Equation (3).

Set 1
we get Set 2 we get Set 3 we get Set 4 we get we obtain Set 6 we get Set 7 we get

CONCLUSION
In this article, we have successfully used the extended sinh-Gordon expansion method to solve the problem for the non-linear cubic-quartic Schrödinger equations involving fractional derivatives with the Parabolic law. A traveling wave transforms in the sense of the comfortable derivative has been used to convert the governing equation into a NODE. The various optical solutions of the studied model have been constructed, for example, the singular soliton solutions as shown in Figures 1-6, and the dark-bright soliton solution as seen in Figure 7. Comparing our solutions to the results obtained in references [16][17][18], our findings solutions are new and different.
To better analyze the dynamic attitude, and the characteristics of these solutions, the 2D, 3D and counter-surface of all obtained solutions are plotted. The study shows that this method is the effective and appropriate technique for finding the exact solution of the model considered in the paper.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/supplementary material.

AUTHOR CONTRIBUTIONS
HD contributed in developing the proofs and edited the article for possible improvement. HG, KA, and RY contributed in developing the main results and proofs. All authors read the final version and approved it.