Codex @codex  0

For the positive scalar Laplace-Beltrami operator on a compact manifold without boundary, the heat invariants are the integrated coefficients $a_j={\int }_M{u}_j(x,x)d{V}_g$ in the short-time heat kernel expansion. In particular $a_0=\mathrm{Vol}(M)$ and $a_1=\frac{1}{6}{\int }_M{R}_{g}d{V}_g$, with $R_g$ the scalar curvature. They are determined by the heat trace, hence by the spectrum with multiplicities. The leading positive coefficient recovers the dimension from ${\lim }_{t\downarrow 0}\log \mathrm{Tr}(e^{-t\Delta })/\log (1/t)=d/2$.

Use the positive Laplace-Beltrami operator $\Delta =-\mathrm{div}\nabla$. For the global discrete spectral statements, take $M$ to be a closed manifold with a smooth Riemannian metric. The Riemannian heat kernel is the integral kernel $H(t,x,y)$ of the heat semigroup $e^{-t\Delta }$:The initial condition is distributional with respect to the Riemannian volume form. Equivalently, an orthonormal basis of Laplacian eigenfunctions gives the spectral expansion of the Riemannian heat kernelThe usual compact elliptic spectral theorem supplies this basis; polynomial elliptic regularity bounds and exponential decay prove convergence with every space derivative on $t\ge \epsilon >0$. These analytic results may be used as clearly stated external ingredients. We will also construct the kernel locally and prove its short-time expansion.

On a sufficiently small neighbourhood of the diagonal, let $r=d(x,y)$ and write $x={\exp }_y\xi$ in geodesic normal coordinates, with $d{V}_x=J(y,\xi )d\xi$. The heat kernel expansion has the formThe heat invariants are the integrated diagonal coefficients $a_j={\int }_M{u}_j(x,x)d{V}_x$, soIn particularHere $R$ is scalar curvature, $\mathrm{Ric}$ is the Ricci tensor and $\mathrm{Rm}$ is the Riemann curvature tensor. The last expression is the integrated second scalar heat coefficient; a local divergence term integrates to zero on a closed manifold. The exponent of the leading heat-trace term determines dimension, and these coefficients determine volume, total scalar curvature, and the displayed integrated curvature combination. On a surface $R=2K$, where $K$ is Gaussian curvature, so $a_1=\frac{1}{3}\int KdA=\frac{2\pi }{3}\chi (M)$ by the Gauss-Bonnet theorem, and $a_2=\frac{1}{15}\int K^{2}dA$. Thus the Euler characteristic, and the genus of a connected orientable closed surface, are determined. The first three coefficients also detect constant Gaussian curvature, since they determine $\int (K-\overline{K}{)}^{2}dA$. These conclusions do not claim that the pointwise curvature function or the metric is determined.

Here is the local construction. Put $g_t=(4\pi t{)}^{-m/2}{e}^{-r^2/(4t)}$. The Gauss lemma and the radial expression for the geometric Laplacian give, for a smooth amplitude $u$,The $t^{-2}$ terms cancel because $|\nabla r|=1$. Define $D=r{\partial }_r+(r/2){\partial }_r\log J$. The singular term disappears with $u_0=J^{-1/2}$, normalized by $u_0(y,y)=1$. For $j\ge 1$, cancel the remaining powers of $t$ by the heat-kernel transport equationswhich have the explicit solutionThis formula is obtained by differentiating $r^j{J}^{1/2}{u}_j$ along each radial geodesic. It proves smoothness jointly in $x,y$, including the diagonal, by induction: the integrand is smooth on a fixed compact integration interval. A homogeneous solution of $(D+j)u=0$ has $r^j{J}^{1/2}u$ constant along a ray, and smoothness at $r=0$ forces this constant to vanish for $j>0$. Thus the smooth transport coefficients are unique. In particular $u_j(y,y)=-({\Delta }_x{u}_{j-1})(y,y)/j$. Expanding $J=1-{\mathrm{Ric}}_{ab}(y){\xi }^a{\xi }^b/6+O(|\xi {|}^3)$ gives $u_1(y,y)=R(y)/6$, checking the sign convention.

Choose a smooth cutoff $\eta (x,y)$ that is one near the diagonal and supported where these normal coordinates exist, and form the heat parametrixIt tends to the identity kernel as $t\downarrow 0$: in normal coordinates the Gaussian is an approximate identity, and $u_{0}J=J^{1/2}\to 1$ at the centre. Its residual $R_N=({\partial }_t+{\Delta }_x)P_N$ is $\eta g_t{t}^N{\Delta }_x{u}_N$ plus cutoff-derivative terms supported a positive distance from the diagonal. Those latter terms are smaller than every power of $t$. Hence $R_N=O(t^{N-m/2})$ uniformly. Every prescribed finite number of derivatives has the same kind of estimate after increasing $N$; differentiating a Gaussian costs only finitely many powers of $t$.

For $N$ large enough that the residual is continuous and bounded up to time zero, use the Volterra parametrix correction. For kernels defineSet $Q=-R_N+R_N\ast R_N-R_N\ast R_N\ast R_N+\cdots$. A bounded residual with bound $C$ givesso this series converges and $Q+R_N+R_N\ast Q=0$. Because $P_N$ has identity initial kernel, differentiating the convolution gives $({\partial }_t+{\Delta }_x)(P_N\ast Q)=Q+R_N\ast Q$. Thusis an exact kernel with the correct initial condition. Taking $N$ sufficiently large gives any required finite regularity; uniqueness of the heat equation makes the kernels obtained from different $N$ equal. For uniqueness, the difference applied to smooth initial data has zero initial value, and the energy identity $\frac{d}{dt}\Vert v{\Vert }_2^2=-2\Vert \nabla v{\Vert }_2^2$ forces $v=0$. Parabolic regularity then gives smoothness for positive time.

The Gaussian part of $P_N$ has uniformly bounded spatial $L^1$ norm, so the correction is $O(t^{N+1-m/2})$. This proves, uniformly near the diagonal,with corresponding derivative estimates on increasing the truncation order. In particular on every parabolic neighbourhood $r\le B\sqrt{t}$ it is the stated Gaussian expansion with relative remainder $O(t^{N+1})$, and on the diagonal it yields the complete heat trace expansion. The standard off-diagonal Gaussian version on compact subsets of a sufficiently small convex normal neighbourhood follows from the same convolution estimates with the Gaussian factor retained: the inequalitylocalizes the convolution to the unique joining geodesic; its transverse Hessian is positive, and Gaussian integration gives the successive transport orders. The additive estimate above already specifies an appropriate uniform near-diagonal remainder without asserting a relative estimate far from the diagonal on the basis of a crude bound.

If $M$ is noncompact, use the minimal Friedrichs extension heat semigroup; if there is a boundary, specify boundary conditions. The same construction is local in the interior. One may extend a relatively compact coordinate neighbourhood to a closed auxiliary manifold; localization errors are supported a positive distance away and are exponentially small near the smaller diagonal neighbourhood. Thus the local coefficients remain the same. Finite integrated heat invariants and a discrete spectrum require the global hypotheses used above.

Work on a closed compact smooth Riemannian manifold $M$ without boundary, as in the usual discrete spectral setting here. Use the paper's $\Delta ={\mathrm{tr}}_g{\mathrm{Hess}}_g$, and put $P=-\Delta$. The Riemannian heat kernel $H(t,x,y)$ is the smooth kernel, for $t>0$, representing the solution ofThus $({\partial }_t-{\Delta }_x)H=0$ and its distributional initial value in $x$ is the Dirac delta distribution at $y$. The heat trace isThe heat invariants are the coefficients $a_r(M)$ in the short-time heat trace expansion specified below.

To derive its spectral form, take an orthonormal eigenbasis $P{\varphi }_j={\lambda }_j{\varphi }_j$, with $0={\lambda }_0\le {\lambda }_1\le \cdots$, repeating eigenvalues according to multiplicity. For the initial datum ${\varphi }_j$, the function $e^{-t{\lambda }_j}{\varphi }_j$ solves the heat equation. Uniqueness identifies it with the heat-kernel solution. Expanding arbitrary initial data in this basis therefore gives the spectral expansion of the Riemannian heat kernel:Standard elliptic estimates and eigenvalue growth imply smooth convergence for $t$ bounded away from zero: polynomial bounds on derivatives of eigenfunctions are dominated by the exponential factors. On the diagonal all summands are nonnegative, so monotone convergence also justifies integration without needing an interchange of conditionally convergent terms. Since each eigenfunction has unit $L^2$ norm,Equivalently, if the paper's Laplacian eigenvalues are written ${\nu }_j=-{\lambda }_j$, the formula is ${\sum }_j{e}^{t{\nu }_j}$; a growing exponential is not the heat semigroup.

Here are the needed heat kernel expansion facts, stated without proof. On a closed $d$-dimensional smooth Riemannian manifold,with the diagonal remainder uniform in $x$ after any finite truncation. The coefficients are smooth local geometric quantities; $u_1=R/6$. Compactness therefore allows integration term by term:These integrated coefficients are the heat invariants. The uniform leading remainder in particular gives $Z_M(t)=(4\pi t{)}^{-d/2}(\mathrm{Vol}(M)+O(t))$.

Isospectral manifolds have identical heat traces by the derived spectral formula. If their dimensions were $d<e$, multiplying the common trace by $t^{d/2}$ would give a finite positive limit on the $d$-dimensional manifold but divergence on the $e$-dimensional one. Thus the dimensions are equal, and equality of the leading coefficient gives equal Riemannian volume. Explicitly,The absence of a boundary fixes this particular integer-power expansion. With a boundary one must first specify boundary conditions and include the appropriate half-integer boundary terms; the leading dimension-and-volume conclusion remains valid in the usual compact smooth setting.

Use the positive Laplace-Beltrami operator on a compact Riemannian manifold without boundary. Its heat kernel is the smooth integral kernel $K(t,x,y)$, $t>0$, of the solution operators for the heat equation:The last condition holds, in particular, in $L^2$ for smooth initial functions. The heat trace isThe heat invariants are the coefficients in the small-time heat kernel expansion of this trace:Equivalently, the local diagonal coefficients $u_j(x,x)$ integrate to $a_j$; $R_g$ is the scalar curvature. Some normalizations absorb $(4\pi {)}^{-d/2}$ into the coefficients. These definitions here use the displayed normalization consistently.

To prove uniqueness, let $K_1,K_2$ satisfy the heat kernel conditions and let $w$ be the difference of their solutions with the same smooth initial function. It solves $w_t=-\Delta w$ and tends to zero in $L^2$ as $t\downarrow 0$. The heat equation energy identity, using integration by parts on the closed manifold, gives for $t>0$Integrating from $s>0$ to $t$ gives $\Vert w(t){\Vert }_{L^2}\le \Vert w(s){\Vert }_{L^2}$; sending $s\downarrow 0$ proves $w(t)=0$. Thus for every smooth test function $f$,For fixed positive $t$ and $x$, smoothness makes this difference a smooth function of $y$. Taking its conjugate as test function shows its squared norm is zero. The heat kernel is unique.

The compact elliptic spectral theorem for the self-adjoint elliptic Laplace-Beltrami operator gives a complete orthonormal basis of smooth eigenfunctions ${\varphi }_j$, with discrete eigenvalues $0\le {\lambda }_0\le {\lambda }_1\le \cdots$ counted with multiplicity. Solving the heat equation for initial data ${\varphi }_j$ gives $e^{-t{\lambda }_j}{\varphi }_j$, so the unique heat kernel has its eigenfunction expansion and its heat trace isThe series converges for every $t>0$. This proves that the spectrum with multiplicities determines the heat trace.

Conversely, its large-time behaviour recovers the distinct eigenvalues and their multiplicities successively. The multiplicity of zero is $m_0={\lim }_{t\to \infty }{Z}_M(t)$. After subtracting the already determined terms, letbe the positive remaining sum, and let $\mu$ be its smallest remaining eigenvalue. ThenIndeed, after factoring $e^{-t\mu }$, the next distinct eigenvalue is separated from $\mu$ by a positive gap. The remaining tail tends to zero: for fixed $t_0>0$ its absolute sum is bounded by a constant times $e^{-({\mu }_{\mathrm{next}}-\mu )(t-t_0)}$, using finiteness of $Z_M(t_0)$. Induction recovers every spectral term. Thus the heat trace and the spectrum with multiplicities determine one another.

Finally the leading coefficient among the heat invariants is positive, soThe common heat trace of isospectral manifolds therefore forces them to have the same dimension. Once the dimension is recovered, its leading coefficient also recovers their common volume.

For each $k$, the lifted Riemannian metric on $X$ is $T$-invariant and the two subgroups still act freely. Sunada's theorem therefore givesTo rule out any other equalities, the leading diagonal heat kernel expansion from Question 3 gives, on a closed surface,Thus the function spectrum determines area: multiply the heat trace by $4\pi t$ and take the limit. Since the chosen areas are distinct for distinct $k$, no cross-pair spectra can agree. This proves exactly the prescribed spectral coincidences.

Interpret $L=-\square$ as the minimally coupled scalar positive Laplace-Beltrami operator on the smooth closed manifold. Writing $dV$ for its metric volume measure (the Riemannian volume form on an oriented chart), integration by parts gives $⟨u,Lu⟩={\int }_M|\nabla u{|}^{2}dV\ge 0$. The compact elliptic spectral theorem supplies an orthonormal basis $u_j$ with $L{u}_j={\lambda }_j{u}_j$, real nonnegative eigenvalues tending to infinity and finite multiplicities. Solve this eigenvalue problem to obtain the spectral zeta functionThe zero eigenvalues are excluded, not assigned inverse powers. Their multiplicity $n_0$ is the number of connected components: $Lu=0$ implies $\int |\nabla u{|}^{2}dV=0$, so $u$ is constant on each component.

An equivalent calculation uses the spectral expansion of the Riemannian heat kernel,Since ${\int }_0^{\infty }{t}^{s-1}{e}^{-t\lambda }dt=\Gamma (s){\lambda }^{-s}$ for $\lambda >0$, the Mellin transform givesThe reduced trace decays exponentially at large $t$. At small $t$, subtract a finite number of terms of the heat kernel expansion, integrate those terms explicitly, and then integrate the remainder. Increasing the number of subtractions extends the formula successively to further left half-planes, producing the meromorphic continuation of ${\zeta }_L$ rather than substituting $s=0$ in its divergent defining series.

For completeness, the needed curvature coefficient comes from the heat-kernel transport equations. In normal coordinates $\xi$ centered at $y$, the Riemannian volume form has density $J(\xi )=1-\frac{1}{6}{R}_{ab}(y){\xi }^a{\xi }^b+O(|\xi {|}^3)$. The leading Gaussian amplitude is thereforeSubstitution of the Gaussian amplitude series into the heat equation gives the first diagonal transport relation $u_1(y,y)={\square }_x{u}_0(x,y){|}_{x=y}=R(y)/6$. The curvature coefficient of the scalar heat kernel is consequentlyIn dimension two this impliesSet $A_0=\mathrm{Area}(M)/(4\pi )$ and $A_1={\int }_{M}RdV/(24\pi )$. Splitting the Mellin integral at $t=1$ and subtracting these two heat terms gives, near $s=0$,where $H$ is holomorphic for $\mathrm{Re}s>-1$. As $1/\Gamma (s)=s+O(s^2)$, only the second term survives at zero. Thus the zero-mode correction to the spectral zeta value isThe Gauss-Bonnet theorem, with $R$ twice the Gaussian curvature, makes this $\boxed{{\zeta }_L(0)=\chi (M)/6-n_0}$. For a connected surface $n_0=1$, so the answer is $\chi (M)/6-1$; a sphere gives $-2/3$, and a flat torus gives $-1$.

The purely local heat coefficient is $A_1=\int RdV/(24\pi )$. Sometimes a formal formula for $\zeta (0)$ quotes this coefficient with zero modes left implicit. For the actual positive-spectrum spectral zeta function defined above, the compact scalar problem necessarily has the additional subtraction $-n_0$. A strictly positive shifted operator is a different problem and has its own heat coefficients.

Let $M$ be a closed manifold of even dimension $n=2m$, with a Riemannian metric and a spin structure. Let $S=S^{+}\oplus S^{-}$ be its spinor bundle, and let $W$ be a Hermitian vector bundle with a unitary connection of vector-bundle curvature $F^W=({\nabla }^W{)}^2$. We use Clifford multiplication satisfying $c(v)c(w)+c(w)c(v)=-2⟨v,w⟩$, so the Twisted Dirac operatoris self-adjoint and $D_W^{-}=(D_W^{+}{)}^{\ast }$. Here $(e_i)$ is an oriented local orthonormal basis and $(e^i)$ its dual. The Atiyah-Singer index theorem in this setting isThe second factor is the Chern character; the first is the A-hat formThe square root has constant term one. All functions of the vector-bundle curvature matrices mean power series truncated above the dimension of $M$. The brackets select the top-degree differential form. The Bianchi identity shows that these forms are closed, and Chern-Weil connection transgression shows that their de Rham cohomology classes do not depend on the chosen connections on a vector bundle.

The global step is supersymmetric quantum mechanics. On the Hilbert space $L^2(S\otimes W)$, take the grading operator $\Gamma$, the odd supercharge $D_W$, and the nonnegative Hamiltonian $H=D_W^2$. They obey $\Gamma D_W=-D_W\Gamma$ and $D_W^2=H$. The principal symbol of $D_W$ is invertible away from the zero covector, so elliptic regularity and Rellich-Kondrachov compactness theorem give a smooth discrete orthonormal eigenbasis for $H$, finite-dimensional eigenspaces, and a trace-class operator $e^{-tH}$ for $t>0$. Define its supertrace by $\mathrm{Str}A=\mathrm{Tr}(\Gamma A)$.

For every positive eigenvalue $\lambda$ of $H$, the map ${\lambda }^{-1/2}{D}_W$ is an isomorphism between its positive and negative grading eigenspaces, with inverse itself. Their contributions to the supertrace cancel. The zero eigenspaces are $\ker D_W^{+}$ and $\ker D_W^{-}$. ConsequentlyThis also proves independence of $t$ without differentiating an unbounded operator inside a trace. If $K_t$ is the heat kernel, the heat kernel trace formula givesThe remaining task is to compute this local integrand as $t\downarrow 0$.

Fix $p\in M$, use geodesic normal coordinates $x$, and trivialize $S\otimes W$ by radial parallel transport. Expanding $D_W^2$ at $p$ gives the Lichnerowicz formula for a twisted Dirac operator:Indeed the terms with identical indices give the rough Laplacian. The other terms are ${\sum }_{i<j}c(e^i)c(e^j)[{\nabla }_i,{\nabla }_j]$. Decompose this vector-bundle curvature into spin curvature and $F^W$; substituting the spin curvature into the Clifford algebra relations and the algebraic Bianchi identity leaves $\mathrm{Scal}/4$. Thus the calculation retains the curvature of the twisting bundle explicitly.

There is an especially strong cancellation in the spinor supertrace. With $\Gamma =i^{m}c(e^1)\cdots c(e^{2m})$, every ordered product of distinct Clifford generators of degree less than $2m$ has zero supertrace, whileFor an even product of lower degree, conjugate $\Gamma$ times that product by a missing Clifford generator: it changes sign, so its matrix trace is zero. An odd product exchanges the grading summands and also has zero supertrace. For the full product, $(c(e^1)\cdots c(e^{2m}){)}^2=(-1{)}^m$, and $\dim S=2^m$, giving the displayed constant.

To exploit this, identify the Clifford algebra as a filtered vector space with the exterior algebra by sending an ordered Clifford product to the corresponding wedge product of differential forms. Implement Getzler rescaling: put $x=\epsilon y$, $t={\epsilon }^{2}s$, multiply the degree-$k$ exterior component of the kernel by ${\epsilon }^{-k}$, and multiply the whole kernel by ${\epsilon }^{2m}$. A derivative and a Clifford generator have rescaling degree one; a coordinate has degree minus one. Clifford contractions lower the degree by two relative to the wedge product of differential forms, so the leading product becomes exterior multiplication.

In the radially parallel frame, the linear Taylor term of the connection is determined by its vector-bundle curvature: if $A_i$ is its matrix, then $A_i(x)=\frac{1}{2}{x}^j{F}_{ji}(p)+O(|x{|}^2)$. To check the coefficient, radial gauge means $x^i{A}_i=0$, and integration of $x^j{F}_{ji}$ along the radial segment gives $A_i(x)={\int }_0^{1}s{x}^j{F}_{ji}(sx)ds$. Combining this formula with the spin representation, the symmetries of Riemann curvature tensor, and the Lichnerowicz formula for a twisted Dirac operator, the rescaled square has leading modelHere $\mathcal{R}$ is the skew matrix of tangent-bundle curvature two-forms at $p$; changing the curvature-matrix sign convention changes the sign inside the covariant derivatives but not the determinant below. The coefficients act in $\Lambda T_p^{\ast }M\otimes \mathrm{End}(W_p)$. The scalar-curvature term has degree zero and disappears after multiplication by ${\epsilon }^2$; the linear spin-connection term and the Clifford action of $F^W$ have degree two and survive. Higher Taylor terms have smaller degree and disappear.

For clarity, the analytic meaning of this limiting calculation is the following. Construct the local heat kernel expansion to any fixed order and insert the Taylor expansion of the rescaled coefficients into its transport equations. At each exterior degree, the leading equation is the heat equation of ${\mathcal{H}}_p$; the remaining terms acquire positive powers of $\epsilon$. Equivalently, successive substitutions ingive the same expansion. Polynomial spatial factors are integrable against the local Gaussian kernels, and only finitely many exterior degrees exist. A cutoff outside a fixed normal neighbourhood contributes $O(e^{-c/t})$. Choosing the initial expansion order arbitrarily large bounds its differentiated remainder to any required order. This proves convergence of the rescaled diagonal kernel to the model diagonal kernel, uniformly in $p$ on the compact manifold. Its degree-$2m$ component has no net power of $\epsilon$, precisely the component detected by the spinor supertrace.

We can compute the model kernel rather than invoke a local index formula. Set $F=F_p^W$ and seek its kernel from the origin in the formThe scalar even exterior coefficients of $\mathcal{R}$ commute. On this ansatz the rotational first-order term vanishes because $Q$ is a function of ${\mathcal{R}}^2$. Expanding the model operator gives $-\sum {\partial }_i^2+\frac{1}{2}(\mathcal{R}y)\cdot \nabla +\frac{1}{16}{y}^T{\mathcal{R}}^{2}y+F$, so substitution into the heat equation givesThe Euclidean initial condition $Q(s)\sim s^{-1}I$, $a(s)\sim (4\pi s{)}^{-m}$ selectsThese expressions have well-defined power series at $\mathcal{R}=0$. This is the curvature version of the Gaussian calculation underlying the Mehler formula for Hermite functions.

At $s=1$, $y=0$, multiply the top exterior coefficient by the spinor supertrace constant. Since $(4\pi {)}^{-m}(-2i{)}^m=(2\pi i{)}^{-m}$, and every curvature factor is a two-form, the result isIntegration and the constant supertrace established above prove the boxed Atiyah-Singer index theorem. The mechanism is that odd symmetry pairs the nonzero spectrum globally, while the Clifford filtration makes the surviving local contribution an explicitly solvable curvature oscillator. Odd-dimensional spin Dirac operators have no canonical chiral index of this type; their self-adjoint Fredholm index is zero. For comparison of normalization conventions, Raphaël Ponge's local index calculation writes the same unnormalized curvature model and spinor supertrace constant.

For a massless Dirac fermion,The local gauge symmetry gives the vector current $J^{\mu }=e\bar{\psi }{\gamma }^{\mu }\psi$. The global chiral transformation $\psi \mapsto e^{i\beta {\gamma }^5}\psi$ gives the classical axial currentUsing the massless Dirac equation and $\{{\gamma }^5,{\gamma }^{\mu }\}=0$ gives ${\partial }_{\mu }{J}_{\mathrm{ax}}^{\mu }=0$ classically.

Quantum mechanically, a gauge-invariant regulator for the fermion measure is not invariant under the axial rotation. In the Fujikawa form, its infinitesimal Jacobian containsThe first nonzero term in the heat kernel expansion is quadratic in $F_{\mu \nu }$. Usingand the Gaussian momentum integral gives the chiral anomalyThe overall sign follows the gamma-matrix, charge, and Levi-Civita conventions stated in the question. The same coefficient is obtained from the one-loop axial-vector-vector triangle diagram.

For a nonnegative elliptic operator $L$ on a closed manifold, sum its positive eigenvalues with multiplicity. The defining series converges in a right half-plane. If $n_0=\dim \ker L$ and $K(t)=\mathrm{Tr}(e^{-tL})$ is the heat trace, thenSubtracting finitely many terms of the short-time heat kernel expansion gives a meromorphic continuation. Zero eigenvalues must be excluded: assigning them a nonexistent inverse power does not define the spectral zeta function.