Convergence rate for random walk approximations of mean field BSDEs
Boualem Djehiche, Hannah Geiss, Stefan Geiss, Céline Labart, Jani Nykänen
TL;DR
This paper analyzes the convergence rate of scaled random-walk approximations to mean-field backward stochastic differential equations (BSDEs) in Wasserstein-type distances, using a freezing technique instead of particle methods. It extends Donsker-type convergence results to the mean-field setting, where a time-singularity in the generator arises from Hölder-terminal data, and introduces a modified Hölder continuity |r|_{(ε,β)} to recover polynomial rates up to logarithmic factors. The authors establish coupling-based pathwise and Wasserstein bounds, with rates that include a log term and depend on ε, α, and β, and they provide a practical, non-particle algorithm along with numerical validation showing the expected rate behavior. The work advances numerical methods for nonlocal mean-field BSDEs and yields tools that are broadly applicable to the probabilistic analysis of mean-field problems and their associated nonlocal PDEs.
Abstract
We study the rate of convergence w.r.t.~a Wasserstein type distance for random walk approximations of mean field BSDEs. Our method does not use the particle method but instead a freezing technique. We extend results by Briand, Ch. Geiss, S. Geiss, and Labart [Bernoulli, 27(2) 2021] about the rate of convergence of a Donsker-type theorem for BSDEs from the classical setting to the mean field setting. In this connection the mean field setting leads to new phenomena and requires new techniques that should be of independent interest: The Hölder continuous terminal condition causes a singularity in time of the generator when seen as a generator in the non-mean field setting. To handle this singularity we introduce a concept of modified Hölder continuity by which we are able to achieve, up to a logarithmic term, the same polynomial approximation rates as in the classical non-mean field setting (in fact, already when approximating the Brownian motion itself a logarithmic term is necessary). Moreover, the exploited freezing technique of the mean field terms yields to the problem to handle the quantitative behavior of several different generators. Using BMO-techniques we obtain the rate of convergence for the integrated gradient process in the scale of Lorentz (type) spaces of exponential type.
