Lie Symmetry Analysis, Analytical Solution, and Conservation Laws of a Sixth-Order Generalized Time-Fractional Sawada-Kotera Equation

: To discuss the invariance properties of a sixth-order generalized time-fractional Sawada-Kotera equation, on the basis of the Riemann-Liouville derivative, the Lie point symmetry and symmetry reductions are derived. Then the power series theory is used to construct the exact power series solution of the equation. Finally, the conservation laws for a sixth-order generalized time-fractional Sawada-Kotera equation are computed.


Introduction
As is well-known, constructing exact solutions of a partial differential equations (PDEs) is a vital theme in nonlinear science. A considerable number of methods have been developed, such as Lie symmetry analysis, the Painlevé test, inverse scattering transformation (IST), Darboux transformation (DT), and Clarkson-Kruskal (CK) transformation. The Lie symmetry analysis method is generally preferred for obtaining exact solutions of a PDE. This method was introduced by Sophus Lie in the early 19th century.
Recently, the symmetry analysis of fractional partial differential equations (FPDEs) and fractional derivatives were proposed [1]. On this basis, several works have been published on time-fractional differential equations. Baleanua et.al obtained the exact traveling solutions of the time-fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation [2]. Saberi and Hejazi obtained some exact solutions of the time-fractional generalized Hirota-Satsuma-coupled KdV system by means of the invariant subspace method [3]. The coefficient of a fractional-order equation was extended from a constant coefficient to a variable coefficient, and the four generators of the equation were obtained by a classified discussion of the variable coefficient [4]. Two diverse approaches were used to obtain the similarity solutions of the time-fractional Burgers system, and the similarity solutions were approximated by a numerical simulation method [5]. Explicit analytical solutions of several time-fractional equations were obtained by using the singular manifold method [6]. Roul and Prasad Goura designed a high-order numerical method to approximate the solution of time-fractional fourth-order PDEs. The result showed that this new method was more accurate than the previous method [7]. A Crank-Nicolson Legendre spectral method was developed to solve the two-dimensional nonlinear time-fractional diffusion-wave equation [8]. Zou et al. considered the following generalized time-fractional Sawada-Kotera equation [9]: The analytical solution and conservation laws (CLs) were then derived. In this work, we consider a sixth-order generalized time-fractional Sawada-Kotera equation: where u = u(x, t),a,b,c,d are arbitrary constants, and ∂ α u ∂t α is the Riemann-Liouville (RL) fractional derivative of order α about t. Equation (2) is a vital mathematical model that appears in many differential physical environments. It has been widely used in microscopic particle phoronomics, conformal field theory, and nonlinear optics. The long-wave motion in shallow water under the action of gravity and the long-wave motion in a one-dimensional nonlinear lattice can also be described. Equation (2) is a extension of Equation (1), which is more sophisticated and powerful than Equation (1). It is a new equation that has not been studied by other scholars. Therefore, Equation (2) enriches the study of time fractional PDE and has a certain significance of research value.
With α = 1, Equation (2) becomes whose Lie symmetry analysis and exact solutions were studied in [10]. The rest of the article is arranged as follows: in Section 2, some information about the Lie group method is given to analyze the FPDE. In Section 3, symmetry reductions are investigated on the basis of the symmetry groups of the FPDE. In Section 4, by using the power series method, we obtain the explicit analytical power series solution. In Section 5, conservation laws are derived using Ibragimov's nonlocal conservation theorem and the fractional Noether operators. Lastly, several remarks and conclusions are given in Section 6.

Lie Symmetry Analysis
First, let us recall the definition of the RL fractional derivative of order α: Definition 1 ( [11,12]).
Let us also consider the following scalar time FPDE: where α(0 < α < 2) is a parameter . Next, we focus on a one-parameter Lie group of infinitesimal transformations as follows: where D x is the total differential operator and is defined as follows: Then, the associated Lie algebra of symmetries is generated by the vector field below: On the basis of the infinitesimal invariance criterion, we can easily obtain Additionally, the variance condition yields [13] τ(x, t, u)| t=0 = 0.
The αth extended infinitesimal about the RL fractional time derivative with Equation (9) is Let us assume that Equation (2) is invariant under a one-parameter transformation (6); then, the transformed equations are provided that u = u(x, t) satisfies Equation (2). From the invariant Equation (8), we have Inserting Equations (6) and (10) into Equation (13), we obtain the following determining equations: Solving the above equations, we get where c 1 and c 2 are arbitrary constants. Hence, the infinitesimal symmetry group for Equation (2) is generated by the following two vector fields: It is easy to verify that the above generator is closed under the Lie bracket:

Symmetry Reductions
Next, according to the Equation (14), we consider the similarity reductions of Equation (2).
Therefore, we can easily construct the expression u = f (t) by solving the above characteristic equation. If this expression is put into Equation (2), then the fractional differential equation is Consequently, Equation (2) has a group invariant solution of the form where k is an arbitrary constant of integration.
Solving the above equations, we get Theorem 1. Putting Equation (15) into Equation (2), we obtain the following reduced equation: with the Erdélyi-Kober (EK) fractional differential operator where the EK fractional integral operator is (18) Proof. Let n − 1 < α < n, n = 1, 2, 3 . . . . Then, the similarity transformation of RL fractional derivatives is By setting v = t s , we have ds = −( t v 2 )dv. Then, Equation (19) can be written as If we apply the definition of EK integral operator (18) to formula (19), then we obtain the following expression: Let us simplify Equation (21) further: considering ξ = xt − α 6 , φ C 1 (0, ∞), we can obtain Thus, we have If we perform the same step n − 1 times, then the following result can be obtained: Putting operator (17) into Equation (24) leads to the following expression: Then, we can use the form of the EK fractional differential operator to express the time-fractional derivative: Given the above, Equation (2) is transformed into the equation below:

Exact Power Series Solution
Here, on the basis of the power series theory and symbolic computation, a kind of exact power series solution of Equation (27) is obtained with a detailed derivation. If we set then we can obtain Substituting Equations (28) and (29) into Equation (27), we obtain Observing coefficients in Equation (30) carefully, when n = 0, we have When n ≥ 1, we get Then, the explicit solution for Equation (27) can be expressed in the following form: Thus, the exact power series solution for Equation (2) is (32)

Conservation Laws
In this section, we investigate the CLs of Equation (2). The RL left-sided time-fractional derivative is applied as The left-sided time-fractional integral of order n − α, namely, 0 I n−α t u, is defined by In the formula above, Γ(z) is the Gamma function, D t represents the operator of differentiation about t, and n = [α] + 1. A CL for Equation (2) is given as the following formula: which holds for all solutions u(x, t) of Equation (2).
A formal Lagrangian for Equation (2) is written as where v(t, x) is a new dependent variable.The Euler-Lagrange operator is given by [14] δ δu where (D α t ) * is the adjoint operator of (D α t ). The adjoint equation to Equation (2) is defined by [15].
We also have where I represents the identity operator, and N t and N x represent Noether operators. The prolonged vector field X is defined by the expression below.
The Lie characteristic function W is defined by If we use the RL time-fractional derivative in Equation (2), then N t is given by ) [14].
with J given by For Equation (2), the operator N x is given by Any given generator X of Equation (2) has a corresponding invariance condition, and its solution is Consequently, the conservation laws of Equation (2) can be represented as Now, we present the conservation laws for Equation (2) using the basic definitions presented above.
On the basis of different values of the order α, let us reflect on two diverse cases.
Case 1 For the case in which α (0, 1), where i = 1, 2, and the Lie characteristic functions W i have the form Case 2 For the case in which α (1, 2), where i = 1, 2, and the Lie characteristic functions W i have the form

Concluding Remarks
In this paper, we apply the Lie symmetry analysis method to study a sixth-order generalized time-fractional Sawada-Kotera equation. First, we obtain its vector fields and symmetry reductions. Furthermore, we use the power series method to obtain explicit analytical solutions of the equation. Finally, the CLs for the original equation are computed. This paper shows that the Lie symmetry analysis method and the power series method play significant roles in the field of mathematical physics. However, we still have several problems to address: Can we generalize time fractional Sawada-Kotera equations to a higher order? Can we use Lie symmetry analysis method to solve (n + 1)-dimensional time fractional equations? How do we solve time-dependent coefficient time fractional PDEs? How can we study the space fractional equation or time-space fractional equation? These are questions that can be investigated in the future.