Error analysis of finite difference and Markov chain approximations for option pricing

Mijatović and Pistorius proposed an efficient Markov chain approximation method for pricing European and barrier options in general one‐dimensional Markovian models. However, sharp convergence rates of this method for realistic financial payoffs, which are nonsmooth, are rarely available. In this paper, we solve this problem for general one‐dimensional diffusion models, which play a fundamental role in financial applications. For such models, the Markov chain approximation method is equivalent to the method of lines using the central difference. Our analysis is based on the spectral representation of the exact solution and the approximate solution. By establishing the convergence rate for the eigenvalues and the eigenfunctions, we obtain sharp convergence rates for the transition density and the price of options with nonsmooth payoffs. In particular, we show that for call‐/put‐type payoffs, convergence is second order, while for digital‐type payoffs, convergence is generally only first order. Furthermore, we provide theoretical justification for two well‐known smoothing techniques that can restore second‐order convergence for digital‐type payoffs and explain oscillations observed in the convergence for options with nonsmooth payoffs. As an extension, we also establish sharp convergence rates for European options for a rich class of Markovian jump models constructed from diffusions via subordination. The theoretical estimates are confirmed using numerical examples.