Abstract:In the field of computational finance, one is commonly interested in the expected value of a financial derivative whose payoff depends on the solution of stochastic differential equations (SDEs). For multi-dimensional SDEs with non-commutative diffusion coefficients in the globally Lipschitz setting, a kind of one-half order truncated Milstein-type scheme without Lévy areas was recently introduced by Giles and Szpruch (2014), which combined with the antithetic multilevel Monte Carlo (MLMC) gives the optimal overall computational cost O(ε−2) for the required target accuracy ε . Nevertheless, many nonlinear SDEs in applications have non-globally Lipschitz continuous coefficients and the corresponding theoretical guarantees for antithetic MLMC are absent in the literature. In the present work, we aim to fill the gap and analyze antithetic MLMC in a non-globally Lipschitz setting. First, we propose a family of modified Milstein-type schemes without Lévy areas to approximate SDEs with non-globally Lipschitz continuous coefficients. The expected one-half order of strong convergence is recovered in a non-globally Lipschitz setting, where even the diffusion coefficients are allowed to grow superlinearly. This then helps us to analyze the relevant variance of the multilevel estimator and the optimal computational cost is finally achieved for the antithetic MLMC. Since getting rid of the Lévy areas destroys the martingale properties of the scheme, the analysis of both the convergence rate and the desired variance becomes highly non-trivial in the non-globally Lipschitz setting. By introducing an auxiliary approximation process, we develop non-standard arguments to overcome the essential difficulties. Numerical experiments are provided to confirm the theoretical findings.

Multilevel Monte Carlo method for jump-diffusion SDEs

Variable-Step Euler–Maruyama Approximations of Regime-Switching Jump Diffusion Processes

Itô-Taylor Expansion Method of European Spread Option Pricing for Multivariate Diffusions with Jumps

Three Improvements to Multi-Level Monte Carlo Simulation of SDE Systems

A Multi-level Monte Carlo simulation for invariant distribution of Markovian switching Lévy-driven SDEs with super-linearly growth coefficients

Multilevel Monte Carlo for a class of Partially Observed Processes in Neuroscience

Multilevel Monte Carlo Path Simulation

Computing Functionals of Multidimensional Diffusions via Monte Carlo Methods

Local M-Estimation for Jump-Diffusion Processes

Adaptive Weak Approximation of Diffusions with Jumps

Numerical Solutions of Jump Diffusions with Markovian Switching

Adaptive Mesh Construction for the Numerical Solution of Stochastic Differential Equations with Markovian Switching

A micro/macro algorithm to accelerate Monte Carlo simulation of stochastic differential equations

An Antithetic Multilevel Monte Carlo-Milstein Scheme for Stochastic Partial Differential Equations

Strong convergence of a class of adaptive numerical methods for SDEs with jumps

Multilevel Monte Carlo EM scheme for MV-SDEs with small noise

An Adaptive Random Bit Multilevel Algorithm for SDEs

Compensated projected Euler-Maruyama method for stochastic differential equations with superlinear jumps

Jump locations of jump-diffusion processes with state-dependent rates

Antithetic multilevel Monte Carlo method for approximations of SDEs with non-globally Lipschitz continuous coefficients