Topological bifurcations in a mean-field game
Ali Akbar Rezaei Lori, Piyush Grover
TL;DR
This work develops a reduced-order, four-dimensional Hamiltonian BVP for a finite-horizon mean-field game by constraining agent dynamics to the first two moments, the mean $X(t)$ and variance $S(t)$. It leverages phase-space geometry near ergodic equilibria (notably saddle$ imes$saddle and saddle$ imes$center types) and invariant manifolds of unstable periodic orbits to explain the topological structure and multiplicity of solution branches as the horizon $T$ grows. The authors validate the reduced model against full-order MFG PDEs, showing consistent topology (e.g., rotations in the ergodic phase) and density evolution features such as bimodal peaks. The approach highlights a topological bottleneck mechanism that governs transitions between initial and final conditions and suggests natural extensions to higher-dimensional ROMs and richer dynamical regimes for mean-field control and transport problems.
Abstract
Mean-field games (MFG) provide a statistical physics inspired modeling framework for decision making in large-populations of strategic, non-cooperative agents. Mathematically, these systems consist of a forward-backward in time system of two coupled nonlinear partial differential equations (PDEs), namely the Fokker-Plank and the Hamilton-Jacobi-Bellman equations, governing the agent state and control distribution, respectively. In this work, we study a finite-time MFG with a rich global bifurcation structure using a reduced-order model (ROM). The ROM is a 4D two-point boundary value problem obtained by restricting the controlled dynamics to first two moments of the agent state distribution, i.e., the mean and the variance. Phase space analysis of the ROM reveals that the invariant manifolds of periodic orbits around the so-called `ergodic MFG equilibrium' play a crucial role in determining the bifurcation diagram, and impart a topological signature to various solution branches. We show a qualitative agreement of these results with numerical solutions of the full-order MFG PDE system.
