In this article, we first provide a survey of the exponential option pricing models and show that in the framework of the risk-neutral approach, they are governed by the space-fractional diffusion equation. Then, we introduce a more general class of models based on the space-time-fractional diffusion equation and recall some recent results in this field concerning the European option pricing and the risk-neutral parameter. We proceed with an extension of these results to the class of exotic options. In particular, we show that the call and put prices can be expressed in the form of simple power series in terms of the log-forward moneyness and the risk-neutral parameter. Finally, we provide the closed-form formulas for the first and second order risk sensitivities and study the dependencies of the portfolio hedging and profit-and-loss calculations upon the model parameters.
IntroductionFractional Calculus (FC) is nearly as old as conventional calculus. Many prominent mathematicians including Leibniz, Fourier, Laplace, Liouville, Riemann, Weyl, and Riesz suggested their own definitions of the fractional integrals and derivatives and studied their properties. Whereas the mathematical theory of FC was nearly completed a long time ago (see, e.g., [1] or the recent reference books [2,3]), only a few applications of FC outside mathematics were known until recently. During the last three decades, the situation changed completely, and currently, the majority of FC publications is devoted to modeling of a broad class of systems and processes using either the FC operators or the so-called fractional ordinary or partial differential equations. In particular, the FC models were successfully employed in physics [4,5], control theory [6], as well as in engineering, life, and social sciences [7,8], to mention only a few of the many application areas.The two probably most prominent and broadly-recognized FC applications are in linear viscoelasticity [9] and for describing anomalous transport processes [10]. In both cases, the FC models in the form of the fractional ODEs (linear viscoelasticity) and the fractional PDEs (anomalous transport processes), respectively, cannot be derived from "first principles". Instead, they are introduced as interpolations between several models formulated in terms of the ODEs or PDEs, respectively. In the case of linear viscoelasticity, the basic FC model interpolates between the Hooke model (elasticity, ODE Mathematics 2019, 7, 796 2 of 23 of the zero order) and the Newton model (viscosity, ODE of the first order). The FC model for the slow anomalous diffusion interpolates between a stationary state (time-independent Poisson equation) and the diffusion equation (first order equation with respect to time), whereas the fast diffusion is modeled by the fractional diffusion-wave PDE that interpolates between the diffusion equation (first order equation with respect to time) and the wave equation (second order equation with respect to time).Of course, this kind of model needs an additional justification...