Book 6C

© 2026 NeckPinch · www.neckpinch.com · All rights reserved.

Course 6Book 6C: Nonlinear PDE: Variational and Geometric MethodsChapter 5

Fixed Points, Monotonicity and Sub- and Supersolutions

Schauder and Schaefer, Leray–Schauder a priori estimates, Minty–Browder and monotone iteration.

46 min read · Updated Oct 23, 2026

Read with Evans, Partial Differential Equations, sections 9.1 (monotonicity methods), 9.2 (fixed point methods: Banach, Schauder, Schaefer) and 9.3 (the method of subsolutions and supersolutions); Evans section 8.1.4 contains the proof of Brouwer's theorem by null Lagrangians used here. Taylor III, chapter 14, uses the same tools for quasilinear equations. Zeidler, Nonlinear Functional Analysis and its Applications (volumes I and II/B), is the encyclopaedic reference. This chapter proves Brouwer's theorem (by the determinant as a null Lagrangian), Schauder's and Schaefer's theorems, the existence theorem for monotone operators with Minty's trick, and the monotone iteration of sub- and supersolutions.

In this chapter · 8 sections
  1. 5.1When can a population survive?
  2. 5.2Brouwer's fixed point theorem
  3. 5.3Schauder's fixed point theorem
  4. 5.4A priori estimates give existence
  5. 5.5Monotone operators
  6. 5.6Sub- and supersolutions
  7. 5.7History
  8. 5.8Exercises
  9. ↗Free resources

The variational methods of 6C.1 Regularity of Minimisers and Constraints, 6C.2 Critical Points and the Mountain Pass and 6C.4 The Yamabe Problem and Critical Exponents need an energy whose critical points are the solutions. Many equations have none: an equation with a first-order term b(∇u)b(\nabla u), a system that is not a gradient, a nonlinearity depending on xx awkwardly. This chapter gives the three other great sources of existence for nonlinear equations.

Fixed points. Write the equation as u=A(u)u = A(u) for a map AA built by solving a linear problem, and show AA has a fixed point. Banach's contraction principle needs AA to be a contraction, which is how 6A.8 Nonlinear Parabolic Equations solved semilinear heat equations for short times. When AA is merely continuous and compact (it maps bounded sets to precompact ones, which is what elliptic regularity plus Rellich gives), Schauder's theorem still produces a fixed point, and Schaefer's version reduces existence to an a priori estimate: if all solutions of a family of related problems are bounded, a solution exists. This "a priori estimates imply existence" principle, due to Jean Leray and Juliusz Schauder, is the backbone of nonlinear elliptic theory, and it is why most of the work in a nonlinear PDE paper is proving estimates for solutions that are not yet known to exist.

Monotonicity. For divergence-form operators −div⁡a(∇u)-\operatorname{div}a(\nabla u) with aa monotone, as for the pp-Laplacian, a Galerkin scheme converges weakly, and the trick of George Minty identifies the weak limit of the nonlinear term, which weak convergence alone cannot do.

Order. If the equation has a maximum principle, a subsolution below and a supersolution above trap a solution between them, found by a monotone iteration. This is the most practical existence method of all, and it gives the critical patch size of mathematical ecology and the uniformization of negatively curved surfaces.

By the end of this chapter you will be able to:

  • prove Brouwer's fixed point theorem from the fact that ∫det⁡Du\int\det Du depends only on the boundary values of uu;
  • prove Schauder's fixed point theorem by finite-dimensional approximation, and Schaefer's theorem from it;
  • turn an a priori estimate into an existence theorem, for an equation with a gradient nonlinearity;
  • solve −div⁡a(∇u)=f-\operatorname{div}a(\nabla u) = f for monotone aa by Galerkin's method and Minty's trick;
  • prove the method of sub- and supersolutions by monotone iteration, and apply it to the logistic equation and to prescribing negative curvature on a surface.

When can a population survive?

In the world Model Critical patch size

A population of plankton grows logistically in a patch of water, but individuals diffuse out of the patch and are lost (the surrounding water is hostile, or carries them away). With density u(x)u(x) scaled by the carrying capacity, the steady states solve

−DΔu=ru(1−u)  in U,u=0  on ∂U,-D\Delta u = ru(1 - u) \ \text{ in } U, \qquad u = 0 \ \text{ on } \partial U,

with DD the diffusion coefficient and rr the growth rate. Does a positive steady state exist? John Skellam (1951), and independently Henry Kierstead and Lawrence Slobodkin (1953), who were studying why plankton blooms need water masses of a certain size, found the answer: only if the patch is larger than a critical size. In one dimension the patch (0,L)(0, L) supports the population exactly when L>πD/rL > \pi\sqrt{D/r}; in general the condition is rD>λ1(U)\frac rD > \lambda_1(U), the first Dirichlet eigenvalue. Growth must beat the loss rate through the boundary, and that loss rate is the eigenvalue.

The "KISS" model (for Kierstead, Slobodkin and Skellam) is the starting point of the theory of nature reserves and habitat fragmentation: a reserve smaller than the critical size cannot sustain a species however good the habitat inside. The method of sub- and supersolutions proves the existence half in three lines (Example 5.12), and a test against the first eigenfunction proves the converse (Exercise 5.16).

Brouwer's fixed point theorem

Theorem 5.1 Brouwer's fixed point theorem

Every continuous map f:B‾→B‾f : \overline{B} \to \overline{B} of the closed unit ball of Rn\mathbb{R}^n into itself has a fixed point. The same holds for any compact convex set K⊂RnK \subset \mathbb{R}^n in place of B‾\overline B.

The theorem is topological, but there is a short analytic proof, which Evans gives (section 8.1.4) and which uses an idea worth having for its own sake: the determinant of the Jacobian is a null Lagrangian.

Lemma 5.2 The Jacobian integral depends only on boundary values

Let U⊂RnU \subset \mathbb{R}^n be bounded with smooth boundary, and u,v∈C2(U‾;Rn)u, v \in C^2(\overline U; \mathbb{R}^n) with u=vu = v on ∂U\partial U. Then ∫Udet⁡Du dx=∫Udet⁡Dv dx\int_U\det Du\,dx = \int_U\det Dv\,dx.

Proof. In the language of differential forms (8A.8 Differential Forms and Stokes’ Theorem), det⁡Du dx1∧⋯∧dxn=du1∧du2∧⋯∧dun\det Du\,dx^1\wedge\cdots\wedge dx^n = du^1\wedge du^2\wedge\cdots\wedge du^n, the pullback of the volume form by uu. Since d(dui)=0d(du^i) = 0,

du1∧du2∧⋯∧dun=d(u1 du2∧⋯∧dun).du^1\wedge du^2\wedge\cdots\wedge du^n = d\big(u^1\,du^2\wedge\cdots\wedge du^n\big).

By Stokes' theorem, ∫Udet⁡Du dx=∫∂Uu1 du2∧⋯∧dun\int_U\det Du\,dx = \int_{\partial U}u^1\,du^2\wedge\cdots\wedge du^n. The form on the right is restricted to ∂U\partial U, and the restriction of duidu^i to the boundary is d(ui∣∂U)d(u^i|_{\partial U}): it depends only on the values of uu on ∂U\partial U. So the integral is the same for vv.

A retraction of B‾\overline B onto its boundary sphere SS is a continuous map r:B‾→Sr : \overline B \to S with r(x)=xr(x) = x for x∈Sx \in S.

Lemma 5.3 No smooth retraction

There is no C2C^2 map w:B‾→Sw : \overline B \to S with w(x)=xw(x) = x on SS.

Proof. Suppose there were. Apply Lemma 5.2 to ww and the identity map, which agree on ∂B\partial B: ∫Bdet⁡Dw dx=∫Bdet⁡I dx=∣B∣>0\int_B\det Dw\,dx = \int_B\det I\,dx = |B| > 0. But ∣w∣2=1|w|^2 = 1 everywhere, and differentiating, (Dw)Tw=0(Dw)^Tw = 0: the nonzero vector w(x)w(x) is in the kernel of Dw(x)TDw(x)^T, so det⁡Dw(x)=0\det Dw(x) = 0 at every point. Contradiction.

Proof (Brouwer's theorem). Suppose f:B‾→B‾f : \overline B \to \overline B is continuous with no fixed point. Then ∣f(x)−x∣≥3ε>0|f(x) - x| \geq 3\varepsilon > 0 on the compact ball, for some ε\varepsilon. Extend ff continuously to Rn\mathbb{R}^n by f(x)=f(x/∣x∣)f(x) = f(x/|x|) for ∣x∣>1|x| > 1 and mollify it: since ff is uniformly continuous near B‾\overline B, convolution with a narrow mollifier gives a smooth map g0g_0 with ∣g0−f∣<ε|g_0 - f| < \varepsilon on B‾\overline B; then g=g0/(1+ε)g = g_0/(1 + \varepsilon) maps B‾\overline B into B‾\overline B, is smooth, and ∣g−f∣≤∣g0−f∣+ε∣g0∣/(1+ε)<2ε|g - f| \leq |g_0 - f| + \varepsilon|g_0|/(1 + \varepsilon) < 2\varepsilon, so ∣g(x)−x∣>ε|g(x) - x| > \varepsilon: gg has no fixed point either.

Now define w(x)w(x) as the point where the ray from g(x)g(x) through xx leaves the ball: w(x)=x+t(x)(x−g(x))w(x) = x + t(x)\big(x - g(x)\big) with t(x)≥0t(x) \geq 0 the larger root of ∣x+t(x−g(x))∣2=1|x + t(x - g(x))|^2 = 1. Writing e=x−g(x)≠0e = x - g(x) \neq 0, this is the quadratic ∣e∣2t2+2(x⋅e)t+∣x∣2−1=0|e|^2t^2 + 2(x\cdot e)t + |x|^2 - 1 = 0, whose discriminant 4(x⋅e)2+4∣e∣2(1−∣x∣2)4(x\cdot e)^2 + 4|e|^2(1 - |x|^2) is positive except where ∣x∣=1|x| = 1 and x⋅e=0x\cdot e = 0. At such a point x⋅g(x)=∣x∣2=1x\cdot g(x) = |x|^2 = 1, which for ∣g(x)∣≤1=∣x∣|g(x)| \leq 1 = |x| forces g(x)=xg(x) = x, excluded. So the discriminant is positive on B‾\overline B, t(x)t(x) is a smooth function of xx, and ww is smooth with values in SS. For ∣x∣=1|x| = 1, x⋅e=1−x⋅g(x)≥0x\cdot e = 1 - x\cdot g(x) \geq 0, so the larger root is t=0t = 0 and w(x)=xw(x) = x. This is a smooth retraction, contradicting Lemma 5.3.

For a compact convex KK: KK is a retract of a large closed ball BR‾⊇K\overline{B_R} \supseteq K by the nearest-point projection PKP_K, which is continuous (4A.5 Hilbert Spaces and Lax–Milgram; it is 11-Lipschitz). If f:K→Kf : K \to K, then f∘PK:BR‾→K⊆BR‾f\circ P_K : \overline{B_R} \to K \subseteq \overline{B_R} has a fixed point xx, which lies in KK, so PKx=xP_Kx = x and f(x)=xf(x) = x.

Example 5.4 Nash equilibria

In a finite game, each player ii chooses a mixed strategy, a probability vector σi\sigma_i on their finitely many pure strategies, and receives expected payoff πi(σ1,…,σN)\pi_i(\sigma_1, \ldots, \sigma_N), linear in each σj\sigma_j. A Nash equilibrium is a profile from which no player can gain by changing strategy alone. John Nash's 1951 proof of existence (published in the Annals of Mathematics) applied Brouwer's theorem to the map on the product of simplices (compact and convex) that shifts each player's weight towards pure strategies doing better than their current mix: σi(s)↦σi(s)+max⁡(0,πi(s,σ−i)−πi(σ))1+∑s′max⁡(0,πi(s′,σ−i)−πi(σ))\sigma_i(s) \mapsto \frac{\sigma_i(s) + \max(0, \pi_i(s, \sigma_{-i}) - \pi_i(\sigma))}{1 + \sum_{s'}\max(0, \pi_i(s', \sigma_{-i}) - \pi_i(\sigma))}. Its fixed points are exactly the equilibria. The existence of competitive equilibria in economics (Arrow and Debreu, 1954) rests on fixed point theorems in the same way. In every case the theorem proves existence without saying how to find the fixed point, and computing one can be genuinely hard.

Schauder's fixed point theorem

In infinite dimensions Brouwer fails: the closed unit ball of ℓ2\ell^2 has a continuous self-map without fixed points (Exercise 5.14). What restores it is compactness.

Theorem 5.5 Schauder's fixed point theorem

Let KK be a compact convex subset of a Banach space XX and A:K→KA : K \to K continuous. Then AA has a fixed point.

Proof. Fix ε>0\varepsilon > 0. Since KK is compact, it has a finite ε\varepsilon-net u1,…,uN∈Ku_1, \ldots, u_N \in K: every u∈Ku \in K is within ε\varepsilon of some uiu_i. The Schauder projection is

Pε(u)=∑iλi(u) ui∑iλi(u),λi(u)=max⁡(ε−∥u−ui∥,0).P_\varepsilon(u) = \frac{\sum_i\lambda_i(u)\,u_i}{\sum_i\lambda_i(u)}, \qquad \lambda_i(u) = \max\big(\varepsilon - \|u - u_i\|, 0\big).

The denominator is positive (some λi(u)>0\lambda_i(u) > 0), so PεP_\varepsilon is continuous, with values in KεK_\varepsilon, the convex hull of u1,…,uNu_1, \ldots, u_N, which lies in KK (convexity) and is a compact convex subset of a finite-dimensional space. Since Pε(u)P_\varepsilon(u) is an average of points uiu_i with ∥u−ui∥<ε\|u - u_i\| < \varepsilon, ∥Pε(u)−u∥<ε\|P_\varepsilon(u) - u\| < \varepsilon.

The map Pε∘AP_\varepsilon\circ A sends KεK_\varepsilon into KεK_\varepsilon continuously, so by Brouwer's theorem for compact convex sets (in the finite-dimensional span of the uiu_i) it has a fixed point vε=Pε(Avε)v_\varepsilon = P_\varepsilon(Av_\varepsilon). Then ∥Avε−vε∥=∥Avε−Pε(Avε)∥<ε\|Av_\varepsilon - v_\varepsilon\| = \|Av_\varepsilon - P_\varepsilon(Av_\varepsilon)\| < \varepsilon. Take ε=1k\varepsilon = \frac1k; by compactness of KK a subsequence v1/kj→v∈Kv_{1/k_j} \to v \in K, and by continuity ∥Av−v∥=lim⁡∥Av1/kj−v1/kj∥=0\|Av - v\| = \lim\|Av_{1/k_j} - v_{1/k_j}\| = 0.

In PDE, the set KK is rarely compact itself; what is compact is the map. A map A:X→XA : X \to X is compact if it is continuous and maps bounded sets to precompact sets. Elliptic regularity makes solution operators compact: solving −Δw=g-\Delta w = g with g∈L2g \in L^2 gives w∈H2w \in H^2, and H2H^2 embeds compactly in H1H^1 (Rellich, 4A.11 Sobolev Embeddings and Critical Exponents), so g↦wg \mapsto w is compact from L2L^2 to H1H^1.

Theorem 5.6 Schaefer's fixed point theorem

Let A:X→XA : X \to X be compact. Suppose the set

{u∈X:u=λA(u) for some 0≤λ≤1}\{u \in X : u = \lambda A(u) \ \text{for some } 0 \leq \lambda \leq 1\}

is bounded. Then AA has a fixed point.

Proof. Choose MM with ∥u∥<M\|u\| < M for every uu in the set. Define

A~(u)={A(u)if ∥A(u)∥≤M,M A(u)∥A(u)∥if ∥A(u)∥>M,\tilde A(u) = \begin{cases}A(u) & \text{if } \|A(u)\| \leq M,\\ \dfrac{M\,A(u)}{\|A(u)\|} & \text{if } \|A(u)\| > M,\end{cases}

continuous, mapping the closed ball BM‾\overline{B_M} into itself, with A~(BM‾)\tilde A(\overline{B_M}) precompact. Let KK be the closed convex hull of A~(BM‾)\tilde A(\overline{B_M}): it is compact (the closed convex hull of a compact set in a Banach space is compact, Mazur's theorem) and contained in BM‾\overline{B_M}, so A~\tilde A maps KK into KK. By Theorem 5.5 there is u∈Ku \in K with A~(u)=u\tilde A(u) = u. If ∥A(u)∥>M\|A(u)\| > M, then u=λA(u)u = \lambda A(u) with λ=M/∥A(u)∥<1\lambda = M/\|A(u)\| < 1 and ∥u∥=M\|u\| = M, contradicting the choice of MM. So ∥A(u)∥≤M\|A(u)\| \leq M and u=A(u)u = A(u).

The theorem is a version of the Leray–Schauder principle: deform the problem u=A(u)u = A(u) to the trivial problem u=0u = 0 through u=λA(u)u = \lambda A(u); if no solution escapes to infinity during the deformation, a solution survives to λ=1\lambda = 1. Leray and Schauder (1934) made this rigorous with a degree theory for compact perturbations of the identity, which also counts solutions with signs.

A priori estimates give existence

Theorem 5.7 A semilinear equation with a gradient term

Let UU be bounded with smooth boundary, f∈L2(U)f \in L^2(U), and b:Rn→Rb : \mathbb{R}^n \to \mathbb{R} Lipschitz, so that ∣b(p)∣≤C(∣p∣+1)|b(p)| \leq C(|p| + 1). If μ>C2+12\mu > \frac{C^2 + 1}{2}, the problem

−Δu+b(∇u)+μu=f  in U,u=0  on ∂U-\Delta u + b(\nabla u) + \mu u = f \ \text{ in } U, \qquad u = 0 \ \text{ on } \partial U

has a solution u∈H2(U)∩H01(U)u \in H^2(U)\cap H^1_0(U).

Proof. The map. For u∈H01(U)u \in H^1_0(U), let A(u)=wA(u) = w be the solution of the linear problem −Δw+μw=f−b(∇u)-\Delta w + \mu w = f - b(\nabla u), w∈H01w \in H^1_0. The right side is in L2L^2, with ∥f−b(∇u)∥L2≤∥f∥L2+C(∥∇u∥L2+∣U∣1/2)\|f - b(\nabla u)\|_{L^2} \leq \|f\|_{L^2} + C(\|\nabla u\|_{L^2} + |U|^{1/2}), so by Lax–Milgram and boundary regularity (6A.6 Weak Solutions and Elliptic Regularity) w∈H2(U)w \in H^2(U) with ∥w∥H2≤C′(∥f∥L2+∥u∥H1+1)\|w\|_{H^2} \leq C'(\|f\|_{L^2} + \|u\|_{H^1} + 1).

AA is compact. Bounded sets of H01H^1_0 go to bounded sets of H2H^2, which are precompact in H01H^1_0 by Rellich. And AA is continuous: ∥A(u)−A(v)∥H2≤C′∥b(∇u)−b(∇v)∥L2≤C′Lip⁡(b)∥u−v∥H1\|A(u) - A(v)\|_{H^2} \leq C'\|b(\nabla u) - b(\nabla v)\|_{L^2} \leq C'\operatorname{Lip}(b)\|u - v\|_{H^1}.

The a priori estimate. Suppose u=λA(u)u = \lambda A(u), 0≤λ≤10 \leq \lambda \leq 1: then −Δu+μu=λ(f−b(∇u))-\Delta u + \mu u = \lambda\big(f - b(\nabla u)\big). Multiply by uu and integrate:

∫∣∇u∣2+μ∫u2=λ∫(f−b(∇u))u≤∫(∣f∣+C)∣u∣+C∫∣∇u∣∣u∣.\int|\nabla u|^2 + \mu\int u^2 = \lambda\int\big(f - b(\nabla u)\big)u \leq \int\big(|f| + C\big)|u| + C\int|\nabla u||u|.

By Young's inequality, C∫∣∇u∣∣u∣≤12∫∣∇u∣2+C22∫u2C\int|\nabla u||u| \leq \frac12\int|\nabla u|^2 + \frac{C^2}{2}\int u^2 and ∫(∣f∣+C)∣u∣≤12∫u2+12∫(∣f∣+C)2\int(|f| + C)|u| \leq \frac12\int u^2 + \frac12\int(|f| + C)^2. So

12∫∣∇u∣2+(μ−C2+12)∫u2≤12∫(∣f∣+C)2,\frac12\int|\nabla u|^2 + \Big(\mu - \frac{C^2 + 1}{2}\Big)\int u^2 \leq \frac12\int(|f| + C)^2,

a bound on ∥u∥H1\|u\|_{H^1} independent of λ\lambda. By Theorem 5.6, AA has a fixed point, which is the solution.

Notice where the work went: the existence proof is three lines once the estimate for hypothetical solutions is known. For a quasilinear equation the same scheme needs a bound in C1,αC^{1,\alpha} rather than H1H^1, because the coefficients depend on ∇u\nabla u; 6C.6 Quasilinear Elliptic Equations and Minimal Surfaces proves such bounds for the minimal surface equation, and Yau's solution of the Calabi conjecture (6C.10 Fully Nonlinear Equations and Monge–Ampère) is, at its core, a priori estimates fed into the continuity method of 6B.7 Schauder Theory and the Continuity Method, which is the same principle.

Monotone operators

Now consider divergence-form equations

−div⁡a(∇u)=f  in U,u=0  on ∂U,-\operatorname{div}a(\nabla u) = f \ \text{ in } U, \qquad u = 0 \ \text{ on } \partial U,

with a:Rn→Rna : \mathbb{R}^n \to \mathbb{R}^n continuous and

  1. monotone: (a(p)−a(q))⋅(p−q)≥0\big(a(p) - a(q)\big)\cdot(p - q) \geq 0 for all p,qp, q;
  2. linear growth: ∣a(p)∣≤C(1+∣p∣)|a(p)| \leq C(1 + |p|);
  3. coercive: a(p)⋅p≥α∣p∣2−βa(p)\cdot p \geq \alpha|p|^2 - \beta, with α>0\alpha > 0, β≥0\beta \geq 0.

If a=∇Fa = \nabla F for a convex FF, this is the Euler–Lagrange equation of ∫F(∇u)−fu\int F(\nabla u) - fu, and 6A.11 Calculus of Variations and Gradient Flows applies; monotonicity is what convexity of FF becomes when aa is not a gradient. A weak solution is u∈H01u \in H^1_0 with ∫a(∇u)⋅∇v=∫fv\int a(\nabla u)\cdot\nabla v = \int fv for all v∈H01v \in H^1_0.

The difficulty is the limit. A Galerkin scheme gives um⇀uu_m \rightharpoonup u weakly in H01H^1_0, but a(∇um)a(\nabla u_m) need not converge to a(∇u)a(\nabla u): weak convergence does not commute with nonlinear functions (sin⁡(mx)⇀0\sin(mx) \rightharpoonup 0 but sin⁡2(mx)⇀12\sin^2(mx) \rightharpoonup \frac12). Minty's trick uses monotonicity to identify the limit anyway.

Lemma 5.8 Zeros of vector fields pointing outward

Let v:Rm→Rmv : \mathbb{R}^m \to \mathbb{R}^m be continuous with v(x)⋅x≥0v(x)\cdot x \geq 0 whenever ∣x∣=r|x| = r. Then v(x)=0v(x) = 0 for some ∣x∣≤r|x| \leq r.

Proof. If not, w(x)=−rv(x)∣v(x)∣w(x) = -r\frac{v(x)}{|v(x)|} is a continuous map of the closed ball Br‾\overline{B_r} into itself, with a fixed point zz by Brouwer. Then ∣z∣=r|z| = r, and r2=z⋅z=−rv(z)⋅z∣v(z)∣≤0r^2 = z\cdot z = -r\frac{v(z)\cdot z}{|v(z)|} \leq 0, a contradiction.

Theorem 5.9 Existence for monotone operators

Under 1–3 and f∈L2(U)f \in L^2(U), there is a weak solution u∈H01(U)u \in H^1_0(U). If aa is strictly monotone, (a(p)−a(q))⋅(p−q)≥θ∣p−q∣2\big(a(p) - a(q)\big)\cdot(p - q) \geq \theta|p - q|^2 with θ>0\theta > 0, it is unique.

Proof. Galerkin approximations. Let w1,w2,…w_1, w_2, \ldots be the Dirichlet eigenfunctions of −Δ-\Delta (6A.6 Weak Solutions and Elliptic Regularity), an orthogonal basis of H01H^1_0. We look for um=∑k=1mdmkwku_m = \sum_{k=1}^md_m^kw_k with

∫a(∇um)⋅∇wk=∫fwk,k=1,…,m.(G)\int a(\nabla u_m)\cdot\nabla w_k = \int fw_k, \qquad k = 1, \ldots, m. \tag{G}

Define v:Rm→Rmv : \mathbb{R}^m \to \mathbb{R}^m by v(d)k=∫a(∇∑jdjwj)⋅∇wk−∫fwkv(d)^k = \int a\big(\nabla\sum_jd^jw_j\big)\cdot\nabla w_k - \int fw_k, continuous by 2 and dominated convergence. With u=∑djwju = \sum d^jw_j, by 3,

v(d)⋅d=∫a(∇u)⋅∇u−∫fu≥α∥∇u∥L22−β∣U∣−∥f∥L2∥u∥L2.v(d)\cdot d = \int a(\nabla u)\cdot\nabla u - \int fu \geq \alpha\|\nabla u\|_{L^2}^2 - \beta|U| - \|f\|_{L^2}\|u\|_{L^2}.

The norms ∥∇u∥L2\|\nabla u\|_{L^2} and ∣d∣|d| are equivalent on the span of w1,…,wmw_1, \ldots, w_m, so this is positive for ∣d∣=r|d| = r large, and Lemma 5.8 gives a solution of (G).

Bounds. Multiplying (G) by dmkd_m^k and summing, ∫a(∇um)⋅∇um=∫fum\int a(\nabla u_m)\cdot\nabla u_m = \int fu_m, so by 3 and Poincaré, α∥∇um∥2≤β∣U∣+CP∥f∥∥∇um∥\alpha\|\nabla u_m\|^2 \leq \beta|U| + C_P\|f\|\|\nabla u_m\|: the umu_m are bounded in H01H^1_0, and by 2 the a(∇um)a(\nabla u_m) are bounded in L2(U;Rn)L^2(U; \mathbb{R}^n). Pass to a subsequence: um⇀uu_m \rightharpoonup u in H01H^1_0, um→uu_m \to u in L2L^2 (Rellich), and a(∇um)⇀ξa(\nabla u_m) \rightharpoonup \xi in L2L^2. Letting m→∞m \to \infty in (G) with kk fixed, ∫ξ⋅∇wk=∫fwk\int\xi\cdot\nabla w_k = \int fw_k for every kk, hence

∫ξ⋅∇v=∫fvfor all v∈H01.(L)\int\xi\cdot\nabla v = \int fv \quad \text{for all } v \in H^1_0. \tag{L}

Minty's trick: ξ=a(∇u)\xi = a(\nabla u). For any w∈H01w \in H^1_0, monotonicity gives

∫(a(∇um)−a(∇w))⋅(∇um−∇w)≥0.\int\big(a(\nabla u_m) - a(\nabla w)\big)\cdot(\nabla u_m - \nabla w) \geq 0.

Expand. The term ∫a(∇um)⋅∇um=∫fum→∫fu=∫ξ⋅∇u\int a(\nabla u_m)\cdot\nabla u_m = \int fu_m \to \int fu = \int\xi\cdot\nabla u, by the energy identity, strong L2L^2 convergence and (L) with v=uv = u. The other terms are pairings of a weakly convergent sequence with a fixed function. So in the limit

∫(ξ−a(∇w))⋅(∇u−∇w)≥0for all w∈H01.\int\big(\xi - a(\nabla w)\big)\cdot(\nabla u - \nabla w) \geq 0 \quad \text{for all } w \in H^1_0.

Take w=u−λvw = u - \lambda v with λ>0\lambda > 0 and v∈H01v \in H^1_0, and divide by λ\lambda: ∫(ξ−a(∇u−λ∇v))⋅∇v≥0\int\big(\xi - a(\nabla u - \lambda\nabla v)\big)\cdot\nabla v \geq 0. As λ→0\lambda \to 0, a(∇u−λ∇v)→a(∇u)a(\nabla u - \lambda\nabla v) \to a(\nabla u) pointwise, dominated by C(1+∣∇u∣+∣∇v∣)C(1 + |\nabla u| + |\nabla v|), so ∫(ξ−a(∇u))⋅∇v≥0\int\big(\xi - a(\nabla u)\big)\cdot\nabla v \geq 0 for every vv; replacing vv by −v-v gives equality. With (L), ∫a(∇u)⋅∇v=∫fv\int a(\nabla u)\cdot\nabla v = \int fv: uu is a weak solution.

Uniqueness. If u1,u2u_1, u_2 are solutions, subtract the equations and test with u1−u2u_1 - u_2: θ∥∇(u1−u2)∥2≤∫(a(∇u1)−a(∇u2))⋅∇(u1−u2)=0\theta\|\nabla(u_1 - u_2)\|^2 \leq \int\big(a(\nabla u_1) - a(\nabla u_2)\big)\cdot\nabla(u_1 - u_2) = 0.

The abstract form is the Minty–Browder theorem (George Minty, 1962–63, and Felix Browder, 1963): a monotone, coercive, hemicontinuous operator from a reflexive Banach space VV to its dual is surjective. The proof is the one above, with VV in place of H01H^1_0.

The pp-Laplacian. The operator −Δpu=−div⁡(∣∇u∣p−2∇u)-\Delta_pu = -\operatorname{div}\big(|\nabla u|^{p-2}\nabla u\big) comes from the energy 1p∫∣∇u∣p\frac1p\int|\nabla u|^p, and its vector field a(ξ)=∣ξ∣p−2ξa(\xi) = |\xi|^{p-2}\xi is monotone (Exercise 5.15). It has growth ∣ξ∣p−1|\xi|^{p-1}, not linear, so the theorem above applies only for p=2p = 2; the same proof in V=W01,p(U)V = W^{1,p}_0(U), with Lp′L^{p'} in place of L2L^2 for a(∇um)a(\nabla u_m), gives a unique weak solution of −Δpu=f-\Delta_pu = f for f∈W−1,p′f \in W^{-1,p'} and every 1<p<∞1 < p < \infty. The pp-Laplacian is degenerate where ∇u=0\nabla u = 0 (for p>2p > 2) or singular there (for p<2p < 2), and its solutions are only C1,αC^{1,\alpha} in general, a theorem of Nina Ural'tseva, Karen Uhlenbeck, Emmanuele DiBenedetto and Peter Tolksdorf (1968–1984). It models non-Newtonian fluids and, as p→∞p \to \infty, optimal Lipschitz extensions.

Sub- and supersolutions

Consider

−Δu=f(x,u)  in U,u=0  on ∂U,(P)-\Delta u = f(x, u) \ \text{ in } U, \qquad u = 0 \ \text{ on } \partial U, \tag{P}

with ff measurable in xx, and Lipschitz in uu on bounded sets uniformly in xx. A subsolution is u‾∈H1(U)\underline u \in H^1(U) with −Δu‾≤f(x,u‾)-\Delta\underline u \leq f(x, \underline u) weakly, meaning ∫∇u‾⋅∇v≤∫f(x,u‾)v\int\nabla\underline u\cdot\nabla v \leq \int f(x, \underline u)v for all v∈H01v \in H^1_0 with v≥0v \geq 0, and u‾≤0\underline u \leq 0 on ∂U\partial U (in the sense u‾+∈H01\underline u^+ \in H^1_0). A supersolution u‾\overline u satisfies the reverse inequalities, with u‾≥0\overline u \geq 0 on ∂U\partial U.

Lemma 5.10 Weak maximum principle in H1H^1

Let λ≥0\lambda \geq 0 and w∈H1(U)w \in H^1(U) with w+∈H01(U)w^+ \in H^1_0(U) satisfy ∫∇w⋅∇v+λ∫wv≤0\int\nabla w\cdot\nabla v + \lambda\int wv \leq 0 for all v∈H01v \in H^1_0 with v≥0v \geq 0. Then w≤0w \leq 0 almost everywhere.

Proof. Test with v=w+v = w^+, using ∇w+=∇w\nabla w^+ = \nabla w on {w>0}\{w > 0\} and 00 elsewhere (4A.11 Sobolev Embeddings and Critical Exponents): ∫∣∇w+∣2+λ∫(w+)2≤0\int|\nabla w^+|^2 + \lambda\int(w^+)^2 \leq 0. So ∇w+=0\nabla w^+ = 0, and since w+∈H01w^+ \in H^1_0, Poincaré's inequality gives w+=0w^+ = 0.

Theorem 5.11 The method of sub- and supersolutions

If u‾≤u‾\underline u \leq \overline u are a subsolution and a supersolution of (P) and both are bounded, then (P) has a weak solution u∈H01u \in H^1_0 with u‾≤u≤u‾\underline u \leq u \leq \overline u.

Proof. Let m≤u‾≤u‾≤Mm \leq \underline u \leq \overline u \leq M, and choose λ≥0\lambda \geq 0 so that z↦f(x,z)+λzz \mapsto f(x, z) + \lambda z is non-decreasing on [m,M][m, M] (any λ\lambda at least the Lipschitz constant of ff there will do). Set u0=u‾u_0 = \underline u and define uk+1∈H01u_{k+1} \in H^1_0 as the weak solution of the linear problem

−Δuk+1+λuk+1=f(x,uk)+λuk,-\Delta u_{k+1} + \lambda u_{k+1} = f(x, u_k) + \lambda u_k,

which exists by Lax–Milgram whenever uku_k is bounded.

Claim: u‾=u0≤u1≤u2≤⋯≤u‾\underline u = u_0 \leq u_1 \leq u_2 \leq \cdots \leq \overline u. First, w=u0−u1w = u_0 - u_1 satisfies −Δw+λw≤(f(x,u0)+λu0)−(f(x,u0)+λu0)=0-\Delta w + \lambda w \leq \big(f(x, u_0) + \lambda u_0\big) - \big(f(x, u_0) + \lambda u_0\big) = 0 weakly, and w+∈H01w^+ \in H^1_0 because u0≤0u_0 \leq 0 on ∂U\partial U and u1∈H01u_1 \in H^1_0; by Lemma 5.10, u0≤u1u_0 \leq u_1. Next, if uk−1≤uku_{k-1} \leq u_k (both in [m,M][m, M]), then w=uk−uk+1w = u_k - u_{k+1} satisfies

−Δw+λw=(f(x,uk−1)+λuk−1)−(f(x,uk)+λuk)≤0,-\Delta w + \lambda w = \big(f(x, u_{k-1}) + \lambda u_{k-1}\big) - \big(f(x, u_k) + \lambda u_k\big) \leq 0,

by monotonicity of z↦f+λzz \mapsto f + \lambda z, so uk≤uk+1u_k \leq u_{k+1}. Finally, if uk≤u‾u_k \leq \overline u, then w=uk+1−u‾w = u_{k+1} - \overline u satisfies −Δw+λw≤(f(x,uk)+λuk)−(f(x,u‾)+λu‾)≤0-\Delta w + \lambda w \leq \big(f(x, u_k) + \lambda u_k\big) - \big(f(x, \overline u) + \lambda\overline u\big) \leq 0, with w+∈H01w^+ \in H^1_0 since u‾≥0\overline u \geq 0 on ∂U\partial U, so uk+1≤u‾u_{k+1} \leq \overline u. By induction the claim holds, and all uku_k take values in [m,M][m, M].

Convergence. The sequence increases pointwise and is bounded, so uk→uu_k \to u pointwise and in L2L^2 (dominated convergence), and f(x,uk)→f(x,u)f(x, u_k) \to f(x, u) in L2L^2 (Lipschitz). Testing the equation for uk+1u_{k+1} with uk+1u_{k+1} itself shows ∥∇uk+1∥2≤∥f(x,uk)+λuk∥ ∥uk+1∥\|\nabla u_{k+1}\|^2 \leq \|f(x, u_k) + \lambda u_k\|\,\|u_{k+1}\|, bounded, so a subsequence converges weakly in H01H^1_0, to uu. Passing to the limit in ∫∇uk+1⋅∇v+λ∫uk+1v=∫(f(x,uk)+λuk)v\int\nabla u_{k+1}\cdot\nabla v + \lambda\int u_{k+1}v = \int\big(f(x, u_k) + \lambda u_k\big)v gives ∫∇u⋅∇v=∫f(x,u)v\int\nabla u\cdot\nabla v = \int f(x, u)v.

The proof is a discrete version of a parabolic flow: uk+1u_{k+1} is one semi-implicit step of ut=Δu+f(x,u)u_t = \Delta u + f(x, u) started from the subsolution (with time step 1/λ1/\lambda). Indeed the solution of the heat flow started at a subsolution increases in time and converges to a solution, by the parabolic maximum principle (Exercise 5.18). Both constructions give the minimal solution above u‾\underline u: any solution u~\tilde u with u‾≤u~≤u‾\underline u \leq \tilde u \leq \overline u satisfies uk≤u~u_k \leq \tilde u for all kk by the same induction, so u≤u~u \leq \tilde u.

Example 5.12 The logistic equation

For −Δu=μu(1−u)-\Delta u = \mu u(1 - u) in UU, u=0u = 0 on ∂U\partial U, with μ>λ1\mu > \lambda_1 (the first Dirichlet eigenvalue), there is a solution with 0<u<10 < u < 1. Take u‾=1\overline u = 1: then −Δu‾=0=μu‾(1−u‾)-\Delta\overline u = 0 = \mu\overline u(1 - \overline u), and u‾≥0\overline u \geq 0 on the boundary. Take u‾=εϕ1\underline u = \varepsilon\phi_1, with ϕ1>0\phi_1 > 0 the first eigenfunction normalised by max⁡ϕ1=1\max\phi_1 = 1:

−Δu‾=λ1εϕ1≤μεϕ1(1−εϕ1)  ⟺  λ1≤μ(1−εϕ1),-\Delta\underline u = \lambda_1\varepsilon\phi_1 \leq \mu\varepsilon\phi_1(1 - \varepsilon\phi_1) \iff \lambda_1 \leq \mu(1 - \varepsilon\phi_1),

true for all xx once ε≤1−λ1μ\varepsilon \leq 1 - \frac{\lambda_1}{\mu}. And u‾≤u‾\underline u \leq \overline u. Theorem 5.11 gives a solution with εϕ1≤u≤1\varepsilon\phi_1 \leq u \leq 1; it is smooth by bootstrapping (6A.6 Weak Solutions and Elliptic Regularity), and u<1u < 1: w=1−u≥0w = 1 - u \geq 0 satisfies Δw−μuw=0\Delta w - \mu uw = 0, so the strong maximum principle with a zero-order term (6C.3 Symmetry, Noether’s Theorem and Moving Planes) gives w>0w > 0, since w=1w = 1 on the boundary. With μ=rD\mu = \frac rD this is the existence half of the critical patch size.

Where this goes Fixed points, estimates and order in geometric analysis

Hamilton's short-time existence for Ricci flow (1982) used the Nash–Moser implicit function theorem, a fixed-point argument with loss of derivatives, because the Ricci flow is only weakly parabolic; DeTurck's trick (6C.12 Quasilinear Parabolic Systems, 11A.4 Short-Time Existence and Uniqueness) reduces it to a strictly parabolic system where a contraction suffices. Sub- and supersolutions are the classical tool for the prescribed curvature problems of Jerry Kazdan and Frank Warner (1974–75) and for the negative case of the Yamabe problem (Exercise 5.19). The Leray–Schauder philosophy, existence from a priori estimates, is the structure of Yau's proof of the Calabi conjecture (6C.10 Fully Nonlinear Equations and Monge–Ampère) and of most long-time existence theorems for geometric flows: a flow exists as long as its curvature stays bounded, so existence is reduced to estimates.

History

Luitzen Brouwer proved his fixed point theorem in 1911–12 (the case n=3n = 3 had been found by Piers Bohl in 1904). Juliusz Schauder proved his theorem in 1930, and with Jean Leray built the degree theory for compact perturbations of the identity in 1934; Helmut Schaefer gave the form used here in 1955. Schauder was murdered by the Nazis in Lwów in 1943. Monotone operators were introduced by George Minty (1962) and Felix Browder (1963), with roots in electrical network theory, where Minty worked: the current–voltage relation of a nonlinear resistor is a monotone curve. The method of sub- and supersolutions in its monotone-iteration form was developed by David Sattinger (1972) and Herbert Amann in the early 1970s, building on older ideas going back to Perron's method (6B.6 Perron’s Method and Barriers). Nash's equilibrium paper appeared in 1951; he had announced the result, with a proof by Kakutani's fixed point theorem, in 1950.

Recall Where we stand

Brouwer: a continuous self-map of a compact convex set in Rn\mathbb{R}^n has a fixed point (proved by the null Lagrangian det⁡Du\det Du and the no-retraction lemma). Schauder: the same for compact convex sets in Banach spaces, by finite-dimensional projections. Schaefer: a compact map whose solutions of u=λA(u)u = \lambda A(u) are bounded has a fixed point, so a priori estimates give existence, as for −Δu+b(∇u)+μu=f-\Delta u + b(\nabla u) + \mu u = f. For monotone, coercive aa with linear growth, −div⁡a(∇u)=f-\operatorname{div}a(\nabla u) = f is solved by Galerkin and Minty's trick, which identifies weak limits by monotonicity; the pp-Laplacian is the main example. Ordered sub- and supersolutions trap a solution, found by monotone iteration, the minimal one above the subsolution. 6C.6 Quasilinear Elliptic Equations and Minimal Surfaces applies the Leray–Schauder method to quasilinear equations, where the a priori estimate is a gradient bound, and solves the minimal surface equation.

Exercises

Exercise 5.13 Brouwer in one dimension

Prove Brouwer's theorem for n=1n = 1 (f:[−1,1]→[−1,1]f : [-1, 1] \to [-1, 1] continuous) from the intermediate value theorem, and show by example that it fails for the open interval (−1,1)(-1, 1).

Solution

g(x)=f(x)−xg(x) = f(x) - x has g(−1)=f(−1)+1≥0g(-1) = f(-1) + 1 \geq 0 and g(1)=f(1)−1≤0g(1) = f(1) - 1 \leq 0, so gg vanishes somewhere. On (−1,1)(-1, 1), f(x)=x+12f(x) = \frac{x + 1}{2} maps into (0,1)⊂(−1,1)(0, 1) \subset (-1, 1) and its only fixed point would be x=1x = 1, not in the interval.

Exercise 5.14 Brouwer fails in ℓ2\ell^2

On the closed unit ball of ℓ2\ell^2, define T(x)=(1−∥x∥2,x1,x2,…)T(x) = \big(\sqrt{1 - \|x\|^2}, x_1, x_2, \ldots\big). Show that TT is continuous, maps the ball into the unit sphere, and has no fixed point. Which hypothesis of Theorem 5.5 fails?

Solution

∥T(x)∥2=1−∥x∥2+∥x∥2=1\|T(x)\|^2 = 1 - \|x\|^2 + \|x\|^2 = 1, and TT is continuous as the sum of a continuous scalar times e1e_1 and the (isometric) right shift. A fixed point would have ∥x∥=1\|x\| = 1, so x1=1−1=0x_1 = \sqrt{1 - 1} = 0, and then xk+1=xkx_{k+1} = x_k for all kk gives x=0x = 0, contradicting ∥x∥=1\|x\| = 1. The ball is not compact (and TT is not a compact map).

Exercise 5.15 The pp-Laplacian is monotone

For 1<p<∞1 < p < \infty and a(ξ)=∣ξ∣p−2ξa(\xi) = |\xi|^{p-2}\xi, show that (a(ξ)−a(η))⋅(ξ−η)≥(∣ξ∣p−1−∣η∣p−1)(∣ξ∣−∣η∣)≥0\big(a(\xi) - a(\eta)\big)\cdot(\xi - \eta) \geq \big(|\xi|^{p-1} - |\eta|^{p-1}\big)\big(|\xi| - |\eta|\big) \geq 0.

Solution

Expand: (a(ξ)−a(η))⋅(ξ−η)=∣ξ∣p+∣η∣p−(∣ξ∣p−2+∣η∣p−2)ξ⋅η\big(a(\xi) - a(\eta)\big)\cdot(\xi - \eta) = |\xi|^p + |\eta|^p - \big(|\xi|^{p-2} + |\eta|^{p-2}\big)\xi\cdot\eta. By Cauchy–Schwarz, ξ⋅η≤∣ξ∣∣η∣\xi\cdot\eta \leq |\xi||\eta|, so this is at least ∣ξ∣p+∣η∣p−∣ξ∣p−1∣η∣−∣η∣p−1∣ξ∣=(∣ξ∣p−1−∣η∣p−1)(∣ξ∣−∣η∣)|\xi|^p + |\eta|^p - |\xi|^{p-1}|\eta| - |\eta|^{p-1}|\xi| = \big(|\xi|^{p-1} - |\eta|^{p-1}\big)\big(|\xi| - |\eta|\big), which is non-negative because t↦tp−1t \mapsto t^{p-1} is increasing.

Exercise 5.16 The critical patch size

(a) Suppose −Δu=μu(1−u)-\Delta u = \mu u(1 - u) in UU, u=0u = 0 on ∂U\partial U, has a solution with 0<u<10 < u < 1 in UU. Multiply by ϕ1\phi_1 and integrate by parts to show μ>λ1\mu > \lambda_1. (b) In one dimension, for Du′′+ru(1−u)=0Du'' + ru(1 - u) = 0 on (0,L)(0, L) with zero boundary values, deduce that a positive steady state exists if and only if L>πD/rL > \pi\sqrt{D/r}.

Solution

(a) ∫(−Δu)ϕ1=∫u(−Δϕ1)=λ1∫uϕ1\int(-\Delta u)\phi_1 = \int u(-\Delta\phi_1) = \lambda_1\int u\phi_1, so λ1∫uϕ1=μ∫u(1−u)ϕ1<μ∫uϕ1\lambda_1\int u\phi_1 = \mu\int u(1 - u)\phi_1 < \mu\int u\phi_1, using uϕ1>0u\phi_1 > 0 and 1−u<11 - u < 1; since ∫uϕ1>0\int u\phi_1 > 0, λ1<μ\lambda_1 < \mu. (b) The equation is −u′′=rDu(1−u)-u'' = \frac rDu(1 - u) and λ1(0,L)=π2L2\lambda_1(0, L) = \frac{\pi^2}{L^2}. By (a) and Example 5.12, a solution with 0<u<10 < u < 1 exists if and only if rD>π2L2\frac rD > \frac{\pi^2}{L^2}, i.e. L>πD/rL > \pi\sqrt{D/r}. (A positive solution automatically has u<1u < 1: at an interior maximum above 11, −u′′≥0-u'' \geq 0 while rDu(1−u)<0\frac rDu(1 - u) < 0, so u≤1u \leq 1, and then u<1u < 1 as in Example 5.12.)

Exercise 5.17 Where the a priori estimate fails

For −u′′−μu=f-u'' - \mu u = f on (0,π)(0, \pi) with zero boundary values, the equation of Theorem 5.7 with b=0b = 0 and μ\mu replaced by −μ-\mu: show that when μ=1\mu = 1 and f=sin⁡xf = \sin x there is no solution, and identify which step of the proof of Schaefer's method breaks.

Solution

Multiply by sin⁡x\sin x and integrate by parts: ∫0π(−u′′−u)sin⁡x=∫0πu(−sin⁡′′x−sin⁡x)=0\int_0^\pi(-u'' - u)\sin x = \int_0^\pi u(-\sin'' x - \sin x) = 0, while ∫0πsin⁡2x=π2≠0\int_0^\pi\sin^2x = \frac\pi2 \neq 0. So no solution exists. In the method, solutions of u=λA(u)u = \lambda A(u), i.e. −u′′−u=λsin⁡x-u'' - u = \lambda\sin x, would have to stay bounded; but the energy estimate needs the coefficient of ∫u2\int u^2 to be positive, and −1+1=0-1 + 1 = 0 (the Poincaré constant is 11 here), so no bound is available: sin⁡x\sin x is in the kernel of −d2dx2−1-\frac{d^2}{dx^2} - 1, and that is the Fredholm obstruction of 6A.6 Weak Solutions and Elliptic Regularity.

Exercise 5.18 The heat flow from a subsolution increases

Let u‾\underline u be a smooth subsolution of −Δu=f(u)-\Delta u = f(u) on UU with u‾=0\underline u = 0 on ∂U\partial U, and u(x,t)u(x, t) the smooth solution of ut=Δu+f(u)u_t = \Delta u + f(u), u=0u = 0 on ∂U\partial U, u(x,0)=u‾(x)u(x, 0) = \underline u(x), existing for all time. Show that ut≥0u_t \geq 0. (Hint: v=utv = u_t satisfies a linear equation vt=Δv+f′(u)vv_t = \Delta v + f'(u)v; what is its sign at t=0t = 0, and on the boundary?)

Solution

Differentiating the equation in tt, v=utv = u_t solves vt=Δv+f′(u)vv_t = \Delta v + f'(u)v, with v=0v = 0 on ∂U\partial U (the boundary values are constant in time) and v(x,0)=Δu‾+f(u‾)≥0v(x, 0) = \Delta\underline u + f(\underline u) \geq 0, the subsolution inequality. The maximum principle for linear parabolic equations with bounded zero-order coefficient (6A.5 Maximum Principles; substitute v=eKtv~v = e^{Kt}\tilde v with K≥sup⁡∣f′(u)∣K \geq \sup|f'(u)| to make the zero-order term non-positive) gives v≥0v \geq 0 for all tt.

Exercise 5.19 Rehearsal: prescribing curvature −1-1 on a surface

Let (M,g)(M, g) be a closed surface with Gaussian curvature K<0K < 0 everywhere. The metric g~=e2ug\tilde g = e^{2u}g has curvature K~=e−2u(K−Δgu)\tilde K = e^{-2u}(K - \Delta_gu) (5A.5 Uniformization and the Two-Dimensional Ricci Flow). (a) Show that K~≡−1\tilde K \equiv -1 means −Δgu=−K−e2u-\Delta_gu = -K - e^{2u}. (b) Show that the constant cc with e2c=min⁡(−K)e^{2c} = \min(-K) is a subsolution and the constant CC with e2C=max⁡(−K)e^{2C} = \max(-K) is a supersolution. (c) The method of Theorem 5.11 works on a closed manifold, without boundary conditions (use −Δg+λ-\Delta_g + \lambda with λ>0\lambda > 0 in the iteration). Conclude that gg is conformal to a metric of constant curvature −1-1. This is the negative case of the uniformization theorem; Hamilton's Ricci flow on surfaces (11A.12 Ricci Flow on Surfaces) reaches the same metric as the limit of a flow, and in that proof the maximum principle plays the role that the ordering plays here.

Solution

(a) e−2u(K−Δu)=−1  ⟺  −Δu=−K−e2ue^{-2u}(K - \Delta u) = -1 \iff -\Delta u = -K - e^{2u}. (b) For the constant cc: −Δc=0-\Delta c = 0, and f(x,c)=−K(x)−e2c=−K(x)−min⁡(−K)≥0f(x, c) = -K(x) - e^{2c} = -K(x) - \min(-K) \geq 0, so −Δc≤f(x,c)-\Delta c \leq f(x, c): a subsolution. For CC: f(x,C)=−K(x)−max⁡(−K)≤0=−ΔCf(x, C) = -K(x) - \max(-K) \leq 0 = -\Delta C: a supersolution. And c≤Cc \leq C. (c) With f(x,z)=−K(x)−e2zf(x, z) = -K(x) - e^{2z}, Lipschitz on [c,C][c, C] with constant 2e2C2e^{2C}, the iteration −Δuk+1+λuk+1=f(x,uk)+λuk-\Delta u_{k+1} + \lambda u_{k+1} = f(x, u_k) + \lambda u_k is uniquely solvable on the closed manifold for λ>0\lambda > 0, and the weak maximum principle for −Δ+λ-\Delta + \lambda holds there (test with w+w^+, giving ∫∣∇w+∣2+λ∫(w+)2≤0\int|\nabla w^+|^2 + \lambda\int(w^+)^2 \leq 0). The proof goes through and gives a solution uu with c≤u≤Cc \leq u \leq C, smooth by elliptic regularity; e2uge^{2u}g has curvature −1-1.

Free resources for this chapter

Chosen for what this chapter covers, and all free to read. Links checked Oct 20, 2026.

© 2026 NeckPinch (www.neckpinch.com). All content in the guidebook (text, mathematics, figures and exercises) is protected by copyright. All rights reserved. No part may be copied, republished or redistributed without written permission.