Dynamical behaviours and soliton solutions of the conformable fractional Schrödinger–Hirota equation using two different methods

The nonlinear conformable fractional Schrödinger–Hirota (NCFSH) equation, describing the propagation of optical solitons in a dispersive optical fibre, is considered and its exact solutions are investigated. Abundant solutions including dark, bright and singular solutions have been obtained by two distinct methods, namely direct method and trial equation method. These methods provide a direct and efficient strategy for the NCFSH equation. The solutions have been simulated by graphs to check the types of solitons and their dynamical behaviours. With back substitution using Mathematica, all the solutions have been verified. Obtaining optical solutions to the proposed equation will play a fundamental role due to the applications in quantum mechanics, fluid dynamics, nonlinear optics and various other branches of mathematical physics.


Introduction
In recent times, fractional derivatives are more preferred than classical derivatives to model real-world problems. Fractional derivatives are commonly seen in conservative and nonconservative phenomena. To model processes consisting of electromagnetic waves, electro-electrolyte polarization, diffusion wave and heat conduction, etc. fractional derivatives can be used [1].
The class of Schrödinger equations is generally modelled phenomenon, especially in quantum mechanics, energy and energy quantization. The linear Schrödinger equation represents the time evolution of wave function, whereas the nonlinear form is the most common equation for the evolution of slowly changing quasi-monochrome wave packets in weak nonlinear media with dispersion [2]. It models various nonlinearity effects in fibre such as self-phase modulation, second harmonic generation, ultrashort pulses, optical solitons, four-wave mixing, stimulated Raman scattering and so on. Laskin proposed the fractional Schrödinger equation which is an essential equation of fractional quantum mechanics [3].
In this paper, the nonlinear conformable fractional Schrödinger-Hirota (NCFSH) equation is considered as [4][5][6][7] iu (μ) The NCFSH equation describes real-world models in nonlinear optics and dispersive optical fibres [6,7]. Due to the applications in various fields, especially in optical communication systems, this equation and its solutions which enable to interperet the phenomena it represents taken much attention recently. For the fractional operator, many approaches are known based on Gamma function. Therefore, in Equation (1), the fractional operator is considered the conformable derivative operator [8]. Most of the properties of the conformable derivative correspond to the classical derivative [9] and this definition allows to solve nonlinear partial fractional differential equations (NPFDEs) more easily [10].
The conformable derivative of order μ for the function f is given by [8] T where f : [0, ∞) → R, ∀t > 0, μ ∈ (0, 1). All of the classical derivative rules, such as linearity, sum, product, division, etc. are same as the conformable derivative. Additionally, if μ ∈ (0, 1) and f is differentiable, dt can be determined [8]. Some authors [11,12] studied conformable fractional derivative and integral operators and gave geometric meaning of the conformable derivative [13]. The analysis of asymptotical stability for conformable fractional systems is given in recently published paper [14]. Recent studies [15][16][17] show attention to the conformable fractional derivative and its applications.
In this work, in addition to studies in [4][5][6][7], the TEM and the direct method are considered in the view of the definition and rules of the conformable derivative to attain the analytical solutions of the NCFSH equation.

Briefs of methods and results
In this section, briefs of methods and results obtained with these methods are given.

Trial equation method (TEM)
First, a nonlinear PDE, is reduced to an ordinary differential equation, under the transformation Second, consider the trial equation where a i and n are constants, which are derived from the solution of the system and the balancing principle, respectively. Finally, the solution of Equation (5) can be given by the integral form The roots of the denominator can be classified and Equation (6) can be solved by a complete discrimination system for polynomial. Therefore, the exact solution of Equation (3) is hold [32].
Also, the trial equation can be chosen as in the modified trial equation method [33,34], then the same procedure can be applied with the integral form solution of Equation (7) as Similarly, one can solve integral (8) and find the exact travelling wave solutions of Equation (3).

Results for NCFSH Equation via the TEM
Consider the suitable wave transformation where k and η are constants. Using the above transformations and substituting the derivatives in Equation (1), and in point of the real and imaginary parts [7], k = −1/3λ is obtained. Also, z(ξ ) satisfies the following ordinary differential equation: Balancing z and z 3 , we obtain n = 4. So, trial equation is Using the solution procedure of the algorithm yields a system of algebraic equations: As a result of the system, coefficients are determined as Equation (6) is rewritten with Equation (12), If we set a 0 = 0 in Equation (13)  If 2η 3 + 5 81λ 2 > 0, then we have bright and singular solutions, respectively: u(x, t) = ± − 2η + 5 27λ 2 × csch 2η 3 If 2η 3 + 5 81λ 2 < 0, the singular periodic solutions appear as  If we now use the trial Equation (7) in Equation (10) for balancing procedure, we get r − s = 2. Choosing r = 2 and s = 0, the trial equation is The corresponding system is: Solving the corresponding system, we obtain  Using these coefficients in Equation (8), we get When we integrate this equation and use the wave transformation, the exact travelling wave solutions of Equation (1) are obtained as follows: If 5 2 + 27ηλ 2 < 0, then we have dark optical solutions: If 5 2 + 27ηλ 2 > 0, then the solutions purport periodic solutions: u(x, t) = ± i 3λ u(x, t) = ± i 3λ

Direct method
If f is differentiable, to reduce NPFDEs, the following rule of conformable fractional derivative, can be used with the transformation As a result, the reduced equation is nonlinear ordinary differential equation. The reduced equation can be solved via analytical, semi-analytical and numerical methods.

Results for NCFSH equation via the direct method
Now, we consider the direct method to obtain analytical solutions of the NCFSH equation.
the following equations are attained due to imaginary and real parts, respectively: The analytical solution of the first equation is Substituting the analytical solution z(ξ ) into the second equation and setting the coefficients of exponential function to zero gives the following system: Using Maple, the above system has 18 solutions. However, some of the results are trivial solutions and some of them do not match up to our aim. One of them satisfies all conditions, which is Therefore, the exact solution of the considered Equation (1) is the following equations are attained due to imaginary and real parts, respectively. λz Considering σ (ξ) = ξ 6λ + C 1 , the analytical solution of the first equation is Substituting the analytical solution z(ξ ) into the second equation, Therefore, the exact solution of the considered Equation (1) with general transformation is where C 2 is given above.

Conclusion
In this work, TEM and the direct method are considered to investigate exact solutions of the NCFSH equation.
To this purpose, NCFSH equation is converted into an integer order ordinary differential equation by using the rules of the conformable derivative and the considered transformations. Using these methods, Maple and Mathematica, the optical wave solutions including various types of solitons that can be useful in nonlinear optics and telecommunication industry are obtained. Physical meanings and dynamical behaviours of the solutions are exhibited by using figures. Results indicate that these methods are powerful tools for NPDEs. Obtained solutions represent some real-life situations. Additionally, as seen in the proposed results, these suggested methods can be implemented to various fractional models arising in other areas of physics and science.

Disclosure statement
No potential conflict of interest was reported by the author(s).