Divergence of the multilevel Monte Carlo Euler method for nonlinear stochastic differential equations

Abstract

The Euler-Maruyama scheme is known to diverge strongly and numerically weakly when applied to nonlinear stochastic differential equations (SDEs) with superlinearly growing and globally one-sided Lipschitz continuous drift coefficients. Classical Monte Carlo simulations do, however, not suffer from this divergence behavior of Euler's method because this divergence behavior happens on rare events. Indeed, for such nonlinear SDEs the classical Monte Carlo Euler method has been shown to converge by exploiting that the Euler approximations diverge only on events whose probabilities decay to zero very rapidly. Significantly more efficient than the classical Monte Carlo Euler method is the recently introduced multilevel Monte Carlo Euler method. The main observation of this article is that this multilevel Monte Carlo Euler method does - in contrast to classical Monte Carlo methods - not converge in general in the case of such nonlinear SDEs. More precisely, we establish divergence of the multilevel Monte Carlo Euler method for a family of SDEs with superlinearly growing and globally one-sided Lipschitz continuous drift coefficients. In particular, the multilevel Monte Carlo Euler method diverges for these nonlinear SDEs on an event that is not at all rare but has probability one. As a consequence for applications, we recommend not to use the multilevel Monte Carlo Euler method for SDEs with superlinearly growing nonlinearities. Instead we propose to combine the multilevel Monte Carlo method with a slightly modified Euler method. More precisely, we show that the multilevel Monte Carlo method combined with a tamed Euler method converges for nonlinear SDEs with globally one-sided Lipschitz continuous drift coefficients and preserves its strikingly higher order convergence rate from the Lipschitz case.

Published in at http://dx.doi.org/10.1214/12-AAP890 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (5)

• Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients
• Numerical approximations of stochastic differential equations with non-globally Lipschitz continuous coefficients
• Multilevel Monte Carlo algorithms for Lévy-driven SDEs with Gaussian correction
• Statistical Romberg extrapolation: A new variance reduction method and applications to option pricing
• A Semigroup Point Of View On Splitting Schemes For Stochastic (Partial) Differential Equations

Cited by in corpus (36)

• Numerical approximations of stochastic differential equations with non-globally Lipschitz continuous coefficients
• On multilevel Picard numerical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations
• A fundamental mean-square convergence theorem for SDEs with locally Lipschitz coefficients and its applications
• On a perturbation theory and on strong convergence rates for stochastic ordinary and partial differential equations with non-globally monotone coefficients
• Exponential integrability properties of numerical approximation processes for nonlinear stochastic differential equations
• Strong convergence rates for nonlinearity-truncated Euler-type approximations of stochastic Ginzburg-Landau equations
• Strong convergence rates for explicit space-time discrete numerical approximations of stochastic Allen-Cahn equations
• Multilevel Richardson-Romberg extrapolation
• Central limit theorem for the multilevel Monte Carlo Euler method
• Strong convergence of full-discrete nonlinearity-truncated accelerated exponential Euler-type approximations for stochastic Kuramoto-Sivashinsky equations
• On strong $L^2$ convergence of time numerical schemes for the stochastic 2D Navier-Stokes equations
• Semi-Parametric Drift and Diffusion Estimation for Multiscale Diffusions
• Strong convergence rates for an explicit numerical approximation method for stochastic evolution equations with non-globally Lipschitz continuous nonlinearities
• Strong convergence for explicit space-time discrete numerical approximation methods for stochastic Burgers equations
• Convergence in Hölder norms with applications to Monte Carlo methods in infinite dimensions
• Strong and weak divergence of exponential and linear-implicit Euler approximations for stochastic partial differential equations with superlinearly growing nonlinearities
• $V$-Integrability, Asymptotic Stability And Comparison Theorem of Explicit Numerical Schemes for SDEs
• On stochastic differential equations with arbitrarily slow convergence rates for strong approximation in two space dimensions
• Lower and upper bounds for strong approximation errors for numerical approximations of stochastic heat equations
• Strong convergence of the tamed and the semi-tamed Euler schemes for stochastic differential equations with jumps under non-global Lipschitz condition
• A stochastic Gronwall inequality and applications to moments, strong completeness, strong local Lipschitz continuity, and perturbations
• Antithetic multilevel Monte Carlo method for approximations of SDEs with non-globally Lipschitz continuous coefficients
• Weak convergence of Euler-Maruyama's approximation for SDEs under integrability condition
• First order strong approximations of scalar SDEs with values in a domain
• The truncated milstein method for stochastic differential equations
• Convergence of tamed Euler schemes for a class of stochastic evolution equations
• Stability of the semi-tamed and tamed Euler schemes for stochastic differential equations with jumps under non-global Lipschitz condition
• Strongly Asymptotically Optimal Schemes for the Strong Approximation of Stochastic Differential Equations with respect to the Supremum Error
• Computational complexity analysis for Monte Carlo approximations of classically scaled population processes
• Multilevel Monte Carlo theta EM scheme for SDDEs with small noise
• Improving MLMC for SDEs with application to the Langevin equation
• Strong convergence rate of Euler-Maruyama approximations in temporal-spatial Hölder-norms
• Exponential moment bounds and strong convergence rates for tamed-truncated numerical approximations of stochastic convolutions
• Convergence rates of theta-method for neutral SDDEs under non-globally Lipschitz continuous coefficients
• Stability equivalence among stochastic differential equations and stochastic differential equations with piecewise continuous arguments and corresponding Euler-Maruyama methods
• Multi-level Monte Carlo methods with the Truncated Euler-Maruyama Scheme for Stochastic Differential Equations