#### Abstract

The main purpose of this paper is to apply the Lie symmetry analysis method for the two-dimensional time fractional Fokker-Planck (FP) equation in the sense of Riemann–Liouville fractional derivative. The Lie point symmetries are derived to obtain the similarity reductions and explicit solutions of the governing equation. By using the new conservation theorem, the new conserved vectors for the two-dimensional time fractional Fokker-Planck equation have been constructed with a detailed derivation. Finally, we obtain its explicit analytic solutions with the aid of the power series expansion method.

#### 1. Introduction

Fractional calculus has attracted more attention of many researches in various scientific areas including biology, physics, financial theory, gas dynamics, engineering, fluid mechanics, and other areas of science, see for example [1–8]. The theory of fractional calculus is considered as a generalization of classical differential and integral calculus; it is an excellent tool for describing the memory effect and hereditary properties of various processes and viscoelastic materials. Due to its realistic senses, many researchers have tried to look for exact, analytical, and numerical solutions of fractional partial differential equations using different powerful methods such us -expansion method [9], the Variational iteration method [10, 11], functional variable method [12], subequation method [13], Finite difference method [14], Exp function method [15], Homotopy analysis method [16], Adomian decomposition method [17], the First integral method [18], Laplace transform method [19], Sumudu transform method [20], and so many other approaches.

The Lie symmetry method was firstly advocated by the Norwegian mathematician Sophus Lie [21, 22], who has made great achievements in the theories of continuous groups and differential equations. It is an efficient approach and widely employed for solving ordinary differential equations (ODEs), partial differential equations (PDEs), and fractional partial differential equations (FPDEs). This popularity is due to its utility in determining the explicit solutions of both ODEs and PDEs, linearization of some nonlinear equations, reducing the order of independent variables, and so on. Many papers focused on constructing symmetries of different fractional differential equations [23, 24–27]. Furthermore, the concept of conservation laws is fundamental and widely used in the study of the resolution of PDEs. Moreover, they convey a large deal of information about the studied physical system. The new conservation laws were introduced by Ibragimov [28], based on the notion of Lie symmetry generators without Lagrangian for solving FPDEs. Therefore, this new conservation law plays an increasingly important role in solving the conservation laws of FPDEs. More details about conservation laws can be found in [29–32].

In this paper, we consider the following two-dimensional time fractional Fokker-Planck equation It is well known that the Markovian diffusion process can be described with the Fokker-Planck equation. The Fokker-Planck equation is a partial differential equation for the probability density and the transition probability of these stochastic processes. It plays an important role in control theory, fluid mechanics, astrophysics, and quantum [33]. Moreover, it has been applied in various natural science fields such as quantum optics, solid-state physics, chemical physics, theoretical biology, and circuit theory. It is firstly proposed by Fokker and Planck to characterize the Brownian motion of particles [34]. Many researchers have solved the Fokker-Planck equation using various powerful methods, for more details see [35–37].

The fractional derivatives described here are in the Riemann-Liouville sense of order , see [38, 39], which is defined by where is the Gamma function defined by The main motivation behind this article is to make use of the Lie symmetry method to get the infinitesimal generators, group invariant solutions for the time fractional two-dimensional Fokker-Planck equation (FP), and to construct conservation laws given by Ibragimov [28]. Therefore, the new conserved vectors have been obtained using the new conservation theorem. Based on the power series method [32], the explicit power series solutions of the two-dimensional time fractional Fokker-Planck equation are derived.

The rest of this paper is organized as follows: in Section 2, we review some basic definitions of the Lie Symmetry method for fractional partial differential equations (FPDEs) and its properties. By employing the proposed method, Lie point symmetries of the Eq. (1) are obtained; by using similarity variables, the reduced equations are obtained; solving some of them, then the similarity solutions of Eq. (1) are deduced in Section 3. In Section 4, the conservation laws of Eq. (1) are obtained. Section 5 is devoted to constructing the explicit analytical power series solutions. Some conclusions and discussions are given in Section 6.

#### 2. Method of Lie Symmetry Analysis for FPDEs

In this section, we briefly review the main points about Lie symmetry analysis of FPDEs [39–41] of the following form We assume that the Eq. (4) is invariant under a one-parameter Lie group of infinitesimal transformations which are given as where is a group parameter and , , , and are infinitesimals and is extended infinitesimal. The explicit expressions of , , , are given by

, , and are the total derivatives with respect to , , and , respectively, which are defined as where , , and so on.

The infinitesimal generator is given by the following expression The infinitesimal generator satisfy the following invariance condition of Eq. (4): where The structure of the Riemann-Liouville derivative must be invariant under transformations (5), because the lower limit of the integral (2) is fixed. The invariance condition yields

The explicit form of the extended infinitesimal can be obtained as follows: where It is worth noting that if the infinitesimal is linear in , due to the presence of .

*Definition 1 (see [42]). *The function is an invariant solution of Eq. (4) if and only if
(i) is an invariant surface, that is to say(ii) satisfies Eq. (4)

#### 3. Symmetry Analysis and Similarity Reductions of the Two-Dimensional Time Fractional Fokker-Planck Equation

In the present section, the Lie symmetry analysis method has been applied for deriving the infinitesimal generators of the two-dimensional time fractional Fokker-Planck (1). By using the third prolongation [43, 44], the symmetry determining equation for Eq. (1) has been obtained as

By substituting the expressions given in Eq. (6) and Eq. (12) into Eq. (15), and equating various powers of derivatives of to zero, we obtain an overdetermined system of linear equations; by solving this system, we obtain the following infinitesimals

where , are arbitrary constants. So, the associated vector fields of Eq. (1) are given by

*Case 2. *For , the characteristic equation is

By solving the above characteristic equation, we obtain the solution . Substituting it into Eq. (1), we derive the following reduced fractional ordinary Fokker-Planck equation:

Using the symmetry method, we obtain the following infinitesimals: where , are arbitrary constants. Then, the Lie algebra of infinitesimal symmetries of Eq. (19) is given by

*Case 3. *The characteristic equation for the infinitesimal generator can be expressed symbolically as follows:

By solving the above characteristic equation, we obtain the solution . Substituting it into Eq. (19), we derive the reduced fractional ordinary equation: The above can be solved through the Laplace transform method Since the Laplace transform of the Riemann–Liouville derivative is defined by the following form then, According to Eq. (25), we have By using the inverse Laplace transform, it gives where is the Mittag-Leffler function.

*Case 4. *For , the similarity transformation corresponding to this generator can be derived by solving the associated characteristic equation
which take the form
replacing it in Eq. (19) yields the following reduced FODE:
By using the Laplace transform, we obtain the following solution

#### 4. Conservation Laws of the Two-Dimensional Time Fractional Fokker-Planck Equation

In this section, the conservation laws of the two-dimensional time fractional Fokker-Planck equation have been investigated by using a new conservation theorem [28]. The conserved vectors have been obtained, and it satisfies the following conservation equation:

The formal Lagrangian for Eq. (1) can be written as follows:

where is the new dependent variable. Based on the definition of the Lagrangian, the action integral of Eq. (35) is given by

The Euler–Lagrange operator is defined as The adjoint operator of is defined by where is the right-sided operator of fractional integration of order that is defined by So, the adjoint equation of Eq. (1) as the Euler–Lagrange equation, given by

For the case of three independent variables and one dependent variable , we get

where is the identity operator and is denoted as the Euler-Lagrange operator. So is presented as and the Lie characteristic function is given by

Using the Lie symmetries , we have

Based on the fractional generalizations of the Noether operators, the components of conserved vectors can be presented as follows:

where is defined by And the other components are defined as where , , . Using Eqs. (45) and (47), we obtain the following components of conserved vectors.

*Case 5. *For , we have

*Case 7. *For , we have

*Case 8. *For , we have

*Case 9. *For , we have

*Case 10. *For , we have

*Case 11. *For , we have

#### 5. Power Series and Analytical Solutions for Eq. (1)

In this section, based on the power series method [45], the exact analytic solutions are a kind of exact power series solutions for Eq. (1), constructed with a detailed derivation. where and are arbitrary. The time fractional Fokker-Planck Eq. (1) is reduced to the following ODE We assume that the solution of Eq. (1) has the following form: where are constants to be determined later. According to Eq. (57), we get

Substituting (57) and (58) into (56), we obtain Observing coefficients in Eq. (59), when , we have By comparing coefficients of , we get When , we have

The power series solution for Eq. (5) can be rewritten as follows:

#### 6. Conclusion

In this paper, the invariance properties of the two-dimensional time fractional Fokker-Planck equation with the Riemann-Liouville fractional derivative have been investigated in the sense of Lie point symmetries. Then, the power series method has been applied to get an explicit solution for the two-dimensional time fractional Fokker-Planck equation. For obtaining new components of conserved vectors, a new theorem of conservation law has been employed along with the formal Lagrangian, which allows us to construct conservation laws for the two-dimensional time fractional Fokker-Planck equation. Our results show that the extended Lie group analysis approach and the power series method provide powerful mathematical tools to investigate other FDEs in different fields of applied mathematics. In addition, it shows that the proposed analysis is very efficient to construct conservation laws of the two-dimensional time fractional Fokker-Planck equation. Moreover, we can employ symmetry analysis to the time-space fractional Fokker-Planck equation; it will be valuable as future subject works.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.