Option pricing in subdiffusive Bachelier model

The earliest model of stock prices based on Brownian diffusion is the Bachelier model. In this paper we propose an extension of the Bachelier model, which reflects the subdiffusive nature of the underlying asset dynamics. The subdiffusive property is manifested by the random (infinitely divisible) periods of time, during which the asset price does not change. We introduce a subdiffusive arithmetic Brownian motion as a model of stock prices with such characteristics. The structure of this process agrees with two-stage scenario underlying the anomalous diffusion mechanism, in which trapping random events are superimposed on the Langevin dynamics.We find the corresponding fractional Fokker-Planck equation governing the probability density function of the introduced process. We construct the corresponding martingale measure and show that the model is incomplete. We derive the formulas for European put and call option prices. We describe explicit algorithms and present some Monte-Carlo simulations for the particular cases of alpha-stable and tempered alpha-stable distributions of waiting times.