# Heat kernel

Last updated

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.

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

${\displaystyle K(t,x,y)={\frac {1}{(4\pi t)^{d/2}}}e^{-|x-y|^{2}/4t}\,}$

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  Rd, where Δ is the Laplacian 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 Rd, 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 L2(U). The inverse Dirichlet Laplacian Δ−1 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.

## Related Research Articles

In mathematics and physics, the heat equation is a certain partial differential equation. Solutions of the heat equation are sometimes known as caloric functions. The theory of the heat equation was first developed by Joseph Fourier in 1822 for the purpose of modeling how a quantity such as heat diffuses through a given region.

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a function on Euclidean space. It is usually denoted by the symbols , , or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf(p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f(p).

In mathematics, the discrete Laplace operator is an analog of the continuous Laplace operator, defined so that it has meaning on a graph or a discrete grid. For the case of a finite-dimensional graph, the discrete Laplace operator is more commonly called the Laplacian matrix.

In the mathematical field of graph theory, the Laplacian matrix, also called the graph Laplacian, admittance matrix, Kirchhoff matrix or discrete Laplacian, is a matrix representation of a graph. The Laplacian matrix can be used to find many useful properties of a graph. Together with Kirchhoff's theorem, it can be used to calculate the number of spanning trees for a given graph. The sparsest cut of a graph can be approximated through the second smallest eigenvalue of its Laplacian by Cheeger's inequality. It can also be used to construct low dimensional embeddings, which can be useful for a variety of machine learning applications.

In differential geometry, the Laplace–Beltrami operator is a generalization of the Laplace operator to functions defined on submanifolds in Euclidean space and, even more generally, on Riemannian and pseudo-Riemannian manifolds. It is named after Pierre-Simon Laplace and Eugenio Beltrami.

In mathematics, the viscosity solution concept was introduced in the early 1980s by Pierre-Louis Lions and Michael G. Crandall as a generalization of the classical concept of what is meant by a 'solution' to a partial differential equation (PDE). It has been found that the viscosity solution is the natural solution concept to use in many applications of PDE's, including for example first order equations arising in dynamic programming, differential games or front evolution problems, as well as second-order equations such as the ones arising in stochastic optimal control or stochastic differential games.

In mathematics, a flow formalizes the idea of the motion of particles in a fluid. Flows are ubiquitous in science, including engineering and physics. The notion of flow is basic to the study of ordinary differential equations. Informally, a flow may be viewed as a continuous motion of points over time. More formally, a flow is a group action of the real numbers on a set.

In mathematics, the Fredholm alternative, named after Ivar Fredholm, is one of Fredholm's theorems and is a result in Fredholm theory. It may be expressed in several ways, as a theorem of linear algebra, a theorem of integral equations, or as a theorem on Fredholm operators. Part of the result states that a non-zero complex number in the spectrum of a compact operator is an eigenvalue.

In mathematics, Fredholm theory is a theory of integral equations. In the narrowest sense, Fredholm theory concerns itself with the solution of the Fredholm integral equation. In a broader sense, the abstract structure of Fredholm's theory is given in terms of the spectral theory of Fredholm operators and Fredholm kernels on Hilbert space. The theory is named in honour of Erik Ivar Fredholm.

The Minakshisundaram–Pleijel zeta function is a zeta function encoding the eigenvalues of the Laplacian of a compact Riemannian manifold. It was introduced by Subbaramiah Minakshisundaram and Åke Pleijel (1949). The case of a compact region of the plane was treated earlier by Torsten Carleman (1935).

In mathematics, an elliptic boundary value problem is a special kind of boundary value problem which can be thought of as the stable state of an evolution problem. For example, the Dirichlet problem for the Laplacian gives the eventual distribution of heat in a room several hours after the heating is turned on.

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

In mathematics, the Dirichlet eigenvalues are the fundamental modes of vibration of an idealized drum with a given shape. The problem of whether one can hear the shape of a drum is: given the Dirichlet eigenvalues, what features of the shape of the drum can one deduce. Here a "drum" is thought of as an elastic membrane Ω, which is represented as a planar domain whose boundary is fixed. The Dirichlet eigenvalues are found by solving the following problem for an unknown function u ≠ 0 and eigenvalue λ

In mathematics, especially spectral theory, Weyl's law describes the asymptotic behavior of eigenvalues of the Laplace–Beltrami operator. This description was discovered in 1911 by Hermann Weyl for eigenvalues for the Laplace–Beltrami operator acting on functions that vanish at the boundary of a bounded domain . In particular, he proved that the number, , of Dirichlet eigenvalues less than or equal to satisfies

A heat kernel signature (HKS) is a feature descriptor for use in deformable shape analysis and belongs to the group of spectral shape analysis methods. For each point in the shape, HKS defines its feature vector representing the point's local and global geometric properties. Applications include segmentation, classification, structure discovery, shape matching and shape retrieval.

In scientific computation and simulation, the method of fundamental solutions (MFS) is a technique for solving partial differential equations based on using the fundamental solution as a basis function. The MFS was developed to overcome the major drawbacks in the boundary element method (BEM) which also uses the fundamental solution to satisfy the governing equation. Consequently, both the MFS and the BEM are of a boundary discretization numerical technique and reduce the computational complexity by one dimensionality and have particular edge over the domain-type numerical techniques such as the finite element and finite volume methods on the solution of infinite domain, thin-walled structures, and inverse problems.

In mathematics, Sobolev spaces for planar domains are one of the principal techniques used in the theory of partial differential equations for solving the Dirichlet and Neumann boundary value problems for the Laplacian in a bounded domain in the plane with smooth boundary. The methods use the theory of bounded operators on Hilbert space. They can be used to deduce regularity properties of solutions and to solve the corresponding eigenvalue problems.

In mathematics, the Neumann–Poincaré operator or Poincaré–Neumann operator, named after Carl Neumann and Henri Poincaré, is a non-self-adjoint compact operator introduced by Poincaré to solve boundary value problems for the Laplacian on bounded domains in Euclidean space. Within the language of potential theory it reduces the partial differential equation to an integral equation on the boundary to which the theory of Fredholm operators can be applied. The theory is particularly simple in two dimensions—the case treated in detail in this article—where it is related to complex function theory, the conjugate Beurling transform or complex Hilbert transform and the Fredholm eigenvalues of bounded planar domains.

In mathematics, differential forms on a Riemann surface are an important special case of the general theory of differential forms on smooth manifolds, distinguished by the fact that the conformal structure on the Riemann surface intrinsically defines a Hodge star operator on 1-forms without specifying a Riemannian metric. This allows the use of Hilbert space techniques for studying function theory on the Riemann surface and in particular for the construction of harmonic and holomorphic differentials with prescribed singularities. These methods were first used by Hilbert (1909) in his variational approach to the Dirichlet principle, making rigorous the arguments proposed by Riemann. Later Weyl (1940) found a direct approach using his method of orthogonal projection, a precursor of the modern theory of elliptic differential operators and Sobolev spaces. These techniques were originally applied to prove the uniformization theorem and its generalization to planar Riemann surfaces. Later they supplied the analytic foundations for the harmonic integrals of Hodge (1940). This article covers general results on differential forms on a Riemann surface that do not rely on any choice of Riemannian structure.