The Power Fractional Calculus: First Deﬁnitions and Properties with Applications to Power Fractional Differential Equations

: Using the Laplace transform method and the convolution theorem, we introduce new and more general deﬁnitions for fractional operators with non-singular kernels, extending well-known concepts existing in the literature. The new operators are based on a generalization of the Mittag– Lefﬂer function, characterized by the presence of a key parameter p . This power parameter p is important to enable researchers to choose an adequate notion of the derivative that properly represents the reality under study, to provide good mathematical models, and to predict future dynamic behaviors. The fundamental properties of the new operators are investigated and rigorously proved. As an application, we solve a Caputo and a Riemann–Liouville fractional differential equation.


Introduction
Fractional calculus theory plays a crucial role in bridging the gap on the modeling of many neglected phenomena with memory effects. Unlike Markov-chain processes, where the current value of the function under consideration depends only on that of the recent past, long-range memory is naturally included under fractional modeling [1,2].
An in-depth examination of the literature of fractional calculus confirms that the modeling of memory effects has undergone several transformations in recent years, namely by considering the exponential effect under the Caputo-Fabrizio derivative [3], the Mittag-Leffler effect with Atangana-Baleanu and Al-Refai operators [4,5], and the new generalized fractional operator of Hattaf [6]. Here, we propose new, and more general, fractional operators based on a generalized Mittag-Leffler function, which we call the "power Mittag-Leffler function". Our new mathematical concept allows us to unify and extend the fractional literature by developing a family of power fractional operators (PFOs) that expand the existing generalized fractional operators and their many consequences [3][4][5][6]. Broadly speaking, the exponential function is converted to the expanded power function, and the generalized Mittag-Leffler function is transformed into the power Mittag-Leffler counterpart that we propose here.
Advanced mathematical results have recently been proved in the framework of fractional calculus: see, e.g., [7][8][9][10][11] and the references therein. However, to effectively describe realistic phenomena, all available definitions suffer from some limitations, depending on the application at hand, which has motivated us to propose here new, more general, notions, containing the key power parameter p. The currently introduced power fractional calculus enables the generalization and unification of many of the cited results, allowing engineers, researchers, and scientists to select the appropriate fractional derivative with respect to the phenomenon under study in a natural way via the presence of the parameter p in our new definitions.
The action of the parameter p on a system is illustrated in the numerical simulation phase, where it is essential to find the appropriate value of p to describe real data with the adopted model, to describe the current trajectories to correctly predict the asymptotic behavior in the future: see our section devoted to the resolution of some power fractional differential equations (PFDEs). Furthermore, the defined power fractional derivative derives its legitimacy from the construction of its inverse power fractional integral operator (PFIO), using the Laplace transform and the convolution theorem. Finally, we claim that our PFOs have considerable potential, both for the development of mathematical modeling, in various fields, and in the mathematics discipline itself. All these reasons support the originality, importance, relevance and robustness of our definitions and results.
The paper is organized as follows. Section 2 is devoted to the introduction of the new power Mittag-Leffler function (Definition 1) accompanied with its convergence (Theorem 1). Section 3 contains novel definitions of the PFOs in both Caputo (Definition 2) and Riemann-Liouville senses (Definition 3), as well as establishing the connection between them (Theorem 3). Section 4 is dedicated to the discovery of the corresponding PFIO (Definition 4). To show the significance and usefulness of our PFOs, the resolution of two PFDEs is performed in Section 5. Section 6 concludes the paper and highlights some directions for future research.

The Power Mittag-Leffler Function
In this section, we introduce a new generalization of the Mittag-Leffler function, which we call the power Mittag-Leffler function. , where p ∈ R * + , and min α, β > 0.
It follows, from Cauchy's criterion, that as n −→ ∞ when α > 0, which leads to the absolute convergence for all values of z ∈ C with the radius of convergence of the power series being infinite.

The Power Fractional Derivatives
In this section, we present a new fractional derivative. Along the text, , which is a Hilbert space. In addition, we adopt the following notations: In applications, the choice of a suitable and concrete normalization function N may depend on the phenomenon under study. Along the paper, we denote Definition 2 (The power fractional derivative of order α in the Caputo sense). Let 0 ≤ α < 1 and min(β, p) > 0. The power fractional derivative of order α in the Caputo sense, of a function f ∈ H 1 (a, b) with respect to the weight function w(t), is defined as where We note that "p C " in the operator p C D α,β,p a,t,w stands for "power Caputo".

Remark 2.
Our power fractional derivative in the Caputo sense given by Definition 2 generalizes many existing notions found in the literature: if w(t) ≡ 1, p = e, and β = 1, then we obtain the Caputo-Fabrizio fractional derivative [3] given by if w(t) ≡ 1, p = e, and β = α, then we get the Atangana-Baleanu fractional derivative [4] given by 3. if p = e and β = α, then we obtain the weighted Atangana-Baleanu fractional derivative defined by Al-Refai in [5], given by if p = e, then we obtain the weighted generalized fractional derivative introduced by Hattaf [6], which is given by Remark 3. It is worth observing that the power fractional derivative in the Caputo sense satisfies the following two properties: and p C D α,β,p a,t,1 f (t) = 0 for any constant function f (t).
Definition 3 (The power fractional derivative of order α in the Riemann-Liouville sense). Let 0 ≤ α < 1 and min(p, β) > 0. The power fractional derivative of order α in the Riemann-Liouville sense, of a function f ∈ H 1 (a, b) with respect to the weight function w(t), is defined as where Remark 4. The statements of Remark 2 are also verified in the Riemann-Liouville sense.

Remark 5.
The following property of the power fractional derivative in the Riemann-Liouville sense is satisfied: Theorem 2. The power fractional derivatives in the Caputo and Riemann-Liouville senses are linear operators.
Proof. We easily see that and for all scalars c 1 and c 2 and all functions f , g ∈ H 1 (a, b).

Theorem 3.
Let w f be an analytic function. Then, Proof. Because of the analyticity of the function w f , we get n!Γ(βk + 1)(βk + n + 1) The proof is complete.

The Power Fractional Integral
With the intention of establishing the associated power fractional integral, we begin this section by computing the Laplace transform of the power fractional derivatives in Caputo and Riemann-Liouville senses multiplied by w(t).
The following equalities hold: Proof. We begin by proving the first statement of Lemma 1: To prove the second statement, we get: The result is proved.
Theorem 4. The fractional differential equation has a unique solution given by where RL I β 0,w is the standard weighted Riemann-Liouville fractional integral of order β given by Proof. The equality (13) is equivalent to Using Lemma 1, we conclude that Moreover, the action of the inverse Laplace transform yields which completes the proof.
Our Theorem 4 allows us to define an appropriate inverse operator for power fractional differentiation.
Definition 4 (The power fractional integral of order α). Let 0 ≤ α < 1 and min(p, β) > 0. The power fractional integral of order α, of a function f ∈ H 1 (a, b) with respect to the weight function w(t), is defined by where

Examples of Power Fractional Differential Equations
In this section, we treat two examples of power fractional differential equations (PFDEs). Our first example considers a non-autonomous PFDE in the Riemann-Liouville sense. Example 1. Consider the following non-autonomous PFDE on [0, 100]: Using Theorem 4, we obtain that The action of the parameter p on the obtained solution is shown in Figure 1. We now consider an autonomous PFDE in the Caputo sense.

Example 2.
Consider the following autonomous PFDE: The action of the Laplace transform on both sides of Equation (20) yields: Using Lemma 1, we obtain that The effect of the inverse Laplace transform operator yields Applying the integration by parts formula on d dt , we can state that the solution to problem (20) is given by t β * w (t).

Conclusions
In this paper, some new mathematical concepts, enabling the introduction of a new extended fractional calculus, are provided. The new approach allows choice of the most appropriate notion of differentiation to suitably describe real dynamic phenomena under study, describing the observed trajectories and correctly predicting future behaviors. More precisely, we introduce the "power Mittag-Leffler function" p E α,β (·) that extends several important functions: the Mittag-Leffler function e E α,1 (·), first introduced by Mittag-Leffler in [12]; the function e E α,β (·) of Wiman [13]; and those introduced by Prabhakar [14] and Salim [16]. With the help of the new power Mittag-Leffler function, we then introduce the new power fractional derivatives p C D α,β,p a,t,w (·) and p RL D α,β,p a,t,w (·), which generalize those available in the literature, namely the Caputo-Fabrizio [3], Atangana-Baleanu [4], weighted Atangana-Baleanu [5], and weighted generalized fractional derivatives [6]. Moreover, an appropriate power fractional integral operator (PFIO) p I α,β,p a,t,w (·) is introduced, which is an important tool for the solution of power fractional differential equations (PFDEs). As examples, we investigated two PFDEs. The first is a non-autonomous PFDE: using our PFIO, we compute its solution and illustrate, numerically, the impact of the parameter p on the solution. The second example considered is an autonomous PFDE and its solution is obtained using the Laplace transform operator.
Here, we have only introduced the power fractional calculus and provided the most fundamental results with some applications to power fractional differential equations. In future work, several investigations may be designed to develop the new fractional calculus, which will enable the setting up of numerous applications on many parallel domains, e.g., in the fractional neural networks framework, analogously to what is done in [20][21][22].