Simultaneously Solving Infinitely Many LQ Mean Field Games In Hilbert Spaces: The Power of Neural Operators
Dena Firoozi, Anastasis Kratsios, Xuwei Yang
TL;DR
The paper develops a principled framework for solving infinite-dimensional LQ mean field games by learning the rules-to-equilibrium operator with neural operators. It establishes that regularized rules-to-equilibrium maps are locally Lipschitz and PAC-learnable by Lipschitz Residually-Guided Neural Operators in Hilbert spaces, with explicit finite-sample guarantees. A universal Lipschitz NO approximation theorem, together with optimal-transport–based generalization bounds, yields practical guarantees on small NOs (depth O(1), width O(ε^{-c}) up to log factors) to approximate the operator across a broad class of problem variations. Under structured input sets (exponentially ellipsoidal, tempered sampling), the results provide faster rates and scalable learning for infinite families of LQ MFGs, enabling robust, data-driven solution descriptions across perturbations and continuum-parameterized agents.
Abstract
Traditional mean-field game (MFG) solvers operate on an instance-by-instance basis, which becomes infeasible when many related problems must be solved (e.g., for seeking a robust description of the solution under perturbations of the dynamics or utilities, or in settings involving continuum-parameterized agents.). We overcome this by training neural operators (NOs) to learn the rules-to-equilibrium map from the problem data (``rules'': dynamics and cost functionals) of LQ MFGs defined on separable Hilbert spaces to the corresponding equilibrium strategy. Our main result is a statistical guarantee: an NO trained on a small number of randomly sampled rules reliably solves unseen LQ MFG variants, even in infinite-dimensional settings. The number of NO parameters needed remains controlled under appropriate rule sampling during training. Our guarantee follows from three results: (i) local-Lipschitz estimates for the highly nonlinear rules-to-equilibrium map; (ii) a universal approximation theorem using NOs with a prespecified Lipschitz regularity (unlike traditional NO results where the NO's Lipschitz constant can diverge as the approximation error vanishes); and (iii) new sample-complexity bounds for $L$-Lipschitz learners in infinite dimensions, directly applicable as the Lipschitz constants of our approximating NOs are controlled in (ii).
