Picture of a black hole

Strathclyde Open Access research that creates ripples...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of research papers by University of Strathclyde researchers, including by Strathclyde physicists involved in observing gravitational waves and black hole mergers as part of the Laser Interferometer Gravitational-Wave Observatory (LIGO) - but also other internationally significant research from the Department of Physics. Discover why Strathclyde's physics research is making ripples...

Strathprints also exposes world leading research from the Faculties of Science, Engineering, Humanities & Social Sciences, and from the Strathclyde Business School.

Discover more...

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.