Optical solitons of space-time fractional Fokas–Lenells equation with two versatile integration architectures

Nonlinear Schrödinger’s equation and its variation structures assume a significant job in soliton dynamics. The soliton solutions of space-time fractional Fokas–Lenells equation with a relatively new definition of local M-derivative have been recovered by utilizing improved tan( 2 )-expansion method and generalized projective Riccati equation method. The obtained solutions are periodic, dark, bright, singular, rational, along with few forms of combo-soliton solutions. These solutions are given under constraints conditions which ensure their existence. The impact of local fractional parameter is featured by its graphical portrayal. 2D and 3D diagrams are drawn to illustrate the efficacy of the conformable fractional order on the behavior of some of those solutions. The secured solutions of this model have dynamic and significant justifications for some real-world physical occurrences. Our study shows that the suggested schemes are effective, reliable, and simple for solving different types of nonlinear differential equations.


Introduction
From the past three decades, optical solitons emerge as a fast growing area of research due to their use in transmission technology, through different forms of wave-guides. Solitons are utilized to represent the particle-like properties of nonlinear pulses. The importance of solitons is due to their presence in a variety of nonlinear differential equations portraying many complex nonlinear phenomena, including acoustics, nonlinear optics, telecommunication industry, convictive fluids, plasma physics, condensed matter, and solid-state physics. Nonlinear Schrödinger's equation and its variant forms are used in dispersive mediums in different fields of mathematical physics and have been studied mathematically in recent years [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15].
Solitons exist due to an accurate balance among nonlinearity and group velocity dispersion (GVD) in the area. If the value of GVD is small, this balance may be at risk. Therefore, to keep the balance among the two, expression terms with dispersive effects need to be in-vestigated. One of the known models which is relevant is the Fokas-Lenells (FL) equation [16][17][18]. FL equation is one of the known forms of nonlinear Schrödinger's equation, introduced about a decade ago. Due to its vast applications in fiber-optics communication, this sort of equation is crucial in research and thus the search for different forms of wave solutions is very significant. Many computational methodologies have been developed for constructing wave solutions for such equations, including the improved tan(φ(η)/2)expansion method, the trial equation method, improved Bernoulli subequation function method, the extended Fan subequation method, extended and the modified simple equation methods, Riccati-Bernoulli's sub-ODE method, the Lie group analysis, extended Jacobi's elliptic function approach, and many others [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35].
Nowadays, fractional science is a flourishing area of mathematical analysis along with fractional operators, such as Caputo, Grunwald-Letnikov, and Riemann-Liouville [36][37][38][39][40][41][42][43][44]. In 1695, in a letter to Leibniz, l'Hospital asked him about the detection of expanding the sense of an integer-order derivative d λ y dx λ to the case of a fraction of the order. This problem started the development of a modern calculus that was named the calculus of arbitrary order and is now commonly called the calculus of fractions. Many forms of fractional derivatives, including Riemann-Liouville, Caputo, Hadamard, Caputo-Hadamard, and Riesz [40,41,43], have been developed to date. Almost all of these derivatives are described in the Riemann-Liouville sense based on the corresponding fractional integral.
In 2017, Sousa introduced a new fractional derivative that generalizes the so-called alternative fractional derivative [45]. This new differential operator is denoted by D λ,μ M , where λ is the order, such that 0 < λ ≤ 1, μ > 0, and M denotes that the derived function includes a Mittag-Leffler function along with one parameter. This new type of derivative is known as a local M-derivative, it fulfills certain characteristics of integer-order calculus, e.g., linearity, quotient rule, product rule, chain rule, and function composition. Furthermore, the local M-derivative of a constant is zero. Since the Mittag-Leffler function is the generalization of the exponential function, some of the classical outcomes of calculus of the integer-order can be extended, namely the mean value theorem, Rolle's theorem, and its extensions. Moreover, when the derivative order is λ = 1 and the Mittag-Leffler function parameter is also unitary, our specification is analogous to that of the ordinary derivative of order one.

Governing model
Using the definition of the local M-derivative and its properties, the space-time fractional FL equation is introduced as follows: ) is a complex-valued wave function. The first term of Eq. (1) gives the fractional temporal evolution of the pulse; a 1 , a 2 are the spatio-temporal dispersion (STD), group velocity dispersion (GVD) coefficients, while ρ, δ, and γ represent the self-steepening, inter-modal dispersion (IMD), and nonlinear dispersion (ND) coefficients respectively [16][17][18].
When α = 1, Eq. (1) is converted to the original FL equation [16][17][18]. The FL equation emerges as a model equation that defines the nonlinear pulse propagation in optical fibers by maintaining terms up to the next leading asymptotic order (the nonlinear Schrödinger equation (NLSE) results in the leading asymptotic order). In the context of nonlinear optics, this equation sculpts the promulgation of nonlinear light pulses in monomode optical fibers as assumed nonlinear higher-order effects are captured in the elaboration [51]. It is worth noticing that the FL equation is a fully integrable equation in nonlinear PDEs which has been developed as an integrable generalization of the NLSE using bi-Hamiltonian techniques [52].

Local M-derivative
Consider the function f : [0, ∞) → R where t > 0. For 0 < μ < 1, let us define the local M-derivative of order μ for the function f , denoted by D μ;δ M [53][54][55][56], by where E δ (·) is the Mittag-Leffler function with one parameter. Here The local M-derivative possess the following properties: This local fractional-order M-derivative also has the following chain rule property: Using Eqs. (4)- (6), we obtain the following expression: with where m is a constant. The last property is given as

Traveling wave hypothesis
Consider the following complex traveling wave transformation: Substituting Eq. (10) into Eq. (1) yields from the real part, and from the imaginary part.

Soliton solutions
In this section soliton solutions are extracted for Eq. (13) with the help of two different integration schemes, namely the improved tan( φ(η) 2 )-expansion method and generalized projective Riccati equation method. In order to obtain these solutions, it is enough to solve the real part of Eq. (13), which is given in Eq. (14).
Substituting these values in Eq. (19) and using the relation in Eq. (18) yields the following soliton solutions for Eq. (13).

Generalized projective Riccati equation method
Consider the following initial solution in order to solve Eq. (14) to find soliton solution of Eq. (13) with the aid of generalized projective Riccati equation method [49,50]: In Eq. (14), the homogeneous balance principle gives N = 1. Equation (34) becomes where (η) and τ (η) satisfy the projective Riccati system (PRS) and PRS first integral is expressed as Substituting Eqs. (35), (36), and (38) into Eq. (14) provides a polynomial in ( i (η), τ j (η)), setting whose coefficients to zero yields a set of algebraic equations. When solving this system of equation with the aid of Maple, we will get Set 1.

Results and discussion
This section deals with graphical demonstration of the obtained results and provides a brief discussion on the effect of fractional parameter α. Figure 1 ter. A graphical illustration of the dark soliton solution | 2 (x, t)| can be viewed in Fig. 2(a), and its 2D fractional parameter effects are depicted in Fig. 2(b).

Conclusion
An M-fractional FL equation representing the propagation of short light pulses in the monomode optical fibers is investigated using the improved tan( φ(η) 2 )-expansion method and GPREM. The FL model is a higher-order nonlinear Schrödinger shape equation that gives bright soliton solutions with internal freedom. Furthermore, the dark soliton solutions for the equation with the M-fractional effect, which have no internal freedom and exist for both focusing and lake of focusing equations, are investigated. The improved

tan( φ(η)
2 )-expansion method is used to extract dark, singular, and rational soliton solutions, and GPREM provides dark, singular, periodic, and some forms of combo soliton solutions. These solutions are also demonstrated graphically. The fractional parameter effect on the dispersion is also highlighted through 2D graphical representation. The reported outcomes are useful in the empirical application of fiber optics. The essential advantages of the proposed schemes over all the other methods are that these methods provide new explicit analytic wave solutions including many real free parameters. The closed-form wave solutions of the nonlinear PDEs have their significant meaning to reveal the interior device of the complex physical phenomena. More problems in applied mathematics, mathematical physics, and engineering might be solved through the presented methods. In the future work, we will find the multisoliton solutions for the FL equation by the aid of the Hirota method.