Heat kernel

In the mathematical study of heat conduction and diffusion, a heat kernel is the fundamental solution to the heat equation on a specified domain with appropriate boundary conditions. It is also one of the main tools in the study of the spectrum of the Laplace operator, and is thus of some auxiliary importance throughout mathematical physics. The heat kernel represents the evolution of temperature in a region whose boundary is held fixed at a particular temperature (typically zero), such that an initial unit of heat energy is placed at a point at time t = 0.

Fundamental solution of the one-dimensional heat equation. Red: time course of ${\displaystyle \Phi (x,t)}$. Blue: time courses of ${\displaystyle \Phi (x_{0},t)}$ for two selected points. Interactive version.

The most well-known heat kernel is the heat kernel of d-dimensional Euclidean space R^d, which has the form of a time-varying Gaussian function,

${\displaystyle K(t,x,y)=\exp \left(t\Delta \right)(x,y)={\frac {1}{(4\pi t)^{d/2}}}e^{-\|x-y\|^{2}/4t}\qquad (x,y\in \mathbb {R} ^{d},t>0)\,}$

This solves the heat equation

${\displaystyle {\frac {\partial K}{\partial t}}(t,x,y)=\Delta _{x}K(t,x,y)\,}$

for all t > 0 and x,y ∈ R^d, where Δ is the Laplace operator, with the initial condition

${\displaystyle \lim _{t\to 0}K(t,x,y)=\delta (x-y)=\delta _{x}(y)}$

where δ is a Dirac delta distribution and the limit is taken in the sense of distributions. To wit, for every smooth function φ of compact support,

${\displaystyle \lim _{t\to 0}\int _{\mathbf {R} ^{d}}K(t,x,y)\phi (y)\,dy=\phi (x).}$

On a more general domain Ω in R^d, such an explicit formula is not generally possible. The next simplest cases of a disc or square involve, respectively, Bessel functions and Jacobi theta functions. Nevertheless, the heat kernel (for, say, the Dirichlet problem) still exists and is smooth for t > 0 on arbitrary domains and indeed on any Riemannian manifold with boundary, provided the boundary is sufficiently regular. More precisely, in these more general domains, the heat kernel for the Dirichlet problem is the solution of the initial boundary value problem

${\displaystyle {\begin{aligned}&{\frac {\partial K}{\partial t}}(t,x,y)=\Delta K(t,x,y){\text{ for all }}t>0{\text{ and }}x,y\in \Omega \\[6pt]&\lim _{t\to 0}K(t,x,y)=\delta _{x}(y){\text{ for all }}x,y\in \Omega \\[6pt]&K(t,x,y)=0,\quad x\in \partial \Omega {\text{ or }}y\in \partial \Omega .\end{aligned}}}$

It is not difficult to derive a formal expression for the heat kernel on an arbitrary domain. Consider the Dirichlet problem in a connected domain (or manifold with boundary) U. Let λ_n be the eigenvalues for the Dirichlet problem of the Laplacian

${\displaystyle \left\{{\begin{array}{ll}\Delta \phi +\lambda \phi =0&{\text{in }}U\\\phi =0&{\text{on }}\ \partial U.\end{array}}\right.}$

Let φ_n denote the associated eigenfunctions, normalized to be orthonormal in L²(U). The inverse Dirichlet Laplacian Δ⁻¹ is a compact and selfadjoint operator, and so the spectral theorem implies that the eigenvalues satisfy

${\displaystyle 0<\lambda _{1}<\lambda _{2}\leq \lambda _{3}\leq \cdots ,\quad \lambda _{n}\to \infty .}$

The heat kernel has the following expression:

${\displaystyle K(t,x,y)=\sum _{n=0}^{\infty }e^{-\lambda _{n}t}\phi _{n}(x)\phi _{n}(y).}$

(1)

Formally differentiating the series under the sign of the summation shows that this should satisfy the heat equation. However, convergence and regularity of the series are quite delicate.

The heat kernel is also sometimes identified with the associated integral transform, defined for compactly supported smooth φ by

${\displaystyle T\phi =\int _{\Omega }K(t,x,y)\phi (y)\,dy.}$

The spectral mapping theorem gives a representation of T in the form

${\displaystyle T=e^{t\Delta }.}$

There are several geometric results on heat kernels on manifolds; say, short-time asymptotics, long-time asymptotics, and upper/lower bounds of Gaussian type.