Operator and coefficient

Let DD be the untwisted spin Dirac operator on a smooth Riemannian spatial manifold, away from any boundary, with spinor rank rsr_s. With positive connection Laplacian, the Lichnerowicz identity and local expansion give

D2=∇∗∇+R4I,KD2(t;x,x)∼1(4πt)d/2[I−tR(x)12I+O(t2)]. D^2=\nabla^*\nabla+\frac R4 I,\qquad K_{D^2}(t;x,x)\sim\frac1{(4\pi t)^{d/2}} \left[I-\frac{tR(x)}{12}I+O(t^2)\right].

Indeed a Laplace-type operator L=∇∗∇+VL=\nabla^*\nabla+V has first local correction RI/6−VRI/6-V. Substitution gives −RI/12-RI/12. Tracing multiplies the scalar terms by rsr_s. See Vassilevich, equations (3.27) and (4.27). Twisting connections or additional potentials change the endomorphism; this formula is not a substitute for calculating the specified V2 operator.

Here tt is a heat parameter with length-squared units when DD has inverse-length units. It is not the shared physical time coordinate.

Independent sphere benchmark

For the round two-sphere of radius aa, the untwisted spin Dirac eigenvalues are ±j/a\pm j/a, each of multiplicity 2j2j, j=1,2,…j=1,2,\ldots. Thus

Tr⁡e−tD2=4∑j=1∞je−tj2/a2=2a2t−13+O(t/a2). \operatorname{Tr}e^{-tD^2}=4\sum_{j=1}^\infty j e^{-tj^2/a^2} =\frac{2a^2}{t}-\frac13+O(t/a^2).

Euler–Maclaurin gives the constant term from −B2f′(0)/2!=−1/12-B_2 f'(0)/2!=-1/12 for f(s)=se−(t/a2)s2f(s)=s e^{-(t/a^2)s^2}, before multiplying by four. The local coefficient agrees: integrating rank two over area 4πa24\pi a^2 with R=2/a2R=2/a^2 gives exactly the constant −1/3-1/3.

Curvature is not an Euler characteristic

Rescaling a metric by a2a^2 changes its scalar curvature by a−2a^{-2}, whereas the Euler characteristic is unchanged. On a closed surface,

∫R dμ=4πχ. \int R\,d\mu=4\pi\chi.

This is an integrated identity. It does not make R(x)=χR(x)=\chi. For the round sphere, their numerical equality at radius one is a choice of units and metric. In particular, the round three-sphere has χ(S3)=0\chi(S^3)=0 but R=6/a2R=6/a^2, so its local Dirac heat correction is nonzero.

Applications to finite mode bands or mixing vertices require the associated operator, physical scale, and normalized matrix element. Numerical use at a particular nonzero tt additionally requires control of the asymptotic remainder.

interaction_operator_checks.py compares the sphere heat sum with its constant term at several radii and heat parameters, and bounds the omitted positive tail by an integral.