On the asymptotic nature of first order mean field games

TL;DR

This paper proves that symmetric N-player Nash equilibria in finite horizon first order mean field games converge to solutions of the mean field game, linking Lagrangian and PDE formulations, and establishing convergence in Markov strategies.

Contribution

It provides a simple proof of convergence of Nash equilibria to mean field game solutions and connects Lagrangian and PDE approaches.

Findings

Nash equilibria converge to mean field game solutions

Lagrangian solutions relate to PDE formulations

Convergence established for distributed Markov strategies

Abstract

For a class of finite horizon first order mean field games and associated N-player games, we give a simple proof of convergence of symmetric N-player Nash equilibria in distributed open-loop strategies to solutions of the mean field game in Lagrangian form. Lagrangian solutions are then connected with those determined by the usual mean field game system of two coupled first order PDEs, and convergence of Nash equilibria in distributed Markov strategies is established.

Equations

$$\phi ♯m(A):=m({\phi }^{-1}(A))\forall A\in \mathcal{B}(Y),$$

$$d_p({\mu }_1,{\mu }_2):=\inf \left\{{\left({\int }_{X\times X}{d}_X^p(x,y)\mathrm{d}\gamma (x,y)\right)}^{\frac{1}{p}}|\gamma \in \mathcal{P}(X\times X),{\pi }_1♯\gamma ={\mu }_1,{\pi }_2♯\gamma ={\mu }_2\right\},$$

$$d_1({\mu }_1,{\mu }_2)=\sup \left\{{\int }_{X}f(x)\mathrm{d}\left({\mu }_1-{\mu }_2\right)(x)|f\in \text{mbox}{Lip}_1(X)\right\},$$

$${\underline{c}}_{\ell }|\alpha {|}^q-C_{\ell }\le \ell (\alpha ,x,\mu )\le {\overline{c}}_{\ell }|\alpha {|}^q+C_{\ell }\forall \alpha \in {\mathbb{R}}^d,x\in {\mathbb{R}}^d,\mu \in {\mathcal{P}}_1({\mathbb{R}}^d),$$

$$\mathcal{A}(x):=\left\{\gamma\in\Gamma\;\big{|}\;\gamma\in W^{1,q}([0,T];\mathbb{R}^{d})\;\;\mbox{and }\;\gamma(0)=x\right\},\quad x\in\mathbb{R}^{d}.$$

$$j_i^N({\gamma }_1,\dots ,{\gamma }_N):={\int }_0^T\ell \left({\dot{\gamma }}_i(t),{\gamma }_i(t),({\gamma }_j(t){)}_{j\ne i}\right)\mathrm{d}t+\Phi \left({\gamma }_i(T),({\gamma }_j(T){)}_{j\ne i}\right),$$

$$\ell\left(\dot{\gamma}_{i}(t),\gamma_{i}(t),(\gamma_{j}(t))_{j\neq i}\right):=\ell\left(\dot{\gamma}_{i}(t),\gamma_{i}(t),\frac{1}{N-1}\sum_{{}_{j\neq i}}\delta_{\gamma_{j}(t)}\right),\hskip 5.69046pt\mbox{with the same convention for $\Phi$.}$$

$$j^N({\gamma }_1;({\gamma }_j{)}_{j=2}^N):={\int }_0^T\ell \left({\dot{\gamma }}_1(t),{\gamma }_1(t),({\gamma }_j(t){)}_{j=2}^N\right)\mathrm{d}t+\Phi \left({\gamma }_1(T),({\gamma }_j(T){)}_{j=2}^N\right),$$

$$\mathcal{A}:=\left\{\gamma:\mathbb{R}^{d}\to\Gamma\;\big{|}\;\gamma\text{ is Borel measurable and }\gamma(x)\in\mathcal{A}(x)\text{ for all }x\in\mbox{supp}(m_{0})\right\}.$$

$$J_i^N({\gamma }_1,\dots ,{\gamma }_N):={\int }_{({\mathbb{R}}^d{)}^N}{j}_i^N({\gamma }_1^{x_1},\dots ,{\gamma }_N^{x_N}){\otimes }_{j=1}^N\mathrm{d}{m}_0(x_j).$$

$$J^N({\gamma }_1;({\gamma }_j{)}_{j=2}^N):={\int }_{({\mathbb{R}}^d{)}^N}{j}^N({\gamma }_1^{x_1};({\gamma }_j^{x_j}{)}_{j=2}^N){\otimes }_{j=1}^N\mathrm{d}{m}_0(x_j),$$

$$J_i^N({\gamma }_1,\dots ,{\gamma }_N)=J_{\sigma (i)}^N({\gamma }_{\sigma (1)},\dots ,{\gamma }_{\sigma (N)}).$$

$$J^N({\gamma }_{i,N};({\gamma }_{j,N}{)}_{j\ne i})\le J^N(\gamma ;({\gamma }_{j,N}{)}_{j\ne i})\forall \gamma \in \mathcal{A},i=1,\dots ,N.$$

$$\int_{0}^{T}|\dot{\gamma}_{i,N}^{x}(t)|^{q}{\rm d}t\leq C\hskip 11.38092pt\mbox{{\rm for $m_{0}$-a.e.} $x\in\mathbb{R}^{d}$},\;\;i=1,\ldots,N.$$

$$j_i(\gamma )={\int }_0^T{\ell }_i(\dot{\gamma }(t),\gamma (t),t)\mathrm{d}t+{\Phi }_i(\gamma (T)),$$

$$\begin{array}[]{rcl}\ell_{i}(\alpha,x,t)&:=&\int_{\mathbb{R}^{(N-1)\times d}}\ell\left(\alpha,x,(\gamma_{j,N}^{x_{j}}(t))_{j\neq i}\right)\otimes_{j\neq i}{\rm d}m_{0}(x_{j}),\\[6.0pt] \Phi_{i}(x)&:=&\int_{\mathbb{R}^{(N-1)\times d}}\Phi\left(x,(\gamma_{j,N}^{x_{j}}(T))_{j\neq i}\right)\otimes_{j\neq i}{\rm d}m_{0}(x_{j}).\end{array}$$

$$J^N({\gamma }_{i,N};({\gamma }_{j,N}{)}_{j\ne i})={\int }_{{\mathbb{R}}^d}{j}_i({\gamma }_{i,N}^x)\mathrm{d}{m}_0(x)\le J^N(\gamma ;({\gamma }_{j,N}{)}_{j\ne i})={\int }_{{\mathbb{R}}^d}{j}_i({\gamma }^x)\mathrm{d}{m}_0(x).$$

$${\text{SS}}_i^N(x):=\text{mbox}argmin\left\{j_i(\gamma )|\gamma \in \mathcal{A}(x)\right\}\subseteq \Gamma .$$

$$j_i({\gamma }_n^{x_n})\le j_i({\gamma }^{x_n})\forall {\gamma }^{x_n}\in \mathcal{A}(x_n),$$

$$Q_C:=\left\{\gamma \in W^{1,q}([0,T];{\mathbb{R}}^d)|{\int }_0^T|\dot{\gamma }(t){|}^q\mathrm{d}t\le C,\gamma (0)\in \text{mbox}supp(m_0)\right\}.$$

$$\gamma (t)\in K_C\text{for all }t\in [0,T],\gamma \in Q_C.$$

$$J^N({\gamma }_i;({\gamma }_j{)}_{j\ne i})={\int }_{{\Gamma }^N}{j}^N({\tilde{\gamma }}_i;({\tilde{\gamma }}_j{)}_{j\ne i}){\otimes }_{j=1}^N\mathrm{d}{m}_j({\tilde{\gamma }}_j).$$

$${\mathcal{A}}_{rel}:=\left\{m\in \mathcal{P}(\Gamma )|e_0♯m=m_0,\text{mbox}supp(m)\subseteq Q_C\right\},$$

$$d_1(e_t♯m,e_s♯m)\le C^{\prime }|t-s{|}^{\frac{1}{q^{\prime }}}\forall s,t\in [0,T].$$

$$J_{rel}^N(m_i;(m_j{)}_{j\ne i}):={\int }_{Q_C^N}{j}^N({\gamma }_i;({\gamma }_j{)}_{j\ne i}){\otimes }_{j=1}^N\mathrm{d}{m}_j({\gamma }_j).$$

$$J_{rel}^N(m_{i,N};(m_{j,N}{)}_{j\ne i})\le J_{rel}^N(m;(m_{j,N}{)}_{j\ne i})\forall m\in {\mathcal{A}}_{rel},i=1,\dots ,N.$$

$$\mathcal{A}_{rel}\ni\mu\rightrightarrows\SS^{N}(\mu):=\mbox{argmin}\left\{J^{N}_{rel}(\mu^{\prime},\mu,\ldots,\mu)\;\big{|}\;\mu^{\prime}\in\mathcal{A}_{rel}\right\}\subseteq\mathcal{A}_{rel}.$$

$$J_{rel}^N({\nu }_n;{\mu }_n,\dots ,{\mu }_n)\le J_{rel}^N({\mu }^{\prime };{\mu }_n,\dots ,{\mu }_n)\forall {\mu }^{\prime }\in {\mathcal{A}}_{rel}.$$

$$J_{rel}^N(\bar{\nu };\bar{\mu },\dots ,\bar{\mu })\le J_{rel}^N({\mu }^{\prime };\bar{\mu },\dots ,\bar{\mu })\forall {\mu }^{\prime }\in {\mathcal{A}}_{rel},$$

$$\begin{array}[]{rcl}\hat{\ell}(\alpha,x,t)&:=&\int_{Q_{C}^{N-1}}\ell\left(\alpha,x,(\gamma_{j}(t))_{j=2}^{N}\right)\otimes_{j=2}^{N}{\rm d}m_{N}(\gamma_{j}),\\[6.0pt] \hat{\Phi}(x)&:=&\int_{Q_{C}^{N-1}}\Phi\left(x,(\gamma_{j}(T))_{j=2}^{N}\right)\otimes_{j=2}^{N}{\rm d}m_{N}(\gamma_{j}),\end{array}$$

Full text

Dipartimento di Matematica “Tullio Levi-Civita”, Università degli Studi di Padova, via Trieste, 63, 35121 Padova, Italia (fischer@math.unipd.it)

Institut de recherche XLIM-DMI, UMR-CNRS 7252, Faculté des Sciences et Techniques, Université de Limoges, 87060 Limoges, France (francisco.silva@unilim.fr)

AMS-Subject Classification: 49N70, 60B10, 91A06, 91A13.

Keywords: Mean field games, Lagrangian form, deterministic dynamics, Nash equilibrium, distributed strategies.

1. Introduction

The purpose of this article is to illustrate a simple way of establishing convergence of open-loop Nash equilibria in the case of first-order non-stationary Mean Field Games (MFGs). Introduced by J.-M. Lasry and P.-L. Lions and, independently, by M. Huang, R.P. Malhamé and P.E. Caines about fifteen years ago (cf. [34, 31]), mean field games are limit models for non-cooperative symmetric $N$-player differential games as the number of players $N$ tends to infinity; see, for instance, the lecture notes [13] and the recent two-volume work [18]. The notion of solution usually adopted for the prelimit models is that of a Nash equilibrium. A standard way of making the connection with the limit model rigorous is to show that a solution of the mean field game yields approximate Nash equilibria for the $N$-player games, with approximation error vanishing as $N\to\infty$. In the opposite direction, one aims to prove that a sequence of $N$-player Nash equilibria converges, as $N$ tends to infinity, to the mean field game limit.

When Nash equilibria are considered in stochastic open-loop strategies, then their convergence is well understood and can be established under mild conditions; see [28] and [32], both for finite horizon games with general, possibly degenerate, Brownian dynamics. The convergence analysis is much harder when Nash equilibria are defined over Markov feedback strategies with full state information.

A first result in this setting was given by Gomes, Mohr, and Souza [29] for continuous time games with finite state space. There, convergence of Markovian Nash equilibria is proved, but only if the time horizon is small enough. A breakthrough was achieved by Cardaliaguet, Delarue, Lasry, and Lions in [15]. In the setting of games with non-degenerate Brownian dynamics, possibly including common noise, convergence to the mean field game limit is established there for arbitrary time horizon provided the so-called master equation associated with the mean field game possesses a unique sufficiently regular solution. In this case, the convergence analysis can be refined, yielding not only convergence of minimal costs and propagation of chaos for the Nash equilibrium state trajectories, but also fluctuation and large deviations results for the associated empirical measures; see Cecchin and Pelino [19] and, independently, Bayraktar and Cohen [7] for finite state games and the papers by Delarue, Lacker, and Ramanan [23, 24] for Brownian dynamics without or with common noise.

Well-posedness of the master equation implies uniqueness of solutions for the mean field game. But also the situation where the mean field game possesses multiple solutions, while the $N$-player Nash equilibria in full Markov feedback strategies are still uniquely determined, occurs. In this case, the convergence problem is in part open. The most general result appears to be the recent preprint [33] by Lacker. There, for a class of games with non-degenerate Brownian dynamics, it is shown that all limit points of the $N$-player Nash equilibria are concentrated on weak solutions of the mean field game; these are more general than randomizations of ordinary (“strong”) solutions of the mean field game. In the three recent works [36], [25] and [20], the authors present case studies, giving criteria for characterizing those mean field game solutions that can be attained as limits of feedback Nash equilibria with full state information.

Here, we consider a much simpler situation: The underlying dynamics are deterministic with direct control of players’ states; randomness enters only through the players’ initial positions, which are assumed to be independently and identically distributed. Thanks to the deterministic dynamics without explicit interaction, players can directly control their entire state trajectories. Thus, the set of strategies (or actions) of each player consists of trajectory-valued functions depending on time and player’s own position. We call these actions distributed open-loop strategies. By considering randomizations of these strategies, obtaining what we call randomized distributed open-loop strategies, we obtain existence of symmetric Nash equilibria in this new class of actions (through a standard fixed point argument), as well as precompactness of sequences of such equilibria. Convergence to the mean field game equilibrium in Lagrangian form (see e.g. [16, 8, 17, 10])) along weakly converging subsequences of Nash equilibria is then established in Theorem 3.1 by showing that the variance of the associated empirical measures (evaluated by integrating over test functions from a countable family) tends to zero as $N\to\infty$. Thus, to obtain convergence, we essentially prove a law of large numbers for a triangular array $(Y^{N}_{i})_{i\in\{1,\ldots,N\},N\in\mathbb{N}}$ where $Y^{N}_{1},\ldots,Y^{N}_{N}$ are independent and identically distributed with common marginal law that however depends on $N$. Let us point out that after obtaining a suitable compactness property for the set of admissible strategies, the convergence result in Theorem 3.1 can actually be deduced from the more general results in [28] or [32]. Thus, the main purpose of this result is therefore to give a simple proof of convergence, which takes advantage of the structure of the dynamics and the deterministic nature of the underlying differential games. We believe that the techniques considered here could be useful in order to justify the asymptotic nature of more sophisticated deterministic MFGs (see e.g. [10, 11] dealing with state constrained problems).

In the second part of this article, we consider additional second order assumptions on the data and we assume that the initial distribution is absolutely continuous with respect to the Lebesgue measure. In this framework, and following [16], we link the notion of Lagrangian MFG equilibrium to the well known PDE characterization of mean field games in terms of two coupled first order partial differential equations, namely a backward first order Hamilton-Jacobi-Bellman equation and a forward continuity equation; see equation $(MFG)$ in Subsection 4.1 below. Under the stronger assumptions mentioned above, for each $N\in\mathbb{N}$, any symmetric randomized distributed open-loop Nash equilibrium of the $N$-player game can be identified with a symmetric distributed open-loop Nash equilibrium (non-randomized). Moreover, the strategies associated to such equilibria can be described by controls which are feedback with respect to the individual states. We call these actions distributed Markov strategies. The latter are determined by the solutions of a coupled system of first order PDEs analogous to the mean field game PDE system; see Eq. $(MFG_{N})$ in Subsection 4.2. This is in contrast with $N$-player Nash equilibria defined over Markov feedback strategies with full state information, which are determined through a system of $N$ coupled PDEs. Building on the equivalence of characterizations and the convergence result in Theorem 3.1, we also establish in Theorem 4.1 the convergence of solutions of $(MFG_{N})$ to solutions of $(MFG)$, as well as the convergence of the Nash equilibria in distributed Markov strategies.

The rest of this paper is organized as follows. In Section 2, we introduce the $N$-player games together with some notation and our standing assumptions. Existence of Nash equilibria in randomized distributed open-loop strategies is verified in Proposition 2.1. In Section 3, the associated mean field game is introduced in Lagrangian form; see Definition 3.1 and Proposition 3.1. We show convergence of symmetric $N$-player randomized distributed open-loop Nash equilibria to the mean field game limit in Theorem 3.1. Section 4 links, under additional assumptions, the mean field game in Lagrangian form with the mean field game PDE system. Similarly, $N$-player Nash equilibria in distributed Markov strategies are described in terms of a PDE system analogous to the limit system, but dependent on the number of players $N$. The corresponding convergence results are given in Theorem 4.1 and Corollary 4.1, respectively.

2. The $N$-player game

Before introducing the $N$-player deterministic differential game that we are interested in, let us first fix some notations. Given a Polish space $(X,d_{X})$, we denote by $\mathcal{P}(X)$ the set of probability measures over $X$. If $(Y,d_{Y})$ is another Polish space and $\varphi:X\to Y$ a Borel measurable function, then the push-forward $\varphi\sharp m\in\mathcal{P}(Y)$ of a measure $m\in\mathcal{P}(X)$ by $\varphi$ is defined by

[TABLE]

where $\mathcal{B}(Y)$ denotes the $\sigma$-algebra of Borel sets in $Y$. Given $p\in[1,+\infty)$ we denote by $\mathcal{P}_{p}(X)$ the set of probability measures $\mu$ over $X$ such that $\int_{X}d_{X}(x,x_{0})^{p}{\rm d}\mu(x)<+\infty$ for some $x_{0}\in X$. The set $\mathcal{P}_{p}(X)$ is endowed with the Monge-Kantorovich distance $d_{p}:\mathcal{P}_{p}(X)\times\mathcal{P}_{p}(X)\to[0,\infty)$ defined by

[TABLE]

where $\pi_{i}:X\times X\to X$ ($i=1$, $2$) is the projection on the $i$-th coordinate, that is, $\pi_{i}(x_{1},x_{2}):=x_{i}$. It is well known (see e.g. [37]) that in the particular case $p=1$ we have

[TABLE]

where $\mbox{Lip}_{1}(X)$ denotes the set of Lipschitz functions on $X$ with Lipschitz constant equal to one.

We will fix as data some functions $\ell:\mathbb{R}^{d}\times\mathbb{R}^{d}\times\mathcal{P}_{1}(\mathbb{R}^{d})\to\mathbb{R}$, $\Phi:\mathbb{R}^{d}\times\mathcal{P}_{1}(\mathbb{R}^{d})\to\mathbb{R}$ and a probability measure $m_{0}\in\mathcal{P}(\mathbb{R}^{d})$. The functions $\ell$, $\Phi$ will represent running and terminal costs, respectively, while $m_{0}$ will be an initial state distribution. We will assume that:

(A1) (i) The functions $\ell$ and $\Phi$ are continuous. Moreover, the following properties hold true:

• (i.1)

For every $(x,\mu)\in\mathbb{R}^{d}\times\mathcal{P}_{1}(\mathbb{R}^{d})$, the function $\ell(\cdot,x,\mu)$ is convex.

• (i.2)

There exist $q>1$, $\underline{c}_{\ell}>0$, $\overline{c}_{\ell}>0$ and $C_{\ell}>0$ such that

[TABLE]

where we denote by $|\cdot|$ the Eulidean norm in $\mathbb{R}^{d}$.

• (i.3)

The function $\Phi$ is bounded from below.

(ii) The support of $m_{0}$, denoted by $\mbox{supp}(m_{0})$, is a compact subset of $\mathbb{R}^{d}$.

Choose $T>0$, the finite time horizon, and set $\Gamma:=C([0,T];\mathbb{R}^{d})$, the space of continuous $\mathbb{R}^{d}$-valued trajectories on $[0,T]$. The space $\Gamma$ is naturally endowed with the topology of uniform convergence. Let $W^{1,q}([0,T];\mathbb{R}^{d})$ denote the Sobolev space of $\mathbb{R}^{d}$-valued absolutely continuous functions that possess first order weak sense partial derivatives in $L^{q}((0,T);\mathbb{R}^{d})$.

In order to introduce the game that we will consider, assume first that there are $N$ individuals ($N\geq 2$), which, from now on, will be called players, positioned at $x_{1},\ldots,x_{N}\in\mathbb{R}^{d}$ at time $t=0$. For $i\in\{1,\ldots,N\}$, player $i$ chooses a trajectory $\gamma_{i}\in\mathcal{A}(x_{i})$, where

[TABLE]

Given these initial positions, the cost $j_{i}^{N}:\prod_{j=1}^{N}\mathcal{A}(x_{j})\to\mathbb{R}$ for player $i$ is defined by

[TABLE]

where, for notational convenience, we have denoted

[TABLE]

Assumption (A1) implies that $j_{i}^{N}$ is well-defined. Note that defining

[TABLE]

we have that $j_{i}^{N}(\gamma_{1},\ldots,\gamma_{N})=j^{N}(\gamma_{i};(\gamma_{j})_{j\neq i})$.

In the game that we will consider, the initial position of each player is random, independent of the initial positions of the other players and with the same law $m_{0}$. In this new context, we define the set of distributed open-loop strategies of the players as

[TABLE]

For notational simplicity, for $\gamma\in\mathcal{A}$, we will write $\gamma^{x}:=\gamma(x)$. Given a profile of actions $(\gamma_{1},\ldots,\gamma_{N})\in\mathcal{A}^{N}$, it is natural to define the cost that it induces on player $i$ as the mean with respect to the initial conditions of costs having the form (2.4). Namely, the cost function $J_{i}^{N}:\mathcal{A}^{N}\to\mathbb{R}$ for player $i$ is defined as

[TABLE]

Recalling (2.5) and defining

[TABLE]

we have $J_{i}^{N}(\gamma_{1},\ldots,\gamma_{N})=J^{N}(\gamma_{i};(\gamma_{j})_{j\neq i})$ for all $i=1,\ldots,N$. In particular, the costs are symmetric, i.e. for every permutation $\displaystyle\sigma:\{1,\ldots,N\}\to\{1,\ldots,N\}$, we have that

[TABLE]

Let us recall the classical notion of Nash equilibrium when applied to the game defined by the $N$ players, the action set $\mathcal{A}$ (which is the same for all players) and the individual costs $(J_{i})_{i=1}^{N}$.

Definition 2.1**.**

We say that $(\gamma_{1,N},\ldots,\gamma_{N,N})\in\mathcal{A}^{N}$ is a Nash equilibrium in distributed open-loop strategies if

[TABLE]

The existence of a Nash equilibrium for this symmetric game is not immediate because of the lack of compactness of the set $\mathcal{A}$. The following simple result shows that the time derivatives of the strategies in a Nash equilibrium configuration (provided that it exists) enjoy a uniform boundedness property in $L^{q}((0,T);\mathbb{R}^{d})$. This fact will allow us to reduce the set of action strategies in Definition 2.1.

Lemma 2.1**.**

Assume that (A1) holds and that $(\gamma_{1,N},\ldots,\gamma_{N,N})\in\mathcal{A}^{N}$ is a Nash equilibrium. Then, there exists a constant $C>0$, independent of $N$, such that

[TABLE]

Proof.

Fix $i\in\{1,\ldots,N\}$ and define $j_{i}:W^{1,q}([0,T];\mathbb{R}^{d})\to\mathbb{R}$ as

[TABLE]

where $\ell_{i}:\mathbb{R}^{d}\times\mathbb{R}^{d}\times[0,T]\to\mathbb{R}$ and $\Phi_{i}:\mathbb{R}^{d}\to\mathbb{R}$ are defined by

[TABLE]

By Fubini’s theorem, and (2.7), for all $\gamma\in\mathcal{A}$ we have

[TABLE]

For $x\in\mathbb{R}^{d}$ define

[TABLE]

Assumption (A1) implies that $j_{i}$ is lower semicontinuous, w.r.t. the weak topology in $W^{1,q}([0,T];\mathbb{R}^{d})$ (see e.g. [22, Corollary 3.24]). Using this fact, the direct method in the Calculus of Variations and the first inequality in (2.3), we obtain that $\SS_{i}^{N}(x)\neq\emptyset$ for all $x\in\mathbb{R}^{d}$.

*Claim: * The set-valued map $\mathbb{R}^{d}\ni x\rightrightarrows\SS_{i}^{N}(x)\in 2^{\Gamma}$ takes closed values and is upper semicontinuous, i.e. for all closed sets $M\subseteq\Gamma$ we have that $\{x\in\mathbb{R}^{d}\;|\;\SS_{i}^{N}(x)\cap M\neq\emptyset\}$ is closed.

Let us assume for a moment that the claim is true. Then, as a consequence of its second statement, the set valued map $\SS_{i}^{N}$ is Borel measurable and, hence, by the Kuratowski-Ryll-Nardzewski Selection Theorem (see e.g. [1, Theorem 18.13]) we have the existence of a Borel measurable function $\mathbb{R}^{d}\ni x\mapsto\hat{\gamma}(x)\in\Gamma$ such that $\hat{\gamma}(x)\in\SS_{i}^{N}(x)$ for all $x\in\mathbb{R}^{d}$. Since $\hat{\gamma}\in\mathcal{A}$, relation (2.10) implies that $\gamma_{i,N}^{x}\in\SS_{i}^{N}(x)$ for $m_{0}$-almost every $x\in\mathbb{R}^{d}$. Taking $x\in\mathbb{R}^{d}$ such that $\gamma_{i,N}^{x}\in\SS_{i}^{N}(x)$, we have that $j_{i}(\gamma_{i,N}^{x})\leq j_{i}(\bar{\gamma}^{x})$, where $\bar{\gamma}^{x}(t):=x$ for all $t\in[0,T]$. Assumption ${\bf(A1)}$ implies that $j_{1}(\bar{\gamma}^{x})$ is bounded by a constant, which is uniform for $x\in\mbox{supp}(m_{0})$ and independent of $N$. Using this fact, the first inequality in (2.3) easily yields (2.8).

It remains to prove the claim. First note that if $(\gamma_{n}^{x})$ is a sequence in $\SS_{i}^{N}(x)$ converging to $\gamma^{x}$ uniformly in $[0,T]$, then, by the first inequality in (2.3), the sequence $(\dot{\gamma}_{n}^{x})$ is bounded in $L^{q}([0,T];\mathbb{R}^{d})$. If $g$ is a weak limit point of $(\dot{\gamma}_{n}^{x})$ in $L^{q}([0,T];\mathbb{R}^{d})$, then, passing to the limit along a subsequence in the relation $\gamma_{n}^{x}(t)=x+\int_{0}^{t}\dot{\gamma}_{n}^{x}(s){\rm d}s$ for all $t\in[0,T]$, we get that $\gamma^{x}\in W^{1,q}([0,T];\mathbb{R}^{d})$, with $\dot{\gamma}^{x}=g$, and, hence, the whole sequence $(\dot{\gamma}_{n}^{x})$ converges weakly to $\dot{\gamma}^{x}$ in $L^{q}([0,T];\mathbb{R}^{d})$. Using this fact and the weak lower semicontinuity of $j_{i}$ in $W^{1,q}([0,T];\mathbb{R}^{d})$, we obtain that $\gamma^{x}\in\SS_{i}^{N}(x)$ and, hence, $\SS_{i}^{N}(x)$ is closed in $\Gamma$. In order to show that $\SS_{i}^{N}$ is upper semicontinuous, let $M$ be a closed subset of $\Gamma$ and let $(x_{n})$ be a sequence in $\{x\in\mathbb{R}^{d}\;|\;\SS_{i}^{N}(x)\cap M\neq\emptyset\}$ converging to some $\bar{x}\in\mathbb{R}^{d}$. Then, by definition, there exists $\gamma_{n}^{x_{n}}\in\SS_{i}^{N}(x_{n})\cap M$. Arguing as before, using the first inequality in (2.3), we obtain that $(\dot{\gamma}_{n}^{x_{n}})$ is bounded in $L^{q}([0,T];\mathbb{R}^{d})$. This implies that, up to some subsequence, $\gamma_{n}^{x_{n}}$ converge uniformly to some $\gamma^{\bar{x}}\in W^{1,q}([0,T];\mathbb{R}^{d})$ and $\dot{\gamma}^{x_{n}}\to\dot{\gamma}^{\bar{x}}$ weakly in $L^{q}([0,T];\mathbb{R}^{d})$. Using that $M$ is closed, w.r.t. the uniform convergence, we get that $\gamma^{\bar{x}}\in M$. On the other hand, noticing that

[TABLE]

and the fact that any $\gamma\in\mathcal{A}(\bar{x})$ satisfies $\gamma-\bar{x}+x_{n}\in\mathcal{A}(x_{n})$, the weak lower semicontinuity of $j_{i}$ yields that $j_{i}(\gamma^{\bar{x}})\leq j_{i}(\gamma)\hskip 8.5359pt\forall\;\gamma\in\mathcal{A}(\bar{x})$, i.e, $\gamma^{\bar{x}}\in\SS_{i}^{N}(\bar{x})\cap M$, and, hence, $\bar{x}\in\{x\in\mathbb{R}^{d}\;|\;\SS_{i}^{N}(x)\cap M\neq\emptyset\}$, which implies that the latter set is closed.

∎

Now, we focus our attention on the existence of Nash equilibria for the described game. Using Lemma 2.1, a reformulation of the set of admissible strategies and cost functionals of the $N$-players game will be useful. Let $C>0$ be given by Lemma 2.1 and define the set

[TABLE]

Since $\mbox{supp}(m_{0})$ is compact, Hölder’s inequality yields the existence of a compact set $K_{C}\subseteq\mathbb{R}^{d}$ such that

[TABLE]

Using this fact and arguing as in the proof of Lemma 2.1, we have that $Q_{C}$ is compact as a subset of $\Gamma$, that is, when it is endowed with the topology of uniform convergence.

Given $\gamma_{i}\in\mathcal{A}$, let us set $m_{i}:=\gamma_{i}\sharp m_{0}\in\mathcal{P}(\Gamma)$. By (2.1), for any profile of strategies $(\gamma_{i})_{i=1}^{N}\in\mathcal{A}^{N}$, the cost for player $i$ is given by

[TABLE]

This expression for the cost motivates a relaxation of the game with strategies in $\mathcal{A}$ by considering strategies taking values in $\mathcal{P}(\Gamma)$. Let us define the set $\mathcal{A}_{rel}$ of randomized distributed open-loop strategies by

[TABLE]

where, for each $t\in[0,T]$, $e_{t}\!:\Gamma\rightarrow\mathbb{R}^{d}$ is defined by $e_{t}(\gamma):=\gamma(t)$.

Lemma 2.2**.**

The set $\mathcal{A}_{rel}$ is convex and compact as a subset of $\mathcal{P}(\Gamma)$.

Proof.

Convexity follows directly from the definition. On the other hand, since $\mathcal{A}_{rel}\subseteq\mathcal{P}(Q_{C})$ and $\mathcal{P}(Q_{C})$ is compact as a subset of $\mathcal{P}(\Gamma)$ (because $Q_{C}$ is a compact subset of $\Gamma$), it suffices to check that $\mathcal{A}_{rel}$ is closed in $\mathcal{P}(\Gamma)$, but this follows directly from [3, Proposition 5.1.8] and the fact that $Q_{C}$ is closed.

∎

Remark 2.1**.**

For later use, note that if $m\in\mathcal{A}_{rel}$, then $[0,T]\ni t\to e_{t}\sharp m\in\mathcal{P}_{1}(\mathbb{R}^{d})$ is well-defined and, by (2.2), belongs to $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$. Moreover, since $\mbox{{\rm supp}}(m)\subseteq Q_{C}$, we easily check that there exists $C^{\prime}>0$, independent of $m\in\mathcal{A}_{rel}$, such that

[TABLE]

Therefore, by (2.11), [3, Proposition 7.1.5] and Lemma 2.2, the set $\{[0,T]\ni t\to e_{t}\sharp m\in\mathcal{P}_{1}(\mathbb{R}^{d})\;|\;m\in\mathcal{A}_{rel}\}$ is compact in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$.

Motivated by (2.12), we introduce the new relaxed game which has $\mathcal{A}_{rel}$ as set of strategies for each player and, given a strategy profile $(m_{j})_{j=1}^{N}\subseteq\mathcal{A}_{rel}^{N}$, the cost for player $i$ is given by

[TABLE]

Note that the this game is still symmetric. In this framework, a profile of strategies $(m_{1,N},\ldots,m_{N,N})\in\mathcal{A}_{rel}^{N}$ is called a Nash equilibrium in randomized distributed open-loop strategies if

[TABLE]

Proposition 2.1**.**

Under assumption (A1), there exists at least one Nash equilibrium in randomized distributed open-loop strategies having the form $(m_{N},\ldots,m_{N})\in\mathcal{A}_{rel}^{N}$.

Proof.

It suffices to show the existence of a fixed point of the following set-valued map

[TABLE]

First note that, as a consequence of [3, Lemma 5.1.7], for every $\mu\in\mathcal{A}_{rel}$ the map $\mathcal{A}_{rel}\ni\mu^{\prime}\to J_{rel}^{N}(\mu^{\prime},\mu,\ldots,\mu)\in\mathbb{R}$ is lower semicontinuous. Thus, the compactness of $\mathcal{A}_{rel}$ yields that $\SS^{N}(\mu)\neq\emptyset$ for all $\mu\in\mathcal{A}_{rel}$. A similar argument implies that $\SS^{N}(\mu)$ is closed for all $\mu\in\mathcal{A}_{rel}$. Notice also that $\SS^{N}(\mu)$ is convex for all $\mu\in\mathcal{A}_{rel}$. Let us show that $\SS^{N}$ is upper semicontinuous. Consider a closed set $M\subseteq\mathcal{A}_{rel}$ and a sequence $(\mu_{n})$ in $\{\mu\in\mathcal{A}_{rel}\;|\;\SS^{N}(\mu)\cap M\neq\emptyset\}$ converging to some $\bar{\mu}\in\mathcal{A}_{rel}$. By definition, there exists $\nu_{n}\in M$ such that

[TABLE]

Since $M$ is compact (because $\mathcal{A}_{rel}$ is compact), there exists $\bar{\nu}\in M$ such that, up to some subsequence, $\nu_{n}\to\bar{\nu}$ narrowly. By [9, Theorem 3.2] and (A1) for all $\mu^{\prime}\in\mathcal{A}_{rel}$ we have that $J_{rel}^{N}(\mu^{\prime};\mu_{n},\ldots,\mu_{n})$ converges to $J_{rel}^{N}(\mu^{\prime};\bar{\mu},\ldots,\bar{\mu})$. On the other hand, by [22, Theorem 3.23] and (A1), the function $Q_{C}^{N}\ni(\gamma_{1},\ldots,\gamma_{N})\to j^{N}(\gamma_{1};(\gamma_{j})_{j\geq 2})\in\mathbb{R}$ is lower semi-continuous, which, using [3, Lemma 5.1.7] again, implies that $J_{rel}^{N}(\bar{\nu};\bar{\mu},\ldots,\bar{\mu})\leq\liminf_{n\to\infty}J_{rel}^{N}(\nu_{n};\mu_{n},\ldots,\mu_{n})$. Therefore, passing to the limit in (2.15) we obtain that

[TABLE]

i.e. $\bar{\nu}\in\SS^{N}(\bar{\mu})\cap M$, which implies the closedness $\{\mu\in\mathcal{A}_{rel}\;|\;\SS^{N}(\mu)\cap M\neq\emptyset\}$ and the upper semicontinuity of $\SS^{N}$. Using the properties above, the existence of a fixed point for $\SS^{N}$ follows from the Kakutani-Fan-Glicksberg fixed-point theorem (see e.g. [1, Corollary 17.55]). ∎

Corollary 2.1**.**

Let $(m_{N},\ldots,m_{N})\in\mathcal{A}_{rel}^{N}$ be a Nash equilibrium for the game defined by the cost $J_{rel}^{N}$ and the set of strategies $\mathcal{A}_{rel}$. Define $\hat{\ell}:\mathbb{R}^{d}\times\mathbb{R}^{d}\times[0,T]\to\mathbb{R}$ and $\hat{\Phi}:\mathbb{R}^{d}\to\mathbb{R}$ by

[TABLE]

and assume that for $m_{0}$-almost every $x\in\mathbb{R}^{d}$ the optimization problem

[TABLE]

admits a unique solution. Then, there exists $\gamma_{N}\in\mathcal{A}$ such that $m_{N}=\gamma_{N}\sharp m_{0}$. Moreover, $\gamma_{N}$ is $m_{0}$-uniquely determined.

In particular, $(\gamma_{N},\ldots,\gamma_{N})\in\mathcal{A}^{N}$ is a Nash equilibrium in distributed open-loop strategies.

Proof.

Define $\widehat{j}:Q_{C}\to\mathbb{R}$ by

[TABLE]

By definition of Nash equilibrium and Fubini’s theorem, for all $m\in\mathcal{A}_{rel}$ we have that

[TABLE]

By the desintegration theorem (see e.g. [3, Theorem 5.3.1]), there exists a Borel family $\{m_{N}^{x}\;|\;x\in\mathbb{R}^{d}\}\subseteq\mathcal{P}(Q_{C})$, such that $m_{N}^{x}\left(\{\gamma\in Q_{C}\;|\;\gamma(0)=x\}\right)=1$, for $m_{0}$-almost every $x\in\mathbb{R}^{d}$, and

[TABLE]

Define the set-valued function $\hat{\SS}^{N}:\mathbb{R}^{d}\to 2^{\Gamma}$ by

[TABLE]

Arguing as in the proof of Lemma 2.1, we have the existence of a Borel measurable selection $\mathbb{R}^{d}\ni x\mapsto\gamma_{N}^{x}\in\hat{\SS}^{N}(x)$, which, by assumption, is $m_{0}$-uniquely determined. Moreover, by definition, $\gamma_{N}\in\mathcal{A}$. Now, let us define $\hat{m}:=\gamma_{N}\sharp m_{0}\in\mathcal{A}_{rel}$. Then, by (2.19), taking $m=\hat{m}$, and (2.20), we have that

[TABLE]

Since the integrand in the expression above is non-negative, by definition of $\gamma_{N}^{x}$, we deduce that $\widehat{j}(\gamma_{N}^{x})=\int_{Q_{C}}\widehat{j}(\gamma){\rm d}m_{N}^{x}(\gamma)$ for $m_{0}$-a.e. $x\in\mathbb{R}^{d}$, and, hence, $\widehat{j}(\gamma)=\widehat{j}(\gamma_{N}^{x})$ for $m_{N}^{x}$-a.e. $\gamma\in Q_{C}$. Since, by assumption, $\hat{\SS}^{N}(x)=\{\gamma_{N}^{x}\}$ for $m_{0}$-a.e. $x\in\mathbb{R}^{d}$, we deduce that $m^{x}_{N}=\delta_{\gamma_{N}^{x}}$ for $m_{0}$-a.e. $x\in\mathbb{R}^{d}$, hence $m_{N}=\gamma_{N}\sharp m_{0}$. The result follows. ∎

Remark 2.2**.**

An example of application of Corollary 2.1 is provided in Section 4.2 below.

3. Convergence to a mean field game equilibrium

In this section, we study the limit behavior, as $N\to\infty$, of symmetric Nash equilibria in randomized distributed open-loop strategies. The existence of such Nash equilibria is ensured by Proposition 2.1. We begin by defining the limit object, i.e. the MFG equilibrium. Then we will prove that any cluster point of the sequence $(m_{N})$ is a MFG equilibrium.

Let us define $J:W^{1,q}([0,T];\mathbb{R}^{d})\times\mathcal{P}_{1}(\Gamma)\to\mathbb{R}$ by

[TABLE]

It is straightforward to check that $m\in\mathcal{P}_{1}(\Gamma)$ implies that $e_{t}\sharp m\in\mathcal{P}_{1}(\mathbb{R}^{d})$ for all $t\in[0,T]$, which implies that $J$ is well-defined.

Following the terminology in [35], we consider next the notion of Lagrangian MFG equilibrium (see e.g. [16, 8, 17, 10]).

Definition 3.1**.**

We say that $m_{\ast}\in\mathcal{P}_{1}(\Gamma)$ is a Lagrangian MFG equilibrium if $e_{0}\sharp m_{\ast}=m_{0}$ and

[TABLE]

Reasoning as in the proof of Lemma 2.1, assumption (A1) implies that if $m_{\ast}$ is a Lagrangian MFG equilibrium, then $\mbox{supp}(m_{\ast})\subseteq Q_{C}$ and, hence, $m_{\ast}\in\mathcal{A}_{rel}$. Thus, $m_{\ast}$ is a Lagrangian MFG if and only if $m_{\ast}\in\mathcal{A}_{rel}$, $e_{0}\sharp m_{\ast}=m_{0}$ and

[TABLE]

We still denote by $J$ the restriction of $J$ to $Q_{C}\times\mathcal{A}_{rel}$ and recall that $Q_{C}$, endowed with the topology of uniform convergence, is a compact set. For later use, let us state the following simple results.

Lemma 3.1**.**

*Suppose that (A1) holds. Then, the following assertions hold true:

(i) The relative topology on $\mathcal{A}_{rel}$, as a subset of $\mathcal{P}_{1}(\Gamma)$, coincides with the topology induced by the narrow convergence.

(ii) The function $J$ is lower semicontinuous in $Q_{C}\times\mathcal{A}_{rel}$.

(iii) For all $\gamma\in Q_{C}$, the function $J(\gamma,\cdot)$ is bounded, uniformly in $\gamma$, and continuous in $\mathcal{A}_{rel}$.*

Proof.

Assertion (i) follows from the fact that both topologies coincide on $\mathcal{P}(Q_{C})$, because $Q_{C}$ is a compact subset of $\Gamma$. In particular, if $m_{n}\to m$ narrowly in $\mathcal{A}_{rel}$, then for all $t\in[0,T]$ we have that

[TABLE]

Thus, assertion (ii) follows from (3.3), (A1) and the proof of [22, Theorem 3.23]. Assertion (iii) follows directly from (2.3), (3.3) and dominated convergence. ∎

If $m_{\ast}$ is a Lagrangian MFG equilibrium, we will denote by $\{m_{\ast}^{x}\;|\;x\in\mathbb{R}^{d}\}$ the $m_{0}$-uniquely determined Borel family of probability measures on $Q_{C}$ satisfying that $m_{\ast}^{x}(Q_{C}\setminus\mathcal{A}(x))=0$ and ${\rm d}m_{\ast}(\gamma)={\rm d}m_{\ast}^{x}(\gamma)\otimes{\rm d}m_{0}(x)$. The existence of such a family is ensured by the disintegration theorem (see e.g. [3, Theorem 5.3.1]). We have the following equivalent characterization of a Lagrangian MFG equilibrium.

Proposition 3.1**.**

*The following assertions are equivalent:

(i) The measure $m_{\ast}\in\mathcal{A}_{rel}$ is a Lagrangian MFG equilibrium.

(ii) For $m_{0}$-a.e. $x\in\mathbb{R}^{d}$ we have that*

[TABLE]

(iii)* The following inequality holds true:*

[TABLE]

Proof.

Let us prove the equivalence between (i) and (ii). Let $m_{\ast}\in\mathcal{A}_{rel}$ be a Lagrangian MFG equilibrium. If (3.4) does not hold, there exists $A\in\mathcal{B}(\mathbb{R}^{d})$, with $m_{0}(A)>0$, such that $m_{\ast}^{x}(\SS(x)^{c})>0$ for all $x\in A$. Define the set $E:=\left\{\gamma\in Q_{C}\;|\;\gamma(0)\in A,\;\gamma\in\SS(\gamma(0))^{c}\right\}=e_{0}^{-1}(A)\cap\left\{\gamma\in Q_{C}\;|\;\gamma\in\SS(\gamma(0))^{c}\right\}$. Arguing as in the proof of the claim in Lemma 2.1, the set $\left\{\gamma\in Q_{C}\;|\;\gamma\in\SS(\gamma(0))\right\}$ is closed in $Q_{C}$, which implies that $E\in\mathcal{B}(Q_{C})$. Since $m_{\ast}(E)=\int_{A}m_{\ast}^{x}(\SS(x)^{c}){\rm d}m_{0}(x)>0$, we obtain a contradiction with (3.2). Conversely, suppose that (ii) holds and that $m_{\ast}$ is not a Lagrangian MFG equilibrium. Then, defining $E^{\prime}:=\{\gamma\in Q_{C}\;|\;\gamma\in\SS(\gamma(0))^{c}\}$, which is an open set and, hence, belongs to $\mathcal{B}(Q_{C})$, we have that

[TABLE]

which is impossible because (3.4) implies that the r.h.s. above is equal to [math].

Let us now prove that (ii) $\Leftrightarrow$ (iii). Notice that

[TABLE]

Analogously, given $m\in\mathcal{A}_{rel}$, we disintegrate it w.r.t. $m_{0}$ and write ${\rm d}m(\gamma)={\rm d}m^{x}(\gamma)\otimes{\rm d}m_{0}(x)$, where $\{m^{x}\;|\;x\in\mathbb{R}^{d}\}$ is a $m_{0}$-uniquely determined Borel family of probability measures on $Q_{C}$ such that $m^{x}(Q_{C}\setminus\mathcal{A}(x))=0$ for $m_{0}$-a.e. $x\in\mathbb{R}^{d}$. Thus,

[TABLE]

If (ii) holds, then for $m_{0}$-a.e. $x\in\mathbb{R}^{d}$ and $m_{\ast}^{x}$-a.e. $\gamma\in Q_{C}$ we have

[TABLE]

Integrating both sides of (3.8), first with respect to ${\rm d}m^{x}(\gamma^{\prime})$ and then with respect to ${\rm d}m_{\ast}^{x}(\gamma)$, and using (3.6)-(3.7) we obtain (3.5). Conversely, using the notations introduced above, suppose that (3.5) holds and let $\hat{\gamma}\in\mathcal{A}$ be a Borel measurable selection of $\SS$ (the existence of such selection can be justified arguing exactly as in the proof of Lemma 2.1). Then, taking the measure $m\in\mathcal{A}_{rel}$ defined by ${\rm d}m(\gamma)={\rm d}\delta_{\hat{\gamma}^{x}}(\gamma)\otimes{\rm d}m_{0}(x)$ in (3.5) and using (3.6), we deduce that

[TABLE]

and, hence, $m^{x}_{\ast}$-almost every $\gamma$ belongs to $\SS(x)$. The conclusion follows. ∎

Now, consider the symmetric $N$-player game defined in Section 2, with randomized distributed open-loop strategies, and let $(m_{N},\ldots,m_{N})\in\mathcal{A}_{rel}^{N}$ be a symmetric equilibrium. Our main result in this section, stated in the next theorem, shows that any limit point of this sequence is a Lagrangian MFG equilibrium.

Theorem 3.1**.**

Suppose that (A1) holds. Then any limit point $m_{\ast}$ of $(m_{N})_{N\in\mathbb{N}}$ (there exist at least one) is a Lagrangian MFG equilibrium. Moreover, if $(m_{N_{k}})_{k\in\mathbb{N}}$ is a subsequence of $(m_{N})_{N\in\mathbb{N}}$ converging to $m_{\ast}$, then $\sup_{t\in[0,T]}d_{1}\left(e_{t}\sharp m_{N_{k}},e_{t}\sharp m_{\ast}\right)\to 0$ as $k\to\infty$.

Proof.

We only prove the first assertion, since, having this result, the second one follows directly from Remark 2.1. For $N\geq 2$, let $Y^{N}_{1},\ldots,Y^{N}_{N}$ be independent and identically distributed (i.i.d.) $Q_{C}$-valued random variables with common distribution $m_{N}$ defined on some probability space $(\Omega_{N},\mathcal{F}_{N},\mathbf{P}_{N})$. We denote by $\mathbf{E}_{N}$ the expectation with respect to $\mathbf{P}_{N}$. For $i\in\{1,\ldots,N\}$, let $\mu^{N,i}$ denote the (random) empirical measure of $Y^{N}_{1},\ldots,Y^{N}_{N}$ excluding $Y^{N}_{i}$, that is,

[TABLE]

Notice that $Y^{N}_{i}$ and $\mu^{N,i}$ are independent for every $i$, while $\mu^{N,1},\ldots,\mu^{N,N}$ are identically distributed (not independent in general) with common distribution depending on $N$. Moreover, for every $i\in\{1,\ldots,N\}$,

[TABLE]

Let $(m_{N_{k}})_{k\in\mathbb{N}}$ be a converging subsequence of $(m_{N})_{N\in\mathbb{N}}$ with limit $m_{\ast}$ for some $m_{\ast}\in\mathcal{A}_{rel}$. The existence of such a subsequence follows from the compactness of $\mathcal{A}_{rel}$. Let us prove that $(\mu^{N_{k},1})_{k\in\mathbb{N}}$ converges in distribution to the deterministic limit $m_{\ast}$. Let $\mathcal{T}\subset C_{b}(Q_{C})$ be countable and measure determining (or separating). Thus, $\mathcal{T}$ is a countable collection of bounded continuous functions on $Q_{C}$ such that two probability measures $\nu,\tilde{\nu}\in\mathcal{P}(Q_{C})$ are equal whenever $\int\psi\,d\nu=\int\psi\,d\tilde{\nu}$ for all $\psi\in\mathcal{T}$. Observe that $\mathcal{T}$ can be chosen countable since $Q_{C}$ is a Polish space, hence separable, under the supremum norm topology.

For $\psi\in\mathcal{T}$ set

[TABLE]

By construction and symmetry, for every $i\in\{1,\ldots,N\}$,

[TABLE]

Since $Y^{N}_{1},\ldots,Y^{N}_{N}$ are independent and the functions in $\mathcal{T}$ are bounded, it follows that

[TABLE]

Moreover, considering the subsequence $(m_{N_{k}})_{k\in\mathbb{N}}$, we have that

[TABLE]

Since $\mathcal{P}(Q_{C})$ is compact, the set $\mathcal{P}(\mathcal{P}(Q_{C}))$ is also compact. Thus, there exists a sub-subsequence $(N_{k_{l}})_{l\in\mathbb{N}}$ such that $(\mu^{N_{k_{l}},1})_{l\in\mathbb{N}}$ converges in distribution to $\mu$ for some $\mathcal{P}(Q_{C})$-valued random variable $\mu$ defined on some probability space $(\Omega,\mathcal{F},\mathbf{P})$. We denote by $\mathbf{E}$ the expectation under $\mathbf{P}$. Then, by the continuous mapping theorem (see e.g. [9, Theorem 5.1]) and (3.11), we have that

[TABLE]

On the other hand, thanks to (3.10), $\lim_{l\to\infty}v^{N_{k_{l}}}_{\psi}=0$. It follows that for every $\psi\in\mathcal{T}$,

[TABLE]

Since $\mathcal{T}$ is countable, we have the existence of $A\in\mathcal{F}$ such that $\mathbf{P}(A)=1$ and for every $\omega\in A$,

[TABLE]

Since $\mathcal{T}$ is measure determining, we find that

[TABLE]

As we can always choose converging (sub-)subsequences, we deduce that $(\mu^{N_{k},1})_{k\in\mathbb{N}}$ converges in distribution to $m_{\ast}$. By independence of $Y^{N_{k}}_{1}$ and $\mu^{N_{k},1}$, and using [9, Theorem 3.2], we find that

[TABLE]

This implies, thanks to (3.9), Lemma 3.1(ii) and [3, Lemma 5.1.7], that

[TABLE]

Let $m\in\mathcal{A}_{rel}$. By Lemma 3.1(iii) and dominated convergence, it follows that

[TABLE]

Passing to the limit in the Nash equilibrium inequality

[TABLE]

and using (3.12)-(3.13), we obtain that (3.5) holds. The result now follows from Proposition 3.1. ∎

Remark 3.1**.**

In particular, under assumption (A1), Theorem 3.1 ensures the existence of at least one Lagrangian MFG equilibrium.

4. The first order Mean Field Game system

In this section, we first discuss, following [16], the relation between the notion of Lagrangian MFG equilibrium and the first order PDE system introduced by Lasry and Lions in [34, Section 2.5] for some particular data. Next, in Section 4.2, we consider symmetric randomized distributed open-loop Nash equilibria for the $N$-player game, which, thanks to Corollary 2.1 and Assumption (A2) below, can be identified with symmetric distributed open-loop Nash equilibria (non-randomized). Arguing as in the MFG limit, we connect these equilibria with a first order PDE system which is similar to the one appearing in the limit case. Consequently, analytic techniques can also be used in order study the limit behavior of these equilibria as the number of players tends to infinity.

Let $L:\mathbb{R}^{d}\times\mathbb{R}^{d}\to\mathbb{R}$, $f$, $g:\mathbb{R}^{d}\times\mathcal{P}_{1}(\mathbb{R}^{d})\to\mathbb{R}$ and $m_{0}\in\mathcal{P}_{1}(\mathbb{R}^{d})$ satisfying that

(A2)(i) The function $L$ belongs to $C^{2}(\mathbb{R}^{d}\times\mathbb{R}^{d})$, is bounded from below and

(i.1) there exist $C_{L}>0$, $\overline{L}>0$ such that

[TABLE]

(i.2) There exist $c_{L},\;c_{L}^{\prime}>0$ such that

[TABLE]

(ii) The functions $f$ and $g$ are continuous. Moreover, for every $m\in\mathcal{P}_{1}(\mathbb{R}^{d})$ the functions $f(\cdot,m)$ and $g(\cdot,m)$ belong to $C^{2}(\mathbb{R}^{d})$ and there exists a constant $C_{f,g}>0$ such that

[TABLE]

where, for $h=f$, $g$, we have set

[TABLE]

(iii) The measure $m_{0}$ is absolutely continuous w.r.t. to the Lebesgue measure $\mathcal{L}^{d}$, with density still denoted by $m_{0}$, and has a compact support.

A typical example of function $L$ satisfying (A2)(i) is given by $\mathbb{R}^{d}\times\mathbb{R}^{d}\ni(\alpha,x)\to L(\alpha,x):=b_{1}(x)|\alpha|^{2}+b_{2}(x)$, where, for $i=1$, $2$, $b_{i}\in C^{2}(\mathbb{R}^{d})$, $b_{i}$ is Lipschitz, and there exist constants $\underline{b}_{1}>0$, $\underline{b}_{2}\in\mathbb{R}$ and $\overline{b}_{i}>0$ such that $\underline{b}_{i}\leq b_{i}\leq\overline{b}_{i}$.

Remark 4.1**.**

(i)* Assumption (A2)(i) above implies the convexity of $L(\cdot,x)$ and the existence of $\underline{L}$, $C_{L}^{\prime}$, $c^{{}^{\prime\prime}}_{L}>0$ such that*

[TABLE]

(ii)* For $x\in\mathbb{R}^{d}$ let us denote by $H(\cdot,x)$ the convex conjugate of $L(\cdot,x)$. The bound in ${\rm(i.1)}$ and the first bound in (4.3) imply the exitence of constants $\underline{H}$, $\overline{H}$ and $C_{H}>0$ such that*

[TABLE]

(iii)* By the first estimate in (4.2) we have that $\partial_{\xi}H(\xi,x)$ is characterized as the unique solution $\alpha(\xi,x)$ of the optimization problem $\max_{\alpha\in\mathbb{R}^{d}}\{\xi\cdot\alpha-L(\alpha,x)\}$. As a consequence of this fact and the first relation in (4.3), we obtain the existence of $c_{H}>0$ such that*

[TABLE]

Moreover, from the convexity of $L(x,\cdot)$, for all $\xi$, $x\in\mathbb{R}^{d}$ we have that $\partial_{\xi}H(\xi,x)$ is the unique solution to

[TABLE]

Using the relation above, the regularity $L\in C^{2}(\mathbb{R}^{d}\times\mathbb{R}^{d})$, the first estimate in (4.2) and the implicit function theorem, we obtain that $\mathbb{R}^{d}\times\mathbb{R}^{d}\ni(\xi,x)\to\partial_{\xi}H(\xi,x)\in\mathbb{R}^{d}$ belongs to $C^{1}(\mathbb{R}^{d}\times\mathbb{R}^{d};\mathbb{R}^{d})$. Using this fact, we get that $\mathbb{R}^{d}\times\mathbb{R}^{d}\ni(\xi,x)\mapsto\partial_{x}H(\xi,x)=-\partial_{x}L(\partial_{\xi}H(\xi,x),x)\in\mathbb{R}^{d}$ is also of class $C^{1}$. As a consequence $H$ is of class $C^{2}$.

Let us define $\ell:\mathbb{R}^{d}\times\mathbb{R}^{d}\times\mathcal{P}_{1}(\mathbb{R}^{d})\to\mathbb{R}$ and $\Phi:\mathbb{R}^{d}\times\mathcal{P}_{1}(\mathbb{R}^{d})\to\mathbb{R}$ by

[TABLE]

Clearly, Assumption (A2) implies that $\ell$, $\Phi$ and $m_{0}$ satisfy (A1).

4.1. Lagrangian MFG equilibria and the MFG PDE system

As pointed out in [16], under (A2) the existence of a Lagrangian equilibrium for the MFG problem defined by $\ell$, $\Phi$ and $m_{0}$, is equivalent the existence of a solution $(u,\rho)$ of the following PDE system, which was first introduced in [34],

[TABLE]

In the system above, $u:\mathbb{R}^{d}\times[0,T]\to\mathbb{R}$ is a solution to the first equation, with the associated terminal condition, if it is globally Lipschitz, locally semi-concave with respect to its first argument (see [12, Section 2]), uniformly in $t\in[0,T]$, and the equation holds in the viscosity sense (see e.g. [4, Chapter III, Section 3]). In $(MFG)$, $\nabla u$ denotes a Borel measurable selection of the set-valued map

[TABLE]

The existence of such measurable selection follows from the fact that the above set-valued map has a closed graph (thanks to the semi-concavity property of $u$, see e.g. [12, Proposition 3.3.4]). Moreover, since $u$ is Lipschitz, $\nabla u(x,t)$ is uniformly bounded in $(x,t)\in\mathbb{R}^{d}\times[0,T]$.

The function $\rho:[0,T]\to\mathcal{P}(\mathbb{R}^{d})$ is a solution to the second equation, with the associated initial condition, if $\rho\in C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$, and the equation is satisfied in the sense of distributions, i.e. for all $\phi\in C^{\infty}(\mathbb{R}^{d})$ with compact support we have

[TABLE]

Note that by the previous considerations, the second term in the right-hand-side of (4.8) is well-defined.

Any pair $(u,\rho)$ satisfying $(MFG)$ is called an equilibrium of the first order MFG problem.

For the sake of completeness, let us provide the main arguments that justify the equivalence between both notions of equilibria. Let $m_{\ast}\in\mathcal{A}_{rel}$ be a Lagrangian MFG equilibrium. Let us define $\rho:[0,T]\to\mathcal{P}(\mathbb{R}^{d})$ by $\rho(t):=e_{t}\sharp m_{\ast}$ for all $t\in[0,T]$, and the value function $u:\mathbb{R}^{d}\times[0,T]\to\mathbb{R}$ by

[TABLE]

for all $(x,t)\in\mathbb{R}^{d}\times[0,T]$. Since $m_{\ast}\in\mathcal{A}_{rel}$, we have that $\rho\in C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$. Using this fact, assumption (A2) and [5, Proposition 1.1 and Remark 1.1], we obtain that $u$ is a viscosity solution of

[TABLE]

Moreover, by [21, Theorem 2.1], the Hamilton-Jacobi-Bellman equation above admits a comparison principle, which implies that $u$ is its unique viscosity solution. We will need the following result, whose proof follows from standard arguments.

Lemma 4.1**.**

*Under (A2) we have:

(i) For every $(x,t)\in\mathbb{R}^{d}\times[0,T]$ the set $\SS(x,t)$ of paths $\gamma_{\ast}\in W^{1,2}([t,T];\mathbb{R}^{d})$ such that $\gamma_{\ast}(t)=x$ and*

[TABLE]

is non-empty. Moreover, there exists a constant $C>0$, independent of $(x,t)$, such that

[TABLE]

(ii)* The value function $u$ is globally Lipschitz.

(iii) The value function $u$ is locally semi-concave w.r.t. to the space variable, uniformly in $t\in[0,T]$. More precisely, for any compact set $K\subseteq\mathbb{R}^{d}$ there exists a constant $C_{K}$, independent of $t$, such that for every $\lambda\in[0,1]$, $x$, $y\in K$, such that the segment $[x,y]$ is contained in $K$, the following inequality holds*

[TABLE]

Proof.

The proof being standard, we only sketch the main ideas. The fact that $\SS(x,t)$ is non-empty follows directly from (A2) and the direct method in the Calculus of Variations. Moreover, by (A2) and arguing as in the proof of Lemma 2.1 we obtain the existence of $c>0$, independent of $(x,t,\rho)$, such that

[TABLE]

Using this fact, the Euler-Lagrange equation associated to any element $\gamma_{\ast}\in\SS(x,t)$, the second estimate in (4.2) and arguing as in the proof of [12, Theorem 6.2.5], we easily obtain (4.12), which proves assertion (i). In order to prove (ii), notice that (i) implies that the value function can also be written as

[TABLE]

where $\hat{A}_{C}:=\left\{\alpha\in L^{\infty}([0,T];\mathbb{R}^{d})\;|\;\|\alpha\|_{L^{\infty}}\leq C\right\}.$ Using the estimate $|\inf_{\alpha\in\hat{A}_{C}}A(\alpha)-\inf_{\alpha\in\hat{A}_{C}}B|\leq\sup_{\alpha\in\hat{A}_{C}}|A(\alpha)-B(\alpha)|$ for any functions $A$, $B:L^{\infty}([0,T];\mathbb{R}^{d})\to\mathbb{R}$, expression (4.7), the uniform Lipschitz property for $f$ and $g$ in (A2)(ii), and the second estimate in (4.2), we easily obtain that $u(\cdot,t)$ is globally Lipschitz, with a Lipschitz constant which is independent of $t\in[0,T]$. Similarly, using (4.15) and the estimate (4.1), we get that $u(x,\cdot)$ is globally Lipschitz, with a Lipschitz constant which is independent of $x\in\mathbb{R}^{d}$. Assertion (ii) follows. Finally, assertion (iii) follows directly from [12, Theorem 6.4.1]. ∎

Now, let us consider the set-valued map

[TABLE]

Since Lemma 4.1(ii) implies that $u(\cdot,0)$ is a.e. differentiable, [12, Corollary 6.4.10] yields that for a.e. $x\in\mathbb{R}^{d}$ we have $\SS(x)=\{\tilde{\gamma}^{x}\}$ for some $\tilde{\gamma}^{x}\in\mathcal{A}(x)$. Now, as in the proof of Proposition 3.1, let $\gamma_{\ast}\in\mathcal{A}$ be a Borel measurable selection of $\SS$. Then, for a.e. $x\in\mathbb{R}^{d}$ we have that $\gamma_{\ast}^{x}=\tilde{\gamma}^{x}$. Thus, Proposition 3.1(ii) yields $m_{\ast}^{x}=\delta_{\gamma_{\ast}^{x}}$ for a.e. $x\in\mbox{supp}(m_{0})$ and, hence, $m_{\ast}=\gamma_{\ast}\sharp m_{0}$. In particular, $\rho(t)=\gamma_{\ast}^{(\cdot)}(t)\sharp m_{0}$ for all $t\in[0,T]$.

On the other hand, by [12, Theorem 6.4.9, Theorem 6.3.3 and Theorem 6.4.8], for a.e. $x\in\mathbb{R}^{d}$, we have

[TABLE]

where we underline that $u$ is differentiable w.r.t. to its first argument at the point $(\gamma_{\ast}^{x}(t),t)$ if $t\in(0,T)$ (see [12, Theorem 6.4.7]). Denoting by still by $\nabla u$ a measurable selection of $(x,t)\mapsto D_{x}^{+}u(x,t)$, for every $\phi\in C^{\infty}(\mathbb{R}^{d})$ with compact support and $t\in[0,T]$, we have

[TABLE]

which implies that $\rho$ satisfies (4.8) and, hence, the couple $(u,\rho)$ solves $(MFG)$. Notice that under (A2) a Lagrangian MFG equilibrium $m_{\ast}$ exists (see Remark 3.1) and, hence, the previous arguments show, in particular, the existence of at least one solution $(u,\rho)$ to $(MFG)$.

Conversely, if $(u,\rho)$ solves $(MFG)$, then the first equation therein implies that $u$ and $\rho$ are still related by (4.9). By the second equation in $(MFG)$ and [3, Theorem 8.2.1], there exists a probability measure $m_{\ast}\in\mathcal{P}(\Gamma)$ such that $\rho(t)=e_{t}\sharp m_{\ast}$ for all $t\in[0,T]$ and, considering the disintegration ${\rm d}m_{\ast}(\gamma)={\rm d}m_{\ast}^{x}(\gamma)\otimes{\rm d}m_{0}(x)$, for a.e. $x\in\mbox{supp}(m_{0})$ the support of the measure $m_{\ast}^{x}$ is contained in the set of solutions of (4.16). By Lemma 4.1 and arguing as in the proof of [13, Lemma 4.11], we have that every solution to (4.16) solves the optimization problem in the r.h.s. of (4.9) with $t=0$. Thus, by Proposition 3.1(ii) we obtain that $m_{\ast}$ is a Lagrangian MFG equilibrium. Notice also that $\SS(x)$ being a singleton for a.e. $x\in\mathbb{R}^{d}$, the previous argument shows, in particular, that $[0,T]\ni t\mapsto\gamma_{\ast}^{(\cdot)}(t)\sharp m_{0}\in\mathcal{P}_{1}(\mathbb{R}^{d})$ is the unique solution in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$ of the continuity equation

[TABLE]

In addition to the relation between Lagrangian MFG equilibria and the solutions of $(MFG)$, assumption (A2) has also consequences on the regularity of the time marginals $\{\rho(t)\;|\;t\in[0,T]\}$ as the following result shows.

Proposition 4.1**.**

*In addition to (A2)(iii), assume that $m_{0}\in L^{p}(\mathbb{R}^{d})$ for some $p\in(1,+\infty]$ and let $(u,\rho)$ be a solution to $(MFG)$. Then, the following assertions hold true:

(i) There exists $c_{1}>0$, independent of $t\in[0,T]$, such that $\mathop{\rm supp}(\rho(t))\subseteq B(0,c_{1})$ for all $t\in[0,T]$.

(ii) For all $t\in[0,T]$ the measure $\rho(t)$ is absolutely continuous w.r.t. the Lebesgue measure. Moreover, the density of $\rho(t)$, that we will still denote by $\rho(t)$, belongs to $L^{p}(\mathbb{R}^{d})$ and there exists a constant $c_{2}>0$, independent of $p\in(1,+\infty]$ and $t\in[0,T]$, such that*

[TABLE]

Proof.

Assertion (i) follows directly from the formula $\rho(t)=\gamma_{\ast}^{(\cdot)}(t)\sharp m_{0}$, where $\gamma_{\ast}^{x}\in\SS(x)$ for all $x\in\mathbb{R}^{d}$, Lemma 4.1(i) and the fact that $\mathop{\rm supp}(m_{0})$ is compact. In order to prove (ii), let $\beta\in C^{\infty}(\mathbb{R}^{d})$, non-negative, with support contained in the unit ball and such that $\int_{\mathbb{R}^{d}}\beta(x){\rm d}x=1$. For $\varepsilon>0$, let us define $\beta_{\varepsilon}(x):=\varepsilon^{-d}\beta(x/\varepsilon)$, $u_{\varepsilon}(x,t):=\left(\beta_{\varepsilon}\ast u(\cdot,t)\right)(x)$ and consider the equation

[TABLE]

For every $x\in\mathbb{R}^{d}$, let us define $\gamma_{\varepsilon}^{x}\in C^{1}([0,T];\mathbb{R}^{d})$ as the unique solution to

[TABLE]

By [3, Proposition 8.1.8], equation (4.19) admits a unique solution in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$, which is given by $\rho_{\varepsilon}(t):=\gamma_{\varepsilon}^{(\cdot)}(t)\sharp m_{0}$ for all $t\in[0,T]$. Moreover, by a standard change of variable argument (see e.g. [2, Section 2]), for every $t\in[0,T]$ we have that $\rho_{\varepsilon}(t)$ is absolutely continuous, with density given by

[TABLE]

where, for each $y\in\mathbb{R}^{d}$, $Y(y,\cdot)$ is defined as the unique solution to

[TABLE]

with $\mathbb{R}^{d}\times[0,T]\ni(x,t)\mapsto L_{\varepsilon}(x,t)\in\mathbb{R}^{d\times d}$ being given by

[TABLE]

Let us assume that $p\in(1,+\infty)$. By a change of variable again, we obtain that

[TABLE]

Now, for all $x\in\mathbb{R}^{d}$ and $t\in[0,T]$, we have (see [2, Section 2, estimate $(2.4)$])

[TABLE]

where $[a]_{-}:=\max\{0,-a\}$, $[a]_{+}=a+[a]_{-}$ and for any $a\in\mathbb{R}$, and $\mbox{Tr}\left(L_{\varepsilon}(x,t)\right)$ denotes the trace of the matrix $L_{\varepsilon}(x,t)$. On the other hand, Lemma 4.1(ii), and the definition of $u_{\varepsilon}$, imply that $\nabla u_{\varepsilon}(x,t)$ is bounded, uniformly in $(x,t)\in\mathbb{R}^{d}\times[0,T]$ and $\varepsilon>0$. Moreover, by Lemma 4.1(iii), the compactness of $\mathop{\rm supp}(m_{0})$ and the definition of $u_{\varepsilon}(\cdot,t)$ again, we can assume that $u_{\varepsilon}(\cdot,t)$ is uniformly semiconcave in a bounded open set $\mathcal{O}$ containing $\mathop{\rm supp}(m_{0})$, i.e. $u_{\varepsilon}(\cdot,t)$ satisfies (4.13) for all $x$, $y\in\mathcal{O}$, with $C_{K}$ replaced by $\tilde{c}$, for some $\tilde{c}$ independent of $t$ and $\varepsilon$ small enough. By [12, Proposition 1.1.3], we have that $\partial^{2}_{xx}u_{\varepsilon}(x,t)-\tilde{c}I_{d}$ negative semidefinite for all $(x,t)\in\mathcal{O}\times[0,T]$ and, hence, using that $\partial_{\xi,\xi}H(\xi,x)$ is positive semidefinite for all $\xi\in\mathbb{R}^{d}$ and $x\in\mathbb{R}^{d}$, there exists a constant $\hat{c}>0$, independent of $\varepsilon$ and $t$, such that

[TABLE]

is negative semidefinite for all $x\in\mathcal{O}$. As a consequence, $\mbox{Tr}(L_{\varepsilon}(x,t))$ is bounded from above by a constant which is independent of $\varepsilon>0$ small enough, $x\in\mathcal{O}$, and $t\in[0,T]$. Thus, by (4.23) and taking the power $1/p$ in (4.22), there exists $c_{2}>0$, independent of $\varepsilon$, $t$ and $p$, such that

[TABLE]

The previous estimate shows the existence of $\tilde{\rho}\in L^{\infty}([0,T];L^{p}(\mathbb{R}^{d}))$ and a sequence $(\rho_{\varepsilon_{n}})_{n\in\mathbb{N}}$ such that, as $n\to\infty$, $\varepsilon_{n}\to 0$ and $\rho_{\varepsilon_{n}}\to\tilde{\rho}\in L^{\infty}([0,T];L^{p}(\mathbb{R}^{d}))$ in the weak* topology. By dominated convergence, we have that $\partial_{\xi}H(\nabla u_{\varepsilon_{n}},\cdot)\to\partial_{\xi}H(\nabla u,\cdot)$ in $L^{1}([0,T];L^{s}(\mathbb{R}^{d}))$ for any $s\in[1,+\infty)$. As a consequence, $\tilde{\rho}$ satisfies estimate (4.24) and, passing to the limit in (4.19), we get that the measure $[0,T]\ni t\to\tilde{\rho}(t)\mathcal{L}^{d}\in L^{p}(\mathbb{R}^{d})$ satisfies (4.8). Using that $[0,T]\ni t\mapsto\rho(t)\in\mathcal{P}_{1}(\mathbb{R}^{d})$ is the unique solution to (4.17) in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$, [3, Lemma 8.1.2] implies that $\tilde{\rho}(t)\mathcal{L}^{d}=\rho(t)$ for a.e. $t\in[0,T]$. Thus, for a.e. $t\in[0,T]$, $\rho(t)$ is absolutely continuous w.r.t. to the Lebesgue measure and estimate (4.18) holds for its density. Using this fact and that $\rho\in C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$, the previous statement is valid in the whole time interval $[0,T]$, which proves (ii) when $p<+\infty$. Since $c_{2}$ does not depend on $p$, assertion (ii) for $p=\infty$ follows by taking the limit in (4.18) when $p\to\infty$. ∎

Remark 4.2**.**

Similar regularization techniques have been recently employed in [26], in order to establish $L^{p}$-estimates for the time evolving distributions describing equilibria in optimal-exit MFGs.

4.2. The $N$-player equilibria: associated time marginals and value functions

Let us consider the game with $N$ players defined in Section 2 with $\ell$ and $\Phi$ given by (4.7). Let $(m_{N},\ldots,m_{N})\in\mathcal{A}_{rel}^{N}$ be a symmetric equilibrium in randomized distributed open-loop strategies for the $N$-player game. Note that if for $h=f$, $g$ we define

[TABLE]

we have that $f_{N}$ and $g_{N}$ satisfy the assumptions for $f$ and $g$ in (A2)(ii) (with the same constant $C_{f,g}$). As a consequence of this fact, the results in [12, Chapter 6] and Corollary 2.1 we obtain the existence of $\gamma_{N}\in\mathcal{A}$ such that $m_{N}=\gamma_{N}\sharp m_{0}$, i.e. $(m_{N},\ldots,m_{N})$ can be identified with the non-randomized symmetric equilibrium in distributed open-loop strategies given by $(\gamma_{N},\ldots,\gamma_{N})\in\mathcal{A}^{N}$. Furthermore, setting $\rho_{N}(t):=e_{t}\sharp m_{N}=\gamma_{N}^{(\cdot)}(t)\sharp m_{N}$ for all $t\in[0,T]$, we have that $\gamma_{N}^{x}\in\SS^{N}(x):=\SS^{N}(x,0)$, where

[TABLE]

Remark 4.3**.**

(i)* Recall that the representation $m_{N}=\gamma_{N}\sharp m_{0}$ is only $m_{0}$-uniquely determined. In particular, if $\gamma_{N}^{\prime}\in\mathcal{A}$ is different from $\gamma_{N}$ but coincides with it on a set $A$ such that $m_{0}(A)=1$, then we also have that $m_{N}=\gamma_{N}^{\prime}\sharp m_{0}$. For the sake of simplicity, we have chosen to represent always $m_{N}$ via a measurable selection $\gamma_{N}$ of the set-valued map $\SS^{N}$. Notice that (A2) and the results in [12, Chapter 6] imply that $\gamma_{N}^{x}$ is uniquely defined for a.e. $x\in\mathbb{R}^{d}$.

(ii) Exactly as in the limit case (see Lemma 4.1(i)), we have the existence of a constant $C>0$, independent of $(x,t)$ and $N\in\mathbb{N}$, such that*

[TABLE]

As a consequence, there exists a compact set $K_{C}\subseteq\mathbb{R}^{d}$ such that $\gamma_{N}^{x}(t)\in K_{C}$ for all $N\in\mathbb{N}$, $x\in\mathop{\rm supp}(m_{0})$ and $t\in[0,T]$. In particular, the representation $\rho_{N}(t)=e_{t}\sharp m_{N}=\gamma_{N}^{(\cdot)}(t)$ implies that $\mathop{\rm supp}(\rho_{N}(t))\subseteq K_{C}$ for all $N\in\mathbb{N}$ and $t\in[0,T]$.

Let us define $u_{N}:\mathbb{R}^{d}\times[0,T]\to\mathbb{R}$ by

[TABLE]

Remark 4.4**.**

Mimicking the proofs of Lemma 4.1 and [12, Theorem 6.4.1] we obtain that $u_{N}$ is globally Lipschitz and locally semi-concave. Moreover, the Lipschitz and local semi-concavity constants are independent of $N$.

Arguing as in the previous subsection, the pair $(u_{N},\rho_{N})$ solves

[TABLE]

Conversely, associated to any solution $(u_{N},\rho_{N})$ we have the existence of a symmetric equilibrium in distributed open-loop strategies $(\gamma_{N},\ldots,\gamma_{N})\in\mathcal{A}^{N}$ for the $N$-player game. Moreover, using again the results in [12, Chapter 6], any $\gamma_{N}\in\mathcal{A}$ defining such equilibrium satisfies

[TABLE]

for a.e. $x\in\mbox{supp}(m_{0})$. Thus, we can think of the r.h.s. above as an optimal control which is feedback with respect to the individual states. We call $(\gamma_{N},\ldots,\gamma_{N})$ a Nash equilibrium in distributed Markov strategies for the $N$-player game. As a consequence of the previous discussion, such equilibria exist for all $N\in\mathbb{N}$ provided that (A2) holds true.

Let us consider a sequence $(\gamma_{N})_{N\in\mathbb{N}}$ of elements in $\mathcal{A}$ defining Nash equilibria in distributed Markov strategies for the $N$-player games. Theorem 3.1 yields the existence of a Lagrangian equilibrium $m_{\ast}\in\mathcal{A}_{rel}$ and a subsequence $(\gamma_{N_{k}})_{k\in\mathbb{N}}$ such that $\gamma_{N_{k}}\sharp m_{0}\to m_{\ast}$ as $k\to\infty$. Let $\gamma_{\ast}\in\mathcal{A}$ be such that $m_{\ast}=\gamma_{\ast}\sharp m_{0}$ and $\gamma_{\ast}^{x}\in\SS(x)$ for all $x\in\mathbb{R}^{d}$, i.e. an equilibrium in distributed Markov strategies for the MFG.

Our aim now is to study the convergence of the associated time marginals $\rho_{N_{k}}$ to $\rho$, the convergence of the associated value functions $u_{N_{k}}$ to $u$ and, finally, the convergence of $\gamma_{N_{k}}^{x}$ to $\gamma_{\ast}^{x}$ for a.e. $x\in\mathbb{R}^{d}$.

We will need the following preliminary result.

Lemma 4.2**.**

Assume (A2)(ii) and let $K\subseteq\mathbb{R}^{d}$ be a nonempty compact set. Consider a sequence of measures $(\mu_{k})_{k\in\mathbb{N}}\subseteq\mathcal{P}_{1}(\mathbb{R}^{d})$ such that $\mbox{{\rm supp}}(\mu_{k})\subseteq K$, for all $k\in\mathbb{N}$, and, as $k\to\infty$, $\mu_{k}\to\mu$ for some $\mu\in\mathcal{P}_{1}(\mathbb{R}^{d})$. Then, for any sequence $(x_{k})_{k\in\mathbb{N}}$ and $x\in\mathbb{R}^{d}$ such that $x_{k}\to x$, we have

[TABLE]

where $h=f$, $g$ and $h_{k}$ is defined by (4.25).

Proof.

Notice that (A2)(ii) implies that

[TABLE]

Now, let $Y^{k}_{1},\ldots,Y^{k}_{k}$ be independent and identically distributed $K$-valued random variables, defined on some probability space $(\Omega_{k},\mathcal{F}_{k},\mathbf{P}_{k})$, with common distribution $\mu_{k}$. Using that $\mathcal{P}(\mathcal{P}(K))$ is compact and arguing as in the proof of Theorem 3.1, we obtain that, as $k\to\infty$, the $\mathcal{P}_{1}(\mathbb{R}^{d})$-valued random sequence $\left(\frac{1}{k-1}\sum_{j=2}^{k}\delta_{Y^{k}_{j}}\right)_{k\in\mathbb{N}}$ converges in distribution to the deterministic measure $\mu$. Since (4.30) can be written as

[TABLE]

relation (4.29) follows by letting $k\to\infty$. ∎

Theorem 4.1**.**

*Assume that (A2) holds. Then, the following assertions hold true:

(i) The sequence $(\rho_{N_{k}})_{k\in\mathbb{N}}$ converges to $\rho$ in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$.

(ii) The sequence $(u_{N_{k}})_{k\in\mathbb{N}}$ converges to $u$ uniformly on compact subsets of $\mathbb{R}^{d}\times[0,T]$.

(iii) For a.e. $x\in\mathbb{R}^{d}$, the sequence $(\gamma^{x}_{N_{k}})_{k\in\mathbb{N}}$ converges to $\gamma^{x}_{\ast}$ uniformly in $[0,T]$ and $(\dot{\gamma}^{x}_{N_{k}})_{k\in\mathbb{N}}$ converges to $\dot{\gamma}^{x}_{\ast}$ in the weak* topology in $L^{\infty}([0,T];\mathbb{R}^{d})$.*

Proof.

Assertion (i) follows directly from Theorem 3.1. Note that (A2) implies that $(u_{N_{k}})_{k\in\mathbb{N}}$ is a sequence of uniformly bounded functions on $\mathbb{R}^{d}\times[0,T]$. Let us fix $(x,t)\in\mathbb{R}^{d}\times[0,T]$. The definition of $u_{N_{k}}$, Remark 4.3(ii) and Lemma 4.2 imply that

[TABLE]

Let $\gamma_{N_{k}}^{x,t}\in\SS^{N}(x,t)$ and $\gamma^{x,t}\in C([t,T];\mathbb{R}^{d})$ be a cluster point of $(\gamma_{N_{k}}^{x,t})_{k\in\mathbb{N}}$, with respect to the uniform convergence. The existence of $\gamma^{x,t}$ is ensured by (4.26) and the Arzelà-Ascoli theorem. Up to the extraction of a subsequence, we can assume that $\liminf_{k\to\infty}u_{N_{k}}(x,t)=\lim_{k\to\infty}u_{N_{k}}(x,t)$ and $\lim_{k\to\infty}\gamma_{N_{k}}^{x,t}=\gamma^{x,t}$ in $C([t,T];\mathbb{R}^{d})$. Using estimate (4.26) again, we get that $\dot{\gamma}^{x,t}$ exists and $\dot{\gamma}_{N_{k}}^{x,t}\to\dot{\gamma}^{x,t}$ in the weak* topology in $L^{\infty}([0,T];\mathbb{R}^{d})$. By the weak lower semi-continuity of the cost functional we obtain

[TABLE]

Thus, by (4.31) we get the pointwise convergence

[TABLE]

and hence, using that $u_{N_{k}}$ is Lipschitz continuous, with a Lipschitz constant which is independent of $k$, assertion (ii) follows from the Arzelà-Ascoli theorem. Finally, (iii) is a consequence of the previous analysis with $t=0$ and the fact that [12, Corollary 6.4.10] implies that for a.e. $x\in\mathbb{R}^{d}$ the optimization problem associated with $u(x,0)$ admits a unique solution. ∎

Recall that, as in the case of $(MFG)$, to each solution $(u_{N},\rho_{N})$ of $(MFG_{N})$ we can associate a symmetric equilibrium $(m_{N},\ldots,m_{N})\in\mathcal{A}_{rel}^{N}$ of the $N$-player game. As a consequence of this fact, Theorem 3.1 and Theorem 4.1, we have the following result.

Corollary 4.1**.**

Let $((u_{N},\rho_{N}))_{N\in\mathbb{N}}$ be a sequence of solutions to $(MFG_{N})$ ($N\in\mathbb{N}$). Then, there exists a solution $(u,\rho)\in C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$ to $(MFG)$ such that, up to some subsequence, $u_{N}\to u$ uniformly over compact subsets of $\mathbb{R}^{d}\times[0,T]$ and $\rho_{N}\to\rho$ in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$.

Remark 4.5**.**

If $h=f$, $g$ satisfies

[TABLE]

then the solution $(u,\rho)$ to $(MFG)$ is unique (see [34] and [30, Corollary 5.2]). Since any Lagrangian equilibrium $m_{\ast}$ can be represented by $\gamma_{\ast}\sharp m_{0}$, where $\gamma_{\ast}^{x}\in\SS(x)$ is uniquely determined for a.e. $x\in\mathbb{R}^{d}$, the Lagrangian equilibrium must also be unique. In this case, the results in Theorem 4.1 hold for the entire sequence $(u_{N},\rho_{N},\gamma_{N})$ and the result in Corollary 4.1 holds for the entire sequence $(u_{N},\rho_{N})$.

Finally, let us point out that the convergence result in Corollary 4.1 can also be established directly, without appealing to Theorem 3.1, under a stronger regularity assumption than (A2)(ii). Indeed, assume that, in addition to (A2), $m_{0}\in L^{p}(\mathbb{R}^{d})$ for some $p\in(1,+\infty]$. If $(u_{N},\rho_{N})$ solves $(MFG_{N})$, then by [3, Theorem 8.2.1] and the results in [12, Chapter 6], we must have that $\rho_{N}(t)=\gamma_{N}^{(\cdot)}(t)\sharp m_{0}$, for some $\gamma_{N}\in\mathcal{A}$ such that $\gamma_{N}^{x}\in\SS^{N}(x)$ for all $x\in\mathop{\rm supp}(m_{0})$. Arguing as in the proof of Lemma 4.1 we get the existence of $C^{\prime}>0$ such that

[TABLE]

Therefore, there exists $C>0$ and a compact set $K\subseteq\mathbb{R}^{d}$, both independent of $N$, such that $\gamma_{N}^{x}(t)\in K$ for a.e. $x\in\mathop{\rm supp}(m_{0})$ and all $t\in[0,T]$. In particular, $\mathop{\rm supp}(\rho_{N}(t))\subseteq K$ and $d_{1}(\rho_{N}(s),\rho_{N}(t))\leq C|s-t|$ for all $s$, $t\in[0,T]$ and $N\in\mathbb{N}$. This implies the existence of $\rho\in C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$ such that, up to some subsequence, $\rho_{N}\to\rho$ in $C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$ as $N\to\infty$. Since Lemma 4.2 implies that $f_{N}(\cdot,\rho_{N}(\cdot))$ converges uniformly to $f(\cdot,\rho(\cdot))$ on compact subsets of $\mathbb{R}^{d}\times[0,T]$, standard stability results for viscosity solutions of Hamilton-Jacobi-Bellman equations imply that, up to some subsequence, $u_{N}\to u$ uniformly on compact subsets of $\mathbb{R}^{d}\times[0,T]$, $u$ being the unique viscosity solution to (4.10). In particular, $u_{N}$ being locally semi-concave with respect to the space variable, uniformly in $N$, for all $t\in[0,T]$ we have that $\nabla u_{N}(x,t)\to\nabla u(x,t)$ for a.e. $x\in\mathbb{R}^{d}$. Using that $u_{N}$ is globally Lipschitz, uniformly in $N$, by dominated convergence we deduce that

[TABLE]

On the other hand, using again the uniform local semiconcavity of $u_{N}(\cdot,t)$ and arguing as in the proof Proposition 4.1 we get the existence of $c_{3}>0$, independent of $N$, such that

[TABLE]

Using this bound, we obtain that $\rho(t)$ is absolutely continuous w.r.t. the Lebesgue measure for a.e. $t\in[0,T]$, and its density, denoted likewise by $\rho(t)$, satisfies $\|\rho(t)\|_{L^{p}}\leq c_{3}\|m_{0}\|_{L^{p}}$ for a.e. $t\in[0,T]$. Since $\rho\in C([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}))$, the previous bound implies that for all $t\in[0,T]$, the measure $\rho(t)$ is absolutely continuous w.r.t. the Lebesgue measure and the estimate $\|\rho(t)\|_{L^{p}}\leq c_{3}\|m_{0}\|_{L^{p}}$ holds. Moreover, using (4.32)-(4.33) we can pass to the limit in the second equation $(MFG_{N})$ to obtain that the pair $(u,\rho)$ solves $(MFG)$.