The Annals of Applied Probability

Cubature on Wiener space for McKean–Vlasov SDEs with smooth scalar interaction

Dan Crisan and Eamon McMurray

Full-text: Access denied (no subscription detected)

We're sorry, but we are unable to provide you with the full text of this article because we are not able to identify you as a subscriber. If you have a personal subscription to this journal, then please login. If you are already logged in, then you may need to update your profile to register your subscription. Read more about accessing full-text


We present two cubature on Wiener space algorithms for the numerical solution of McKean–Vlasov SDEs with smooth scalar interaction. First, we consider a method introduced in [Stochastic Process. Appl. 125 (2015) 2206–2255] under a uniformly elliptic assumption and extend the analysis to a uniform strong Hörmander assumption. Then we introduce a new method based on Lagrange polynomial interpolation. The analysis hinges on sharp gradient to time-inhomogeneous parabolic PDEs bounds. These bounds may be of independent interest. They extend the classical results of Kusuoka and Stroock [J. Fac. Sci., Univ. Tokyo, Sect. 1A, Math. 32 (1985) 1–76] and Kusuoka [J. Math. Sci. Univ. Tokyo 10 (2003) 261–277] further developed in [J. Funct. Anal. 263 (2012) 3024–3101; J. Funct. Anal. 268 (2015) 1928–1971; Cubature Methods and Applications (2013), Springer, Cham] and, more recently, in [Probab. Theory Related Fields 171 (2016) 97–148]. Both algorithms are tested through two numerical examples.

Article information

Ann. Appl. Probab., Volume 29, Number 1 (2019), 130-177.

Received: August 2016
Revised: January 2018
First available in Project Euclid: 5 December 2018

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Primary: 60H30: Applications of stochastic analysis (to PDE, etc.) 60H10: Stochastic ordinary differential equations [See also 34F05] 60H07: Stochastic calculus of variations and the Malliavin calculus
Secondary: 34F05: Equations and systems with randomness [See also 34K50, 60H10, 93E03] 35K55: Nonlinear parabolic equations 35K65: Degenerate parabolic equations

Cubature on Wiener space Lagrange polynomial interpolation McKean–Vlasov SDEs Kusuoka–Stroock functions


Crisan, Dan; McMurray, Eamon. Cubature on Wiener space for McKean–Vlasov SDEs with smooth scalar interaction. Ann. Appl. Probab. 29 (2019), no. 1, 130--177. doi:10.1214/18-AAP1407.

Export citation


  • [1] Abramowitz, M. and Stegun, I. A. (1972). Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables.
  • [2] Antonelli, F. and Kohatsu-Higa, A. (2002). Rate of convergence of a particle method to the solution of the McKean–Vlasov equation. Ann. Appl. Probab. 12 423–476.
  • [3] Bayer, C. and Friz, P. K. (2013). Cubature on Wiener space: Pathwise convergence. Appl. Math. Optim. 67 261–278.
  • [4] Bossy, M. (2005). Some stochastic particle methods for nonlinear parabolic PDEs. In ESAIM: Proceedings 15 18–57. EDP Sciences.
  • [5] Bossy, M. and Talay, D. (1997). A stochastic particle method for the McKean–Vlasov and the Burgers equation. Math. Comp. 66 157–192.
  • [6] Cardaliaguet, P. (2012). Notes on mean field games (from P.-L. Lyons’ lectures at Collège de France).
  • [7] Cattiaux, P. and Mesnager, L. (2002). Hypoelliptic non-homogeneous diffusions. Probab. Theory Related Fields 123 453–483.
  • [8] Chassagneux, J.-F., Crisan, D. and Delarue, F. (2014). A Probabilistic approach to classical solutions of the master equation for large population equilibria. ArXiv e-prints. Available at arXiv:1411.3009.
  • [9] Crisan, D. and Delarue, F. (2012). Sharp derivative bounds for solutions of degenerate semi-linear partial differential equations. J. Funct. Anal. 263 3024–3101.
  • [10] Crisan, D. and Ghazali, S. (2007). On the convergence rates of a general class of weak approximations of SDEs. In Stochastic Differential Equations: Theory and Applications. Interdiscip. Math. Sci. 2 221–248. World Sci. Publ., Hackensack, NJ.
  • [11] Crisan, D., Kurtz, T. G. and Lee, Y. (2014). Conditional distributions, exchangeable particle systems, and stochastic partial differential equations. Ann. Inst. Henri Poincaré Probab. Stat. 50 946–974.
  • [12] Crisan, D., Litterer, C. and Lyons, T. (2015). Kusuoka–Stroock gradient bounds for the solution of the filtering equation. J. Funct. Anal. 268 1928–1971.
  • [13] Crisan, D. and Manolarakis, K. (2012). Solving backward stochastic differential equations using the cubature method: Application to nonlinear pricing. SIAM J. Financial Math. 3 534–571.
  • [14] Crisan, D. and Manolarakis, K. (2014). Second order discretization of backward SDEs and simulation with the cubature method. Ann. Appl. Probab. 24 652–678.
  • [15] Crisan, D., Manolarakis, K. and Nee, C. (2013). Cubature methods and applications. In Paris–Princeton Lectures on Mathematical Finance 2013. Lecture Notes in Math. 2081 203–316. Springer, Cham.
  • [16] Crisan, D. and McMurray, E. (2016). Smoothing properties of McKean–Vlasov SDEs. Probab. Theory Related Fields 171 97–148.
  • [17] Crisan, D. and Ortiz-Latorre, S. (2013). A Kusuoka–Lyons–Victoir particle filter. In Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 469 20130076.
  • [18] de Raynal, P. C. and Trillos, C. G. (2015). A cubature based algorithm to solve decoupled Mckean–Vlasov forward–backward stochastic differential equations. Stochastic Process. Appl. 125 2206–2255.
  • [19] Gyurkó, L. G. and Lyons, T. J. (2010). Efficient and practical implementations of cubature on Wiener space. In Stochastic Analysis 2010. 73–111. Springer, Berlin.
  • [20] Jourdain, B., Méléard, S. and Woyczynski, W. A. (2008). Nonlinear SDEs driven by Lévy processes and related PDEs. ALEA Lat. Am. J. Probab. Math. Stat. 4 1–29.
  • [21] Jourdain, B., Méléard, S. and Woyczynski, W. A. (2008). Nonlinear SDEs driven by Lévy processes and related PDEs. ALEA Lat. Am. J. Probab. Math. Stat. 4 1–29.
  • [22] Kloeden, P. E. and Platen, E. (1992). Numerical Solution of Stochastic Differential Equations. Applications of Mathematics (New York) 23. Springer, Berlin.
  • [23] Kohatsu-Higa, A. and Ogawa, S. (1997). Weak rate of convergence for an Euler scheme of nonlinear SDE’s. Monte Carlo Methods Appl. 3 327–345.
  • [24] Kolokoltsov, V. N. (2004). Hydrodynamic limit of coagulation-fragmentation type models of $k$-nary interacting particles. J. Stat. Phys. 115 1621–1653.
  • [25] Kolokoltsov, V. N. (2010). Nonlinear Markov Processes and Kinetic Equations. Cambridge Tracts in Mathematics 182. Cambridge Univ. Press, Cambridge.
  • [26] Kotelenez, P. M. and Kurtz, T. G. (2010). Macroscopic limits for stochastic partial differential equations of McKean–Vlasov type. Probab. Theory Related Fields 146 189–222.
  • [27] Kunita, H. (1990). Stochastic Flows and Stochastic Differential Equations. Cambridge Studies in Advanced Mathematics 24. Cambridge Univ. Press, Cambridge.
  • [28] Kurtz, T. G. and Protter, P. E. (1996). Weak convergence of stochastic integrals and differential equations. In Probabilistic Models for Nonlinear Partial Differential Equations. Lecture Notes in Math. 1627 1–41. Springer, Berlin.
  • [29] Kusuoka, S. (2001). Approximation of expectation of diffusion process and mathematical finance. In Taniguchi Conference on Mathematics Nara ’98. Adv. Stud. Pure Math. 31 147–165. Math. Soc. Japan, Tokyo.
  • [30] Kusuoka, S. (2003). Malliavin calculus revisited. J. Math. Sci. Univ. Tokyo 10 261–277.
  • [31] Kusuoka, S. and Stroock, D. (1985). Applications of the Malliavin calculus. II. J. Fac. Sci., Univ. Tokyo, Sect. 1A, Math. 32 1–76.
  • [32] Kusuoka, S. and Stroock, D. (1987). Applications of the Malliavin calculus. III. J. Fac. Sci., Univ. Tokyo, Sect. 1A, Math. 34 391–442.
  • [33] Lee, W. and Lyons, T. (2015). The adaptive patched cubature filter and its implementation. arXiv preprint. Available at arXiv:1509.04239.
  • [34] Lions, P.-L. (2014) Cours au collège de france. Available at
  • [35] Litterer, C. and Lyons, T. (2011). Introducing cubature to filtering. In The Oxford Handbook of Nonlinear Filtering 768–796. Oxford Univ. Press, Oxford.
  • [36] Lyons, T. and Victoir, N. (2004). Cubature on Wiener space. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 460 169–198.
  • [37] Morale, D., Capasso, V. and Oelschläger, K. (2005). An interacting particle system modelling aggregation behavior: From individuals to populations. J. Math. Biol. 50 49–66.
  • [38] Ninomyia, S. and Victoir, N. (2008). Weak approximation scheme of stochastic differential equations and applications to derivatives pricing. Appl. Math. Finance. 15 107–121.
  • [39] Ogawa, S. (1995). Some problems in the simulation of nonlinear diffusion processes. Math. Comput. Simulation 38 217–223.
  • [40] Ricketson, L. F. (2015). A multilevel Monte Carlo method for a class of McKean–Vlasov processes. Available at
  • [41] Tachet, R. (2011). Des Combes. Calibration non paramétriques de modèles en finance. Theses, Ecole Centrale Paris.
  • [42] Teichmann, J. (2006). Calculating the Greeks by cubature formulae. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 462 647–670.