Fractional Partial Differential Equations Associated with L ˆ e vy Stable Process

: In this study, we ﬁrst present a time-fractional Lˆ e vy diffusion equation of the exponential option pricing models of European option pricing and the risk-neutral parameter. Then, we modify a particular Lˆ e vy-time fractional diffusion equation of European-style options. Further, we introduce a more general model based on the Lˆ e vy-time fractional diffusion equation and review some recent ﬁndings associated with risk-neutral free European option pricing.


Introduction
One of the significant problems in finance is to derive value from financially traded assets that is also known as the pricing of financial instruments, for example, stocks, and it is a very interesting problem. In the literature, Merton (1990, [1]) was among the first researchers who gave the systematic solution for this problem, and proposed the Black-Scholes (BS) model where the model rests on the assumption that the natural logarithm of the stock price S t defined as follows: where µ > 0 is the average compounded growth rate of the stock S t , and dB t is the increment of Brownian motion which assumed to have the Normal or Gaussian distribution, and σ ≥ 0 represents the volatility of the returns from holding S t . The Equation (1) is also known as componential equation. On either side a Lêvy process is a stochastic process with independent, stationary increments that represents the motion of a point whose successive displacements are random and independent, and statistically identical over different time intervals of the same length. The mathematical theory of Lêvy process can be found in Bertoin (1996) [2] or Sato (1999) [3]. An example of a Lêvy process that is well-known from, for instance, the Black-Scholes-Merton option pricing theory is the Brownian motion (or Wiener process), where the increments are normally distributed.
Thus if we substitute the Lêvy process by the Brownian motion in componential Equation (1), the pricing partial differential equation becomes a partial integro-differential equation; further details related to the partial integro-differential equations (PIDEs) can be found in [4,5]. The partial integro-differential equations are also studied in order to understand the non-locality phenomena which are produced by the jumps in the Lêvy process.
One of the methods to solve PIDEs was the numerical method which was proposed by (Cont,[6]) and that was the finite-difference method for option pricing, having jump-diffusion as well as exponential Lêvy process models, see Lewis [7]. On the other hand, the second method was the fast Fourier transform of European-style options (see [8]). Further similar strategies were also proposed, for example, ref. ([9]) proposed a model to use fractional calculus.
In this article, we modify European-style options under a risk-neutral probability condition for the stock-price assets, followed by liquidity market in the financial literature. We also consider to generate some partial integro-differential equation for possible application to less-studied issues such as barrier options for finite moment having log-stable (FMLS) processes in the future.
The article is based on the following: Section 2 reviews the basic concepts of Lêvy operations and applications in financial modeling. Section 3 introduces the concepts of fractional calculus and how to solve fractional differential equations, and reviews the main concepts of Lêvy process. Section 4 introduce the main result. Finally, Section 5 will conclude and discuss some applications.

Fractional Diffusion Model and Option Pricing
In a fully liquid market, regardless of the trading size, the options trader cannot influence the price of the underlying asset in the trading of the asset in order to duplicate the option. In the literature (Chen, et al. (2014), [10]) studied this model, where L α,β t is a Lêvy α-stable process with skew parameter β. Before viewing the idea of that research we will define α-stable distribution, that is, the distribution is said to be stable if the location and scale parameters have the same distributions of linear combination of two independent random variables with respect to this distribution. Similarly, a random variable is said to be stable if its distribution is stable. The stable dissemination family at times indicates as the Lêvy alpha-stable distribution.

Definition 1.
Any random variable X is s-stable if for each n ∈ N with X 1 , X 2 , . . . , X n being infinitely divisible copies of X X 1 + X 2 + . . . + X n = bX + c or some constants b = b(n) > 0 and c = c(n) ∈ R d . It is called strictly stable for any n ∈ N if c(n) = 0 .
For an infinitely divisible random vector X * t define the alpha-stable as follows.

Definition 2.
A stable X is called alpha-stable, whenever X * t = t 1 α X + c or some constants c = c(r) ∈ R d , t > 0, and 0 < α ≤ 2. When c(t) = 0, for t > 0, then X is called strictly alpha-stable . Now, consider the following dynamic under a risk neutral probability measure for the stock price S t for time 0 < t < T, where index α of stability satisfies 1 < α < 2, and volatility σ > 0. When σ = 0, we will get the original BS model. Moreover where r and q respectively denote deterministic parameters corresponding to the risk-free rate and dividend yield. We restrict our selves to the case where β = −1 to obtain finite moments of S t and negative skewness in the return density distribution. In particular for n > 0, then The model in the Equation (2) is known as Finite Moment Log Stable (FMLS) for short model. Under the risk-neutral measure the log price satisfies the following SDE: where v = − 1 2 σ α sec πα 2 represents the convexity adjustment. Let u(t, x) be the price of the European call option with x = x t := ln (S t ). (Chen et al. (2014), [10]) In order to find FPDE let u(t, x) satisfy under FMLS the following fractional PDEs where K is the strike price and and Γ(.) is the gamma function defined by:

The Model
In this research, we incorporate L α,β t as a Lêvy α-stable process with skew parameter β. Consider the following dynamic under a risk-neutral probability measure for the stock price S t ; the goal is to consider a modified model to Equation (2) that consists on an illiquid market with impact additional term that for 0 ≤ t < T, 0 < γ < 1 and 1 < α ≤ 2,with boundary condition where λ(t, S t ) ≥ 0 is the price impact function of the trader and β t denotes the number of shares that the trader has in the stock at time t. The term λ(t, S t )dβ t represents the price impact of the investor's trading is additional term of Chen model (2) and similarly, the Caputo fractional partial derivative of u defined by the expression ∂y .
In this work, we consider trading strategies written in the following form for some processes (η t ) t≥0 and (ζ t ) t≥0 to be determined endogenously and β 0 is the initial number of shares in the stock. Next we consider the wealth process (V t ) t≥0 corresponds to a self-financing strategy (θ t , β t ) t≥0 for the trader and given by To find the fractional partial differential equations that satisfy our model in Equation (4) we need a method to solve fractional equation. In this way we follow the (Demirci and Ozalp, 2012) [11] as an example where they solved the fractional differential equation for the initial value problem in the sense of Caputo type FDE given by D γ x(t) = f (t, x(t)), x(0) = x 0 which has a solution

) is a solution for an equation having integer order differentials.
We also need the method of the from literature (Jumarie, [12]) of the equation, where 0 < α ≤ 1, has a solution defined by the equality Furthermore, we will use in our model the following formula The general Fourier transform is defined bŷ (8) and the inverse Fourier transformation is defined by

Lêvy Process
The distribution of Lêvy process is characterized by the Lêvy-Khintchine formula and it is considered as a modified model that characterize the Lêvy process in a very compact way. Today it is known as the Lêvy-Khintchine process. More definitely, a time-dependent random variable X t is a Lêvy process if and only if it has independent and stationary increments having log-characteristic function given by the Lêvy-Khintchine theorem: Theorem 1 (Lêvy-Khintchine presentation theorem). Let (X t ) t≥0 Lêvy process on R with characteristic triplet (m, σ, w), then E[e izX t ] = e tΨ(z) , t ∈ R, with characteristic exponent of the Lêvy process where, To accommodate how the Lêvy processes being incorporated in the derivatives pricing models, we recall the standard Black-Scholes framework and see how it was built by Gaussian shocks. To find the fair or arbitrage-free prices of a financial instrument whose value are deriven from the underlying share price S t , it is also necessary to express the dynamics of S t under what is known as a neutral risk measure or the equivalent martingale scale. In the price, the European option may be expressed as the neutral condition for a risk as Fourier transform of European option can be written as (Du,[13]) where and the indicator function of set A where I A : A ⊂ X → {0, 1} and defined by

Lêvy Stable Processes
Let w(x) = w LS (x) be Lêvy density function and given by where D > 0, q + p = 1 and α ∈ (0, 2). Then by using the Equation (9) we obtain the characteristic exponent of an LS process in the parameters as follows: σ, α, β and m: An equivalent form can be written as where β = p − q. If β = −1, then p = 0 and q = 1, that is (Alvaro Cartea et al. [9])

Main Results
Consider the fractional differential Lêvy equation That can be rewritten as λ t = λ(t, S(t)) and let x t = ln S t Next, we derive revised and updated FPDEs for options which are written on assets and follow the Lêvy operations that were mentioned in the previous section. In order to find the relation between the fractional price equations and LP process, [14], we make use of Fourier transform, as in of the value of European style option priceÛ (ξ, T), written on S t , and satisfies Let ζ T (v) denote the Fourier transform of the time value, where Let U (ξ, T) denote the Fourier transform of a European-style option and defined by with boundary condition U (ξ, T) = Π(ξ, T). Now substitute the Equation (15) in Equation (19) and taking the inverse Fourier transform we reach then taking the inverse Fourier transform delivered to To prove Equation (17) Take the integral for the above equation and using method (6) we get By the same way and using method of (Demirci and Ozalp (2012) To prove Equation (17) That isÛ is a solution of the Equation (19).

Conclusions
In this paper, we modified the particular Lêvy-time fractional diffusion equation, and applied to the price of fractional financial derivatives of European-style options. A more general class of model based on the fractional diffusion equation as Lêvy process was also presented in the form of fractional partial differential equation (FPDE) then the solution was obtained and applied to the European option pricing having risk-free parameters.