Take the 2-minute tour ×

The $\varepsilon$-regularity phenomenon shows up in several different contexts. I try to describe it focussing on the harmonic map situation, but I really would like to understand the situation in general. The following is the Schoen-Uhlenbeck $\varepsilon$-regularity lemma, extracted from Tobias H. Colding, William P. Minicozzi II, An excursion into geometric analysis.

Let $N$ be a Riemannian manifold and $B_{r}$ be the ball of radius $r$ centred at the origin in $\mathbf{R}^k$. Then there exists $\varepsilon(k,N)$ such that if $u:B_{r}\in\mathbf{R}^k\rightarrow N$ is an energy minimizing map and $$\frac{\int_{B_{r}}|\nabla u|^2}{r^{k-2}}<\varepsilon,$$ then $u$ is smooth in a neighborhood of $0$ and $$|\nabla u|^2(0)\leq \frac{C}{r}.$$

Thus if a (conformally invariant) rescaling of the energy that $u$ minimizes is small (I suppose $u$ should be in a suitable Sobolev space), then $u$ is automatically smooth in some smaller ball. This rescaling is monotonically increasing thanks to a monotonicity lemma. I am not sure how to interpret the bound on the derivative at zero, though. The $\varepsilon$-regularity lemma quickly implies that the singular set $S$ of $u$ has $(k-2)$-dimensional Hausdorff measure zero.

My questions are:

  1. What are the basic ingredients (I suppose I am talking about the properties of the energy functional here) that guarantee that such a lemma holds?
  2. What is the meaning of the supremum of the set of all $\varepsilon$ such that the energy bound holds, and how can it be computed?
  3. Is there a simple intuitive picture that I am missing that explains the situation?
  4. Is there an instance of this phenomenon that predates the Schoen-Uhlenbeck paper?

Many thanks.

share|improve this question
This isn't an answer, but one thing that may be helpful is to not consider the $\epsilon$-regularity lemma for classical solutions $u_i$. In this case the smoothness is not an issue (as you get it "standard elliptic estimates" i.e. Schauder theory) but rather the point is the uniformity (independent of the $u_i$) of the pointwise estimate. In particular, let $u_i$ have uniformly bounded energy. A subsequence will converge weakly to a weak solution $u$. $\epsilon$-regularity quantifies how the sequence can fail to smoothly converge. Namely, energy concentrating on small scales. –  Rbega Sep 19 '10 at 1:40
I should add another advantage of considering classical solutions is that there are then really slick proofs of $\epsilon$-regularity type theorems. A good example is the proof of the Choi-Schoen theorem or of the smooth version of Allard regularity both of which are (I believe) in Colding and Minicozzi's book in not in the "excursion". These sorts of proofs might may also provide you with some intuition as they strip out all the technicalities and really get at the essence. –  Rbega Sep 19 '10 at 1:45
Thanks! Regarding your first comment: is the pointwise estimate you refer to the estimate on the derivative at zero? –  hce Sep 19 '10 at 8:37
When you said '$\varepsilon$-regularity quantifies how the sequence can fail to smoothly converge' you reminded me of Tristan Rivière ( arxiv.org/abs/math/0304396 ), who calls the Schoen-Uhlenbeck lemma 'the earliest example of energy quantization in non-linear analysis'. I guess this means that once we know that the conformally invariant energy is less than $\Lambda$, but $\Lambda$ is greater than $\varepsilon'$ (the sup of all $\varepsilon$ for which regularity is guaranteed), then we can expect to have a singular set (I admit this is quite tautological). –  hce Sep 19 '10 at 8:39
This is why I find it interesting to understand what exactly $\varepsilon'$ is. –  hce Sep 19 '10 at 8:40

3 Answers 3

The way I think of it is to view semilinear PDEs, such as the harmonic map equation, as a contest between the linear portion of the equation ($\Delta u$ in this case) and the nonlinear portions (which, in the case of harmonic maps, are roughly of the shape $|\nabla u|^2$). Intuitively, if the nonlinear part is small compared to the linear part then we expect the linear behaviour to dominate. In the case of harmonic maps, this means that we expect the solutions to behave like solutions to Laplace's equation $\Delta u = 0$, which are known to be regular.

A bit of dimensional analysis then tells us that the condition $\frac{\int_{B_r} |\nabla u|^2}{r^{k-2}} < \varepsilon$ has the right scale-invariance properties to have a chance of making the nonlinear term smaller than the linear term. (To make this rigorous, one of course needs to deploy various harmonic analysis estimates in well-chosen function space norms, such as Sobolev embedding.)

I discuss these heuristics (though more for dispersive equations than for elliptic ones) a bit at


The question of what happens at the critical value of epsilon is an interesting one. Often, the limiting non-regular solutions at that value of epsilon, after rescaling and taking limits, tend to be quite symmetric and smooth, away from a very simple singular set (e.g. a subspace). I don't know the elliptic case too well, but one obvious candidate for such a solution would be a singular 2D harmonic map (such as the map from C -> S^1 given by x -> x/|x|) extended to k dimensions by adding k-2 dummy variables. In the dispersive case, the analogous concept is that of the minimal energy blowup solutions, and these tend to be soliton solutions (so, typically, they obey a time translation invariance symmetry), associated to the ground state solution of the associated time-independent equation.

share|improve this answer

Below is a rather longwinded description of the special case when the singularity is at worst an isolated point. I suspect you know all this already. The magic comes at the very end (see paragraph that starts with "Here's the critical trick"). I don't know if this is the same thing that gives Schoen-Uhlenbeck the extra oomph or not.

There are three applications I know of: Minimal hypersurfaces, self-dual Yang-Mills connections, and Einstein manifolds. The regularity theory described below is used for both a convergence theorem except possibly a finite number of points and a removable singularity theorem. These theorems are then used to establish the so-called bubbling phenomenon. The story below applies to the latter two applications; the details for minimal hypersurfaces are slightly different.

Assume for convenience that we're on a smooth $n$-dimensional complete Riemannian manifold, where $n > 2$. Denote the Laplacian on both functions and tensors by $\Delta = g^{ij}\nabla_i\nabla_j$.

Denote the $L_p$ norm of a function or tensor $u$ with respect to the Riemannian metric by $\|u\|_p$.

Throughout the discussion below we will restrict to a geodesic ball $B(x, r)$ and assume that the following Sobolev inequality holds for a fixed constant $C_S$ and any smooth function $u$ compactly supported in $B$: $$ \|\nabla u\|_2 \ge C_S\|u\|_{2n/(n-2)}. $$

First, you consider the scalar elliptic inequality $-\Delta u \le bu$, where $b$ can be viewed as a given potential function. Using Moser iteration, you show that if

$$ \|b\|_{q/2}, \|u\|_p < C, $$

where $q > n$, for some $p > 1$ on $B(x,r)$, then there is a bound on $\|u\|_\infty$ on, say, $B(x,r/2)$.

Second, you use Moser iteration to show that if $ \|b\|_{n/2} $ is sufficiently small (depends on $C_S$) on $B(x,r)$, then there is a bound for $\|u\|_{q/2}$ for some $q > n$ on $B(x,r/2)$.

Combining the first two shows that if $u$ satisfies $-\Delta u \le cu^2$ and $ \|u\|_ {n/2} $ is sufficiently small on $ B(x,r) $, then there is a bound on $ \|u\|_\infty $ on $ B(x,r/2) $.

In each application there is a curvature tensor $F$ that satisfies a PDE of the form $$ -\Delta F = Q(F), $$ where $Q$ depends quadratically on $F$. Moreover, there is a convergence theorem when there is a uniform pointwise bound on $ F $ (for Einstein manifolds you use the Cheeger-Gromov convergence theorem). Applying the results above to $u = |F|$ using coverings with smaller and smaller balls leads to a convergence theorem when there is a uniform bound on $ \|F\|_ {n/2} $ where the convergence can fail at only a finite number of points (where in the limit there is too much of $ \|F\|_{n/2} $ for the estimates above to hold).

Now you want to study the limit object near each point singularity. If you keep close track of the dependence on $r$ in the estimates above, the best you can do is a bound on $F$ that blows up like $r^{-2}$, where $r$ is the distance to the singularity. This is not enough to remove the singularity, so you need to use more than the elliptic PDE above.

Here's the critical trick: When doing Moser iteration on $u = |F|$, you use the standard Cauchy-Schwarz inequality to obtain the following pointwise inequality: $$ |F\cdot\nabla F| \le |F||\nabla F| $$ But in all of the applications, you have extra information about $F$ and its covariant derivative. In particular, $F$ and/or its covariant derivative have certain symmetries, which allow you to prove a pointwise bound of the form $$ |F\cdot\nabla F| \le c|F||\nabla F|, $$ where $c < 1$. This improvement when used with Moser iteration allows you to show that $F$ blows up more slowly than $r^{-2}$. Iterating this improvement leads to a uniform pointwise bound on $F$, which in turn allows the singularity to be removed using a straightforward geometric ODE argument.

The removable singularity theorem allows you to analyze both the limiting object with the bubbles removed as well as the bubbles themselves.

ADDED: I can't resist adding an anecdote to this: Right after I learned the trick in the paragraph above from a paper of Schoen-Simon-Yau, I went to a colleague's office to show it to him. As it happens, Eli Stein was there, and he exclaimed, "But it's in my book!" And indeed it is. You will find it presented very nicely in VII.3.1 "A subharmonic property of the gradient" of Stein's 1970 book, "Singular Integrals and Differentiability Properties of Functions". It is obvious that S-S-Y did not know this or forgot, because their proof is much messier than Stein's.

share|improve this answer
Might just be me but some of your text is coming out as latex and vice-versa. Is anyone else getting that? I can't edit. –  Spencer Sep 20 '10 at 20:32
Sorry, but I can't figure out how to fix it. The LaTeX input appears fine to me and in fact compiles fine using LaTeX itself. But for some reason it's not working here. –  Deane Yang Sep 20 '10 at 20:46
Bizarrely, the solution seems to be to put an extra space between the underscore and the bracketed n/2 in the affected region. I was inspired by Harald's comments in this meta thread tea.mathoverflow.net/discussion/637/… –  j.c. Sep 21 '10 at 3:26
Many thanks for this answer! I still haven't digested it properly, but I understand the critical trick. I will probably ask for details soon as I still find it difficult to see what's happening geometrically. –  hce Sep 21 '10 at 9:34
jc, thank you! I had learned the trick of adding spaces but it never occurred to me to do it there! –  Deane Yang Sep 21 '10 at 15:27

I can comment on the $\epsilon$-regularity lemma for 4-dimensional Einstein manifolds. Namely, there is an $\epsilon$ depending on the dimension and Sobolev constant so that $\int_{B_x(r)} |Rm|^2 dV_g < \epsilon$ implies that $\sup_{B_x(r / 2)} |Rm|^2 \leq Cr^{-2} \int_{B_x(r)} |Rm|^2 dV_g$.

The key ingredient that makes this lemma work is that for Einstein manifolds, the function $|Rm|$ satisfies an elliptic inequality (it's ``subharmonic" for some elliptic operator). From there a standard PDE argument using Moser iteration gives the $\epsilon$-regularity. (It's like a non-linear version of the mean value inequality for subharmonic functions in Euclidean space).

share|improve this answer
Thanks! Just to clarify: does the elliptic inequality lead to regularity once the inequality you give is satisfied, or does it lead to said inequality? Also, is there a known computation of the sup of all $\varepsilon$ in this case? –  hce Sep 19 '10 at 8:50

Your Answer


