arxiv: v1 [math.ap] 31 Jul 2014

Transcription

1 EXISTENCE FOR STATIONARY MEAN FIELD GAMES WITH QUADRATIC HAMILTONIANS WITH CONGESTION arxiv: v1 [math.ap] 31 Jul 2014 DIOGO A. GOMES AND HIROYOSHI MITAKE Abstract. In this paper, we investigate the existence and uniqueness of solutions to a stationary mean field game model introduced by J.-M. Lasry and P.-L. Lions. This model features a quadratic Hamiltonian with possibly singular congestion effects. Thanks to a new class of a-priori bounds, combined with the continuation method, we prove the existence of smooth solutions in arbitrary dimensions. 1. Introduction Since the seminal papers [HMC06, HCM07, LL06a, LL06b, LL07a, LL07b], research on mean field games has been extremely active see, for instance, the recent surveys [LLG10, Car11, Ach13, GS14] and the references therein). Nevertheless, several fundamental questions have not yet been answered. In this paper, we address one of those, and prove existence and uniqueness of smooth solutions for stationary mean field games with congestion and quadratic Hamiltonian. Mean field games model large populations of rational agents who move according to certain stochastic optimal control goals. To simplify the presentation, we will work in the periodic setting, that is in the d dimensional standard torus, d 1. We consider a large stationary population of agents, whose statistical information is encoded in an unknown probability density m : [0,+ ) R. Each individual agent wants to minimize an infinite horizon discounted cost given by ux,t) = inf { E [ + t e s mxs),s) α vs) bxs)) 2 2 where the infimum is taken over all progressively measurable controls v, dx = vdt+ 2dW t with X0) = x, ) ]} +VXs),mXs),s)) ds, W t and E denote a standard d-dimensional Brownian motion and the expected value, respectively. The constant 0 < α < 1 determines the strength of congestion effects in the term m α v bx) 2, and makes it costly to move in areas of high density with a drift v substantially different from a reference vector field b : R d. The function V : R + R accounts for additional spatial preferences of the agents. We assume b and V to be smooth functions. Date: August 1, Mathematics Subject Classification. 35J47, 35A01. Key words and phrases. Mean Field Game; Quadratic Hamiltonians; Congestion. DG was partially supported by KAUST SRI, Uncertainty Quantification Center in Computational Science and Engineering, CAMGSD-LARSys through FCT-Portugal and by grants PTDC/MAT- CAL/0749/2012. HM was partially supported by JST program to disseminate tenure tracking system. 1

2 2 DIOGO A. GOMES AND HIROYOSHI MITAKE Under standard assumptions of rationality and symmetry, one can derive a mean field problem which models this setup. A detailed discussion can be found in [Lio11, LL07a], where the following problem, consisting of a viscous Hamilton-Jacobi equation for u, coupled with a Fokker-Planck equation for m, was introduced: u t +u u+ Du 2 +bx) Du = Vx,m) 2mα m t +m m div m 1 α Du ) divmb) = 0, together with initial conditions for mx, 0) and suitable asymptotic behavior for u. In the present paper, we consider a stationary version of this problem, which for u,m : R, m > 0, is given by the system u u+ Du 2 2m +bx) Du = Vx,m) in α Td, 1.1) m m div m 1 α Du ) divmb) = 1, in. 1.2) where the right hand side of the second equation is an additional source term for m to avoid the trivial solution m = 0). In [Lio11], only the uniqueness of smooth solutions was proven. However, existence of solutions, the main result of this paper, was not yet known for both stationary and time-dependent problems. The fundamental difficulty lies in the possibly singular behavior due to congestion. The dependence on m in the optimal control problem causes the singularity in the equation 1.1), for which we had to develop a new class of estimates. Thanks to those we obtain our main result, which is Theorem 1.1. Assume the following: A1) 0 α < 1; A2) V : R + R, Vx,m) C R + ) is globally bounded with bounded derivatives and non-decreasing with respect to m; A3) b : R d, b C ). Then there exists a solution u,m) C ) C ) to 1.1)-1.2) with m > 0. Furthermore, if V is strictly increasing with respect to m, then a solution is unique. Numerous a-priori bounds for mean field games have been proved by various authors, including the first authorsee, for instance[ll06a], [LL06b], [LL07a], [GSM14], [GISMY10], [CLLP12], [GPSM12], [GPSM14], [GPSM13] [GPV14], [Por14], [Por13]). However, these bounds were designed to address a different coupling, namely mean field games where the local dependence on m is not singular when m = 0. A typical example is the following system u u+ Du 2 = m α in 2 1.3) m m divmdu) = 1 in. In 1.3), the main difficulties are caused by the growth of the nonlinearity m, especially for large α > 0, rather than singularities caused by m vanishing. Besides, 1.3) can be regarded as an Euler-Lagrange equation of a suitable functional, whereas 1.1)-1.2) does not have this structure. In Section 2, we start by exploring the special form of 1.1) and 1.2) to obtain a bound for m 1 L T n ). This estimate, combined with the techniques from [AC78], yields

3 STATIONARY MEAN FIELD GAMES WITH CONGESTION AND QUADRATIC HAMILTONIANS 3 a-prioriregularityinw 2,p )foranyp 1. Fromthis, asimpleargument showsthatany solution to 1.1)-1.2) is bounded in any Sobolev space W k,p ). Then, in Section 3, we prove the existence of solutions to 1.1)-1.2) by using the continuation method together with the aforementioned a-priori estimates. In Appendix A, for completeness, we present the uniqueness proof for solutions to 1.1)-1.2), based upon the ideas in [Lio11] see also [Gue14]). In a forthcoming paper [GM13], we will study general mean field games with congestion, for which the techniques of the present paper cannot be applied directly, as discussed in Remark 1 at the end of the next section. 2. A-priori Estimates In this section, we obtain a-priori bounds for solutions of 1.1)-1.2). In particular, we prove a L bound for m 1. From this, we derive estimates for u,m in W 2,p ) for any p 1. Then, by a bootstrapping argument, we establish smoothness of solutions. Proposition 2.1. There exists a constant C := C V ) 0 such that for any classical solution u,m) of 1.1)-1.2) we have u L ) C. Furthermore, m 0 on Td, and m L 1 ) = 1. Proof. The L bound is obtained by evaluating the equation at points of maximum of u resp., minimum) and using the facts that at those points Du = 0, u 0 resp., 0), and V is bounded on [0, ). We then observe that m is non-negative by the maximum principle. Moreover, it has total mass 1 by integrating 1.2). Proposition 2.2. There exists a constant C := C b, V ) 0 such that for any classical solution u, m) of 1.1)-1.2) we have 1 m C. L ) Proof. Let r > α. Subtract equation 1.2) divided by r +1 α)m r+1 α from equation 1.1) divided by rm r. Then, ] [u u+ Du 2 T 2m +b Du V 1 dx 2.1) d α rmr [ m m div m 1 α Du ) ] 1 divmb) dx T r +1 α)m d r+1 α 1 = dx. T r +1 α)m d r+1 α Next, observe that u T rm dx = Du Dm dx, d r T m d r+1 and div m 1 α Du ) Du Dm dx = dx. T r +1 α)m d r+1 α T m d r+1 Hence u T rm dx div m 1 α Du ) dx = ) d r r +1 α)mr+1 α

4 4 DIOGO A. GOMES AND HIROYOSHI MITAKE Also, note the identity divbm) m r 1+α = r +1 α = ) m r+α b D r α m r 1+α b Dm divb)m r+α. = 1 r α Then 2.1) is reduced to 1 Td Du 2 Td Dm 2 dx+ dx+ dx T r +1 α)m d r+1 α 2rmr+α mr+2 α [ = V T rm u d r rm + 1 b Du 1 ] r r +1 α)mr α rm r r α divb)m r+α C rm dx+ C Td b 2 Du 2 dx+ + dx r r α)mr α rmr α 4rmr+α in view of Proposition 2.1. Consequently, r +1 α)m C rm dx+ r 1 r+1 α dx+ C dx. r α)mr α By Young s inequality, for α [0,1), we have and with Therefore, C 1 r Td Du 2 Td Dm 2 dx+ dx 4rmr+α mr+2 α C rm 1 r 4r +1 α)m r+1 α +C1 r, C r α)m r α 1 4r α)m r+1 α +C2 r, r r+1 α 1 α)41 α C 1 α :=, Cr 2 rr +1 α) := 4r α C r+1 α r α) r α 1. r +1 α) r+1 α 1 1 r+1 α T m d r α+1 2C1 r +Cr). 2 dx Thus, we get 1 m L r+1 α ) [2r+1 α)c 1r +C2r ) ] 1 r+1 α We can easily check that, for any r 0 > α there exists C α for which C α r) C α for all r [r 0, ). =: C α r). Proposition 2.3. For any p 1 there exists a constant C := C p b, V ) > 0 such that for any classical solution u,m) of 1.1)-1.2), we have u W 2,p )+ m W 2,p ) C.

5 STATIONARY MEAN FIELD GAMES WITH CONGESTION AND QUADRATIC HAMILTONIANS 5 Proof. Let u,m) be a classical solution u,m) to 1.1)-1.2). In view of Lemma [AC78, Lemma 4], combined with Proposition 2.2 we conclude that for all p [1, ) there exists C = C V, b,p) such that u W 2,p ) C. In light of the Sobolev embedding theorem, we get u C 1,γ ) C u W 2,p ) C. 2.3) Then, multiplying 1.2) by m p and using Young s inequality yield m p+1 dx+p m p 1 Dm 2 dx = m p dx p m p g +b) Dmdx T [ ] [ d ] 1 p m p+1 dx+c + m p 1 Dm 2 dx+cp g 2 + b 2 )m p+1 dx, 2 T 2 d whereg := Du/m α. Notingthat g C inviewof 2.3)andm p 1 Dm 2 = C p Dm p+1)/2 2, we get m p+1 dx+c p Dm p+1)/2 2 dx C +C p m p+1 dx. 2.4) Using Hölder s inequality we have ) 1/p+1) ) 2/2+dp) m p+1 m m 2 p+1)/2 )dp/2+dp) 2 p+1)/2. Moreover, using the Sobolev embedding theorem, we get )dp/2+dp) m p+1 2 m 2 p+1)/2 /2 T d ) dp/2+dp) C m p+1 dx+ Dm p+1)/2 2 dx. Then, using the previous estimate and the fact that dp/2 + dp) < 1 in the right-hand side of 2.4), we conclude that m p+1 dx+ Dm p+1)/2 2 dx C. 2.5) Note now that if m W 1,q ), we have m m = m 1 α u+1 α)m α Du Dm+divmb)+1 L q ). 2.6) Thus, by standard elliptic regularity m W 2,q ). Consequently m W 1,q ). In light of 2.5) for p = 1 we have m W 1,2 ). Thus we obtain m W 2,2 ) and m W 1,2 ). By iterating this argument, we finally get m W 2,q ) for any q <. Proposition 2.4. For any integer k 0 there exists a constant C := C b, V,k) > 0 suchthat anyclassicalsolution u,m)of 1.1)-1.2) satisfies u W k, )+ m W k, ) C.

6 6 DIOGO A. GOMES AND HIROYOSHI MITAKE Proof. Note that Dm α ) = m 1+α) Dm L p ) for large p > 1 in view of Propositions 2.2, 2.3. This implies m α W 1,p ). Hence, by Morrey s theorem, we have m α C γ ) for some γ 0,1). Also note that Du 2 C γ ). Therefore, going back to the equation 1.1) we have u u = Du 2 2m b Du+V α Cγ ). 2.7) Then, in view of the elliptic regularity theory we get u C 2,γ ). Note that the norm in C γ ) of the right hand side of 2.7) is estimated by a constant which only depends on b, V. Thus, u C b, V ). Next, going back to the equation 2.6) for m and noting that the right hand side is C γ ) now, we get m C 2,γ ) with m C b, V ). Once we know that u,m C 2,γ ), 2.7) and 2.6) imply u,m C 3,γ ). By continuing this so-called bootstrap argument, we get the conclusion. Remark 1. The methods in this section rely strongly on the particular structure of the Hamiltonian which allows for the cancellation in 2.2). In the forthcoming paper [GM13], we address mean field games which generalize 1.1)-1.2), namely of the form u u+m α Hx, Du mα) = Vx,m) in Td m m div md p Hx, Du m α)) = 1, in, but forwhich the cancellation in 2.2) no longer holds. Such problems arenatural in many applications, for instance 2.2) is not valid for the anisotropic quadratic Hamiltonian Hx,p) = 2 1 Ax)p + bx) 2, where A is a strictly positive definite matrix. This is addressed in [GM13] with a different approach. 3. Existence by Continuation Method In this section we prove the existence of a unique classical solution to 1.1)-1.2) by using the continuation method. We work under the assumptions of Theorem 1.1. For 0 λ 1 we consider the problem u λ u λ + Du λ 2 +λbx) Du 2m α λ λvx,m) 1 λ)v 0 m) = 0 in, λ m λ m λ div ) 3.1) m 1 α λ Du λ λdivbmλ ) = 1 in, where V 0 m) := arctanm). Let E k := H k ) H k ) for k N, and E 0 := L 2 ) L 2 ). For any k 0 N with k 0 > d/2, we define the map F : [0,1] E k0+2 E k 0 by Fλ,u,m) := Then, we can rewrite 3.1) as u u+ Du 2 2m α +λbx) Du λvx,m) 1 λ)v 0m) m m div m 1 α Du ) λdivbm) 1. Fλ,u λ,m λ ) = 0..

7 STATIONARY MEAN FIELD GAMES WITH CONGESTION AND QUADRATIC HAMILTONIANS 7 Note that for any γ > 0, the map F is C in the set {u,m) E k 0+2 ),m > γ}. This is because for k 0 > d/2, the Sobolev space H k 0 ) is an algebra. Moreover, if k 0 is large enough, then any solution u λ,m λ ) in E k 0+2 is, in fact, in E k+2 for all k N, by the a-priori bounds in Section 2. We define the set Λ by Λ := { λ [0,1] 3.1) has a classical solution u,m) E k 0+2 }. When λ = 0 we have an explicit solution, namely u 0,m 0 ) = π/4,1). Therefore, Λ. The main goal of this section is to prove Λ = [0,1]. To prove this, we show that Λ is relatively closed and open on [0,1]. The closeness of Λ is a straightforward consequence of the estimates in Section 2: Proposition 3.1. The set Λ is closed. Proof. To prove that Λ is closed, we need to check that for any sequence λ n Λ such that λ n as n, we have Λ. Fix such a sequence and corresponding solutions u λn,m λn )to3.1)withλ = λ n. Sincethea-prioriboundinProposition2.4isindependent of n N, by taking a subsequence, if necessary, we may assume that u λn,m λn ) u,m) in E k0+2. Moreover, m 1 λ n m 1 in C ). Therefore, we can take the limit in 3.1), and we get that u,m) is solution to 3.1) with λ =. This implies Λ. To prove that Λ is relatively open in [0,1], we need to check that for any Λ there exists a neighborhood of contained in Λ. To do so, we will use the implicit function theorem see, forexample, [Die69], chapter X). Forafixed Λ, we consider the Fréchet derivative L λ0 : E k 0+2 E k 0 of u,m) F,u,m) at the point u λ0,m λ0 ), which is given by L λ0 v,f) = v v + Du Dv m α α Du 2 f 2m α+1 f f div m 1 α Dv ) 1 α)div m α + b Dv D m V +1 )D m V 0 ) f fdu λ0 ) λ0 divbf). 3.2) Because of the a-priori bounds for u and m in Section 2, we can extend the domain of L λ0 by continuity to E k+2 for any k k 0. We will prove that L λ0 is an isomorphism from E k+2 to E k for any k 0. Define the bilinear mapping B λ0 [w 1,w 2 ] : E 1 R by B λ0 [w 1,w 2 ] [ := v 1 + Du Dv 1 α Du 2 f 1 T m α d 2m α+1 + b Dv 1 ] ) D m V +1 )D m V 0 f1 +Dv 1 Df 2 m 1 α Dv 1 Dv 2 + [ f 1 1 α)div m α f 1 Du λ0 ) λ0 divbf 1 ) ] v 2 ) Df 1 Dv 2 dx. f 2

8 8 DIOGO A. GOMES AND HIROYOSHI MITAKE We set Pw := f, v) for w = v,f), and we observe that if w 1 E k with k 2, then B λ0 [w 1,w 2 ] = L λ0 w 1 ) Pw 2 dx. The boundedness of B λ0 is a straightforward result of Proposition 2.4: Lemma 3.2. There exists a constant C > 0 such that for any w 1,w 2 E 1. B λ0 [w 1,w 2 ] C w 1 E 1 w 2 E 1. Thus, in view of the Riesz representation theorem for Hilbert spaces, there exists a linear mapping A : E 1 E 1 such that Lemma 3.3. The operator A is injective. B λ0 [w 1,w 2 ] = Aw 1,w 2 ) E 1. Proof. Let w = v,f). By Young s inequality we have αdu λ0 Dv B λ0 [w,w] = f α Du 2 T m α d 2m α+1 f 2 ) D m V +1 )D m V 0 f 2 m 1 α Dv 2 dx ) D m V +1 )D m V 0 f 2 + α 2)m 1 α Dv 2 T 2 d C λ0 Dv 2 L 2 ) + f 2 L 2 ) ) for a constant C λ0 which depends on bounds for m λ0, and Du λ0, but it is strictly positive for any solution to 3.1) since 0 α < 1. We have used Assumption A1) and the strict monotonicity of V on m. This implies that if Aw = 0 we have w = µ,0), for some constant µ. Then, by computing 0 = Aw,0,µ)) = B[µ,0),0,µ)] = µ 2, we conclude that µ = 0. Remark 2. Note that the injectivity of the operator A holds for all 0 α < 2. However, the a-priori estimates of the previous section are only valid for 0 α < 1. Lemma 3.4. The range RA) is closed, and RA) = E 1. Proof. Take a Cauchy sequence z n in the range of A, that is z n = Aw n, for some sequence w n = v n,f n ). We claim that w n is a Cauchy sequence. We have z n z m,w n w m ) E 1 = Aw n w m ),w n w m ) E 1 C λ0 Dv n v m ) 2 L 2 ) + f n f m 2 L 2 ) ).

9 STATIONARY MEAN FIELD GAMES WITH CONGESTION AND QUADRATIC HAMILTONIANS 9 Note that z n z m,w n w m ) E 1 z n z m E 0 w n w m E 0 + Dz n z m ) E 0 Dw n w m ) E 0 z n z m E 0 vn v m L 2 ) + f n f m L 2 T )) d + Dz n z m ) E 0 Dvn v m ) L 2 ) + Df n f m ) L 2 T )) ) d ε f n f m 2 L 2 ) + Dv n v m ) 2 L 2 ) +C ε z n z m 2 E 1 + z n z m E 1 vn v m L 2 ) + Df n f m ) L 2 T )) d for ε > 0. If we fix a suitable small ε and combine the inequalities above, then we obtain We have Dv n v m ) 2 L 2 ) + f n f m 2 L 2 ) C z n z m 2 E 1 +C z n z m E 1 v n v m L 2 ) + Df n f m ) L 2 )) C z n z m 2 E 1 +C z n z m E 1 w n w m E ) B[w n w m, f n +f m,v n v m,)] = w n w m 2 E 1 +E nm, where, using 3.3), E nm satisfies E nm C v n v m L 2 ) Dv n v m ) L 2 ) +C f n f m L 2 ) v n v m L 2 ) +C Df n f m ) L 2 ) Dv n v m ) L 2 ) +C f n f m L 2 ) Df n f m ) L 2 ) C w n w m E 1 zn z m 2 E 1 + z n z m E 1 w n w m E 1) 1/2. 3.4) On the other hand, by Lemma 3.2 we have B[w n w m, f n +f m,v n v m,)] C z n z m E 1 w n w m E ) Combining 3.4) and 3.5) we deduce w n w m 2 E 1 C z n z m E 1 w n w m E 1 +C w n w m E 1 zn z m 2 E 1 + z n z m E 1 w n w m E 1) 1/2. By using Young s inequality we conclude w n w m 2 E 1 C z n z m 2 E 1. From this we get convergence in E 1. Finally, we prove that RA) = E 1. Suppose that RA) E 1. Since RA) is closed, there would exist z RA) with z 0 such that B λ0 [z,z] = 0. The argument in the proof of Lemma 3.3 implies z = 0 which is a contradiction. Lemma 3.5. The operator L λ0 : E k+2 E k is an isomorphism for all k N with k 2.

10 10 DIOGO A. GOMES AND HIROYOSHI MITAKE Proof. Since L λ0 is injective, it suffices to prove that it is surjective. To do so, fix w 0 E k with w 0 = v 0,f 0 ). We claim there exists a solution w 1 E k+2 to L λ0 w 1 = w 0. Consider the bounded linear functional w w 0,w) E 0 in E 1. By the Riesz representation theorem, there exists w E 1 such that w 0,w) E 0 = w,w) E 1 for any w E 1. In light of Lemmas 3.3 and 3.4, there exists the inverse of A. We define w 1 := A 1 w, and write w 1 = v,f). Set g1 [v,f] g 2 [v,f] ) Du λ0 Dv := α Du 2 f m α 2m α+1 + b Dv ) D m V +1 )D m V 0 f div m 1 α Dv ) 1 α)div ). m α fdu λ0 λ0 divbf) Then, the identity Aw 1,w) E 1 = w,w) E 1 = w 0,w) E 0 for any w E 1, means that v is a weak H 1 ) solution to v v = g 1 [v,f]+v 0, and that f H 1 ) is also a weak solution to f f = g 2 [v,f]+f 0. Observe that if v,f H j+1 ) then g 1,g 2 H j ). Additionally, elliptic regularity yields, from g i [v,f] H j ), that v,j H j+2 ). Since we have v,f H 1 ), we conclude by induction that v,f H j+2 ), for all j k. A straightforward result of Lemma 3.5 and the implicit function theorem in Banach space is Proposition 3.6. The set Λ is relatively open in [0,1]. We finally address the existence of solutions to 1.1), 1.2), and complete the proof of Theorem 1.1. Proof of Theorem 1.1. If V is strictly increasing on m, the existence of a classical solution to 1.1), 1.2) is a straightforward result of Proposition 3.6. Uniqueness of solution is discussed in the next appendix, Proposition A.1. If we only assume V to be nondecreasing on m, existence can be obtained by using a perturbation argument similar to the one in 3.1). More precisely, we add a small perturbation ε arctanm) to V so that we make the potential term strictly monotone. This problem admits a unique classical solution u ǫ,m ǫ ). Because the a-priori bounds in the previous section do not depend on the strict monotonicity of V, u ǫ,m ǫ ) satisfy uniform bounds in any Sobolev space. Thus, by compactness, we can extract a convergent subsequence to a limit u, m) which solves 1.1)-1.2). Remark 3. In this paper, we focus on the case where V is bounded on m, since our main concern is a lower bound on m. In principle, unbounded potentials V can be studied by adapting the techniques in [LL07a, CLLP12, GPSM12, GPSM14], for instance.

11 STATIONARY MEAN FIELD GAMES WITH CONGESTION AND QUADRATIC HAMILTONIANS 11 Appendix A. Uniqueness Uniqueness of solutions of 1.1)-1.2) is well understood see [LL06b, Gue14] for a related problem). However, to make this paper self-contained, we give a proof based on Lions ideas in [Lio11]. Proposition A.1. The system 1.1)-1.2) admits at most one classical solution u, m). Proof. Let u 0,m 0 ) and u 1,m 1 ) be classical solutions to 1.1)-1.2). Subtract 1.1) for u 1,m 1 ) from 1.1) for u 0,m 0 ), and 1.2) for u 1,m 1 ) from 1.2) for u 0,m 0 ), respectively, and then u 0 u 1 = u 0 u 1 )+ Du 1 2 2m α 1 Du 0 2 2m α 0 +b Du 1 u 0 )+Vx,m 0 ) Vx,m 1 ), A.1) m 0 m 1 = m 0 m 1 )+divm 1 α 0 Du 0 ) divm 1 α 1 Du 1 )+divbm 0 m 1 )). A.2) Subtract A.2) multiplied by u 0 u 1 from A.1) multiplied by m 0 m 1, and then Du1 2 Du ) 0 2 m T 2m α d 1 2m α 0 m 1 )dx+ m 1 α 0 Du 0 m 1 α 1 Du 1 ) Du0 u 1 )dx 0 T d = Vx,m 1 ) Vx,m 0 ))m 0 m 1 )dx. A.3) We prove that the left-hand side of A.3) is non-negative if α [0,2] following the technique in [Gue14]. Set u θ := u 0 +θu 1 u 0 ) and m θ := m 0 +θm 1 m 0 ) for θ [0,1]. Define [ Iθ) := Then Duθ 2 2m α θ Du ) 0 2 2m α 0 m 1 m 0 ) + m 1 α θ d Du θ Du 1 u 0 )m 1 m 0 ) Iθ) = α dt T m α d θ Td Du θ 2 m 1 m 0 ) 2 + α 2 1 α ) 2 m 1+α θ dx+ ] Du θ m 1 α 0 Du 0 ) Du1 u 0 )dx. dx m 1 α θ m 1 α θ Du 1 u 0 ) 2 dx 0 Du 1 u 0 ) 2 dx for α [0,2]. Noting that I0) = 0, we conclude that I1) 0 which is the claim. Thus, by A.3) and the assumption that V is strictly increasing on m, we get the conclusion. Acknowledgment. The authors would like to thank Hung V. Tran for suggestions and comments.

12 12 DIOGO A. GOMES AND HIROYOSHI MITAKE References [AC78] H. Amann and M. G. Crandall. On some existence theorems for semi-linear elliptic equations. Indiana Univ. Math. J., 275): , [Ach13] Y. Achdou. Finite difference methods for mean field games. In Hamilton-Jacobi Equations: Approximations, Numerical Analysis and Applications, pages Springer, [Car11] P. Cardaliaguet. Notes on mean-field games [CLLP12] P. Cardaliaguet, J.-M. Lasry, P.-L. Lions, and A. Porretta. Long time average of mean field games. Netw. Heterog. Media, 72): , [Die69] J. Dieudonné. Foundations of modern analysis. Academic Press, New York-London, Enlarged and corrected printing, Pure and Applied Mathematics, Vol. 10-I. [GISMY10] D. A. Gomes, R. Iturriaga, H. Sánchez-Morgado, and Y. Yu. Mather measures selected by an approximation scheme. Proc. Amer. Math. Soc., 13810): , [GM13] D. A. Gomes and H. Mitake. On mean field games with congestion. Work in progress, [GPSM12] D. A. Gomes, G. E. Pires, and H. Sánchez-Morgado. A-priori estimates for stationary meanfield games. Netw. Heterog. Media, 72): , [GPSM13] D. A. Gomes, E. Pimentel, and H Sanchez-Morgado. Time dependent mean-field games in the superquadratic case. arxiv: , [GPSM14] D. A. Gomes, E. Pimentel, and H Sanchez-Morgado. Time dependent mean-field games in the subquadratic case. To appear in Comm. Partial Differential Equations, [GPV14] D. A. Gomes, S. Patrizi, and V. Voskanyan. On the existence of classical solutions for stationary extended mean field games. Nonlinear Anal., 99:49 79, [GS14] D. A. Gomes and J. Saúde. Mean Field Games Models A Brief Survey. Dyn. Games Appl., 42): , [GSM14] D. A. Gomes and H. Sánchez-Morgado. A stochastic Evans-Aronsson problem. Trans. Amer. [Gue14] Math. Soc., 3662): , O. Gueant. Existence and uniqueness result for mean field games with congestion effect on graphs. Preprint. [HCM07] M. Huang, P. E. Caines, and R. P. Malhamé. Large-population cost-coupled LQG problems with nonuniform agents: individual-mass behavior and decentralized ǫ-nash equilibria. IEEE Trans. Automat. Control, 529): , [HMC06] M. Huang, R. P. Malhamé, and P. E. Caines. Large population stochastic dynamic games: closed-loop McKean-Vlasov systems and the Nash certainty equivalence principle. Commun. Inf. Syst., 63): , [LL06a] J.-M. Lasry and P.-L. Lions. Jeux à champ moyen. I. Le cas stationnaire. C. R. Math. Acad. Sci. Paris, 3439): , [LL06b] J.-M. Lasry and P.-L. Lions. Jeux à champ moyen. II. Horizon fini et contrôle optimal. C. R. Math. Acad. Sci. Paris, 34310): , [LL07a] J.-M. Lasry and P.-L. Lions. Mean field games. Jpn. J. Math., 21): , [LL07b] J.-M. Lasry and P.-L. Lions. Mean field games. Cahiers de la Chaire Finance et Développement Durable, [LLG10] J.-M. Lasry, P.-L. Lions, and O. Gueant. Mean field games and applications. Paris-Princeton lectures on Mathematical Finance, [Lio11] P.-L. Lions. College de France course on mean-field games [Por13] [Por14] A. Porretta. Weak solutions to Fokker-Planck equations and mean field games. Preprint, A. Porretta. On the Planning Problem for the Mean Field Games System. Dyn. Games Appl., 42): , 2014.

13 STATIONARY MEAN FIELD GAMES WITH CONGESTION AND QUADRATIC HAMILTONIANS 13 D. A. Gomes) King Abdullah University of Science and Technology KAUST), CSMSE Division, Thuwal Saudi Arabia, and KAUST SRI, Uncertainty Quantification Center in Computational Science and Engineering. address: H. Mitake) Institute for Sustainable Sciences and Development, Hiroshima University Kagamiyama, Higashi-Hiroshima-shi , Japan address: