Antithetic Multilevel Monte Carlo Estimation for Multidimensional SDEs

Abstract

In this paper we develop antithetic multilevel Monte Carlo (MLMC) estimators for multidimensional SDEs driven by Brownian motion. Giles has previously shown that if we combine a numerical approximation with strong order of convergence O(Δt) with MLMC we can reduce the computational complexity to estimate expected values of Lipschitz functionals of SDE solutions with a root-mean-square error of ∈ from O(∈-3)/ to O(∈-2). However, in general, to obtain a rate of strong convergence higher thnan O(Δt1/2) requires simulation, or approximation, of Lévy areas. Recently, Giles and Szpruch [5] constructed an antithetic multilevel estimator thnnat avoids thnne simulation of Lévy areas and still achieves an MLMC correction variance which is O(Δt2) for smooth payoffs and almost O(Δt3/2) for piecewise smooth payoffs, even though there is only O(&Deltat1/2) strong convergence. This results in an O(∈-2) complexity for estimating the value of financial European and Asian put and call options. In this paper, we extend these results to more complex payoffs based on the path minimum. To achieve this, an approximation of the Lévy areas is needed, resulting in O(Δt3/4) strong convergence. By modifying the antithetic MLMC estimator we are able to obtain O(∈-2 log(∈)2) complexity for estimating financial barrier and lookback options. © Springer-Verlag Berlin Heidelberg 2013.