Markov chain approximations to scale functions of Lévy processes

Aleksandar Mijatović, Matija Vidmar, Saul Jacka

Published 2013-10-07, updated 2014-09-11Version 3

We introduce a general algorithm for the computation of the scale functions of a spectrally negative L\'evy process $X$, based on a natural weak approximation of $X$ via upwards skip-free continuous-time Markov chains with stationary independent increments. The algorithm consists of evaluating a finite linear recursion with, what are nonnegative, coefficients given explicitly in terms of the L\'evy triplet of $X$. Thus it is easy to implement and numerically stable. Our main result establishes sharp rates of convergence of this algorithm providing an explicit link between the semimartingale characteristics of $X$ and its scale functions, not unlike the one-dimensional It\^o diffusion setting, where scale functions are expressed in terms of certain integrals of the coefficients of the governing SDE.