research33 min studyLesson 815 in Mathematics

Heat-kernel asymptotics and geometric coefficients

Read volume and curvature from the short-time trace while keeping boundary and zero-mode corrections explicit.

On this page

What you’ll understand

  • Relate the heat kernel diagonal to the multiplicity-weighted heat trace.
  • Identify the scalar volume and integrated-curvature coefficients on a closed manifold.
  • Calculate interval boundary corrections and avoid applying closed formulas to boundary problems.

Heat flow suppresses large eigenvalues and turns the full spectrum into a smooth function of positive time. At very short times, its trace has an expansion with geometric coefficients. We use the positive scalar Laplacian and heat operator e−tΔe^{-t\Delta}. The closed-manifold coefficients below assume a smooth compact manifold without boundary; standard boundary problems have additional half-integer terms, and singular boundaries can introduce further changes. This separation is essential before interpreting a measured constant term as curvature.

The diagonal kernel records all modes

The heat kernel propagates an initial temperature by integration against volume. Its spectral expansion weights each eigenfunction pair by the exponential of minus time times its eigenvalue. Integrating the diagonal uses orthonormality to produce the sum of those exponentials, counting multiplicity. The sum converges for every positive time on a compact smooth manifold. A constant zero mode contributes one at every time, even though high-frequency behavior dominates the short-time limit.

The full trace approaches the number of connected components as time tends to infinity, because all positive modes decay. At short time, by contrast, many high modes contribute, and a finite spectral truncation cannot reproduce the complete asymptotic uniformly down to zero time. One must distinguish a mode experiment at a fixed positive time from an infinite-spectrum short-time expansion. Their limiting operations occur in a particular order and cannot be exchanged without estimates.

K(t,x,y)=∑je−tλjϕj(x)ϕj(y)‾,Z(t)=∫MK(t,x,x)dV=∑je−tλjK(t,x,y)=\sum_je^{-t\lambda_j}\phi_j(x)\overline{\phi_j(y)},\qquad Z(t)=\int_MK(t,x,x)dV=\sum_je^{-t\lambda_j}
This is the full trace, including all zero modes and multiplicities.↔ Scroll wide equations horizontally.

Closed scalar coefficients

For a closed d-dimensional manifold, the leading bracket coefficient is volume and the coefficient multiplying time is one sixth of the integrated scalar curvature. Higher terms involve curvature contractions and derivatives. These are scalar-function coefficients: changing to a form Laplacian, adding a potential, or changing an operator connection changes the formulas. Integrated coefficients are spectral invariants, but a single integrated curvature number does not determine the pointwise curvature function.

An asymptotic expansion means that subtracting finitely many terms leaves an error of the stated smaller order as time approaches zero. It need not be a convergent power series for a fixed positive time. On a closed surface the integrated scalar curvature equals four pi times the Euler characteristic, so the constant heat coefficient is the Euler characteristic divided by six. The ever-present zero mode is already included in the full spectral sum; it should not be added again to that geometric constant.

Z(t)∼(4πt)−d/2[Vol⁡(M)+t6∫MScal⁡ dV+O(t2)]Z(t)\sim(4\pi t)^{-d/2}\left[\operatorname{Vol}(M)+\frac t6\int_M\operatorname{Scal}\,dV+O(t^2)\right]
This smooth closed scalar-Laplacian formula has no boundary half-powers.↔ Scroll wide equations horizontally.

Boundary terms change the interpretation

For a smooth compact manifold with smooth boundary, Dirichlet or Neumann conditions add a term inside the bracket equal to minus or plus one half of square root pi times boundary area times square root time. Subsequent coefficients include boundary curvature contributions as well as interior curvature. In dimension one, an interval has two boundary points, so this correction gives a constant of negative one half for Dirichlet and positive one half for Neumann. The different constants are not evidence for different interior scalar curvature.

For polygonal domains, corner angles contribute to later terms and a smooth-boundary formula cannot simply be reused. Likewise, removing the zero eigenspace changes the trace by its finite dimension and therefore changes its constant term, although the leading volume singularity is unchanged. When comparing heat data, identify the boundary realization and whether the trace is full or reduced. Those choices determine which coefficient can legitimately be interpreted as an audible geometric quantity.

The first two sphere heat coefficients

For the round two-sphere of radius two, calculate the leading and constant terms of the full scalar heat trace. Use its area and scalar curvature, then check how removing the single zero mode changes the result.

  1. The area is 4π⋅22=16π4\pi\cdot2^2=16\pi. In dimension two the scalar curvature is twice Gaussian curvature, so it is 2/22=1/22/2^2=1/2. Its integral is therefore 8π8\pi.

  2. Substitution gives Z(t)=16π/(4πt)+(8π/6)/(4π)+O(t)=4/t+1/3+O(t)Z(t)=16\pi/(4\pi t)+(8\pi/6)/(4\pi)+O(t)=4/t+1/3+O(t). The constant agrees with Euler characteristic two divided by six.

  3. The sphere is connected, so the full trace includes one constant eigenfunction. The reduced trace is Z(t)−1=4/t−2/3+O(t)Z(t)-1=4/t-2/3+O(t). Removing that mode changes the constant term and does not alter the volume coefficient.

The result

The full constant term is one third. Adding another one for the constant eigenfunction would double-count a mode already present in the trace.

Explore the model

Interactive lab

Interactive lab · heat

A mode is one term, not the whole trace

Pause/Reset, then set both mode to 1 and diffusivity to 0.5. Predict amplitude exp(−0.25) after one Step. Reset and set both mode to 3 and diffusivity to 0.5; predict exp(−2.25). Compare how a high mode is more strongly suppressed at the same positive time.

A mode is one term, not the whole traceThe exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Current readings: Current amplitude: 1; Spatial L² norm: 1.772; Decay rate · κn²: 3.6 /s. All parameters are adjustable with the labeled controls below.Periodic heat equation mode 3 at time 0: the amplitude is 1.0π/2π3π/22π−1−0.500.51x · periodic intervalTemperature perturbation uExact mode · n = 3
Current temperature modeInitial sineCurrent amplitude
Current amplitude
1
Spatial L² norm
1.772
Decay rate · κn²
3.6 /s
t = 0 s

The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe exact periodic mode illustrates individual exponential weights. It cannot recover the full diagonal heat kernel or curvature coefficients from one sine wave.

Check your understanding

For a length-L Dirichlet interval, what is the constant term after its leading L divided by square root of 4 pi t?

Try it yourself

Compare the full Dirichlet and Neumann heat traces of a length-three interval at short time. State their leading and constant terms and determine their exact difference for every positive time.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Their positive spectra agree term by term, while the Neumann realization includes one extra zero eigenvalue.

Compare with the solution

The Dirichlet trace is 3/4πt−1/2+o(1)3/\sqrt{4\pi t}-1/2+o(1) and the Neumann trace is 3/4πt+1/2+o(1)3/\sqrt{4\pi t}+1/2+o(1). Their difference is exactly one for every positive time because the Neumann spectrum adds only the constant zero mode to the same positive eigenvalue sequence. The difference of the asymptotic constants matches this exact identity, showing how boundary conditions and zero-mode bookkeeping fit together.

Take it with you

  • The heat trace sums every eigenmode, including zeros unless reduced explicitly.
  • Closed scalar heat coefficients encode volume and integrated scalar curvature.
  • Boundary terms and zero-mode removal alter constant-term interpretations.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.