Picture of a sphere with binary code

Making Strathclyde research discoverable to the world...

The Strathprints institutional repository is a digital archive of University of Strathclyde research outputs. It exposes Strathclyde's world leading Open Access research to many of the world's leading resource discovery tools, and from there onto the screens of researchers around the world.

Explore Strathclyde Open Access research content

Convergence of Monte Carlo simulations involving the mean-reverting square root process

Higham, D.J. and Mao, X. (2005) Convergence of Monte Carlo simulations involving the mean-reverting square root process. Journal of Computational Finance, 8 (3). pp. 35-61. ISSN 1460-1559

[img]
Preview
PDF (strathprints000160.pdf)
strathprints000160.pdf

Download (269kB) | Preview

Abstract

The mean-reverting square root process is a stochastic differential equation (SDE) that has found considerable use as a model for volatility, interest rate, and other financial quantities. The equation has no general, explicit solution, although its transition density can be characterized. For valuing path-dependent options under this model, it is typically quicker and simpler to simulate the SDE directly than to compute with the exact transition density. Because the diffusion coefficient does not satisfy a global Lipschitz condition, there is currently a lack of theory to justify such simulations. We begin by showing that a natural Euler-Maruyama discretization provides qualitatively correct approximations to the first and second moments. We then derive explicitly computable bounds on the strong (pathwise) error over finite time intervals. These bounds imply strong convergence in the limit of the timestep tending to zero. The strong convergence result can be used to justify the method within Monte Carlo simulations that compute the expected payoff of financial products. We spell this out for a bond with interest rate given by the mean-reverting square root process, and for an up-and-out barrier option with asset price governed by the mean-reverting square root process. We also prove convergence for European and up-and-out barrier options under Heston's stochastic volatility model - here the mean-reverting square root process feeds into the asset price dynamics as the squared volatility.