ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game (MFG) theory is that the coupling between the Hamilton-Jacobi equation and the transport equation is monotonically non-decreasing in the density of the population. In many cases, this assumption implies the existence and uniqueness of solutions. Here, we drop that assumption and construct explicit solutions for one-dimensional MFGs. hese solutions exhibit phenomena not present in monotonically increasing MFGs: lowregularity, non-uniqueness, and the formation of regions with no agents. 1. Introduction Mean-field game (MFG) theory [3, 17, 1, 3] describes on-cooperative differential games with infinitely many identical players. hese games were introduced by Lasry and Lions [9, 3, 31] and, independently around the same time, by Huang, Caines and Malhamé [7, 8]. Often, MFGs are given by a Hamilton-Jacobi equation coupled with a Fokker-Planck equation. A standard example is the stationary, one-dimensional, first-order MFG: { (ux+p) + V (x) = g(m) + H, (1.1) (m(u x + p)) x =, with its elliptic regularization, { ɛu xx + (p+ux) + V (x) = H + g(m) ɛm xx ((p + u x )m) x =, (1.) where ɛ >. Here, p is a fixed real number and the unknowns are the constant H and the functions u and m. he function g is C on R +. o simplify the presentation, we consider the periodic case and work in the one-dimensional torus,. Accordingly, V : R is a C potential. We search for periodic solutions, u, m : R. Here, we examine this problem and attempt to understand its features in terms of the monotonicity properties of g. A standard assumption in MFGs is that g is increasing. Heuristically, this assumption means that agents prefer sparsely populated areas. In this case, the existence and uniqueness of smooth solutions to (1.1) is well understood for stationary problems [18, 19,, 33], weakly coupled MFG systems [11], the obstacle MFG problem [1] and extended MFGs [13]. In the time-dependent setting, similar results are obtained in [14, 15, 4] for standard MFGs and in [16, 3] for forward-forward problems. he theory of weak solutions is also well developed for first-order and second-order problems (see [4, 5, 7] and [6, 8, 34, 35], respectively). Congestion problems, see [9,, 5], are also of interest and our results extend straightforwardly [1]. Date: July 8, 18. 1 Mathematics Subect Classification. 35J47, 35A1. Key words and phrases. Mean-field games; stationary problems; dynamic games. D. Gomes, L. Nurbekyan and M. Prazeres were partially supported by KAUS baseline and start-up funds. 1
DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES he case of a non-monotonically increasing g is relevant: if g is decreasing, agents prefer clustering in high-density areas. he case where g first decreases and then increases is also natural; here, agents have a preferred density given by the minimum of g. However, little is known about the properties of (1.1) when g is not increasing. One of the few known cases is a second-order MFG with g(m) = ln m and a quadratic cost. In this case, due to the particular structure of the equations, there are explicit solutions, see [6, 3]. A triplet, (u, m, H), solves (1.1) if i. u is a Lipschitz viscosity solution of the first equation in (1.1); ii. m is a probability density; that is, m, m = 1; iii. m is a weak (distributional) solution of the second equation in (1.1). Because (1.1) is invariant under addition of constants to u, we assume that u() =. Here, u is Lipschitz continuous. However, m can be discontinuous. In this case, viscosity solutions of the first equation in (1.1) are interpreted as discontinuous viscosity solutions; see, for example, [1] and the discussion in Section 6. Our problem is one-dimensional and the Hamiltonian is convex. If u is a piecewise C 1 function and m is continuous, then u is a viscosity solution if the following conditions hold: a. u solves the equation at the points where it is C 1 and m is continuous; b. lim x x u x (x) lim x x + u x (x) at points of discontinuity of u x. When g is not increasing, (1.1) may not admit m continuous. Solutions must, therefore, be considered in the framework of discontinuous viscosity solutions. In this case, the above characterization of one-dimensional viscosity solutions is not valid, and (1.1) admits a large family of discontinuous viscosity solutions (see Section 6). On the other hand, solutions that satisfy the above conditions (a. and b.) have nice structural properties that we discuss in this paper. Furthermore, in their analysis we see the appearance of discontinuities in m, which in turn motivates the study of discontinuous viscosity solutions. Overall, these conditions seem to be good selection criteria for discontinuous solutions of (1.1). We call solutions that satisfy conditions a. and b. regular (they can still be discontinuous). In this paper, we always consider regular solutions except in Section 6, where we discuss general discontinuous viscosity solutions. Furthermore, when m is continuous the term regular is superfluous. hus, except in Section 6, we refer to regular solutions. Our goal is to solve (1.1) explicitly and to understand the qualitative behavior of solutions. For that, in Section, we reformulate (1.1) in terms of the current, = m(u x + p). (1.3) From the second equation in (1.1), is constant. hus, the current becomes the main parameter in our analysis. While we focus our attention into non-increasing MFGs, our methods are also valid for increasing MFGs. o illustrate and contrast these two cases, we begin our analysis in Section 3 by addressing the latter. For >, we show the existence of a unique smooth solution. However, for =, we uncover new phenomena: the existence of non-smooth solutions and the lack of uniqueness.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 3 In Section 4, we consider the elliptic regularization of monotone MFGs. We establish a new variational principle that gives the existence and uniqueness of smooth solutions. Moreover, we address the vanishing viscosity problem using Γ-convergence. In Section 5, we study regular solutions of (1.1) for non-increasing g. In this case, if, m >. However, for certain values of, (1.1) does not have continuous solutions. In contrast, if is large enough, (1.1) has a unique smooth solution. Moreover, if V has a single point of maximum, there exists a unique solution of (1.1) for each >. If V has multiple maxima, there are multiple solutions. If =, the behavior of (1.1) is more complex and m can be discontinuous or vanish. Next, in Section 6, we consider MFGs with a decreasing nonlinearity, g, and discuss the properties of discontinuous viscosity solutions. Subsequently, in Section 7, we study the elliptic regularization of anti-monotone MFGs. here, we use calculus of variations methods to prove the existence of a solution. In Section 8, we examine the regularity of solutions as a function of the current and, in Section 9, we study the asymptotic behavior of solutions of (1.1) as converges to and. Finally, in Sections 1 and 11, we analyze the regularity of H in terms of and p.. he current formulation and regularization Here, we discuss the current formulation of (1.1) and (1.). After some elementary computations, we show that the current formulation of (1.) is the Euler-Lagrange equation of a suitable functional..1. Current formulation. Let be given by (1.3). From the second equation in (1.1), is constant. We split our analysis into the cases, and =. If, m(x) for all x and u x + p = /m. hus, (1.1) can be written as F (m) = H V (x), m >, mdx = 1, (.1) 1 m dx = p, where F (m) = g(m). For each x, the first equation in (.1) is an algebraic m equation for m. If g is increasing, for each x and H R, there exists a unique solution. In contrast, if g is not increasing, there may exist multiple solutions, as we discuss later. For =, (1.1) gives (u x+p) g(m) = H V (x), m, mdx = 1, (.) m(u x + p) =. From the last equation in (.), either m =, in which case u solves (u x + p) g() = H V (x), or m > and g(m) + H V (x) =. Hence, if g is increasing or decreasing, m(x) is determined in a unique way; otherwise, multiple solutions can occur.
4 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES.. Elliptic regularization. Now, we consider the elliptic MFG (1.). From the second equation in that system, we conclude that = ɛm x + m(p + u x ) is constant. hus, we solve for u x and replace it in the first equation. Accordingly, we get ( ) ɛmx ɛ + ( ɛm x) m m + V (x) = H + g(m). (.3) hen, using the identity x ɛ ( ɛm x)m x m + ( ɛm x) m we obtain the following equation for m: ɛ m xx m Now, let Φ be such that Φ (m) = F (m); that is, = ɛ m x m, ɛ m x m + F (m) = H V (x). (.4) Φ (m) = m G(m), where G (m) = g(m). hen, (.4) is the Euler-Lagrange equation of the functional ɛ m x m Φ (m) V (x)m dx (.5) under the constraint m = 1; the constant H is the Lagrangian multiplier for the preceding constraint. 3. First-order monotone MFGs We continue our analysis by considering monotonically increasing nonlinearities, g. In the case of a non-vanishing current, solutions are smooth. However, if the current vanishes, solutions can fail to be smooth, m can vanish, and u may not be unique. he non-smooth behavior for a generic non-decreasing nonlinearity, g, was observed in heorem.8 in [31] where the authors find limits of smooth solutions of second-order MFGs as the viscosity coefficient converges to. 3.1., g increasing. Here, in contrast to the case =, examined later, the solutions are smooth. Elementary computations give the following result. Proposition 3.1. Let g be monotonically increasing. hen, for every >, (1.1) has a unique smooth solution, (u, m, H ), with current. his solution is given by where p = x m (x) = F 1 (H V (x)), u (x) = m (y) dy p x, m (y) dy, F (t) = t g(t), and H is such that m (x)dx = 1.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 5 3.. =, g increasing. o simplify the discussion and illustrate our methods, we consider (.) with g(m) = m. he analysis is similar for other choices of an increasing function, g. Accordingly, we have (u x+p) m = H V (x); m, mdx = 1; (3.1) m(u x + p) =. It is easy to see that m(x) = (V (x) H) + for x. he map H (V (x) H) + dx is decreasing (strictly decreasing at its positive values). Hence, there exists a unique number, H, such that m(x)dx = 1. Moreover, H < max V and H V 1. If min V < H < max V, m is non-smooth and there are regions where it vanishes. In contrast, if H < min V, m is always positive. In this case, u x + p =, and, by periodicity of u, p =. Furthermore, from the first equation in (3.1), we have H = V (x)dx 1. Given H, we find from (3.1) that Hence, for u x + p = u ± (x) = ± (H V (x)) +, x. x (H V (y)) + dy px, where p = ± (H V (y)) + dy, the triplets (u ±, m, H) solve (3.1). However, there are also solutions with a discontinuous derivative, u x. For that, let x be such that V (x ) < H. Such a point always exists if H > min V or, equivalently, when V (x)dx 1 < min V. Let (u x (x)) x = (H V (x)) + χ x<x (H V (x)) + χ x>x p x, where p x = (H V (y)) + dy (H V (y)) + dy and χ denotes the y<x y>x characteristic function. herefore, u x solves the first equation of (3.1) almost everywhere, and u x x has only negative umps. Since m is continuous, u x is a viscosity solution of that equation. Consequently, (u x, m, H) solves (3.1). o summarize, (3.1) has a unique, smooth solution if and only if u x + p or, equivalently, m(x) = V (x) H. he latter holds if and only if V (x)dx 1 + min V. (3.) his is the case for small perturbations V ; that is, oscv 1. For A R, set V A (x) = A sin(π(x + 1 4 )) and let m(x, A), H(A) solve (3.1) for V = V A. In Fig. 1, we plot m(x, A) for A. We observe that m(x, A) is smooth for small values of A and becomes non-differentiable for large A, as expected from our analysis. If A =, (3.) does not hold. hus, m(x, ) is singular and we have multiple solutions, u(x, ). In Fig., we plot m(x, ) and two distinct solutions, u(x, ).
6 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES Fig. 1. m(x, A). m(x,).5 u(x,).5..4 1.5.3 1...5.1.5 1. 1.5. x.5 1. 1.5. x Fig.. m(x, ) (left) and two distinct solutions u(x, ) (right). 4. Monotone elliptic mean-field games o study (1.), we examine the variational problem determined by (.5). As before, for concreteness, we consider the case g(m) = m. In this case, (.5) becomes ( ) J ɛ [m] = ɛ m x m + m + m V (x)m dx. (4.1) he preceding functional is convex and, as we prove next, the direct method in the calculus of variations gives the existence of a minimizer on the set { } A = m W 1, () : m m = 1. Proposition 4.1. For each R, there exists a unique minimizer, m, of J ɛ [m] in A. Moreover, m > and solves ( ɛ m ) x m x m + m + H V (x) = for some constant H R. Proof. he uniqueness of a positive minimizer is a consequence of the strict convexity of J ɛ. he existence of a non-negative minimizer requires separate arguments for the cases and =. We first examine the case. We begin by taking a minimizing sequence, m n A. hen, there exists a constant, C >, such that (m n ) x + 1 dx C. m n m n hus, by Morrey s theorem, the functions m n are equi-hölder continuous of exponent 1. herefore, because m n = 1, this sequence is equibounded and, through
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 7 m(x). 1.5 1..5...4.6.8 1. x Fig. 3. Solution m of (1.) when g(m) = m, = 1, V (x) = sin(π(x + 1/4)) for ɛ =.1 (dashed) and for ɛ = (solid). some subsequence, m n m for some function m. Moreover, by Fatou s lemma, 1 dx C. m Suppose that min m = m(x ) =. hen, because m is Hölder continuous, we have m(x) C x x. However, 1 x x dx is not finite, which is a contradiction. hus, m is a strictly positive minimizer. Moreover, it solves the corresponding Euler-Lagrange equation. For =, we rewrite the Euler-Lagrange equation as ɛ (ln m) xx + m V (x) = H. (4.) Let P be the set of non-negative functions in L ( d ) and consider the map Ξ : P P defined as follows. Given η P, we solve the PDE ɛ w xx + η V (x) = H, where H satisfies the compatibility condition H = V dx 1, and w : R is such that e w dx = 1. An elementary argument shows that w is uniformly bounded from above and from below. Next, we set Ξ(η) = e w. he mapping Ξ is continuous and compact. Accordingly, by Schauder s Fixed Point heorem, there is a fixed point, m, that solves (4.). By the convexity of the variational problem (4.1), this fixed point is the unique solution of the Euler-Lagrange equation. Next, to study the convergence as ɛ, we investigate the Γ-convergence as ɛ of J ɛ. A simple modification of the arguments in [], Chapter 6, shows that where J[m] = J ɛ Γ J, ( ) m + m V (x)m dx, if m and m = 1. In Fig. 3, we observe numerical evidence for this Γ- convergence.
8 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES 5. Regular viscosity solutions in anti-monotone mean-field games Here, we investigate MFGs with decreasing g. o simplify, we assume that g(m) = m. However, our arguments are valid for a general decreasing g. In contrast with the monotone case, m may not be unique. Furthermore, m can be discontinuous and, thus, viscosity solutions of the Hamilton-Jacobi equation in (1.1) should be interpreted in the discontinuous sense. In this section, we are interested in regular discontinuous viscosity solutions; that is, solutions satisfying conditions a. and b. stated in the Introduction. Here, we examine existence, uniqueness, and additional properties of such solutions. In Section 6, we prove that these solutions are indeed discontinuous viscosity solutions. 5.1., g decreasing. o simplify the presentation, we consider >. With g(m) = m, (.1) becomes + m = H V (x); m m >, mdx = 1; (5.1) 1 m dx = p. he minimum of t /t +t is attained at t min = /3. hus, /t +t 3 /3 / for t >. herefore, a lower bound for H is H H cr = max V + 3/3, (5.) where the superscript cr stands for critical. he function t /t + t is decreasing on the interval (, t min ) and increasing on the interval (t min, + ). For any H satisfying (5.), let m and H m+ be the H solutions of (m ± + m± (x) = H V (x), H (x)) H with m (x) t H min m + (x). Due to (5.), H m and H m+ are well defined. H Furthermore, if (u, m, H) solves (1.1), then m(x) agrees with either m + (x) or H m (x), H almost everywhere in. Let m := m H cr and m + := m+ H cr. Note that m (x) m+ (x) for all x, and the equality holds only at the maximum points of V. Hence, m (x) < m+ (x) on a set of positive Lebesgue measure unless V is constant. he two fundamental quantities for our analysis are α + 1 () = m + (x)dx, α 1 (5.3) () = m (x)dx. If V is not constant, we have for >. α () < α + () Proposition 5.1. Suppose that x = is the single maximum of V. hen, for every >, there exists a unique number, p, such that (1.1) has a regular solution with a current level,. Moreover, the solution of (5.1), (u, m, H ), is unique and given as follows.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 9 i. If α + () 1, m (x) = m + H (x), u (x) = where p = dy m (y) and H is such that m (x)dx = 1. ii. If α () 1, m (x) = m H (x), u (x) = where p = dy m (y) and H is such that m (x)dx = 1. iii. If α () < 1 < α + (), we have that H = H cr, and m (x) = m (x)χ [,d ) + m + (x)χ [d,1), u (x) = where p = dy m (y) and d is such that x x dy m (y) p x, (5.4) dy m (y) p x, (5.5) x dy m (y) p x, (5.6) m (x)dx = d 1 m (x)dx + d m + (x)dx = 1. Proof. Case i. he function /t + t is increasing on the interval (t min, + ). herefore, H m + (x) is increasing for all x. Hence, the mapping H H m + (x)dx, H is increasing. By assumption, m + H cr (x)dx = m + (x)dx 1. herefore, there exists a unique H H cr such that m + (x)dx = 1. hus, (u H, m ) given by (5.4) is the unique solution of (1.1) with H = H and p = p. Case ii. he function /t +t is decreasing on the interval (, t min ). herefore, m (x) is decreasing in H for all x. Hence, the mapping H H m (x)dx H is decreasing. By assumption, m H cr (x)dx = m (x)dx 1. hus, there exists a unique number, H H cr, such that m (x)dx = 1. Hence, (u H, m ) given by (5.5) is the unique solution of (1.1) with H = H and p = p. Case iii. We first show that (1.1) does not have regular solutions for H > H cr. By contradiction, suppose that (1.1) has a regular solution, (u, m, H), for some H > H cr and p R. Evidently, m(x) = m + (x)χ H E + m (x)χ H \E for some subset E. Furthermore, inf (m+ (x) H m (x)) > (5.7) H
1 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES because H > H cr. Moreover, m(x)dx = E m + H (x)dx + \E m H (x)dx and m H (x)dx < m (x)dx < 1 < m + (x)dx < m + H (x)dx. herefore, neither E nor \E can be empty or have zero Lebesgue measure. Because E and \ E are not negligible, there exists a real number, e, such that for every ε >, (e ε, e) E and (e, e + ε) E c. According to (5.7), m has a negative ump, m(e ) m(e + ) <, at x = e. Hence, u x = /m p has a positive ump, m (e) >, at x = e. However, derivatives of regular solutions can only have negative umps and, thus, this contradiction m + (e) implies H = H cr. Next, we construct m and u and determine p. We look for a function m of the form { m (x), x [, d), m (x) = m + (x), x [d, 1). (5.8) Note that (5.8) is the only possibility for m because m can switch from m + to m only if there is no ump at the switching point; that is, m + and m are equal at that point, which only holds at maximum of V. hus, by periodicity, m can switch to m from m + only at x = and x = 1. It remains to choose d (, 1) such that m (x)dx = 1. Let φ(d) = 1 m (x)dx = d 1 m (x)dx + d m + (x)dx. Because φ() > 1 and φ(1) < 1 and because φ (d) = m (d) m+ (d) < for d (, 1), there exists a unique d (, 1) such that φ(d ) = 1. he triplet defined by (5.6), (u, m, H ), solves (1.1). By the previous proposition, if V has a single maximum point then, for every current, >, there exists a unique p and a unique triplet, (u, m, H ), that solves (5.1) for p = p. In contrast, as we show next, if V has multiple maxima and > is such that Case iii in Proposition 5.1 holds, there exist infinitely many solutions. Proposition 5.. Suppose that V attains a maximum at x = and at x = x (, 1). Let be such that α () < 1 < α + (). hen, there exist infinitely many numbers, p, and pairs, (u, m), such that (u, m, H cr ) is a regular solution of (1.1). Proof. We look for solutions of the form m d 1,d (x) = { m (x), x [, d 1) [x, d 1 ), m + (x), x [d 1, x ) [d, 1),
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 11 m(x). 1.5 1. m +.5...4.6.8 1. x Fig. 4. Solution m for =.1 and V (x) = 1 sin(π(x + 1/4)). m(x) 1.4 1. m - 1..998.996..4.6.8 1. x Fig. 5. Solution m for = 1 and V (x) = 1 sin(π(x + 1/4)). where < d 1 < x and x < d < 1. Note that m d 1,d has two discontinuity points. At these points, m d 1,d has positive umps. Hence, if we define x u d 1,d (x) = dy m d 1,d (y) pd1,d x, x, where p d 1,d = dy m d1,d (y), the triplet (ud 1,d, m d 1,d, H cr ) is a regular solution of (1.1) if m d 1,d (x)dx = 1. o determine d 1 and d, we consider the function φ(d 1, d ) = 1 We have that φ(, x ) = m d 1,d (x)dx = d 1 + m (x)dx + x 1 d 1 x m (x)dx + d 1 m + (x)dx d m + (x)dx, (d 1, d ) (, x ) (x, 1). m + (x)dx > 1 and φ(x, 1) = 1 m (x)dx < 1. Because φ is continuous, there exists a pair, (d 1, d ) (, x ) (x, 1), such that φ(d 1, d ) = 1. In fact, there are infinitely many such pairs. For arbitrary continuous curve γ [, x ] [x, 1] connecting the points (, x ) and (x, 1), there exists at least one pair, (d 1, d ) γ, such that φ(d 1, d ) = 1. o each such pair corresponds a triplet (u d 1,d, m d 1,d, H cr ) that is a regular solution of (1.1). Let V (x) = 1 sin(π(x + 1 4 )). Because V has a single maximum, Proposition 5.1 gives that (1.1) admits a unique regular solution for all values of >. In Figs. 4, 5, 6, we plot m for different values of. In Fig. 4, we plot m in the low-current
1 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES m(x) 1.5 1. m +.5 m -...4.6.8 1. x Fig. 6. Solution m for =.5 and V (x) = 1 sin(π(x + 1/4)). m(x). m(x) 1.5 m + m + 1.5 m + m + 1. 1..5 m - m -..4.6.8 1. x.5 m - m -..4.6.8 1. x Fig. 7. wo distinct solutions for =.5 and V (x) = 1 sin(4π(x + 1/8)). regime, =.1; that is, Case i in Proposition 5.1. As we can see, m is smooth as predicted by the proposition. In Fig. 5, we plot m in the high-current regime, = 1; that is, Case ii in Proposition 5.1. As before, we observe that m is smooth. Finally, in Fig. 6, we plot m for the intermediate-current regime, =.5; that is, Case iii in Proposition 5.1. As we can see, m is discontinuous. Next, we consider the potential V (x) = 1 sin(4π(x + 1 8 )) that has two maxima. By Proposition 5., we have infinitely many two-ump solutions. In Fig. 7, we plot two such solutions. 5.. =, g decreasing. Now, we examine the case when the current vanishes, and, thus, we consider the system (u x+p) + m = H V (x); m, mdx = 1; (5.9) m(u x + p) =. Suppose that (5.9) has a solution. Because m, we have H V (x) for x. hus, H max V. On the other hand, ( ) H V (x) dx mdx = 1. Consequently, H 1 + V. herefore, H max max V, 1 + V =: H. It turns out that H is the only possible value for H as we show next. Proposition 5.3. he MFG (5.9) does not have regular solutions for H > H.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 13 Proof. Suppose that H > H and that the triplet (u, m, H) is a regular solution of (5.9). If m(x) >, then u x (x) + p = and m = H V (x). If m(x) =, then (u x (x) + p) = (H V (x)). hus, on the set Z = {x : m(x) = }, u x (x) + p = (H V (x)) or u x (x) + p = (H V (x)). We have that ( H V (x) ) dx > 1. Hence, the set Z has a positive Lebesgue measure. Otherwise, m(x) = H V (x) everywhere, and thus m(x)dx > 1. Consequently, u x + p is either (H V (x)) or (H V (x)) on Z. Suppose that u x (x) + p takes the value (H V (x)) at some point x. Without loss of generality, we can assume that u x () + p = (H V ()). Let { } e = sup x (, 1) s.t. u x (x) + p = (H V (x)). hen, at x = e, the function u x + p has a ump of size (H V (e)) or (H V (e)). However, this is impossible because u x is a regular solution, and it cannot have positive umps. herefore, u x (x) + p takes only the values (H V (x)) and. But then, u x must have a positive ump from to (H V (x)) at some point, which also contradicts the regularity property. Now, we construct solutions to (5.9) with H = H. It turns out that if V has a large oscillation, then (5.9) has infinitely many regular solutions. Proposition 5.4. i. if 1 + We have that V max V, then the triplet (u, m, H ) with m (x) = H V (x), u (x) =, (5.1) solves (5.9) in the classical sense for p = ; ii. if max V > 1 + V, define m d 1,d (x) = { H V (x), x [d 1, d ],, x \ [d 1, d ], (5.11) and where (u d 1,d ) x (x) = u d 1,d (x) = x (u d 1,d ) x (y)dy, x, (5.1) (H V (x)) p d 1,d, x [, d 1 ), p d 1,d, x [d 1, d ], (H V (x)) p d 1,d, x (d, 1],
14 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES 1 8 6 4..4.6.8 1. Fig. 8. m as defined in (5.11) for V (x) = 5 sin(π(x + 1 4 )) with d =.5 and d 1 such that (5.13) holds. u[x] 1. 1. ux[x] 4.8.6.4..4.6.8 1. x. -..4.6.8 1. x -4 Fig. 9. u (left) and (u ) x (right) as defined in (5.11) for V (x) = 5 sin(π(x + 1 4 )) with d =.5 and d 1 such that (5.13) holds. and p d 1,d = d 1 (d 1, d ), such that (H V (x))dx d 1 d (H V (x))dx. hen, for any pair, d 1 (H V (x))dx = 1, (5.13) the triplet (u d 1,d, m d 1,d, H ) is a regular solution for (5.9) for p = p d 1,d. Furthermore, there exist infinitely many pairs, (d 1, d ), such that (5.13) holds. Proof. Case i. In this case, H = 1 + V (x)dx and straightforward computations show that (5.1) defines a classical solution of (5.9). Case ii. In this case, we have that H = max V and that 1 (H V (x))dx > 1. Without loss of generality, we assume that is a point of maximum for V. Note that (u d 1,d ) x has only negative umps and u d 1,d satisfies (5.9) almost everywhere. hus, the triplet (u d 1,d, m d 1,d, H ) is a regular solution of (5.9) if m d 1,d 1 (x)dx = 1. However, the latter is equivalent to (5.13). Since (H V (x))dx > 1, we can find infinitely many such pairs. We find p d 1,d from the identity (u d 1,d ) x (x)dx =. Figs. 8 and 9 show the solutions of (5.9) for V (x) = 5 sin(π(x + 1 4 )), x. Remark 5.5. If V has multiple maxima and Case ii in Proposition 5.4 holds, there is a larger family of solutions. Let x = x (, 1) be a point of maximum for V.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 15 For fixed real numbers, d 1 < d < e 1 < e, define and m d 1,d,e 1,e (x) = u d 1,d,e 1,e (x) = { H V (x), x [d 1, d ] [e 1, e ],, elsewhere, x (5.14) (u d 1,d,e 1,e ) x (y)dy, x, (5.15) where (H V (x)) p d 1,d,e 1,e, x [, d 1 ) [x, e 1 ), (u d1,d,e 1,e ) x (x) = (H V (x)) p d 1,d,e 1,e, x (d, x ] (e, 1], p d 1,d,e 1,e, elsewhere. Note that (u d1,d,e 1,e ) x (x) is periodic, only has negative umps, and solves (5.9) almost everywhere. Hence, the triplet (u d 1,d,e 1,e, m d 1,d,e 1,e, H ) is a regular solution of (5.9) if for p d 1,d,e 1,e = [,d 1 ) [x,e 1 ) 1 he equality (5.16) is equivalent to Since 1 d m d1,d,e 1,e (x)dx = 1 (5.16) (H V (x))dx d 1 (H V (x))dx + e (d,x ] (e,1] (H V (x))dx. e 1 (H V (x))dx = 1. (5.17) (H V (x))dx > 1, we can find infinitely many quadruples (d 1, d, e 1, e ) such that (5.17) holds. Hence, we can generate infinitely many solutions of the form (5.14), (5.15). From Propositions 5.1 and 5.4, for every regular solution, m, of (1.1) in the lowcurrent regime ( = or Case i in Proposition 5.1), the smaller V (x) is, the larger m(x) is. his is paradoxical because V (x) represents the spatial preference of the agents and preferred regions correspond to high values of V. hus, areas that are less desirable have a high populational density. herefore, it is possible that the most preferred site is empty and agents aggregate at the least preferred site. For example, in (5.11), m vanishes near the maximum of V and is supported in the neighborhood of the minimum of V, as illustrated in Fig. 8. Hence, if agents do not move fast (low current), they prefer staying together rather than being in a better place, see Fig. 4. In the high-current regime (Case ii in Proposition 5.1), the opposite situation occurs: the larger V (x) is, the larger m(x) becomes, see Fig. 5. herefore, preferred areas have a high population density. Hence, if the level of the current is high enough (we give quantitative estimates in the next section), agents are better off at preferred sites and with more agents. Finally, for the intermediate current level (Case iii in Proposition 5.1), we observe a more complex situation. he solution, m, consists of two parts: m (x) and m+ (x). m (x) is larger where V (x) is larger and the opposite holds for m + (x). herefore, in the region where m is m,
16 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES the most preferred sites are more densely populated. In the region where m is m + the less preferred sites are more densely populated. his is illustrated in Fig. 6. 6. Discontinuous viscosity solutions In the anti-monotone case considered in the preceding section, m can be discontinuous. hus, in addition to regular solutions examined before, we need to consider viscosity solutions in the framework of discontinuous Hamiltonians. In what follows, we recall the main definitions in [1]. Given a locally bounded function, F : R R, we define its lower and upper semicontinuous envelopes as F (x, q) = lim inf F (y, r), F (x, q) = lim sup F (y, r) (y,r) (x,q) (y,r) (x,q) for (x, q) R. We say that a locally bounded function, u : R, is a viscosity solution of F (x, Du) = if, for any smooth function, φ : R, we have that and F (x, φ x + p) for all x argmax(u φ) F (x, φ x + p) for all x argmin(u φ). Let m : R, m L (), and set m (x) = lim inf m(y), y x m (x) = lim sup m(y). y x Suppose that V : R is continuous. hen, for our setting, we have Consequently, F (x, q) = q + V (x) + m(x) H. F (x, q) = q + V (x) + m (x) H, F (x, q) = q + V (x) + m (x) H. Here, we look for piecewise smooth solutions of (1.1) for g(m) = m that are not necessarily regular; that is, the condition lim u x(x) lim u x(x) is not necessarily x x x x + satisfied. It turns out that there are infinitely many such solutions for all independent of properties of V, and the ump direction of u x is irrelevant. his contrasts with the fact that for V with a single maximum, there exists ust one regular solution (Proposition 5.1). hus, we select a current level, > ( < is analogous), and fix arbitrary points x < x 1 < < x n 1 and H H cr. We search for solutions (u, m, H) such that m is continuous on the intervals (x i, x i+1 ) for i n 1. From the above discussion, we have: Proposition 6.1. Assume that > and that m > is continuous on (x i, x i+1 ) and m(x) + m(x) = H V (x) for all x x i. H is such that m(x)dx = 1. hen, the triplet (u, m, H) solves (1.1), where u(x) = x dy px, p = m(y) m(y) dy.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 17 Proof. We have that u x + p = m a.e.. hus, the second equation in (1.1) holds in the sense of distributions. Next, we observe that u is differentiable for all x x i and that the first equation in (1.1) is satisfied in the classical sense at those points. hus, we ust need to check the viscosity condition at x = x i. here are two possible cases: 1. m(x i ) > m(x+ i ). In this case, m (x i ) = m(x i ). Moreover, u x (x i ) = /m(x i ) p < /m(x+ i ) p = u x(x + i ). Hence, there is no smooth function touching u from above; it touches only from below. herefore, we need to check that, for any φ touching u from below at x i, we have (φ x (x i ) + p) + V (x i ) + m(x i ) H. Because (1.1) is satisfied at x x i in the classical sense, we have that (u x (x ± i ) + p) + V (x i ) + m(x ± i ) H =. Because φ touches u from below and >, we have Hence, < u x (x i ) + p φ x(x i ) + p u x (x + i ) + p. (φ x (x i ) + p) + V (x i ) + m(x i ) H (u x(x i ) + p). m(x i ) < m(x i +). In this case, m (x i ) = m(x i ) and + V (x i ) + m(x i ) H =. u x (x i ) = /m(x i ) p > /m(x+ i ) p = u x(x + i ). Hence, there is no smooth function touching u from below only from above. herefore, for any φ touching u from above at x i, we have (φ x (x i ) + p) + V (x i ) + m(x i ) H. (6.1) Because (1.1) holds in the classical sense for x x i, we have that (u x (x ± i ) + p) + V (x i ) + m(x ± i ) H =. Because φ touches u from above, we have < u x (x + i ) + p φ x(x i ) + p u x (x i ) + p. Hence, (6.1) holds. Remark 6.. If g is increasing, the construction of piecewise smooth solutions with discontinuous m in the previous proposition fails because m(x i ) = m(x+ i ), necessarily. herefore, the smooth solutions found in Proposition 3.1 are the only possible ones: there are no extra discontinuous solutions as in the case of decreasing g. his is yet another consequence of the regularizing effect of an increasing g.
18 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES 7. Anti-monotone elliptic mean-field games Now, we consider anti-monotone elliptic MFGs and the corresponding variational problem (.5) with g(m) = m. We use the direct method in the calculus of variations to prove the existence of a minimizer of the functional ( ) J ɛ [m] = ɛ m x m + m m V (x)m dx. (7.1) Proposition 7.1. Moreover, m, solves for some H R. For each R, there exists a minimizer, m, of J[m] in { } A = m W 1, () : m > m = 1. ( ɛ m ) x m x m m + H V (x) = Proof. o prove the existence of a positive minimizer, we consider separately the cases and =. Case 1.. We take a minimizing sequence, m n A, and note that there is a constant, C, such that ɛ (m n) x + m n dx C + m n m n. herefore, we seek to control m by the integral expression on the left-hand side. For that, we recall the Gagliardo-Nirenberg inequality, for a 1, with w L p w x a L r w 1 a L q (7.) ( ) 1 1 p = a r 1 + (1 a) 1 q, whenever w =. With p = 4 and r = q =, we obtain a = 1 4. Using these values in (7.), taking into account that m = 1, and choosing w = m, we obtain ( m n C + C (m n ) ) 1 x. m n hus, using a weighted Cauchy inequality, ɛ (m n) x + dx C. m n m n Finally, we argue as in the proof of Proposition 4.1 and show the existence of a minimizer. Case. =. Here, we use a fixed point argument as in the proof of Proposition 4.1. For that, we rewrite the Euler-Lagrange equation as ɛ (ln m) xx m V (x) = H, (7.3) and argue as before. However, because the functional (7.1) is non-convex, uniqueness may fail.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 19 m(x). 1.5 1..5...4.6.8 1. x Fig. 1. Solution m of (1.) when g(m) = m, = 1, V (x) = sin(π(x + 1/4)) for ɛ =.1 (dashed) and for ɛ = (solid). 1 he preceding result does not give a unique minimizer. We note that for large, the functional (7.1) behaves like a convex functional. Finally, we note that as ɛ, numerical evidence suggests that there is no Γ-convergence to a minimizer, see figure 1, where we plot a solution for small ɛ versus the solution with ɛ =. 8. Regularity regimes of the current equation for g(m) = m Now, we analyze the regularity regimes of (5.1); that is, we determine for which values of (5.1) has or fails to have smooth solutions. For simplicity, we assume that is the only point of maximum of V. Moreover, as before, we consider the case, as the case < is analogous. We begin by proving that α +, α, defined in (5.3), are monotone. Proposition 8.1. We have that i. α + and α are increasing on (, ); ii. lim + α+ () = lim + α () = ; iii. lim α + () = max V V (x)dx, lim α () =. Proof. i. First, we prove that m + (x) and m (x) (see Section 5.1 for the definition) are increasing in at every point x. We fix x and set h = (max V ) V (x). If h =, then m (x) = /3, which is an increasing function of. Next, for h > let t() = m (x) < /3. We have that t() + t() 3 /3 = h. By the implicit function theorem, t() is differentiable. Differentiating the previous equation in gives t () = 1/3 t. 1 t 3 Because < t() < /3, t () >. Hence, t() is increasing. he proof for m + (x) is identical. ii. By definition, m + (x) /3. Hence, lim α + () =. On the other hand, for large enough, we have ( /3 /) + /3 3 /3 = /3 > max V V (x) = (m (x)) + m (x) 3 /3.
DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES herefore, m (x) > /3 / and lim α () =. iii. Because m (x) /3 for every x, lim m (x) = for all x. hus, lim α () =. On the other hand, m + (x) /3. hus, herefore, ( ) 3 lim m+ (x) = lim /3 (m + + max (x)) hus, and lim α+ () = max V V (x)dx. (m + (x)) /3. = max V V (x). Next, we define two numbers that characterize regularity regimes of (1.1): Proposition 8.. lower = inf{ > s.t. α + () > 1}, (8.1) upper = inf{ > s.t. α () > 1}. (8.) Let lower and upper be given by (8.1) and (8.). hen i. lower < upper < ; ii. for upper, the system (1.1) has smooth solutions; iii. for lower < < upper, the system (1.1) has only discontinuous solutions; iv. if lower >, the system (1.1) has smooth solutions for < lower. Proof. he proof is a straightforward application of Propositions 5.1 and 8.1. Finally, we characterize the regularity at =. Proposition 8.3. he system (5.9) admits smooth solutions if and only if α + () 1. Proof. he proof follows from iii in Proposition 8.1 and i in Proposition 5.4. Let V (x) = A sin(π(x + 1/4)). In Fig. 11, we plot α + and α for A =.5 and A = 5. From Proposition 8.1, α + () = A. hus, if A =.5, we have α + () < 1 and, for A = 5, we have α + () > 1. herefore, lower > for A =.5 and lower = for A = 5. Hence, if A =.5, (1.1) has smooth solutions for a low enough current level (.18) or for a high enough current level ( 1.75). In contrast, if A = 5, there are no smooth solutions for low currents, only for large currents ( 3.3). We end the section with an a priori estimate for the current level for smooth solutions. Proposition 8.4 (A priori estimate). Suppose that max V > 1 + V (x)dx and let (u, m, H) be a smooth solution of (1.1) with m >. hen, there exists a constant, c(v ) >, such that inf m(u x + p) c(v ). (8.3) Proof. From Proposition 8.1, we have that α + () > 1. hus, by Proposition 8.3, (1.1) does not have smooth solutions for =. Additionally, lower =. Next, take c(v ) = upper and by Proposition 8., we conclude (8.3).
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 1 α A=.5 α A=5 4 α + 15 α + 3 1 α - 1 5 α - 1 3 4 5 5 1 15 Fig. 11. α + and α for V (x) = A sin(π(x + 1/4)). lower =.18, upper = 1.75 (A =.5); lower =, upper = 3.3 (A = 5). he previous Proposition shows that if the potential, V, has a large oscillation (this happens in the example for A = 5, Fig. 11), then only high current solutions are smooth. 9. Asymptotic behavior of solutions as and + In Section 5.1, we studied regular solutions of (1.1) with a current level >. Here, we continue the analysis of the decreasing nonlinearity, g(m) = m, and examine the asymptotic behavior of regular solutions as and. As before, we assume that V has a single maximum at. First, we address the case. Proposition 9.1. For >, let (u, m, H ) solve (5.1). We have that i. lim H = ; ii. For x, lim m (x) = 1, lim u (x) =, and lim p =. Proof. i. According to (5.), we have that H 3/3 +max V. hus, lim H =. ii. For upper, solutions of (5.1) are given by (5.5). Hence, m consists only of the m branch. hus, m (x) /3, which yields m (x) + m (x) m (x) m (x). herefore, 3 m. Consequently, using this inequality in (5.1), we get (x) m (x) (H V (x)) 3. (9.1) (H V (x)) Integrating the previous inequality and taking into account that m (x) = 1, we get dx 1 (H V (x)) 3 dx. (9.) (H V (x)) Because H converges to, and, for every x, y, where V is bounded, we have that (H V (y)) lim = 1. (H V (x)) Hence, for large enough, we have (H V (x)) (H V (y)), x, y. (9.3)
DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES Let x be such that = (H V ( x)) hen, by (9.1), (9.3), and (9.), we get 3 m (x) (H V (x)) Similarly, we have m (x) (H V (x)) Furthermore, we have that hus, m (x) = dx. (H V (x)) 3 (H V ( x)) (H V ( x)) 3. 1 3. m (x) = H V (x) 1 + m3 (x). (9.4) 1 1 + m3 (x). (9.5) (H V (x)) Finally, because m is bounded and its integral is 1, we get from (9.5) that lim (H V (x)) = 1 (9.6) for all x. he preceding limit implies that lim m (x) = 1 for all x. In fact, (9.6) gives precise asymptotics of H, namely H lim = 1. (9.7) Now, we compute the limit of u (x). We have that (u ) x = m (y) dy. From (5.1), we have that m (x) p, where p = m (x) = (H V (x) m (x)). herefore, using (9.4), m (x) m (y) (H = V (x) m (x)) (H V (y) m (y)) m (x) m (y) + V (x) V (y) (H min V ) 3 + oscv (H min V ), as. Hence, ( (u ) x = m (x) p = m (x) ) dy m (y) m (x) m (y) dy, when. Consequently, lim u (x) = lim x (u ) x (y)dy =.
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 3 Next, we study the behavior of solutions as. Proposition 9.. i. lim H = max ii. if 1 + V > max lim m (x) = 1 + We have that ( max V, 1 + V V, then for all x ; iii. if 1 + V max V, then V V (x), ) = H ; lim u (x) =, and lim p = lim m (x) = m d,1 (x), d lim u (x) = u d,1 (x), and lim p = for all x, where m d,1 and u d,1 are given by (5.11) and (5.1). (max V V (x))dx Proof. i. here are two possible cases: lower > and lower =. If lower =, then α + () > 1 for all > and α () < 1 for small enough. Hence, by the results in Section 5.1, we have that H = H cr = 3 /3 + max V. hus, lim H = max V. On the other hand, lower = means that lim α + () 1. Consequently, by Proposition ( ) (8.1), max V V 1. hus, max max V, 1 + V = max V = lim H. If lower >, then α + () < 1 for < lower and solutions (m, u, H ) are given by (5.4). Hence, m (x) /3 and herefore, lim H = lim < V + /3 m (x). (9.8) m + m (x) = V + 1. But max V ( V = lim α + () < 1, so max max V, 1 + ) V = 1 + V = lim H. ii. Since lim α + () = max V V, we have that the condition 1+ V > max V is equivalent to the condition lower >. In this case, we have that lim H = 1 + V. herefore, from (9.8), we have that Furthermore, lim m (x) = lim ( lim u (x) = lim H V (x) x m (y) ) m (x) = 1 + V V (x). m (z) dz dy =.
4 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES iii. he inequality 1 + V max V is equivalent to lower =. Hence, for < < upper solutions are given by (5.6). Because < m (x) /3, lim m (x) =. Furthermore, m+ (x) /3. hus, herefore, lim m+ (x) = lim ( lim (m + =. (x)) H V (x) (m + (x)) ) = max V V (x). Suppose that the ump points, d, of m (x) (see (5.6)) converge to some d [, 1] through a subsequence. hen, through that subsequence lim m (x) = m d,1 (x), where m d,1 is defined in (5.11). Hence, 1 1 = m d,1 (x)dx = (max V V (x)) dx. d Because V has a single maximum, d is defined uniquely by the previous equation. Hence, lim d = d and lim m (x) = m d,1 (x), globally (not only through some subsequence). Consequently, lim u (x) = u d,1 where d is such that 1 d d (x), lim p = p d,1 = (max V V (x)) = 1. (max V V (x))dx, From Proposition 9., we see that we recover only part of the solutions for = as limits of solutions for >. If we consider the solutions of (5.1) for which m takes negative values, we recover all solutions described in Section 5.. Indeed, the first equation in (5.1) is a cubic equation in m(x). hus, for every x, there are three solutions: two positive and one negative. Because we are interested in the MFG interpretation of (5.1), we neglect solutions with negative m. However, we can construct solutions for (5.1) without the constraint m >. As converges to, the negative parts of these solutions converge to, and, in the limit, we obtain all non-negative solutions of (5.9) given in Proposition 5.4. 1. Properties of H In this section, we study various properties of the effective Hamiltonian, H, as a function of. In the following proposition, we collect several properties of H. Proposition 1.1. We have that i. For every R, there exists a unique number, H, such that (1.1) has solutions with a current level ; ii. H is even; that is, H = H ; iii. H is continuous; iv. H increasing on (, ) ( and decreasing on (, ); v. min H = H = max max V, 1 + ) V (x)dx ; R
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 5 H 5 4 3 upper lower.5 1. 1.5..5 3. Fig. 1. H for V (x) = 1 sin(π(x + 1 4 )). vi. lim H / = 1. Proof. i. his follows from Propositions 5.1 and 5.4. ii. his follows from the fact that /t + t is an even function for all t >. iii. Continuity of H follows from the continuity of the mapping (, t) /t +t for, t >. iv. Since H is even, it suffices to show that it is increasing on (, + ). First, we show that H is increasing on ( upper, ). For that, we fix > upper. We have that H H cr. Hence, for any upper < < we have that H H cr > H cr. herefore, the function, m, determined by { ( m + m (x)) (x) = H V (x), (1.1) m (x) /3, is well defined for all upper < <. Next, we show that the mapping m (x) is increasing in ( upper, ) for all x. Indeed, fix x and differentiate (1.1) in to obtain d m (x) d = m (x) m (x) 3 <. Hence, m (x) < m (x), x. Accordingly, m (x)dx < m (x)dx = 1. Finally, the previous inequality implies H < H. he monotonicity of H on (, lower ) (in the case lower > ) can be proven analogously. Next, for lower < < upper, we have that H = H cr = 3 /3 + max V. H is thus evidently monotone. v. his follows from the previous properties of H and Proposition 9.. vi. We have proven this in (9.7). In Fig. 1 we plot H as a function of for V (x) = 1 sin(π(x + 1 4 )). 11. Analysis in terms of p Now, we analyze (1.1) in terms of the variable p. If g(m) is increasing, for every p R, there exists a unique number, H(p), for which (1.1) has a solution. his solution is unique if m > (see, e.g., [31]). Here, we show that, if g(m) is not increasing,
6 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES p 3..5. 1.5 upper 1..5 lower.5 1. 1.5..5 3. Fig. 13. p for V (x) = 1 sin(π(x + 1 4 )). there may be different values of H(p) for which (1.1) has a regular solution. he uniqueness of H depends both on the monotonicity of g and on the properties of V. For example, if g(m) = m, H is uniquely determined by p if and only if V has a single maximum. Moreover, our prior characterization of regular solutions of (1.1) implies that, for V with a single maximum point, (1.1) admits a unique regular solution for every p R. We start with an auxiliary lemma. Lemma 11.1. Let x = be the single maximum point of V. i. For every, there exists a unique number, p, such that (1.1) has a regular solution. Furthermore, the map p is increasing on (, ) and (, ). ii. If 1 + V (x)dx max V, then p = is the unique number for which (1.1) has a regular solution with =. Moreover, lim p =. iii. If 1 + V (x)dx < max V, then and p > d 1 1 p < d where d 1, d (, 1) are such that d 1 (max V V (x))dx, > (max V V (x))dx, <, 1 (max V V (x))dx = (max V V (x))dx = 1. d Consequently, (1.1) has a regular solution for = if and only if 1 d (max V V (x))dx p d 1 (max V V (x))dx. (11.1) Proof. i. According to Proposition 5.1, for every >, there exists a unique number, p, such that (1.1) has a regular solution with a current level. Let (u, m, H ) be the solution of (1.1) given by (5.4), (5.5) or (5.6). Because p = m (y) dy,
ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 7 to prove that p is increasing it suffices to show that m (x) is increasing for all m (x). x. First, we prove the monotonicity for lower < < upper. Let n (x) := We have that n (x) + n (x) = H V (x), x ; (11.) n (x) = m (x)χ [,d ) + m + (x)χ [d,1). Because the maps m + (x) and m (x) are increasing for all x, the map d is also increasing. Assume that is such that d x. We differentiate in the first equation in (11.) and take into account that H = 3 /3 + max V for lower < < upper to get dn (x) d = H 1 n (x) n (x) = n (x) 1 1 1/3 n (x). n (x) n (x) Let x be such that x = d x. For > x, we have d > x. hus, n (x) = /m (x) > 1/3, which implies dn (x) d >. Similarly, for < x, we have d < x. herefore, n (x) = /m + (x) < 1/3, which implies dn (x) d >. Next, we analyze the behavior of n at x. For > x, n (x) = m (x), and, for < x, n (x) = m + (x). hus, n (x) takes a positive ump, >, m (x) m + (x) at = x. herefore, n (x) has positive derivatives whenever x and a positive ump at = x. It is thus increasing for lower < < upper. Next, we show that n (x) is increasing on ( upper, ). As before, we have dn (x) d = H 1 n (x) n (x). n (x) Because m (x) < /3, we have n (x) > 1/3. herefore, if H 1/n (x), the map n (x) is increasing. Fix and, for >, consider H := H + ( ) min m (x). Define { m + m (x) (x) = H V (x), x ; m (x) /3. Note that H = H and m = m. Now, we compute the derivative of the map m (x) at =. Because m < /3, we have d m (x) d = = H m 1 m 3 hus, d m (x)dx = d = Hence, for small > close to, we get = = m (x)dx > 1. min m m 1 m 3 >. d m (x) dx >. d
8 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES H(p) 5 4 3 p lower 1 p upper.5 1. 1.5..5 3. p Fig. 14. H(p) for V (x) = 1 sin(π(x + 1 4 )). Consequently, for those values of the current, we have that H H. Hence, H H = min m = max 1 n, which completes the monotonicity proof for > upper. he monotonicity for < lower is similar. ii. & iii. hese claims follow from the monotonicity of p and Propositions 5.4 and 9.. In Fig. 13, we plot p as a function of for V (x) = 1 sin(π(x + 1 4 )). Proposition 11.. Let x = be the only point of maximum of V. hen, i. for every p R, there exists a unique number, H(p), for which (1.1) has a regular solution; ii. for every p R, (1.1) has a unique regular solution; iii. if max V > 1 + V (x)dx, H(p) is flat at the origin; iv. H(p) is increasing on (, ) and decreasing on (, ). hus min H(p) = H() = max max V, 1 + V (x)dx ; p R H(p) v. lim p p / = 1. Proof. i & ii. From Lemma 11.1, we have that for every p, there exists a unique such that (1.1) has regular solutions. From Proposition 1.1, we have that for every there exists a unique number, H, such that (1.1) has a regular solution. herefore, for every p, the constant H is determined uniquely. Moreover, from Proposition 5.1, we have that (1.1) has a unique regular solution for this constant. iii. From iii. in Lemma 11.1, we have that if p satisfies (11.1), then H(p) = H. iv. his follows from i in Lemma 11.1 and iv and v in Proposition 1.1. v. his follows from vi in Proposition 1.1, ii in Proposition 9.1, and the formula p = m (y) dy. In Fig. 14, we show H(p) for V (x) = 1 sin(π(x + 1 4 )). Remark 11.3. We conecture that, if V has only one maximum point, H(p) is convex. Let p 1 < p and (u 1, m 1, H(p 1 )) and (u, m, H(p )) solve (1.1) for p = p 1 and p = p, respectively. Consider the traectories, y 1, y, determined by ẏ i (t) = (u i ) x (y i (t)) p i, y i () =, i = 1,.