Heat Kernel Expansion

• Heat kernel expansion is a systematic asymptotic series that expresses the fundamental solution of parabolic equations in terms of local geometric invariants.
• The expansion coefficients, such as the Seeley–DeWitt coefficients, encode detailed curvature, topological, and potential data essential for spectral analysis.
• Derivation methods include transport equations, parametrix construction, and the path integral approach, enabling applications in quantum effective actions and index theorems.

A heat kernel expansion is the systematic asymptotic expansion as $t\to 0^+$ of the fundamental solution $K(t;x,y)$ to the parabolic equation $(\partial_t + D_x)K(t;x,y) = 0$ for an elliptic or hypoelliptic operator $D$ on a (possibly filtered, possibly non-compact) manifold, typically with geometric or analytic structure. These expansions provide precise information about spectral invariants, quantum effective actions, index theorems, and the local and global geometry of the underlying space. Expansion coefficients (“Hadamard–Minakshisundaram–Pleijel” or “Seeley–DeWitt” coefficients) encode detailed geometric, curvature, and topological data, and admit a number of explicit constructions, generalizations, and analytic extensions.

1. General Structure and Canonical Local Expansions

On a compact Riemannian manifold $M^n$ equipped with a Laplace-type operator $D$ acting on sections of a vector bundle, the heat kernel admits the canonical short-time expansion (Minakshisundaram–Pleijel, Hadamard): $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$ where $d(x,y)$ is the Riemannian distance, and $a_k(x,y)$ are smooth bi-tensorial coefficients determined recursively from the geometry and the operator. On the diagonal, this yields

$$K(t;x,x) \sim (4\pi t)^{-n/2} \sum_{k=0}^\infty a_k(x,x) t^k,$$

and integrating over $K(t;x,y)$0 gives the heat trace expansion: $K(t;x,y)$1 where $K(t;x,y)$2 (Hong, 2011, Ludewig, 2016).

The coefficients $K(t;x,y)$3 are local geometric invariants built from the curvature tensors (Riemann, Ricci, scalar curvature) and their derivatives, the structure constants (in the group case), or the potential terms (in Schrödinger or Dirac–type operators). For Laplacians on compact Lie groups with bi-invariant metric, all such tensor invariants reduce to Ad-invariant polynomials of the structure constants and the scalar curvature is constant across $K(t;x,y)$4 (Hong, 2011).

2. Methods of Construction: Transport Equations, Parametrix, and Path Integral

The explicit computation of the $K(t;x,y)$5 proceeds by substituting the expansion into the heat equation $K(t;x,y)$6 and equating powers of $K(t;x,y)$7. This yields a sequence of transport equations along geodesics: $K(t;x,y)$8

$K(t;x,y)$9

with $(\partial_t + D_x)K(t;x,y) = 0$0 the Synge world function and $(\partial_t + D_x)K(t;x,y) = 0$1 the Van Vleck–Morette determinant (Ivanov et al., 2021, Ludewig, 2016).

Equivalently, in the parametrix approach (e.g., for Schrödinger operators), one constructs an approximate solution

$(\partial_t + D_x)K(t;x,y) = 0$2

by solving recursive transport equations for $(\partial_t + D_x)K(t;x,y) = 0$3 along straight lines or geodesics, subject to boundary/gauge data (Bolte et al., 2013).

The path-integral viewpoint gives an alternative derivation, interpreting the fundamental solution as a sum over paths, and yielding, at leading order, expressions for $(\partial_t + D_x)K(t;x,y) = 0$4 in terms of determinants of the Jacobi (second variation) operators along minimizing geodesics. Precisely,

$(\partial_t + D_x)K(t;x,y) = 0$5

where $(\partial_t + D_x)K(t;x,y) = 0$6 is the Van Vleck determinant; off-diagonal behavior and cut-locus effects are controlled via finite-dimensional Laplace integrals over spaces of minimizing geodesics (Ludewig, 2016).

3. Extensions: Hypoelliptic, Rockland, Higher-Order, and Non-local Expansions

Hypoelliptic and Filtered Geometries

For Rockland operators on filtered manifolds (generalizations of elliptic theory to Heisenberg-type and sub-Riemannian structures), a universal heat kernel expansion

$(\partial_t + D_x)K(t;x,y) = 0$7

emerges, where $(\partial_t + D_x)K(t;x,y) = 0$8 is the homogeneous dimension and $(\partial_t + D_x)K(t;x,y) = 0$9 the order (Dave et al., 2017). The relevant "Heisenberg calculus" replaces classical pseudo-differential structure, and coefficients $D$0 reflect the underlying filtration.

Higher-Order and Non-minimal Operators

For operators whose principal symbol is $D$1 with $D$2, the heat kernel admits an expansion in both positive and negative fractional powers of $D$3, with coefficients generated by recursion relations involving curvature and background field data (Barvinsky et al., 2021). The structure is

$D$4

where $D$5 are generalized exponential functions. In the coincidence limit $D$6, negative fractional powers vanish identically, ensuring compatibility with Seeley–Gilkey theory.

Non-local and Covariant Perturbation Expansions

The Barvinsky–Vilkovisky (BV) non-local expansion reorganizes the heat-kernel trace in terms of powers of local curvatures and exact non-local form factors $D$7, resumming infinitely many Seeley–DeWitt contributions at each order (Codello et al., 2012): $D$8 with coefficients for each tensor contraction given by explicit $D$9. These techniques allow systematic construction of one-loop effective actions and partition functions beyond the strictly local regime.

4. Special Geometries, Boundary Conditions, and Asymptotic Types

Compact Lie Groups

When $M^n$0 is a compact Lie group with bi-invariant metric, the Duflo isomorphism relates the heat kernel expansion for the Laplacian on $M^n$1 to the flat Laplacian on the Lie algebra $M^n$2 via the Jacobian of the exponential map and a structure-constant-dependent correction: $M^n$3 with $M^n$4, $M^n$5, $M^n$6, $M^n$7, $M^n$8 constant scalar curvature (Hong, 2011).

Domains with Boundary: Robin, Neumann, Dirichlet

For domains with Robin boundary conditions, the heat kernel converges (in spectral sum) for all real $M^n$9 and admits a local diagonal expansion with coefficients incorporating interior geometry, boundary second fundamental form, and parameter $D$0: $D$1 with explicit leading boundary terms, e.g. $D$2 (Meng et al., 21 May 2025).

Polynomial Potentials and Noncompact Spaces

For Schrödinger operators $D$3 with $D$4 a polynomially confining potential, the trace admits an expansion

$D$5

for confining $D$6 of degree $D$7 (Fucci, 2018, Fucci, 2014). Similar expansions with explicit residue-type computation of the coefficients hold for generic spherically symmetric polynomials.

5. Off-diagonal, Path Space, and Singular Geometric Regimes

Off-diagonal expansions and full two-point asymptotics are central for effective action computations, quantum field theory, and analysis near singularities such as the cut locus. For a pair of points $D$8 and generic second-order operators in flat space $D$9, Gou et al.\ provide an explicit covariant perturbative expansion up to arbitrary order in $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$0: $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$1 with nested Gaussian integrals and explicit analytic control as $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$2, reproducing all Seeley–DeWitt coefficients (Gou et al., 2016).

At the cut locus, or more generally for hypoelliptic/horizontal structures, the expansion takes the form

$$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$3

with $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$4, $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$5 the dimension of the manifold of minimizing controls, and the leading coefficient an oscillatory integral over $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$6: $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$7 (Inahama et al., 2016, Ludewig, 2016, Ludewig, 2016).

6. Noncommutative, Spectral, and Special Function Generalizations

Heat kernel expansions admit formulations in noncommutative geometry, number theory, and analysis on special spaces. For example, for operators whose spectrum is the imaginary parts of nontrivial zeros of the Riemann zeta function (assuming RH): $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$8 with explicit formulas for $$K(t;x,y) \sim (4\pi t)^{-n/2} \exp\left(-\frac{d(x,y)^2}{4t}\right) \sum_{k=0}^\infty t^k a_k(x,y),$$9 in terms of Bernoulli and Euler numbers (Connes, 2024).

Further, expansions involving special function families $d(x,y)$0 and $d(x,y)$1, constructed to diagonalize the Laplacian as shift operator and to expand Green's functions and inverse-heat-kernel transforms, regulate the singularity and provide versatile analytic building blocks (Ivanov et al., 2021).

Convolutions, Borel summation, and analytic continuation techniques enable uniform control and resummation across time scales and even connect solutions between hyperbolic, spherical, and Euclidean settings, as in the Borel summation and gamma-resummed expansions on hyperbolic spaces (Dunne, 2021).

---

References:

• (Hong, 2011): The asymptotic expansion of the heat kernel on a compact Lie group
• (Ludewig, 2016): Heat Kernel Asymptotics, Path Integrals and Infinite-Dimensional Determinants
• (Ludewig, 2016): Strong Short Time Asymptotics and Convolution Approximation of the Heat Kernel
• (Dave et al., 2017): The heat asymptotics on filtered manifolds
• (Barvinsky et al., 2021): Heat kernel expansion for higher order minimal and nonminimal operators
• (Codello et al., 2012): On the non-local heat kernel expansion
• (Fucci, 2018): Asymptotic Expansion of the Heat Kernel Trace of Laplacians with Polynomial Potentials
• (Meng et al., 21 May 2025): The Robin heat kernel and its expansion via Robin eigenfunctions
• (Gou et al., 2016): Covariant perturbation expansion of off-diagonal heat kernel
• (Bolte et al., 2013): Heat kernel asymptotics for magnetic Schrödinger operators
• (Ivanov et al., 2021): Special Functions for Heat Kernel Expansion
• (Connes, 2024): Heat Expansion and Zeta
• (Dunne, 2021): Borel Summation and Analytic Continuation of the Heat Kernel on Hyperbolic Space
• (Fucci, 2014): Heat Kernel Asymptotic Expansion on Unbounded Domains with Polynomially Confining Potentials
• (Dai et al., 2020): Witten Deformation on Non-compact Manifold: Heat Kernel Expansion and Local Index Theorem
• (Alonso-Izquierdo et al., 2013): Gilkey-de Witt heat kernel expansion and zero modes
• (Alonso-Izquierdo et al., 2019): Soliton Fermionic number from the heat kernel expansion