If you want to kill your productivity, move to a new house!

August 3, 2009

I know things have been extremely slow lately, but I was moving last week (and well, that means packing everything, moving things to the new house, cleaning the old house, unpacking things at the new one… well, you get the picture).

My goal for this week is presenting Aleksandrov’s solution to  the Minkowski problem (see an earlier post I did introducing this problem). So I am going to leave you a problem as a preview, it is a sort of discrete version of the Minkowski problem:

Let n_1,...,n_k  be a family of non-coplanar unit vectors in \mathbb{R}^3 and let \alpha_1,...,\alpha_k be positive numbers such that

\sum \limits_{i=1}^k \alpha_in_i=0

Then, show that there exists a convex closed polyhedron with exactly k faces with normal vectors given by n_1,...,n_k and corresponding areas \alpha_1,...,\alpha_k. Plus, this polyhedron is unique up to translation.

This is fact is not surprising (since it is not hard to check that any polyhedron has this property), but the proof is far from trivial. As you may guess, the proof cannot be constructive, it will use a continuity argument to show that there must be at least one such polyhedron. I will present this in my next post.


Books of the trade

July 23, 2009

Pile of math books

Don’t be scared, this has not been turned into a photo blog!  Today I received my copy of Singular Integrals I ordered on-line recently, which is funny given that I have been reading that book for a long time now (at least now I won’t have to borrow a copy from the library or from a friend).

In any case, this reminded me of something I heard once: basically, that throughout the early years of your career as a mathematician, there will be a  list of books (that will depend strongly on your research interests) that you must read completely and in full detail to the point where you are be able to reproduce their contents on command. So, wondering what that list should be for me, I piled up some books from my book shelf and took a picture.

…perhaps I should be reading some of those books instead of writing this blog post


Solving the Monge-Ampere equation (continued… and finished).

July 18, 2009

I have been postponing this post for over a week due to lack of time, but finally here it is. This post ought to finish a series of past posts (here and here) where I have been describing the proof of existence of classical solutions to the Dirichlet problem for the Monge-Ampere equationvia the continuity method. Via the method of continuity  we reduced the question of existence of classical solutions to the problem of proving good a priori estimates for classical solutions, namely, we were trying to prove

Theorem (A priori estimate for the Monge-Ampere equation) Let u be a smooth solution of

det(D^2u)=\psi \mbox{ in } B_1

u= 0 \mbox{ on } \partial B_1

There is a constant C>0 (depending only on n and the C^3 norm of \psi) such that

||u||_{C^{2,\alpha}}\leq C

Last time we got almost there, using the maximum principle and the right barriers we proved the estimate

||u||_{C^{1,1}(B_1)}\leq C

Which is still not strong enough for our needs, so in order to finish the proof of the a priori estimate, we are going to use a powerful interior estimate for concave elliptic equations, proved independently by L.C. Evans and N. Krylov in the 80’s:
Theorem (Evans-Krylov) Let u be a C^{1,1} solution of the elliptic equation

F(D^2u)=f(x) for x \in B_1

if F is concave (or convex), then we have the following interior estimate

||u||_{C^{2,\alpha_0}(B_{\frac{1}{2}})}\leq C||u||_{C^{1,1}(B_1)}

where C depends only on F, and \alpha_0 is a universal constant.

This a well known theorem, a couple of places where one can read it are the book of Gilbarg and Trudinger (last edition), or the book of Caffarelli and Cabre. More recently, Caffarelli and Silvestre have come up with a shorter proof, still based on the original ideas of Evans and Krylov, this proof is available in arxiv. Maybe I will talk about the proof in some other post, but for now I am just going to quote the result.

The Evans-Krylov theorem is an interior result, we need also control at the boundary, that is provided by a result of Krylov

Theorem (Krylov) Let u be a solution to our equation, there is a universal \beta and constant C controlled by |u|_\infty such that  for x \in B_1 and y \in \partial B_1 we have

|D^2u(x)-D^2u(y)|\leq C|x-y|^\beta

This actually a corollary of Krylov’s theorem, which is a more general and remarkable result about equations in non-divergence form with measurable coefficients, but again I want to focus on the Monge-Ampere equation, I will talk about Krylov’s theorem some other time, a good place to read about it is the last chapter of Kazdan’s book. With these two tools its a standard argument to show that for a constant C controlled by the previous two, and for \alpha=min(\alpha_0,\beta), we have

|D^2u(x)-D^2u(y)|\leq C|x-y|^\alpha

I won’t do it in detail, but the proof is not too hard: basically, if the two points are closer to each other than to the boundary, then the Evans-Krylov estimate (properly scaled) gives us the inequality above, otherwise, the two points are closer to the boundary than to each other, so by the estimate of Krylov we get the same inequality in this case, and thats it!

…and that finishes the proof!. I did not present the most general result to simplify the presentation(at least for the weaker a priori estimates, which is where I did most of the details), but one can work in a more general domain (as long as it is convex) and have arbitrary boundary conditions. A much more general result, which includes not only the Monge-Ampere equation but also the k-Hessian equation, was proven  in a  paper by Caffarelli, Nirenberg and Spruck.


Prescribing the Ricci Curvature of Riemannian metrics

July 9, 2009

I would like to spend the next few posts talking about a problem I read about in chapter 5 of the book Einstein Manifolds by Arthur Besse (it turns out that the name Arthur Besse is made up as you can read in the preface of the book). The question that we want to address is the following

Given a compact smooth manifold without boundary M, when is it possible to find a Riemannian metric g in M satisfying Ric(g)=r for a given Ricci candidate r?

It is of course too ambitious to try to answer this question in full generality but we can start by showing some examples of Ricci candidates for which this equation does not have a solution.

Trying to solve Ric(g)= for g amounts to solving a second order, quasilinear PDE on g, however, the main difficulty here is that the operator g\mapsto Ric(g) is not elliptic.

A motivation for considering this problem comes from the question of existence of metrics with constant sectional curvature on 3- manifolds (compact and without boundary). This of course has to do with the celebrated theorem of Richard Hamilton on the description of 3- manifolds with positive Ricci curvature:

Theorem (Hamilton, 1982): Let M be a connected, compact smooth 3 dimensional manifold without boundary and suppose that M admits a metric g such that Ric(g) is positive definite everywhere. Then M also admits a metric with constant sectional curvature.

We will discuss some of the ideas involved in the proof of this theorem in future posts. A consequence of this result is that M is diffeomorphic to the quotient of the 3-sphere S^{3} by a discrete group G.

Back to our original problem, recall that given a Riemannian metric (M,g), the full curvature tensor is defined by

Rm(X,Y,Z,W)=g(R(X,Y)Z,W)=\langle R(X,Y)Z,W\rangle

Where

R(X,Y)Z=\nabla_{X}\nabla_{Y}Z-\nabla_{Y}\nabla_{X}Z-\nabla_{[X,Y]}Z

Here \nabla is the Levi-Civita connection of g.

The Ricci tensor is then defined as

Ric(X,Y)=\mathrm{tr}\left(Z\mapsto R(Z,X)Y\right).

Here are two basic properties of Ric:

1) Ric is symmetric in X and Y, i.e. Ric(X,Y)=Ric(Y,X)

2) In local coordinates Ric looks as follows

R_{ij}=\partial_{l}\Gamma_{ij}^{l}-\partial_{i}\Gamma_{lj}^{l}+\Gamma_{lm}^{l}\Gamma_{ij}^{m}-\Gamma_{im}^{l}\Gamma_{lj}^{m}

We are using the summation convention (i.e we sum over repeated indices). The Christoffel symbols \Gamma_{\alpha\beta}^{\gamma} are defined by

\Gamma_{\alpha\beta}^{\gamma}=\frac{1}{2}g^{\gamma\nu}\left(\partial_{\alpha}g_{\beta\nu}+\partial_{\beta}g_{\alpha\nu}-\partial_{\nu}g_{\alpha\beta}\right)

Where \left(g^{\gamma\nu}\right) are entries of the matrix g^{-1}. This says that in local coordinates we can write schematically Ric(g)=F(g,\partial g,\partial^{2}g) where F is a C^{\infty} function that depends linearly on the entries of \partial^{2}g. If we then want to solve Ric(g)=r locally, property 2) tells us that we have to look at a system of the form

F(g,\partial g,\partial^{2}g)=r

Property 1) tells us that an admissible Ricci candidate r has to be symmetric.

One encounters obstructions for solving this system right away. One of the main difficulties has to do with the fact that the Ricci tensor satisfies the differential Bianchi identity

\delta Ric(g)=-\frac{1}{2}dR

Where \delta is the divergence operator respect to g and R is the scalar curvature of g (the trace of the Ricci tensor). This says that if we define a 1-form Bian(r,g) by Bian(r,g)=\delta r+\frac{1}{2}d\mathrm{tr}(r), then in order for g to satisfy Ric(g)=r we must have

Bian(r,g)=0

To write Bian(r,g) in coordinates, we start with

Bian(r,g)_{k}=g^{ij}\left(\nabla_{i}r_{jk}-\frac{1}{2}\nabla_{k}r_{ij}\right)

From

\nabla_{a}r_{bc}=\partial_{a}r_{bc}-\Gamma_{ab}^{l}r_{lc}-\Gamma_{ac}^{l}r_{bl}

(which is just the definition of \nabla r in coordinates) and from the expression in local coordinates of the symbols \Gamma, we easily see that

Bian(r,g)_{k}=g^{ij}\left[\partial_{i}r_{jk}-\frac{1}{2}\partial_{k}r_{ij}-g^{ls}r_{ks}\left(\partial_{i}g_{js}-\frac{1}{2}\partial_{s}g_{ij}\right)\right]

As discussed in Besse’s book, Dennis DeTurck came up with examples of symmetric tensors that cannot satisfy the Bianchi identity respect to any metric. One of his examples is the following

Consider in \mathbb{R}^{n} a symmetric tensor of the form

r=x_{1}dx_{1}\otimes dx_{1}+\displaystyle{\sum_{1<i,j\le n}q_{ij}(x_{2},\ldots,x_{n})dx_{i}\otimes dx_{j}}

The existence of a metric g such that Ric(g)=r implies as we saw before that Bian(r,g)=0. In particular, from our expression for Bian we must have

0= Bian(r,g)_{1}=g^{ij}\left[\partial_{i}r_{j1}-\frac{1}{2}\partial_{1}r_{ij}-g^{ls}r_{1s}\left(\partial_{i}g_{js}-\frac{1}{2}\partial_{s}g_{ij}\right)\right]=

=\frac{1}{2}g^{11}+x_{1}g^{l1}g^{ij}\left(\partial_{i}g_{js}-\frac{1}{2}\partial_{s}g_{ij}\right)

This implies that on the hyperplane x_{1}=0, the metric g must satisfy g^{11}=0 which is impossible for a Riemannian metric. It follows that for any point p in the hyperplane x_{1}=0 the equation Ric(g)=r has no solution near p. Notice that at these points the tensor r is singular.

In the next post we will interpret the existence of examples like the one we have just discussed as a consequence of the non-ellipticity of the system g\mapsto Ric(g).


Antonio Ache joins!

June 23, 2009

Antonio Ache, a grad student at Wisconsin will be blogging with us starting this week (after all, this is meant to be a group blog!). He is a good friend and we have known each other since our undergraduate days.

Antonio will begin by telling us about the problem of finding Riemannian metrics with given Ricci curvature. Welcome Antonio!

PS: On an unrelated note, I am writing this post from my phone! the wordpress application for the iPhone works perfectly. Yes, this means I can now blog whenever I find myself lost during a math talk…


Solving the Monge-Ampere equation (continued)

June 21, 2009

In a previous post I began the proof of the following theorem:

Let \psi \in C^\alpha be a positive function in B_1, then there exists a unique function u \in C^{2,\alpha}(B_1) such that

det(D^2u(x))=\psi \mbox{ in } B_1

u(x) = 0 \mbox{ on } \partial B_1

To prove the theorem, we looked at the function

\psi_t = (1-t)+t\psi(x).

we noted that it is easy to solve the problem explicitely for \psi=\psi_0.  Then thanks to the inverse function theorem (the Banach space version) we saw that the set A of t’s for which we can solve the equation A is open (recall that A was defined as the set of those t in [0,1] for which we can solve the equation above with right hand side  equal to \psi_t).

Since A is open and non-empty, if we show that its also closed the theorem would be proved.

Part II:  Showing A is closed.

(see previous post of the series for Part I)

That A is closed will be shown to be a consequence of the following a priori estimate:

Theorem (A priori estimate for the Monge-Ampere equation) Let u be a smooth solution of

det(D^2u)=\psi \mbox{ in } B_1

u= 0 \mbox{ on } \partial B_1

There is a constant C>0 (depending only on n and the C^3 norm of \psi) such that

||u||_{C^{2,\alpha}}\leq C

Proving the a priori estimate is usually the hardest step, so first let’s see that the a priori estimate in fact implies that A is closed: let t_k be a sequence in A converging to a number t \in [0,1], by the a priori estimate above the sequence of solutions u^{(t_k)}=u^{(k)} is uniformly bounded in C^{2,\alpha}, thus by the Arzela-Ascoli theorem a subsequence converges to a function u \in C^{2,\alpha}_0, since \psi_{t_k} converges to \psi_t we see this u is a solution of Problem (t), i.e. t \in A. This proves A is closed and therefore the theorem.

Of course, now we “just” have to prove the a priori estimate. ..

The broad outline of its proof is: first we prove a series of a priori estimates, each stronger than the previous one, and at each given step we  use the estimates already achieved, this will bound the first and second derivatives, the main tool here is the maximum principle used as in the classical technique of Bernstein. Then, the Holder regularity of the second derivatives will be consequence of deep regularity results from the theory of fully non-linear equations, among them the Evans-Krylov theorem. Now lets go through the argument in detail, and to try to make this just a bit readable I will number the steps

Step 1) The easiest estimate is the L^\infty bound

(a priori estimate (1))

||u||_{L^\infty(B_1)} \leq C

where C depends only on the L^\infty norm of \psi. To prove it, consider the auxiliary function \phi(x)=\frac{M}{2n}(|x^2|-1), this is a smooth function that vanishes on \partial B_1 and D^2\phi = M^n, picking M large enough (and controlled by the supremum of |\psi|) we have

det(D^2\phi) \geq \psi = det(D^2u)

then the maximum principle implies that

u \geq \phi \geq -\frac{M}{2n} in B_1

plus u is a convex function  vanishing on \partial B_1, so we have u \leq 0 as well. This proves the a priori estimate (1).

Step 2) Now we go for a Lipschitz bound for u, which we will refer to as the

(a priori estimate (2))

||\nabla u||_{L^\infty(B_1)} \leq C

Let’s differentiate the equation det(D^2u)=\psi in the direction of the unit vector \xi, and we obtain

L u_\xi= a_{ij} (u_\xi)_{ij}=\psi_\xi

Here a_{ij}(x) correspond to the entries of the cofactor matrix of D^2u(x), what the left hand side above amounts to is the linearization of det(D^2u) at u, which is the elliptic operator L (we did this before, in the previous post). Notice also that I can talk freely about third derivatives of u because any strictly convex, C^{2,\alpha} solution of the Monge-Ampere equation is actually C^\infty (that this is true is another of the many consequences of  the Schauder theory for second order  linear elliptic equations with Holder continuous coefficients, and its not too hard to check), so the equation above holds classically. In particular, by the maximum principle (again) we conclude that for some C depending on the supremum of |\psi|:

\sup\limits_{B_1} |u_e| \leq \sup \limits_{\partial B_1}|\nabla u|+C

but u \equiv 0 on \partial B_1 thus there \nabla u = u_n n, where n is the outer normal to \partial B_1. Here we bring back our old friend, the function \phi defined in step 1), as we know u\geq \phi on B_1 and they agree on the boundary, therefore

u_n \leq \phi_n on \partial B_1

also, since u \leq 0 and vanishes on \partial B_1 we also have the inequality 0 \leq u_n. Therefore the |u_n| is  controlled by M times a dimensional constants, i.e.  the supremum of |\psi| controls |\nabla u|. Putting the last three identities/inequalities together we get the a priori estimate (2).

Step 3) Notice that up to now we have only used the ellipticity of the equation and not the particular structure of the Monge-Ampere equation (even the differentiation step of Step 2) is a fact that holds for general elliptic equations). This  will change now that we are going to prove the

(a priori estimate (3))

||D^2 u||_{L^\infty(B_1)} \leq C

We start off as in step 2 by differentiating along a fixed, arbitrary direction \xi,  except now we do it twice and thus we get something much more complicated than before. We get around this noting that the map M \to det(M) is a convex function on positive definite matrices, and get an inequality that allows us to drop the lower order terms:

Lu_{\xi \xi} \geq \psi_{\xi\xi} \geq -C

and this says that  \sup \limits_{B_1} u_{\xi\xi} is controlled by the same supremum on \partial B_1 plus a constant C controlled by \sup |\psi|.  Then we only have to estimate things on the boundary (just as for the first derivatives).

So lets pick a point x_0 \in \partial B_1 and estimate the entries of D^2u(x_0). Since on \partial B_1 we have u \equiv 0 all the first, and second order angular derivatives of u at x_0 are zero, and this says a lot about the Hessian of u on the boundary, its just a matter of relating partial derivatives in polar and Cartesian coordinates to see then (after doing an adequate rotation of our coordinate system) that the entries of the Hessian satisfy the relations:

u_{ij}(x_0)=u_n(x_0)\delta_{ij} for i,j<n

here “n” stands both for the number n and the exterior normal direction \partial B_1 since in our rotated coordinate system the vector e_n is normal to \partial B_1 at the point x_0.Since we already have an a priori bound for u_n this gives us an a priori bound for the second partial derivatives on the boundary, except for the partial derivatives of the form u_{in}, 0\leq i \leq n.

Here we use the fact that our equation is rotation invariant (that is, the determinant is invariant under rotations of the coordinate system). Since u_{ij} (i,j<n) is diagonal, we can compute at once det(D^2u(x_0))

det(D^2u(x_0))=\psi(x_0)=u_n^{n-1}u_{nn}

I claim that there exists a \delta_0>0  such that u_n(x_0)>\delta_0 on \partial B_1, and \delta_0 is bounded away from zero by a quantity controlled by \inf \psi. Therefore there exists a  C controlled by \psi, such that

u_{nn}\leq C\psi(x_0)\leq C\sup |\psi|

Plus, u_{nn} is non-negative (since u is convex). Thus we have given a priori estimates for all the pure second order derivatives of u on \partial B_1, and by the remarks made at the beginning of this step we have proved the a priori estimate (3)

Step 4) Putting the 3 previous estimates together, we have:

a priori estimate (4)

||u||_{C^{1,1}(B_1)}\leq C

Next, we have to push this to a C^{2,\alpha} estimate in  B_1. This will be the content of the next post (I am sure we can agree that this post is long enough). In the next post I  will talk about the Evans-Krylov theorem and the boundary behavior of solutions to elliptic equations. Hopefully that will wrap up this series of posts on the Dirichlet problem for Monge-Ampere,  after this I might go back to the Minkowski problem.


Quick post: A theorem and an open problem about Gauss curvature

June 17, 2009

While I finish my next post on the Dirichlet problem for the Monge-Ampere equation, I thought I could mention two neat things I have learned from reading Jerry Kazdan’s survey on Prescribing the curvature of a Riemannian manifold, a short book that I recommend strongly (its a bit out of date, but hey!, geometry has advanced quite a bit in the last couple of years!)

First, a cute solution to an easily-stated problem, as  observed by Wallach and Warner in 1970

Theorem:  Given a compact smooth 2-d manifold and \Omega is a 2-form such that \int_M \Omega = 2\pi \chi(M), then there exists a smooth Riemannian metric g on M such that \Omega = KdA, where K is the Gauss curvature of g and A is the area form of g.

Proof: Pick your favorite metric g_0 on M, we will  prove that the metric we are looking for is in fact conformal to g_0, for any smooth function u, define g_u = e^{2u}g_0. Using the well known formulas for the Gauss curvature and area form of g_u, one arrives at the identity

K_udA_u=K_odA_0-(\Delta_0 u) dA_0

where the sub-indexes 0 and u refer to the object corresponding to the metric g_0 and g_u, in particular, \Delta_0 is the Laplace operator for (M,g_0). Then one wants to pick u such that

K_odA_0-\Delta_0 u dA_0=\Omega

or (\Delta_0 u) dA_0 = K_0dA_0-\Omega

but, using the metric g_0 we can write \Omega = f(x)dA_0 for some f (or equivalently, using the Hodge star operator given by g_0) thus the u we want is the solution of the equation

\Delta_0 u = K_0-f (**)

but since \int_M \Omega =\int_MfdA_0= 2\pi\chi(M)=\int_MK_0dA_0 we have that \int_MK_0-fdA_0=0 thus by standard elliptic theory (or Fredholm theory, etc) there exists a smooth function u solving equation (**). Thus the metric g_u is the one we were looking for and the theorem is proved.

Now, an easily stated problem whose solution is not likely to be as short, and has led to a lot of research in the last 2 decades. I  read about it for the first time in an interview with Louis Nirenberg, which can be found here.

Open problem:  Given a two dimensional Riemannian manifold, can we embed it isometrically (even just locally) in three dimensional Euclidean space?

Recall that Gauss curvature solely determines the local geometry of a 2d manifold, so its not surprising that this problem is equivalent to the following problem involving the  (…tataaaaa! )  Monge-Ampere equation:

Open problem restated: Given a C^\infty function K(x,y), find (even just locally) a function u(x,y) which solves the equation

\frac{u_{xx}u_{yy}-u_{xy}^2}{\left (1+u_x^2+u_y^2 \right)^2} = K(x,y)

but mind you, this is not your grandpa’s Monge Ampere equation, for it is not an elliptic equation unless K is strictly positive. For general K, it is an equation that varies  between hyperbolic or elliptic according to the sign of K, so you run into real problems in the set of points where K vanishes .  If on the other hand K is strictly negative, Kazdan says that the equation above is solved using tools from non-linear hyperbolic equations, of which sadly I know nothing. The case K >0 is then dealt with the techniques from the last two and the next post, which I should finish in a day or two.


Solving the Monge-Ampere equation

June 14, 2009

In my previous post, while talking about the Minkowski problem, I introduced the Monge-Ampere equation. I would like to talk about it for a post or two before continuing with the Minkowski problem, in part because the Dirichlet problem for Monge-Ampere is arguably an easier problem than the Minkowski problem.

To make the presentation cleaner, allow me to work only in the unit ball in Euclidean space, this situation contains all of the important issues one deals with to solve the equation, and it will let us avoid some of the technicalities that give little insight. In any case, all of the arguments can be tweaked without too much effort (but perhaps losing some clarity of exposition) to make them work in any smooth, strictly convex domain \Omega, the same remark goes for dealing with non-zero boundary data,  I will make more detailed remarks about this in the next post. Here we go then.

The Dirichlet problem for the Monge-Ampere equation: Consider a strictly convex domain \Omega \subset \mathbb{R}^n, and a smooth, positive function \psi \bar{\Omega} \to \mathbb{R}. Then find a convex function u twice differentiable in \Omega such that:

det (D^2u(x)) = \psi(x) \mbox{ for all } x \in B_1

u(x)= 0 \mbox{ on } \partial B_1

That a solution exists can be seen via the method of continuity (although, as I pointed out, there are other approaches). The principle behind this method is that if one can show that the set of those positive \psi \in C^\infty(\bar{\Omega}) for which there is a solution is  open, closed and not empty then by connectivity there is a solution for any such positive \psi \in C^\infty . This is a very simple and universal method that has proven to be succesful in many situations in Mathematics (e.g. in many applications of PDEs to differential geometry), but in these post I will talk about its application to the Monge-Ampere equation.

From now on we fix \psi, and define

\psi_t(x)=(1-t)+t\psi(x)

this is a smooth function for every t, and strictly positive for every t \in [0,1]. Thus we can consider the Dirichlet problem above with \psi_t as the right hand side, lets call it “Problem (t)“. Problem (1) is the actual problem we want to solve, and Problem (0) t=0  we can solve explicitely, the solution is given by u_0(x)=\frac{1}{2}(|x|^2-1) which satisfies:

det(D^2u_0(x))=1 = \psi_0(x)\mbox{ in } B_1

u_0(x)= 0 \mbox{ on } \partial B_1

Consider the set

A=\{t \in [0,1] | \mbox{ Problem } (t) \mbox{ can be soved } \}

the example above shows that 0 \in A. Then the strategy is to show that A is at the same time open and closed in [0,1], which would imply the Dirichlet problem has a solution for any \psi. Here is where the real work begins, and we divide it in two parts.

Part I: Showing A is open.

This will be dealt with machinery from classical analysis, namely the inverse and implicit function theorems for Banach spaces. Let us consider the map

T : C^{2,\alpha}_0 \to C^\alpha

given by T(u)=det(D^2u), this is a continuous but non-linear map between these two Banach spaces, we will apply the inverse function theorem to this map. Suppose t \in A and let u^{(t)} be the solution to the correspoding problem, then we can assume without loss of generality that u^{(t)} \in C^{2,\alpha}_0 and that it is a strictly convex function. Then the differential to T at u^{(t)} can be computed by the chain rule, and it turns out that it is just the second order differential operator

L(v)= \sum \limits_{ij}a_{ij}(x)v_{ij}(x)

where a_{ij}(x) are the entries of the cofactor matrix of D^2u^{(t)}(x).  The fact that u^{(t)} is strictly convex means that L is an elliptic operator, moreover it’s coefficients are C^{\alpha}.Now I claim that L is a one to one bounded linear map* from C^{2,\alpha}_0 to C^\alpha.

This means that the differential of T at  u^{(t)} is regular, and by the inverse function theorem T is a diffeomorphism between a neighborhood of u^{(t)} in C^{2,\alpha}_0 and a neighborhood of T(u) in C^{\alpha}. But then we can find \delta>0 such that any s such that |t-s|<\delta the function \psi_s lies in this neighborhood, and thus one can find a u^{(s)} \in C^{2,\alpha}_0 such that det(D^2u^{(s)})=\psi_s,  i.e. Problem (s) has a solution and thus s \in A, which shows A is open.

*Remark: The fact that L is one to one is a consequence of the Schauder theory for linear 2nd order equations, but that is something that could be discussed by itself in several posts! The theory says explicitely that given a differential operator as L, strictly elliptic and with C^\alpha coefficients, then for any f \in C^\alpha(\Omega) there exists a unique u \in C^{2,\alpha}(\Omega) that vanishes on \partial \Omega and such that

Lu=f \mbox{ in } \Omega

That takes care of Part I. Showing that A is closed is usually the more interesting part and it changes from problem to problem (showing A is open is always more or less the same). I will do the second part in the next post


The Minkowski problem and the Monge-Ampere equation

June 4, 2009

So as it has been already noticed I haven’t posted much at all in a couple of months,  Sean has been doing a good job with the posting! Now that is summer and I have more free time I can get back to blogging. I won’t continue with my previous posts on minimal surfaces (at least for now), but rather talk for a couple of posts about something less heavy than geometric measure theory.

I would like to talk about the Minkowski problem. It is not only an  interesting and historically important problem in geometry but it also has deep connections with the theory of non linear elliptic equations.

The Minkowski problem goes as follows:

Given a strictly positive real function f defined on S^2, find a strictly convex compact surface \Sigma \subset \mathbb{R}^3 such that the Gauss curvature of \Sigma at the point x equals f(n(x)), where n(x) denotes the normal to \Sigma at x.

The simplest example would be when f(x) is identically equal to the constant \frac{1}{R^2}, in that case a solution to the Minkowski problem would be the sphere of radius R (it also turns out that the Minkowski problem has a unique solution for a given f, so the sphere is THE solution to the Minkowski problem for f \equiv R^{-2}).

It is my understanding (although I only know very little about the historical background) that this problem was studied intensively (among many, many other people) by A.D. Alexksandrov and  A. Pogorelov. Aleksandrov developed a theory of weak solutions and showed existence of weak solutions. Years later, via the work of many people (particularly people working in non-linear PDE) it was shown that these weak solutions were actually smooth, classical solutions, thus proving that the problem always has a unique solution.

How do PDE’s get in the picture? I will do a computation that will make this matter clear, but first let me point out the analogy with the theory of minimal surfaces and the prescription of the mean curvature of a surface: Minimal surfaces are surfaces for which the mean curvature vanishes, thus if the graph of a function happens to be a minimal surface this function solves a PDE (“the minimal surface equation”), which happens to be a quasilinear PDE. The problem of showing the existence of classical solutions for the minimal surface equation and other problems from the calculus of variations motivated much of the work in PDE through part of the 20th century (read for instance, about the work of Ennio De Giorgi on Hilbert’s 19th problem). For the Minkowski problem what we are prescribing is not the average of the principal curvatures of the surface (mean curvature) but the product of the principal curvatures (Gauss curvature), the corresponding PDE is not a quasilinear PDE, but a fully non linear PDE closely related to the Monge-Ampere equation.

To explain this further, allow me to recall a classical computation for the Gauss curvature.

Let us consider a portion of a surface given by the graph of a function u(x,y) (any small patch can be written this way if we rotate and translate our coordinate system), it is well known that one can express the Gauss curvature of this surface in terms of partial derivatives of u. I am going to assume that if you are reading this you are familiar with geometry of surfaces, if not, there are a couple of good wikipedia articles on this, such as the one on the Weingarten map (aka the Shape operator).

To compute the the Gauss curvature, I am going to use the fact that the Gauss curvature equals the determinant of the differential of the Gauss map. The Gauss map is the map that associates to each point on the surface the unique unit vector parallel to the unit normal at that point), this defines a map from  the surface \Sigma to the sphere S^2. Then its differential its just a map from the tangent plane at a point in \Sigma to a tangent plane to S^2.

In the case when the surface is the graph of u, in the coordinates (x,y) \to (x,y,u(x,y)) the Gauss map is given by

N_u(x)=\left ( \frac{-\nabla u}{\sqrt{1+\nabla u^2}}, \frac{1}{\sqrt{1+\nabla u^2}} \right )

Then the Weingarten map (or second fundamental form, or shape operator) its the differential of this map, its determinant is the Gauss curvature, and computing the differential (I am cheating and won’t write the differential here) we see that the Gauss curvature equals

\frac{det(D^2u)}{\left ( 1+|\nabla u|^2\right )^{\frac{n+2}{2}}}

Then, if the surface u parametrizes is a solution of the Minkowski problem, then

\frac{det(D^2u)}{\left ( 1+|\nabla u|^2\right )^{\frac{n+2}{2}}}=f(N_u(x)) = g(\nabla u)

where g is the obvious function. Thus we have arrived at a PDE for u

det (D^2u) = g(\nabla u) \left ( 1+|\nabla u|^2 \right)^{\frac{n+2}{2}}

this is an equation of Monge-Ampere type.

There are several methods to construct solutions to these equations (weak solutions, approximation via polygonal surfaces, continuity method), some of these were developed independently in works of some people like A. Alexandroff, A. Pogorelov, P.L. Lions and S.Y. Cheng and S.T, Yau. One of these methods consists in building first solutions in a “weak” sense, known as Alexandroff solutions, which a priori might not be twice differentiable, so they are not solutions of the equation in the classical sense, THEN one proves that these weak solutions must actually be smooth (this is is usually called a regularity theory, and its considered the hardest part).  I will talk a bit about the existence of classical solutions on a future post on the continuity method.


Software for Pretty Diagrams

April 30, 2009

I wanted to make some pretty pictures for posts on this blog, and so I had to decide how to go about making them.  I like the hand drawn look, but it’s also nice to be able to edit using the computer.  In the end, I settled on scanning hand-drawn images, editing them in Gimp a bit, and converting them to SVG to edit more in Inkscape.

Since this is a lot of work, I have tried to automate it to the extend possible.  The result is a tiny program which works, somewhat, sometimes.  It makes diagrams that look a bit like whiteboard drawings, so I’m calling it whiteboard.  It’s badly written, badly documented, and it might set your computer on fire.  If that doesn’t scare you away, please use it for whatever you like, and if you do something cool I’d love to hear about it.