Error estimates for finite-dimensional approximations of Hamilton-Jacobi-Bellman equations on the Wasserstein space
Samuel Daudin, Joe Jackson, Benjamin Seeger
TL;DR
This work quantifies the convergence of finite-dimensional HJB equations $V^N$ to the lifted Wasserstein-HJB solution $U$ for problems with purely common noise and nonconvex Hamiltonians. It leverages Lions’ lifting to formulate a Hilbert-space (Lions) viscosity framework and implements a doubling-of-variables strategy augmented by novel simultaneous quantization techniques to control discretization errors. The authors establish dimension-dependent algebraic convergence rates, proving that in 1D the rate is favorable while higher dimensions incur degradation; they also discuss improvements under stronger structural assumptions such as bounded state spaces, periodicity, and smoother data, and they acknowledge a fundamental lower bound $N^{-1/d}$ guiding the rate. Overall, the results provide rigorous, quantitative error bounds for approximating HJB equations on measure spaces by finite-dimensional, tractable formulations, with implications for numerical methods in mean-field-type control under common noise.
Abstract
In this paper, we study a Hamilton-Jacobi-Bellman (HJB) equation set on the Wasserstein space $\mathcal{P}_2(\mathbb{R}^d)$, with a second order term arising from a purely common noise. We do not assume that the Hamiltonian is convex in the momentum variable, which means that we cannot rely on representation formulas coming from mean field control. In this setting, Gangbo, Mayorga, and Święch showed via viscosity solutions methods that the HJB equation on $\mathcal{P}_2(\mathbb{R}^d)$ can be approximated by a sequence of finite-dimensional HJB equations. Our main contribution is to quantify this convergence result. The proof involves a doubling of variables argument, which leverages the Hilbertian approach of P.L. Lions for HJB equations in the Wasserstein space, rather than working with smooth metrics which have been used to obtain similar results in the presence of idiosyncratic noise. In dimension one, our doubling of variables argument is made relatively simply by the rigid structure of one-dimensional optimal transport, but in higher dimension the argument is significantly more complicated, and relies on some estimates concerning the "simultaneous quantization" of probability measures.
