Simulation of McKean Vlasov SDEs with super linear growth

Abstract

We present two fully probabilistic Euler schemes, one explicit and one implicit, for the simulation of McKean-Vlasov Stochastic Differential Equations (MV-SDEs) with drifts of super-linear growth and random initial condition. We provide a pathwise propagation of chaos result and show strong convergence for both schemes on the consequent particle system. The explicit scheme attains the standard $1/2$ rate in stepsize. From a technical point of view, we successfully use stopping times to prove the convergence of the implicit method although we avoid them altogether for the explicit one. The combination of particle interactions and random initial condition makes the proofs technically more involved. Numerical tests recover the theoretical convergence rates and illustrate a computational complexity advantage of the explicit over the implicit scheme. A comparative analysis is carried out on a stylized non-Lipschitz MV-SDE and a mean-field model for FitzHugh-Nagumo neurons. We provide numerical tests illustrating "particle corruption" effect where one single particle diverging can "corrupt" the whole particle system. Moreover, the more particles in the system the more likely this divergence is to occur.

arXiv comment43 pages, 4 figures, 1 Table; Final Author version (to appear in IMA Journal of Numerical Analysis)

References in corpus (6)

• Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients
• Strong and weak divergence in finite time of Euler's method for stochastic differential equations with non-globally Lipschitz continuous coefficients
• A note on tamed Euler approximations
• From the master equation to mean field game limit theory: A central limit theorem
• Time discretization of FBSDE with polynomial growth drivers and reaction-diffusion PDEs
• Adaptive Euler-Maruyama Method for SDEs with Non-globally Lipschitz Drift: Part I, Finite Time Interval

Cited by in corpus (4)

• First order convergence of Milstein schemes for McKean-Vlasov equations and interacting particle systems
• Well-posedness and tamed Euler schemes for McKean-Vlasov equations driven by Lévy noise
• Approximations of Mckean-Vlasov SDEs with Irregular Coefficients
• Explicit Numerical Approximations for McKean-Vlasov Neutral Stochastic Differential Delay Equations