arxiv: v2 [math.ap] 28 Nov 2016

Size: px
Start display at page:

Download "arxiv: v2 [math.ap] 28 Nov 2016"

Transcription

1 ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv: v [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, 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

2 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.

3 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 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.

5 ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES =, 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 )) 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 6 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES Fig. 1. m(x, A). m(x,).5 u(x,) x 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

7 ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 7 m(x) 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 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.

9 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

10 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),

11 ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 11 m(x) m x Fig. 4. Solution m for =.1 and V (x) = 1 sin(π(x + 1/4)). m(x) m 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 )). 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

12 1 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES m(x) m +.5 m x Fig. 6. Solution m for =.5 and V (x) = 1 sin(π(x + 1/4)). m(x). m(x) 1.5 m + m m + m m - m x.5 m - m 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 )) 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.

13 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 14 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES Fig. 8. m as defined in (5.11) for V (x) = 5 sin(π(x )) with d =.5 and d 1 such that (5.13) holds. u[x] ux[x] x x -4 Fig. 9. u (left) and (u ) x (right) as defined in (5.11) for V (x) = 5 sin(π(x )) 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 )), 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.

15 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 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 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.

17 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 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.

19 ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 19 m(x) 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.

20 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).

21 ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 1 α A=.5 α A=5 4 α + 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)

22 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)) m (x) = H V (x) 1 + m3 (x). (9.4) 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 =.

23 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 =.

24 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 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

25 ONE-DIMENSIONAL, FIRS-ORDER, SAIONARY MEAN-FIELD GAMES 5 H upper lower Fig. 1. H for V (x) = 1 sin(π(x )). 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 )). 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,

26 6 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES p upper 1..5 lower Fig. 13. p for V (x) = 1 sin(π(x )). 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 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,

27 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

28 8 DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES H(p) p lower 1 p upper p Fig. 14. H(p) for V (x) = 1 sin(π(x )). 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 )). 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 )). Remark 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,.

Stationary mean-field games Diogo A. Gomes

Stationary mean-field games Diogo A. Gomes Stationary mean-field games Diogo A. Gomes We consider is the periodic stationary MFG, { ɛ u + Du 2 2 + V (x) = g(m) + H ɛ m div(mdu) = 0, (1) where the unknowns are u : T d R, m : T d R, with m 0 and

More information

Price formation models Diogo A. Gomes

Price formation models Diogo A. Gomes models Diogo A. Gomes Here, we are interested in the price formation in electricity markets where: a large number of agents owns storage devices that can be charged and later supply the grid with electricity;

More information

Sobolev Spaces. Chapter Hölder spaces

Sobolev Spaces. Chapter Hölder spaces Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect

More information

1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),

1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989), Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer

More information

Elliptic Kirchhoff equations

Elliptic Kirchhoff equations Elliptic Kirchhoff equations David ARCOYA Universidad de Granada Sevilla, 8-IX-2015 Workshop on Recent Advances in PDEs: Analysis, Numerics and Control In honor of Enrique Fernández-Cara for his 60th birthday

More information

Math212a1413 The Lebesgue integral.

Math212a1413 The Lebesgue integral. Math212a1413 The Lebesgue integral. October 28, 2014 Simple functions. In what follows, (X, F, m) is a space with a σ-field of sets, and m a measure on F. The purpose of today s lecture is to develop the

More information

2 A Model, Harmonic Map, Problem

2 A Model, Harmonic Map, Problem ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or

More information

Global minimization. Chapter Upper and lower bounds

Global minimization. Chapter Upper and lower bounds Chapter 1 Global minimization The issues related to the behavior of global minimization problems along a sequence of functionals F are by now well understood, and mainly rely on the concept of -limit.

More information

New Identities for Weak KAM Theory

New Identities for Weak KAM Theory New Identities for Weak KAM Theory Lawrence C. Evans Department of Mathematics University of California, Berkeley Abstract This paper records for the Hamiltonian H = p + W (x) some old and new identities

More information

Optimal Control. Macroeconomics II SMU. Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112

Optimal Control. Macroeconomics II SMU. Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112 Optimal Control Ömer Özak SMU Macroeconomics II Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112 Review of the Theory of Optimal Control Section 1 Review of the Theory of Optimal Control Ömer

More information

Mean field games and related models

Mean field games and related models Mean field games and related models Fabio Camilli SBAI-Dipartimento di Scienze di Base e Applicate per l Ingegneria Facoltà di Ingegneria Civile ed Industriale Email: Camilli@dmmm.uniroma1.it Web page:

More information

We denote the space of distributions on Ω by D ( Ω) 2.

We denote the space of distributions on Ω by D ( Ω) 2. Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace s Equation. Chapter Mean Value Formulas Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic

More information

Example 1. Hamilton-Jacobi equation. In particular, the eikonal equation. for some n( x) > 0 in Ω. Here 1 / 2

Example 1. Hamilton-Jacobi equation. In particular, the eikonal equation. for some n( x) > 0 in Ω. Here 1 / 2 Oct. 1 0 Viscosity S olutions In this lecture we take a glimpse of the viscosity solution theory for linear and nonlinear PDEs. From our experience we know that even for linear equations, the existence

More information

Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations

Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Course 212: Academic Year Section 1: Metric Spaces

Course 212: Academic Year Section 1: Metric Spaces Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........

More information

REGULARITY FOR INFINITY HARMONIC FUNCTIONS IN TWO DIMENSIONS

REGULARITY FOR INFINITY HARMONIC FUNCTIONS IN TWO DIMENSIONS C,α REGULARITY FOR INFINITY HARMONIC FUNCTIONS IN TWO DIMENSIONS LAWRENCE C. EVANS AND OVIDIU SAVIN Abstract. We propose a new method for showing C,α regularity for solutions of the infinity Laplacian

More information

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological

More information

September Math Course: First Order Derivative

September Math Course: First Order Derivative September Math Course: First Order Derivative Arina Nikandrova Functions Function y = f (x), where x is either be a scalar or a vector of several variables (x,..., x n ), can be thought of as a rule which

More information

Lecture 4 Lebesgue spaces and inequalities

Lecture 4 Lebesgue spaces and inequalities Lecture 4: Lebesgue spaces and inequalities 1 of 10 Course: Theory of Probability I Term: Fall 2013 Instructor: Gordan Zitkovic Lecture 4 Lebesgue spaces and inequalities Lebesgue spaces We have seen how

More information

The Lebesgue Integral

The Lebesgue Integral The Lebesgue Integral Brent Nelson In these notes we give an introduction to the Lebesgue integral, assuming only a knowledge of metric spaces and the iemann integral. For more details see [1, Chapters

More information

SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction

SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.

More information

THE INVERSE FUNCTION THEOREM

THE INVERSE FUNCTION THEOREM THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)

More information

An introduction to Mathematical Theory of Control

An introduction to Mathematical Theory of Control An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018

More information

Continuity. Matt Rosenzweig

Continuity. Matt Rosenzweig Continuity Matt Rosenzweig Contents 1 Continuity 1 1.1 Rudin Chapter 4 Exercises........................................ 1 1.1.1 Exercise 1............................................. 1 1.1.2 Exercise

More information

Math 410 Homework 6 Due Monday, October 26

Math 410 Homework 6 Due Monday, October 26 Math 40 Homework 6 Due Monday, October 26. Let c be any constant and assume that lim s n = s and lim t n = t. Prove that: a) lim c s n = c s We talked about these in class: We want to show that for all

More information

Conservation law equations : problem set

Conservation law equations : problem set Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions

More information

McGill University Math 354: Honors Analysis 3

McGill University Math 354: Honors Analysis 3 Practice problems McGill University Math 354: Honors Analysis 3 not for credit Problem 1. Determine whether the family of F = {f n } functions f n (x) = x n is uniformly equicontinuous. 1st Solution: The

More information

REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS

REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,

More information

ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COUPLED SYSTEMS OF PDE. 1. Introduction

ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COUPLED SYSTEMS OF PDE. 1. Introduction ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COPLED SYSTEMS OF PDE F. CAGNETTI, D. GOMES, AND H.V. TRAN Abstract. The adjoint method, recently introduced by Evans, is used to study obstacle problems,

More information

S chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.

S chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1. Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable

More information

Measure and Integration: Solutions of CW2

Measure and Integration: Solutions of CW2 Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost

More information

Mean Field Games on networks

Mean Field Games on networks Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017

More information

VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS

VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS LUIS SILVESTRE These are the notes from the summer course given in the Second Chicago Summer School In Analysis, in June 2015. We introduce the notion of viscosity

More information

Sobolev spaces. May 18

Sobolev spaces. May 18 Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references

More information

Calculus of Variations. Final Examination

Calculus of Variations. Final Examination Université Paris-Saclay M AMS and Optimization January 18th, 018 Calculus of Variations Final Examination Duration : 3h ; all kind of paper documents (notes, books...) are authorized. The total score of

More information

WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS

WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS ROLAND HERZOG AND FRANK SCHMIDT Abstract. Sufficient conditions ensuring weak lower

More information

g 2 (x) (1/3)M 1 = (1/3)(2/3)M.

g 2 (x) (1/3)M 1 = (1/3)(2/3)M. COMPACTNESS If C R n is closed and bounded, then by B-W it is sequentially compact: any sequence of points in C has a subsequence converging to a point in C Conversely, any sequentially compact C R n is

More information

VISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.

VISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5. VISCOSITY SOLUTIONS PETER HINTZ We follow Han and Lin, Elliptic Partial Differential Equations, 5. 1. Motivation Throughout, we will assume that Ω R n is a bounded and connected domain and that a ij C(Ω)

More information

OPTIMAL POTENTIALS FOR SCHRÖDINGER OPERATORS. 1. Introduction In this paper we consider optimization problems of the form. min F (V ) : V V, (1.

OPTIMAL POTENTIALS FOR SCHRÖDINGER OPERATORS. 1. Introduction In this paper we consider optimization problems of the form. min F (V ) : V V, (1. OPTIMAL POTENTIALS FOR SCHRÖDINGER OPERATORS G. BUTTAZZO, A. GEROLIN, B. RUFFINI, AND B. VELICHKOV Abstract. We consider the Schrödinger operator + V (x) on H 0 (), where is a given domain of R d. Our

More information

3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?

3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure? MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable

More information

JUHA KINNUNEN. Harmonic Analysis

JUHA KINNUNEN. Harmonic Analysis JUHA KINNUNEN Harmonic Analysis Department of Mathematics and Systems Analysis, Aalto University 27 Contents Calderón-Zygmund decomposition. Dyadic subcubes of a cube.........................2 Dyadic cubes

More information

Problem List MATH 5143 Fall, 2013

Problem List MATH 5143 Fall, 2013 Problem List MATH 5143 Fall, 2013 On any problem you may use the result of any previous problem (even if you were not able to do it) and any information given in class up to the moment the problem was

More information

Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations

Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations Journal of Statistical Physics, Vol. 122, No. 2, January 2006 ( C 2006 ) DOI: 10.1007/s10955-005-8006-x Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations Jan

More information

Numerical methods for Mean Field Games: additional material

Numerical methods for Mean Field Games: additional material Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical

More information

1 Independent increments

1 Independent increments Tel Aviv University, 2008 Brownian motion 1 1 Independent increments 1a Three convolution semigroups........... 1 1b Independent increments.............. 2 1c Continuous time................... 3 1d Bad

More information

Mathematical Methods for Physics and Engineering

Mathematical Methods for Physics and Engineering Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory

More information

A metric space X is a non-empty set endowed with a metric ρ (x, y):

A metric space X is a non-empty set endowed with a metric ρ (x, y): Chapter 1 Preliminaries References: Troianiello, G.M., 1987, Elliptic differential equations and obstacle problems, Plenum Press, New York. Friedman, A., 1982, Variational principles and free-boundary

More information

Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation

Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Albert Fathi IAS Princeton March 15, 2016 In this lecture, a singularity for a locally Lipschitz real valued function

More information

Robustness for a Liouville type theorem in exterior domains

Robustness for a Liouville type theorem in exterior domains Robustness for a Liouville type theorem in exterior domains Juliette Bouhours 1 arxiv:1207.0329v3 [math.ap] 24 Oct 2014 1 UPMC Univ Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, Paris,

More information

Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation.

Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation. Lecture 7 Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation div + u = ϕ on ) = 0 in The solution is a critical point or the minimizer

More information

Real Analysis Problems

Real Analysis Problems Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.

More information

Lecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University

Lecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University Lecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University February 7, 2007 2 Contents 1 Metric Spaces 1 1.1 Basic definitions...........................

More information

The Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:

The Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply

More information

REGULARITY RESULTS FOR THE EQUATION u 11 u 22 = Introduction

REGULARITY RESULTS FOR THE EQUATION u 11 u 22 = Introduction REGULARITY RESULTS FOR THE EQUATION u 11 u 22 = 1 CONNOR MOONEY AND OVIDIU SAVIN Abstract. We study the equation u 11 u 22 = 1 in R 2. Our results include an interior C 2 estimate, classical solvability

More information

Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks

Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks C. Imbert and R. Monneau June 24, 2014 Abstract We study Hamilton-Jacobi equations on networks in the case where Hamiltonians

More information

2 Lebesgue integration

2 Lebesgue integration 2 Lebesgue integration 1. Let (, A, µ) be a measure space. We will always assume that µ is complete, otherwise we first take its completion. The example to have in mind is the Lebesgue measure on R n,

More information

Part III. 10 Topological Space Basics. Topological Spaces

Part III. 10 Topological Space Basics. Topological Spaces Part III 10 Topological Space Basics Topological Spaces Using the metric space results above as motivation we will axiomatize the notion of being an open set to more general settings. Definition 10.1.

More information

2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?

2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure? MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is

More information

Non-radial solutions to a bi-harmonic equation with negative exponent

Non-radial solutions to a bi-harmonic equation with negative exponent Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei

More information

Probability and Measure

Probability and Measure Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability

More information

PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION

PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION ALESSIO FIGALLI AND YOUNG-HEON KIM Abstract. Given Ω, Λ R n two bounded open sets, and f and g two probability densities concentrated

More information

L p Spaces and Convexity

L p Spaces and Convexity L p Spaces and Convexity These notes largely follow the treatments in Royden, Real Analysis, and Rudin, Real & Complex Analysis. 1. Convex functions Let I R be an interval. For I open, we say a function

More information

Some Properties of the Augmented Lagrangian in Cone Constrained Optimization

Some Properties of the Augmented Lagrangian in Cone Constrained Optimization MATHEMATICS OF OPERATIONS RESEARCH Vol. 29, No. 3, August 2004, pp. 479 491 issn 0364-765X eissn 1526-5471 04 2903 0479 informs doi 10.1287/moor.1040.0103 2004 INFORMS Some Properties of the Augmented

More information

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively); STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend

More information

Convergence of generalized entropy minimizers in sequences of convex problems

Convergence of generalized entropy minimizers in sequences of convex problems Proceedings IEEE ISIT 206, Barcelona, Spain, 2609 263 Convergence of generalized entropy minimizers in sequences of convex problems Imre Csiszár A Rényi Institute of Mathematics Hungarian Academy of Sciences

More information

VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS, AND ASYMPTOTICS FOR HAMILTONIAN SYSTEMS

VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS, AND ASYMPTOTICS FOR HAMILTONIAN SYSTEMS VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS, AND ASYMPTOTICS FOR HAMILTONIAN SYSTEMS DIOGO AGUIAR GOMES U.C. Berkeley - CA, US and I.S.T. - Lisbon, Portugal email:dgomes@math.ist.utl.pt Abstract.

More information

The effects of a discontinues weight for a problem with a critical nonlinearity

The effects of a discontinues weight for a problem with a critical nonlinearity arxiv:1405.7734v1 [math.ap] 9 May 014 The effects of a discontinues weight for a problem with a critical nonlinearity Rejeb Hadiji and Habib Yazidi Abstract { We study the minimizing problem px) u dx,

More information

Partial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations

Partial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 8 The Calculus of Variations Yung-Hsiang Huang 2018.03.25 Notation: denotes a bounded smooth, open subset of R n. All given functions are

More information

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3) M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation

More information

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov, LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical

More information

On the distributional divergence of vector fields vanishing at infinity

On the distributional divergence of vector fields vanishing at infinity Proceedings of the Royal Society of Edinburgh, 141A, 65 76, 2011 On the distributional divergence of vector fields vanishing at infinity Thierry De Pauw Institut de Recherches en Mathématiques et Physique,

More information

EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS

EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS Electronic Journal of Differential Equations, Vol. 017 (017), No. 15, pp. 1 7. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC

More information

Heat kernels of some Schrödinger operators

Heat kernels of some Schrödinger operators Heat kernels of some Schrödinger operators Alexander Grigor yan Tsinghua University 28 September 2016 Consider an elliptic Schrödinger operator H = Δ + Φ, where Δ = n 2 i=1 is the Laplace operator in R

More information

CHAPTER 6. Differentiation

CHAPTER 6. Differentiation CHPTER 6 Differentiation The generalization from elementary calculus of differentiation in measure theory is less obvious than that of integration, and the methods of treating it are somewhat involved.

More information

The optimal partial transport problem

The optimal partial transport problem The optimal partial transport problem Alessio Figalli Abstract Given two densities f and g, we consider the problem of transporting a fraction m [0, min{ f L 1, g L 1}] of the mass of f onto g minimizing

More information

Overview of normed linear spaces

Overview of normed linear spaces 20 Chapter 2 Overview of normed linear spaces Starting from this chapter, we begin examining linear spaces with at least one extra structure (topology or geometry). We assume linearity; this is a natural

More information

l(y j ) = 0 for all y j (1)

l(y j ) = 0 for all y j (1) Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that

More information

SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction

SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy

More information

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC

More information

On a weighted total variation minimization problem

On a weighted total variation minimization problem On a weighted total variation minimization problem Guillaume Carlier CEREMADE Université Paris Dauphine carlier@ceremade.dauphine.fr Myriam Comte Laboratoire Jacques-Louis Lions, Université Pierre et Marie

More information

NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES

NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES JUHA LEHRBÄCK Abstract. We establish necessary conditions for domains Ω R n which admit the pointwise (p, β)-hardy inequality u(x) Cd Ω(x)

More information

Convex Optimization Notes

Convex Optimization Notes Convex Optimization Notes Jonathan Siegel January 2017 1 Convex Analysis This section is devoted to the study of convex functions f : B R {+ } and convex sets U B, for B a Banach space. The case of B =

More information

u xx + u yy = 0. (5.1)

u xx + u yy = 0. (5.1) Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function

More information

Analysis Finite and Infinite Sets The Real Numbers The Cantor Set

Analysis Finite and Infinite Sets The Real Numbers The Cantor Set Analysis Finite and Infinite Sets Definition. An initial segment is {n N n n 0 }. Definition. A finite set can be put into one-to-one correspondence with an initial segment. The empty set is also considered

More information

MATHS 730 FC Lecture Notes March 5, Introduction

MATHS 730 FC Lecture Notes March 5, Introduction 1 INTRODUCTION MATHS 730 FC Lecture Notes March 5, 2014 1 Introduction Definition. If A, B are sets and there exists a bijection A B, they have the same cardinality, which we write as A, #A. If there exists

More information

A Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth

A Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth A Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth E. DiBenedetto 1 U. Gianazza 2 C. Klaus 1 1 Vanderbilt University, USA 2 Università

More information

BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED

BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis

More information

Maths 212: Homework Solutions

Maths 212: Homework Solutions Maths 212: Homework Solutions 1. The definition of A ensures that x π for all x A, so π is an upper bound of A. To show it is the least upper bound, suppose x < π and consider two cases. If x < 1, then

More information

SOME VIEWS ON GLOBAL REGULARITY OF THE THIN FILM EQUATION

SOME VIEWS ON GLOBAL REGULARITY OF THE THIN FILM EQUATION SOME VIEWS ON GLOBAL REGULARITY OF THE THIN FILM EQUATION STAN PALASEK Abstract. We introduce the thin film equation and the problem of proving positivity and global regularity on a periodic domain with

More information

RIEMANN SURFACES. ω = ( f i (γ(t))γ i (t))dt.

RIEMANN SURFACES. ω = ( f i (γ(t))γ i (t))dt. RIEMANN SURFACES 6. Week 7: Differential forms. De Rham complex 6.1. Introduction. The notion of differential form is important for us for various reasons. First of all, one can integrate a k-form along

More information

Riemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,

Riemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E, Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of

More information

P-adic Functions - Part 1

P-adic Functions - Part 1 P-adic Functions - Part 1 Nicolae Ciocan 22.11.2011 1 Locally constant functions Motivation: Another big difference between p-adic analysis and real analysis is the existence of nontrivial locally constant

More information

i=1 α i. Given an m-times continuously

i=1 α i. Given an m-times continuously 1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable

More information

Stochastic Homogenization for Reaction-Diffusion Equations

Stochastic Homogenization for Reaction-Diffusion Equations Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires

More information

The Relativistic Heat Equation

The Relativistic Heat Equation Maximum Principles and Behavior near Absolute Zero Washington University in St. Louis ARTU meeting March 28, 2014 The Heat Equation The heat equation is the standard model for diffusion and heat flow,

More information

The infinity-laplacian and its properties

The infinity-laplacian and its properties U.U.D.M. Project Report 2017:40 Julia Landström Examensarbete i matematik, 15 hp Handledare: Kaj Nyström Examinator: Martin Herschend December 2017 Department of Mathematics Uppsala University Department

More information

LECTURE 15: COMPLETENESS AND CONVEXITY

LECTURE 15: COMPLETENESS AND CONVEXITY LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other

More information

Minimal time mean field games

Minimal time mean field games based on joint works with Samer Dweik and Filippo Santambrogio PGMO Days 2017 Session Mean Field Games and applications EDF Lab Paris-Saclay November 14th, 2017 LMO, Université Paris-Sud Université Paris-Saclay

More information