I have now discussed what the Laplacian looks like in a general Riemannian manifold and can thus talk about the basic equations of mathematical physics in a more abstract context. Specifically, the key ones are the Laplace equation
for a smooth function on a Riemannian manifold. Since , this often comes up when is the potential energy function of a field which is divergence free, e.g. in electromagnetism. The other major two are the heat equation
for a smooth function on the product manifold for a Riemannian manifold, and the wave equation
in the same setting. (I don’t know the physics behind these at all, but it’s probably in any number of textbooks.) We are often interested in solving these given some kind of boundary data. In the case of the Laplace equation, this is called the Dirichlet problem. In 2-dimensions for data given on a circle, the Dirichlet problem is solved using the Poisson integral, as already discussed. To go further, however, we would need to introduce the general theory of elliptic operators and Sobolev spaces. This will heavily rely on the material discussed earlier on the Fourier transform and distributions, and before plunging into it—if I do decide to plunge into it on this blog—I want to briefly discuss why Fourier transforms are so important in linear PDE. Specifically, I’ll discuss the solution of the heat equation on a half space. So, let’s say that we want to treat the case of . In detail, we have a function , continuous on . We want to extend to a solution to the heat equation which is continuous on and smooth on . To start with, let’s say that . The big idea is that by the Fourier inversion formula, we can get an equivalent equation if we apply the Fourier transform to both sides; this converts the inconvenience of differentiation into much simpler multiplication. When we talk about the Fourier transform, this is as a function of . So, assuming we have a solution as above:
Also, we know what looks like. So this is actually a linear differential equation in for each fixed with initial conditions . The solution is unique, and it is given by
Now recall that multiplication on the Fourier transform level corresponds to conveolution, and the Fourier transform of is . As a result, given a putative solution , we have determined by
So we have a candidate for a solution. Conversely, if the boundary data alone, it is easy to check by differentiation under the integral (justified by the rapid decrease of the exponential) that we have something satisfying the heat equation in the upper half-space. Moreover as by general facts about approximation to the identity and a look at the definition of —note that is just the orthodox version of an approximation to the idnetity. So, we have found a way to solve the heat equation on . It thus seems that the way to solve equations such as the heat equation by convolution with appropriate kernels. In fact, this is more generally true of nonhomogeneous constant-coefficient linear PDE on (we’re forgetting about boundary value problems). Suppose we are given a partial differential operator with constant coefficients, i.e.
where the ‘s are multi-indices. Then it is immediate that extends to an operator on distributions. Moreover,
whenever is a distribution and . (This is clear whenever ; in general any distribution can be approximated in the weak* sense by distributions by convolving with an approximation to the identity.) As a result, if we have a fundamental solution , i.e. one with
we can get a solution to any equation of the form for by taking
which is not only a distribution but also a polynomially increasing one. So we can solve any constant-coefficient PDE given a fundamental solution. There is a big theorem of Malgrange and Ehrenpreis that fundamental solutions always exist to constant-coefficient linear PDE. However, the above statement about solving PDEs can actually be proved in a more elementary fashion; perhaps this will be a future topic. For now, however, I want to show that the Gauss kernel is actually a fundamental solution to the heat equation, once it is extended to with for . (This is no longer smooth, but it is still a distribution.) We need to show that
Let’s take the integral where is integrated over ; then by integration by parts
Since is a solution to the heat equation on (e.g. look at the Fourier transform) it is the second integral that is nonzero. We can write this as
and it is easy to see (the same approximation to the identity argument) that the former term tends to zero as .
So we indeed have a fundamental solution to the heat equation. It thus seems fair that we get solutions to it by convolving with the Gauss kernel.