Abundant Exact Travelling Wave Solutions for a Fractional Massive Thirring Model Using Extended Jacobi Elliptic Function Method

: The fractional massive Thirring model is a coupled system of nonlinear PDEs emerging in the study of the complex ultrashort pulse propagation analysis of nonlinear wave functions. This article considers the NFMT model in terms of a modiﬁed Riemann–Liouville fractional derivative. The novel travelling wave solutions of the considered model are investigated by employing an effective analytic approach based on a complex fractional transformation and Jacobi elliptic functions. The extended Jacobi elliptic function method is a systematic tool for restoring many of the well-known results of complex fractional systems by identifying suitable options for arbitrary elliptic functions. To understand the physical characteristics of NFMT, the 3D graphical representations of the obtained propagation wave solutions for some free physical parameters are randomly drawn for a different order of the fractional derivatives. The results indicate that the proposed method is reliable, simple, and powerful enough to handle more complicated nonlinear fractional partial differential equations in quantum mechanics.


Introduction
Physics can be typically classified into two branches: classical and modern physics. Modern physics can be distinguished by considering spatiotemporal requirements for joint interaction, whereas, in classical physics, we can consider time and space separately because they are independent and absolute. Furthermore, classical physics usually deals with the macroscopic scale, while modern physics deals with microscopic or sub-microscopic scales. Although classical physics has different applications in science and engineering, modern physics can be considered a revolution in applied physics, as it can elucidate many essential phenomena, such as black body radiation, photoelectric effect, Compton's effect and stability of atoms that cannot be explained from a classical physics point of view. However, modern physics focuses on quantum mechanics and the theory of relativity; quantum mechanics considers the physical quantities restricted to be discrete values, where the thinking of the probability is dominant instead of certain measurements, which is represented mathematically by the Schrödinger wave equation. The theory of relativity studies the physical quantities moving at a speed near the speed of light, the time dilation, and the dimensions contraction started to be important concepts, and Einstein's massenergy equation makes a revolution in the science [1][2][3][4].
The contemporary revolution in theoretical and applied physics combines quantum mechanics with the theory of relativity in a multi-body system, which establishes quantum field theory. Quantum field equations represent a general form of the Schrödinger wave equation, where the wavefunction is generalized to an infinite-dimensional space of field configurations [2]. Motivated by this, in this work, we consider the massive Thirring model (MTM) as an important application of the quantum field theory, which was derived by W. Thirring in 1958 [3]. Thereafter, many theoretical and applied studies of such a complex system were conducted. For example, but not limited to, Kondo 1995 studied the bosonization and duality of the MTM with a four-fermion interaction of the current type [4], and the Thirring model was also considered in a separate work as a gauge theory [5]. In 2018, Joshi et al. introduced an integrable semi-discretization of the MTM for the first time in laboratory coordinates [6].
Nevertheless, to find out an alternative methodology for the Schrodinger equation, Dirac discovered the integral path approach, similar to Lagrangian's least-action principle technique in classical mechanics; this approach was developed by Feynman to create Feynman diagrams. Feynman diagrams were modified to Wiener's path integral, which is equivalent to the Brownian path integral in classical mechanics. Recently, the Levy flight random process has been introduced to understand difficult classical and quantum physics phenomena, where the Levy index α is introduced. Now, the consequences of the path integral for the Levy flight paths' studies are an essential issue in fractional quantum mechanics and consequently in fractional quantum field theory [7]. Examining research to obtain novel and additional exact traveling-wave solutions for fractional models is prospering. Indeed, this is not an easy task and is one of the pivotal challenging problems in mathematics and physics. Hence, resorting to sophisticated analytical and digital methods is inevitable. In this direction, many effective and accurate analytical methods for solving these equations have been considered thus far, for example, the Bäcklund transformation method, the Riccati sub-equation method, the extended tanh-function method, the G'/Gexpansion method, the Kudryashov method [8][9][10][11] and so forth.
The analysis in this paper highlights the complex behavior of nonlinear wavefunction, which is notably dependent on the genetic properties and temporal memory that can be explored with great skill using fractional calculus.
In this direction, consider the following semi-discrete nonlinear massive Thirring model (MTM) that can be typically provided by a complex triple system of difference equations: where n denotes the discrete lattice to index iterates, ν denotes the lattice-spacing parameter of aspace discretization, and the symbol i is an imaginary unit. The complex-conjugates of ψ n and φ n are denoted respectively byψ n andφ n . The first equation refers to the case of temporal evolution, while the last two difference equations refer to the semi-discrete massive Thirring equations constrained with the components of {ψ n } n∈Z and {φ n } n∈Z , which can be defined in terms of {χ n } n∈Z in the temporal and spatial coordinates [6]. With the continuity of ν → 0 , the slowly changing solutions between the lattice nodes can be written as: where continuous variables fulfill the following system of partial equations: which leads to an MTM system of two semi-linear equations for (F, G) ∈ C 2 in terms of the variables ψ(x, t) = F(x, t − x) and φ(x, t) = G(x, t − x) in the normalized form: This paper deals with the fractional version of such a system. Therefore, we consider the following nonlinear space-time fractional MTM system: Considerable analytical and numerical investigations of the MTM have been made in the literature using various techniques. The construction of the MTM using the functional integral scheme within quantum field theory was discussed in [12]. The physical states, as well as a solution of the MTM, by means of many-body wave functions, are presented in [13]. In [14], Bethe ansatz solutions of the MTM were tested numerically by solving periodic boundary value problems. Delepine et al. [15] demonstrated that the MTM is equivalent to the quantum sine-Gordon model in quantum field theories at a finite temperature. The white noise of the oscillator MTM was examined in [16] in terms of the phase-space displays. In [17], the non-thermal phase structure of the MTM was studied using ansatz matrix-product states. Using the N-fold Darboux transform, the rogue wave solutions of the MTM equations were derived in [18]. On the other side as well, the fractional versions of the nonlinear complex MTM were numerically solved using advanced semianalytical and approximate methods; for example, the q-HAM was applied in [19] to solve the fractional massive Thirring model in Caputo sense. In [20], the fractional residual power series method was implemented to solve a class of the fractional MTM with conformable derivatives. For more details regarding the numerical and analytical solutions of different fractional models, we refer to [21][22][23][24][25][26][27][28][29][30][31].
Almost all scientific problems can be solved using different fractional calculus techniques, where one or many suitable methods can be chosen for each problem; some problems that are solved using modified Riemann-Liouville fractional calculus techniques were noted as incorrect conditions [32][33][34][35]. Although these cases are not related to this work, this note must be mentioned here. These cases do not affect the Riemann-Liouville fractional calculus technique, which solves a huge number of problems successfully, as do other methods, such as the Mittag-Leffler function, the fractional Riccati method, the fractional double function method, and the fractional Y-function expansion method [36,37].
The novelty of this paper is to explore new travelling wave solutions for fractional MTM equations (5) by employing an effective analytic approach based on a complex fractional transformation and Jacobi elliptic functions. It is worth noting that the previous study of soliton for the fractional MTM equations was performed to provide an approximate solution for or study a special case of the MTM equations [4][5][6]. This paper introduces a general case exact solution for MTM equations for the first time; this study can be considered as a strong motivation to provide the obtained results.
The outline of this analysis has the following sections: In Section 2, some basic definitions and characteristics of the considered fractional operator are presented. In Section 3, the key idea of the proposed method is described. Then, in Section 4, we apply this method to create new sets of exact traveling wave solutions to the fractional massive Thirring model. Finally, a brief conclusion is also provided.

Preliminaries
Recently, many researchers have used various fractional operators to study several models associated with the functions of complex variables, and they proved that these fractional operators are more influential than the classical ones while analyzing the natural behavior of those models. Herein, we introduce the basic definition and some properties of Jumarie's modification of Riemann-Liouville derivative [38][39][40][41] that are very useful for displaying this work in a standardized way.

Definition 1.
Let ω : R → R be a continuous function. Then, the modified Riemann-Liouville derivative of the order α is as follows In this work, if ω(t) has a modified Riemann-Liouville derivative of the order α, it will be defined as D α t -differentiable. Further, it is obvious that the operator D α t of Jumarie's modification satisfies the following interesting properties: t -differentiable and defined in the range of ω 1 . Then, we have:

The Extended Jacobi Elliptic Equation Method
This section presents the definition of Jacobi elliptic functions and reviews some important properties that we will use within the framework of this paper. Then, we introduce the algorithm of the proposed method.

The Jacobi Elliptic Functions
The Jacobi elliptic functions are the standard forms of elliptic functions. There are three double periodic functions, namely the Jacobian elliptic sine function, Jacobian elliptic cosine function, and Jacobian elliptic function of a third kind denoted by sn(u, δ) = sn(u), cn(u, δ) = cn(u) and dn(u, δ) = dn(u), respectively, where δ is the elliptic modulus. In the next segment, we provide the details of the derivation of these functions. To this end, we consider the following nonlinear partial differential equation (PDE): By substituting the linear transformation η = θ(x − µt) into PDE (6), we get the following nonlinear ordinary differential equation (NODE): Then, some simplifications lead to the following equivalent NODE: Letting c = 1, −λ/θ 2 µ = −δ 2 and ω = ϕ/2, then Equation (8) takes the form which is equivalent to 1 where the integral in Equation (10) is called the Legendre elliptic integral of the first kind. Now, we define Provided that u = f (t) so that t = f −1 (u) = sn(u), which is the Jacobi elliptic sine function. Nevertheless, the Jacobi elliptic cosine function can be defined by letting Provided Consequently, one can write the following argument On the other side as well, Jacobi elliptic functions sn(u), cn(u) and dn(u) can be defined respectively as solutions to and possess the following properties in terms of their singular points: when δ → 1 , the Jacobi elliptic functions turn into hyperbolic functions as follows when δ → 0 , they turn into trigonometric functions as follows Furthermore, one can obtain the following identities: The derivatives of the Jacobi elliptic functions are as follows

Extended Jacobi Elliptic Function Expansion Method
Herein, the algorithm of the extended Jacobi elliptic function expansion method will be illustrated to obtain the exact travelling wave solutions of NFPDEs. To this end, let us consider FPDE in the the following form where u = u(t, x, y), P is a polynomial in u, and its partial derivatives, including fractional derivatives, D α t , D β x and D γ y , are a modified Riemann-Liouville derivative of u with respect to the independent variables t, x and y. In the following, the main steps of the proposed algorithm are presented to find out the exact travelling wave solutions of FPDE (16): Step 1. Use the fractional wave transformation where v is the wave velocity that will later be determined. This permits us to reduce FPDE (16) into the following ODE of integer order in terms of ξ: Step 2. Propose that Equation (18) has a solution in the following form where = 1, 2, . . . , 12, in which L is a positive integer,b 0 , b 1 , b 2 , b 3 , and p , r , = 2 , 3, . . . , L and are constants to be determined afterwards. The functions Q (ξ), R (ξ) and S (ξ), = 1, 2, . . . , 12 can be expressed in terms of Jacobi elliptic functions (15) as follows where ρ is an arbitrary constant.
Step 3. Determine the integer L in the predicted solution (19) by balancing the highest order nonlinear terms and the highest-order derivatives O d l 2 dξ l 2 U = L + l 2 for l 2 = 0, 1, 2, . . . .
Step 5. Solve the resulting algebraic system in Step 4 with the aid of Mathematica software to find out the values of ρ, k, b 0 , b 1 , b 2 , b 3 , and p , r , = 2 , 3, . . . , L.
Step 6. Substitute the obtained values in terms of ρ, k, b 0 , b 1 , b 2 , b 3 , and p , r for = 2, 3, . . . , L in the predicted solution (19); new types of abundant traveling wave solutions are provided to FPDEs (16) involving the Jacobi elliptic functions.

Conclusions
In this paper, the fractional massive Thirring model has been considered in the sense of the modified Riemann-Liouville fractional derivative. Based on the nonlinear fractional complex transformation, a series of exact traveling wave solutions for this model has been successfully obtained in terms of Jacobi elliptic functions. With the aid of the Mathematica wolfram computation package, the resulting algebraic system of free parameters was Figure 8. The 3D plots of the real and imaginary parts of F 13 (x, t) and G 13 (x, t) with the parameters r 1 = r 2 = −3, k 1 = k 2 = 1, q 0 = 7, p 2 = 2, q 1 = −5, ρ = 0.5, and τ = 0.1 for the fractional order

Conclusions
In this paper, the fractional massive Thirring model has been considered in the sense of the modified Riemann-Liouville fractional derivative. Based on the nonlinear fractional complex transformation, a series of exact traveling wave solutions for this model has been successfully obtained in terms of Jacobi elliptic functions. With the aid of the Mathematica wolfram computation package, the resulting algebraic system of free parameters was Figure 9. The 3D plots of the real and imaginary parts of F 14 (x, t) and G 14 (x, t) with the parameters r 1 = r 2 = 0.5, k 1 = k 2 = 1, q 0 = 10, p 2 = 0.2, q 1 = 1, ρ = 1 and τ = 0.1 for the fractional order

Conclusions
In this paper, the fractional massive Thirring model has been considered in the sense of the modified Riemann-Liouville fractional derivative. Based on the nonlinear fractional complex transformation, a series of exact traveling wave solutions for this model has been successfully obtained in terms of Jacobi elliptic functions. With the aid of the Mathematica wolfram computation package, the resulting algebraic system of free parameters was solved and graphical representations of some acquired solutions were performed in 3D plots. The proposed method provides a powerful and systematic tool for obtaining novel exact solutions and can be applied to deal with other governing nonlinear fractional evolution equations emerging in mathematical physics.