The Extended Direct Algebraic Method for Extracting Analytical Solitons Solutions to the Cubic Nonlinear Schrödinger Equation Involving Beta Derivatives in Space and Time

: In the ﬁeld of nonlinear optics, quantum mechanics, condensed matter physics, and wave propagation in rigid and other nonlinear instability phenomena, the nonlinear Schrödinger equation has signiﬁcant applications. In this study, the soliton solutions of the space-time fractional cubic nonlinear Schrödinger equation with Kerr law nonlinearity are investigated using an extended direct algebraic method. The solutions are found in the form of hyperbolic, trigonometric, and rational functions. Among the established solutions, some exhibit wide spectral and typical characteristics, while others are standard. Various types of well-known solitons, including kink-shape, periodic, V-shape, and singular kink-shape solitons, have been extracted here. To gain insight into the internal formation of these phenomena, the obtained solutions have been depicted in two-and three-dimensional graphs with different parameter values. The obtained solitons can be employed to explain many complicated phenomena associated with this model.


Introduction
Diverse real-world phenomena have been explained using nonlinear models, leading to the revelation of important information.Fractional nonlinear evolution equations represent an advanced class of differential equations that yield improved results.These equations help to illustrate intricate physical phenomena, attracting many researchers to work in this field due to their significant applications.Within the realm of fractional nonlinear evolution equations, the nonlinear Schrödinger equation plays a crucial role and finds applications in various areas such as quantum mechanics, optical fiber, plasma physics, fluid mechanics, biology, the dispersion of chemically reactive materials, electricity, shallow water wave phenomena, heat flow, finance, and fractal dynamics.
The relationship between the nonlinearity and dispersion components of medium solitons is uncovered, and as they travel through the medium, their undulation structure remains unaltered.The soliton solutions derived from FNLEE have practical and commercial applications in various fields such as optical fiber technology, telecommunications, signal processing, image processing, system identification, water purification, plasma physics, medical device sterilization, chemistry, and other related domains [1,2].Various dynamic approaches have been introduced and implemented in the literature to solve nonlinear fractional differential equations (NFDES) and obtain analytical traveling wave solutions, for example, the exp-function method [3], the Modified Exp-function method [4], the inverse scattering transformation method [5,6], the Bäcklund transformation method [7], the homogenous balance method [8,9], the Jacobi elliptic function method [10], the unified algebraic method [11], the sine-cosine method [12,13], the tanh-coth method [14,15], improved modified extended tanh-function method [16,17], the Lie symmetry analysis method [18], the extended generalized (G /G)-expansion method [19], the modified sim- ple equation method [20], the generalized Kudryashov method [21,22], the sine-Gordon expansion method [23], the Riccati-Bernoulli equation method [24,25], the new extended direct algebraic method [26,27], and the new auxiliary equation method [28].
Fractional derivatives have been widely applied in diverse scientific and engineering fields, including physics, mechanics, signal processing, control systems, biomedical engineering, finance and economics, electromagnetism, and fluid mechanics.For instance, the mathematical modelling of viscoelastic food ingredients experiencing stress and relaxation can be accomplished using fractional calculus [29].These applications showcase the adaptability and practical value of fractional derivatives across a range of scientific and engineering disciplines, enabling improved modelling capabilities and deeper comprehension of intricate phenomena.
In this article, we consider the space and time fractional cubic NLSE with the Kerr law nonlinearity with space and time fraction in the following form [30].
where U(x, t) is a complex-valued wave profile which is related to spatial co-ordinate x, and temporal variable t.In addition, r, s, andz are real coefficients with fractional parameters 0 < α ≤ 1 and 0 < β ≤ 1.The cubic nonlinear Schrödinger equation involving beta derivatives in space and time is used to model certain nonlinear optical phenomena.For example, it can describe the propagation of ultrashort optical pulses in nonlinear media with anomalous dispersion.
Here, by utilizing the equations x U, where i = √ −1, and assuming α = β with beta fractional derivatives, Equation ( 1) is transformed into the following form The model has been investigated in the previous literature using various methods, including Nucci's reduction method and the simplest equation method [31], the fractional Riccati expansion method [32], the fractional mapping expansion method [33], the (G /G)- expansion method [34], and the Adomain decomposition method [35].
To the best of our current knowledge, the extended direct algebraic method has not yet been applied to the model represented by Equation (1) to evaluate soliton solutions.The application of this method extends to various fields of nonlinear sciences, including mathematical physics, quantum physics, and engineering.However, the extended direct algebraic method is modified and implemented on a nonlinear space-time fractional model in Equation (1).By doing so, advanced, fresh, and wide-ranging soliton solutions are obtained.In this study, our primary focus is to establish advanced and widely applicable soliton solutions for the space-time fractional cubic nonlinear Schrödinger equation using the recommended method.The obtained soliton solutions exhibit wave-like behavior and are expressed in trigonometric, hyperbolic, and exponential forms.This research will also provide valuable insight into the internal formation of the travelling wave phenomena by depicting the obtained solutions in two-and three-dimensional graphs with different parameter values.Furthermore, the soliton solutions derived from this study will also contribute to the interpretation of complex phenomena associated with this particular space-time fractional model.
This article organizes its contents as follows: Section 2 presents the properties of the beta derivative.The algorithm of the proposed method is explained in Section 3. Section 4 provides a mathematical analysis.In Section 5, graphical representation and discussion are presented.The comparison scheme is outlined in Section 6, and finally, Section 7 concludes the article.

Definition of Beta Derivative and Its Properties
Several definitions of fractional derivatives, such as Riemann Liouville, the modified Riemann Liouville, the Caputo, the Caputo-Fabrizio, the conformable fractional derivative the Atangana-Baleanu derivatives, have been developed recently by many researchers [36,37].Most of the fractional derivatives do not agree with the well-known properties of classical calculus such as the chain rule, the Leibnitz rule, and the derivative of a constant is zero.Atangana et al. [38] launched a new crucial and progressive definition of fractional derivatives called beta derivative, which follows the fundamental properties of classical calculus.
Definition 1: Let α ∈ R and the function h = h(x) : [α, ∞) → R , then the beta derivative of order α with respect to x is defined as follows [39]: where Γ is the gamma function.D α x h(x) = d dx h(x) for α = 1.Properties: If h(x) and u(x) are α-order differentiable for all x > 0, and d 1 , d 2 are real constants, then the beta derivative encompasses the following properties [39]: dx .By using these properties of beta derivative, fractional differential equations simply turn into ordinary differential equations.As of now, the beta derivative has not been found to have any limitations and it fulfills all the properties associated with integerorder derivatives.Furthermore, it exhibits the property of yielding a derivative of zero for constant functions [40][41][42].The beta derivative is a non-local derivative that exhibits its distinctiveness when applied to functions that embody the entire characteristic of the function itself.It serves as a generalized version of the Caputo and Riemann-Liouville derivatives.In comparison to other derivatives, the beta derivative offers greater flexibility and can accurately model complex systems.Its applications are widespread, ranging from electrochemical systems and complex geometries to modeling electromagnetic waves in dielectric media and cancer treatment [43][44][45].Numerous scientific studies have reported the utilization of the beta derivative in diverse fields, further enhancing its appeal and prompting its application to real-world problems [40,41,46,47].

Algorithm of the Extended Direct Algebraic Method
In this section, we have presented the extended direct algebraic method as an effective technique.This method enables us to obtain fresh and wide-ranging analytical solutions for model (1).By employing this technique, fractional partial differential equations can be transformed into ordinary differential equations, simplifying the calculation process.The algorithm is narrated below: Step 1: Let the general form of the fractional order nonlinear evolution equation be where F is a polynomial of u(x, t), and D α t be fractional derivative of α-order and u(x, t) is the travelling wave variable, where subscripts denote partial derivatives.
Let us hypothesize about the travelling wave solution where u(ξ) is a function of ξ with In Equation ( 5) v, κ are respectively velocity and soliton frequency, ω is wave number and θ is the soliton phase component.
Inserting the above transformation into Equation (3), we obtain the following ordinary differential equation of integer order: where H be the polynomial of function u(ξ), and prime denotes the derivative with respect to ξ.
According to the new extended algebraic method, the solution of Equation ( 6) can be expressed in the form where c j (0 ≤ j ≤ N) are constant coefficients to be evaluated later and H(ξ) satisfies the following ordinary differential equation, where prime denotes derivative with respect to ξ, and µ, γ, λ are constant coefficients.The general solutions of Equation ( 8) (adequate solutions) are given in [27].By substituting Equations ( 7) and ( 8) into Equation ( 6), we obtain a polynomial of H(ξ).By extracting the coefficient terms of different powers of H j (ξ) where (j = 0,1,2,. . . ) and setting them equal to zero, then we obtain a system of algebraic equations with various parametric values such as such as c j (j = 0, 1, 2 . ..), µ, γ, λ, ω, and κ.By solving these algebraic equations, we can determine the values of the unknown parameters.Substituting these values of the parameters along with Equation (8) into Equation (7), as the broadspectrum solutions of Equation ( 8) are known, we obtain new and more general solutions.
For several values of µ, γ, λ and their correlation, Equation (8) gives disparate general solutions of NLSEs.

Mathematical Analysis
In this subsection, we studied the space and time fractional cubic NLSEs to find more general and standard exact wave solutions using an extended direct algebraic method.Furthermore, we discuss the mathematical analysis of the wave solutions.The fractional transformation in Equation (5) converts Equation (2) into the following ordinary differential equation, comprising both real and imaginary parts.
Now balancing between the highest order derivatives and highest power of the nonlinear term in Equation ( 9), we obtain N = 1.Therefore, the solution of Equation ( 9) is of the form By substituting the results from (8) into Equation ( 9) along with Equation ( 10), we obtain a polynomial equation in H(ξ), where (0 ≤ j ≤ N).Taking zero the resemble coefficient power of H j (ξ), we achieve a set of algebraic equations with c 0 , c 1 , µ, γ, λ.Calculating this set of algebraic equations with the software Mathematica, we obtain the values of the parameter as follows: where µ, γ, λ and r, s, k, z are free parameters.Now, embedding the values of ( 12) into (11) and the hypothesis of the auxiliary equation for different conditions, we establish the travelling wave solutions of (1) which are given below.
Case 9: While µ = 0 and γ = 0, we obtain where , with conditions (zs(2 Fractal Fract.2023, 7, 426 8 of 14 The soliton solutions obtained in this study are diverse and novel, originating from the general solutions.

Physical Significance and Explanations
In this section, attained soliton solutions of the space and time fractional cubic NLSEs are presented in Figures 1-5 and discussed the nature of these solitons for several values of unknown parameters through the software Mathematica.
In this section, attained soliton solutions of the space and time fractional cubic NLSEs are presented in Figures 1-5 and discussed the nature of these solitons for several values of unknown parameters through the software Mathematica.
The accomplished solutions are related to two parts including the real part and the imaginary part.The solutions provide various types of solitons such as kink shape soliton, singular kink shape soliton, V shape soliton, periodic soliton, flat kink shape soliton, antisingular kink shape soliton, soliton solutions and in such manners.The wave velocity and wave number have significant effects on the travelling wave profile.
The solution  6 exhibits kink-shaped soliton solution for the modulus part with velocity  = −1.443depicted in Figure 1          Taking other various values of free parameters this model provides the same type of soliton solutions, repeat solitons have been neglected here and solitons profile depend on the value of fractional order, wave velocity, wave number and other wave variables.The accomplished solutions are related to two parts including the real part and the imaginary part.The solutions provide various types of solitons such as kink shape soliton, singular kink shape soliton, V shape soliton, periodic soliton, flat kink shape soliton, antisingular kink shape soliton, soliton solutions and in such manners.The wave velocity and wave number have significant effects on the travelling wave profile.
The solution u 27 exhibits singular bell shape soliton for modulus part with velocity v = 4.1 for the value of α = 0.9, k = 3 depicted in Figure 4, but while decreasing the value α = 0.55, 0.25, k = 1.68 then the wave velocity become v = 2.79 and those graphical representation are provided in Figure 4b,c.Portraits of 3D are shown within the interval 0 ≤ x ≤ 5 and 0 ≤ t ≤ 5, and 2D.Portraits are shown at t = 1 with values of arbitrary parameters s = 2, γ = 0, r = 1.8, µ = 1.
Taking other various values of free parameters this model provides the same type of soliton solutions, repeat solitons have been neglected here and solitons profile depend on the value of fractional order, wave velocity, wave number and other wave variables.

Comparison
We compare the results of the space-time fractional cubic nonlinear Schrodinger equation obtained through the extended algebraic equation method with Abdelwahed et al. [30] solutions.It is noticed that from the attained results few of them are analogous to the results that were established earlier by several approaches and some of them are fresh.

Solution Using the Simplest Equation Method
The Attained Solutions q(x, t) = − In the above table, we discussed the solutions obtained in this paper with the previous study.Hashemi et al. [30] have given more than two solutions which are not homologous with the attained results.

Conclusions
The extended direct algebraic method has been used to derive novel exact analytical soliton solutions of the cubic nonlinear Schrödinger equation with fractional space-time terms.In this study, Kerr's law nonlinearities are utilized, which arise from the nonharmonic motion of bound electrons when light pulses propagate in optical fibers.All solutions are expressed in terms of trigonometric and hyperbolic functions.Computational calculations and graphical representations of the solutions are plotted using the Wolfram Mathematica software.The graphical representations of these solutions help us to visualize and understand the internal features of the system more accurately.Among these solutions, some are new and have not been reported previously in the literature.

Figure 1 .
Figure 1.Three-dimensional and two-dimensional plot of kink shape soliton solution of u 6 .Fractal Fract.2023, 7, x FOR PEER REVIEW 9 of 14

Figure 2 .
Figure 2. The 3D and 2D plots of the periodic travelling wave solution u 6 .

Figure 3 .
Figure 3.The 3D and 2D plots of V shape soliton solution  26 .The solution  27 exhibits singular bell shape soliton for modulus part with velocity  = 4.1 for the value of  = 0.9,  = 3 depicted in Figure 4, but while decreasing the value  = 0.55, 0.25,  = 1.68 then the wave velocity become  = 2.79 and those graphical representation are provided in Figure 4b,c.Portraits of 3D are shown within the interval 0 ≤  ≤ 5 and 0 ≤  ≤ 5, and 2D.Portraits are shown at  = 1 with values of arbitrary parameters  = 2,  = 0,  = 1.8,  = 1,

Figure 4 .
Figure 4.The 3D and 2D plot of singular bell shape soliton solution corresponds to u 27 .

Figure 5 .
Figure 5.The 3D and 2D plot of the soliton solution corresponds to u 35 .