Riemannian heat kernel ($H(t,x,y)$)

For the positive Laplace-Beltrami operator, the kernel solving $({\partial }_t+{\Delta }_x)H=0$ with initial Dirac delta distribution on the diagonal represents the heat semigroup. On a closed manifold with a Riemannian metric it is smooth for $t>0$, symmetric, nonnegative, preserves constants, and obeys the semigroup property.

For a finite free action by Riemannian isometries, the Riemannian heat kernel on the quotient is the image sum shown above. Integrating over a fundamental domain combines all images into the integral on the covering manifold, proving the initial condition and showing that there is no averaging factor in the kernel. The heat trace does have an averaging factor $1/|U|$, because the integral of a quotient function over the cover is $|U|$ times its quotient integral. Mere freeness of a general nondiscrete group action is not enough for this covering formula.

The normalized pullback is a unitary map from $L^2$ of the quotient onto the deck-invariant subspace of $L^2$ of the cover. Local isometry and the sheet-counting integration formula preserve the Dirichlet quadratic form. Therefore the Friedrichs extension and its heat semigroup intertwine with this map. A finite sum of lifted Riemannian heat kernels gives the quotient kernel; independence of lifts follows by reindexing the deck group. This formulation avoids a uniqueness claim for unrestricted heat-equation solutions on a noncompact or incomplete manifold.

On a closed manifold with a Riemannian metric, an orthonormal eigenbasis for the positive Laplace-Beltrami operator gives this formula, with eigenvalues counted with multiplicity. Polynomial elliptic bounds and exponential time decay give convergence of every derivative away from time zero. The conjugate is necessary for a complex basis.

An approximate Riemannian heat kernel with the correct delta initial limit and a residual extending with enough regularity to time zero. A cutoff Gaussian times transport coefficients supplies a local construction. Increasing the expansion order, or summing the full asymptotic series smoothly, makes the residual regular enough for a Volterra parametrix correction.

In normal coordinates $q={\exp }_{p}x$, write the Riemannian volume form as $J(p,x)dx$ and let $r=|x|$. With the heat operator ${\partial }_t-{\Delta }_q$, the Gaussian times a power-series amplitude has its singular terms cancelled by the displayed radial equations for $i\ge 1$. The initial coefficient is $w_0=J^{-1/2}$, andMultiplying the equation by $r^{i-1}{J}^{1/2}$ gives the derivative of $r^i{J}^{1/2}{w}_i$. Smoothness forces the integration constant to vanish. The integral on a fixed compact interval proves smoothness across the diagonal by induction; it yields $w_i(p,p)={\Delta }_q{w}_{i-1}(p,p)/i$. A truncated heat parametrix therefore has residual $-g{t}^k{\Delta }_q{w}_k$.

If $R=({\partial }_t+{\Delta }_x)P$ is smooth and bounded up to time zero on a closed manifold, the Volterra convolution of kernels solves $Q+R+R\ast Q=0$ by $Q={\sum }_{k\ge 1}(-1{)}^k{R}^{\ast k}$. Bounds by $C^k\mathrm{vol}(M{)}^{k-1}{t}^{k-1}/(k-1)!$ prove convergence. Then $H=P+P\ast Q$ is the exact Riemannian heat kernel and has the same initial delta limit.