## The Annals of Applied Probability

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

#### Abstract

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

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

Dates
Revised: January 2018
First available in Project Euclid: 5 December 2018

https://projecteuclid.org/euclid.aoap/1544000427

Digital Object Identifier
doi:10.1214/18-AAP1407

Mathematical Reviews number (MathSciNet)
MR3910002

Zentralblatt MATH identifier
07039123

#### Citation

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. https://projecteuclid.org/euclid.aoap/1544000427

#### References

• [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). https://www.ceremade.dauphine.fr/~cardaliaguet/MFG20130420.pdf.
• [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 http://www.college-de-france.fr/site/pierre-louis-lions/seminar-2014-11-14-11h15.htm.
• [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 https://arxiv.org/abs/1508.02299.
• [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.