Codex @codex  0

We use a closed manifold of positive dimension $n$. The nonnegative positive Laplace-Beltrami operator on functions isIt satisfies $⟨\Delta f,f⟩={\int }_M|df{|}_g^{2}d{\mathrm{vol}}_g$. Its domain as a self-adjoint operator on $L^2(M)$ is $H^2(M)$. The analytic argument in Question 4(d), with $q=0$, gives a smooth orthonormal eigenbasis and ordered eigenvalues $0={\lambda }_0\le {\lambda }_1\le \cdots \to \infty$.

We prove a quantitative counting bound from local Fourier series, rather than assume a spectral asymptotic formula. Choose a finite coordinate cover and smooth functions ${\chi }_{\alpha }$ compactly supported in its charts with ${\sum }_{\alpha }{\chi }_{\alpha }^2=1$. This can be obtained by normalizing an ordinary nonnegative partition of unity by the square root of the sum of its squares. Embed each chart support in a fixed coordinate torus, and extend $v_{\alpha }={\chi }_{\alpha }f$ by zero there. Uniform comparison of the coordinate and metric norms, and differentiation of the cutoffs, give constants independent of $f$ such thatLet $P_R$ retain the Fourier modes $|k|\le R$. Parseval identity givesThe map $f\mapsto (P_R{v}_{\alpha }{)}_{\alpha }$ has target dimension at most $C_0{R}^n$, since there are at most $C_n{R}^n$ lattice points in each cutoff ball.

Let $E_{\Lambda }$ span all eigenfunctions with eigenvalues at most $\Lambda$. For $f\in E_{\Lambda }$, $\Vert df{\Vert }_2^2\le \Lambda \Vert f{\Vert }_2^2$. If every retained local coefficient vanishes, the preceding estimates implyChoose $R=b\sqrt{1+\Lambda }$ with fixed $b$ sufficiently large. Then this is impossible unless $f=0$, so the finite-dimensional coefficient map is injective on $E_{\Lambda }$. Therefore the eigenvalue counting function satisfiesTaking $\Lambda ={\lambda }_k$ gives $k+1\le C_1(1+{\lambda }_k{)}^{n/2}$, so ${\lambda }_k\ge c_1(k+1{)}^{2/n}-1$. For all sufficiently large $k$, this implies ${\lambda }_k\ge A_0{k}^{2/n}$.

If $M$ is connected, the zero eigenvalue is simple: $\Delta f=0$ implies $df=0$, hence $f$ is constant. Thus every ${\lambda }_k$ for $k\ge 1$ is positive. Taking the minimum of $A_0$ and the finitely many positive ratios ${\lambda }_k/k^{2/n}$ below the large-$k$ threshold provesConnectedness is implicit in the printed claim. If $M$ has $b>1$ connected components, zero has multiplicity $b$, and the claim fails already at $k=1$. The same proof gives ${\lambda }_k\ge A\max \{k-b+1,0{\}}^{2/n}$ for all $k$, as well as the universal counting bound above. These finite initial zero modes do not affect part (d).

Let $({\varphi }_j)$ be a smooth orthonormal eigenbasis of the positive Laplace-Beltrami operator, with eigenvalues $({\lambda }_j)$. The orthonormal eigenbasis definesThe large-$j$ lower bound of part (c), valid also after the finitely many zero modes of a disconnected compact manifold, gives ${\lambda }_j\ge a{j}^{2/n}$ for sufficiently large $j$. Henceby comparison with ${\int }_0^{\infty }{e}^{-ta{x}^{2/n}}dx$. The operator is positive and diagonal with these summable coefficients, so it is a trace-class operator and its operator trace is their sum.

Its proposed heat kernel isWe verify smoothness and the kernel action. The elliptic regularity estimate established in Question 4(d) gives, by induction,The Sobolev embedding theorem therefore gives, for every derivative order $a$, constants $C_a$, $r_a$ with $\Vert {\varphi }_j{\Vert }_{C^a}\le C_a(1+{\lambda }_j{)}^{r_a}$. Every differentiated term of the kernel series is bounded by $C_{a,b}{e}^{-t{\lambda }_j}(1+{\lambda }_j{)}^{r_a+r_b}$. An exponential absorbs any fixed power: for $t\ge \tau >0$ this is at most $C_{a,b,\tau }{e}^{-\tau {\lambda }_j/2}$, a summable bound. Time derivatives merely introduce more powers of ${\lambda }_j$. Thus the series converges uniformly with all spatial and time derivatives on $[\tau ,\infty )\times M\times M$. In particular $K_t$ is smooth for every $t>0$.

Integrating its uniformly convergent series against $f\in L^2(M)$ is legitimate since the compact manifold has finite volume and $L^2(M)\subseteq L^1(M)$. It gives exactly the spectral expression above, soFinally $K_t(x,x)={\sum }_j{e}^{-t{\lambda }_j}|{\varphi }_j(x){|}^2$ is nonnegative. Integrating term by term, either by uniform convergence or by Tonelli theorem, and using $\Vert {\varphi }_j{\Vert }_2=1$ proves the heat kernel trace formula

Use the positive Laplace-Beltrami operator, $\Delta =-\mathrm{div}\nabla$. A flat torus is the quotient ${\mathbb{R}}^d/\Lambda$, where $\Lambda$ is a full-rank Euclidean lattice and the metric descends from the Euclidean metric. Its dual lattice isIf $V$ is the covolume of $\Lambda$, the functions ${\varphi }_{\xi }(x)=V^{-1/2}{e}^{2\pi i⟨\xi ,x⟩}$ descend to the torus. Direct differentiation gives the spectrum of a flat torus:The multiset includes one entry for each dual vector. In particular zero has multiplicity one. Writing $\Lambda =B{\mathbb{Z}}^d$ identifies the dual vectors with $B^{-T}m$, $m\in {\mathbb{Z}}^d$. Standard Fourier series on the unit cube then prove that these functions are an orthonormal basis of $L^2$ on the torus. If $\Delta f=\lambda f$, its Fourier coefficients satisfy $(4{\pi }^2|\xi {|}^2-\lambda )\hat{f}(\xi )=0$. Thus every eigenfunction is a linear combination of the indicated shell of frequencies, with no additional eigenvalues. Compactness gives the discrete self-adjoint realization, whose eigenfunctions are smooth. For real functions the two characters at $\xi ,-\xi$ become cosine and sine, giving the same real multiplicities. The opposite Laplacian sign would negate all displayed eigenvalues without changing the rigidity conclusion.

To prove two-dimensional spectral rigidity, let $L={\Lambda }^{\ast }$ and recover its Gram matrix from its vector-length multiset. Let $a$ be the shortest nonzero length and choose $v\in L$ with $|v|=a$. This vector is primitive: a proper integer multiple would have a shorter lattice vector. Hence $L\cap \mathbb{R}v=\mathbb{Z}v$. From the spectral multiset subtract exactly two occurrences of each length $ka$, $k=1,2,\dots$. The remaining multiset consists precisely of the vectors not on this line, including correct multiplicities at coincident lengths. Its least length $b$ is the length of a shortest vector $w$ independent of $v$.

The pair $(v,w)$ is a basis of $L$. Otherwise take a lattice point in a nonzero coset of $\mathbb{Z}v+\mathbb{Z}w$ and reduce its two coefficients into $[-1/2,1/2]$. The resulting nonzero lattice vector $z=\alpha v+\beta w$ is not on $\mathbb{R}v$, by primitivity of $v$. Butwith strict inequality: when both coefficients are nonzero the triangle inequality is strict for independent vectors, and when one vanishes its length is at most $b/2$. This contradicts the definition of $b$. Subtract an integer multiple of $v$ from $w$ and, if necessary, change its sign to arrange $0\le c=⟨v,w⟩\le a^2/2$; minimality of $b$ makes this reduction possible without decreasing its length.

The covolume $A$ of $L$ is also spectral. The count of lattice vectors with length at most $R$ obeysFor example, translate a bounded fundamental parallelogram of diameter bound $D$: the union of tiles with centres in $B_R$ contains $B_{R-D}$ and lies inside $B_{R+D}$. Comparing their areas proves the leading coefficient. The spectrum supplies this count. Since $(v,w)$ is a lattice basis, its parallelogram area satisfies $A^2=a^2{b}^2-c^2$. ThereforeAll three quantities are determined by the spectrum. Equal Gram matrices give an orthogonal map between the dual lattices; taking duals gives an orthogonal map between the original lattices, descending to a Riemannian isometry of the tori. Isospectral flat two-tori are isometric. This is the two-dimensional lattice reconstruction from vector lengths; it uses lengths with multiplicities, not merely the set of distinct lengths.

Take $N$ to be a closed compact Riemannian manifold and use the positive Laplace-Beltrami operator. Its Riemannian heat kernel $K_N(t,x,y)$ is the integral kernel of $e^{-t{\Delta }_N}$:It solves $({\partial }_t+{\Delta }_x)K_N=0$ for $t>0$ and tends to the identity kernel as $t\downarrow 0$. For an orthonormal eigenbasis, the spectral expansion of the Riemannian heat kernel isThe bar is required for a complex basis. With boundary, a specified invariant boundary condition must also be included in this definition.

For the quotient formula use the covering-space interpretation: $U$ acts freely and properly discontinuously by Riemannian isometries. On compact $N$ such a discrete group is finite. Freeness alone, without this covering hypothesis, does not justify the image sum: an infinite dense subgroup of circle rotations, for example, acts freely but has no manifold quotient. Positive-dimensional group actions require a different quotient analysis.

For lifts $x,y$ of $\bar{x},\bar{y}\in M$, the heat kernel on a finite isometric quotient isInvariance of $K_N$ under simultaneous Riemannian isometries and reindexing the sum show that this is independent of both lifts. A local isometry commutes with the Laplacian, so the sum satisfies the quotient heat equation. For its initial condition, integrate over a fundamental domain $F\subset N$ against a lifted function. The terms combine into the integral over all of $N$, whose initial limit is $f(x)$. Uniqueness of the heat evolution proves the formula. There is no factor $1/|U|$ in this kernel formula.

On the diagonal the sum is $U$-invariant. Its integral over $F$ is therefore $1/|U|$ times its integral over $N$, giving the heat traceThus the normalization factor appears in the trace, not the kernel. Put $F_t(g)={\int }_N{K}_N(t,x,gx)dV(x)$. For any finite isometry group $T$, changing variables $x=hy$ proves $F_t(hg{h}^{-1})=F_t(g)$: it is a class function on $T$.

For Gassmann equivalent subgroups $U_1,U_2\le T$, their intersections with each conjugacy class have equal size, and their orders are equal. If both act freely, the heat traces of their quotients satisfyso they agree for every $t>0$. Since $Z_M(t)={\sum }_{\lambda }m(\lambda )e^{-t\lambda }$, equality determines the spectrum with multiplicities: take $t\to \infty$ to recover the smallest eigenvalue and its multiplicity, subtract that term, and repeat. This proves Sunada theorem. Equivalently, projection onto $U_i$-invariant functions averages the group action, and the quotient eigenvalue multiplicity is $|U_i{|}^{-1}{\sum }_{u\in U_i}{\chi }_{\lambda }(u)$, with ${\chi }_{\lambda }$ the eigenspace character of a representation. Almost conjugacy equalizes these averages.

Use the positive Laplace-Beltrami operator $\Delta =-\mathrm{div}\mathrm{grad}$. A Riemannian submersion is a surjective smooth submersion $\pi :(M,g)\to (N,h)$ for which, at every $p$, the restriction of $d{\pi }_p$ to $H_p=(\ker d{\pi }_p{)}^{\perp }$ is a linear isometry onto $T_{\pi (p)}N$. The spaces $V_p=\ker d{\pi }_p$ and $H_p$ are its vertical and horizontal spaces. Its fibres are totally geodesic submanifolds precisely when ${\nabla }_{U}V$ is vertical for vertical vector fields $U,V$: their second fundamental form vanishes. Equivalently, a geodesic initially tangent to a fibre remains in that fibre while defined.

For a smooth $f:N\to \mathbb{R}$, its basic function $F=f\circ \pi$ is constant along each fibre. The Riemannian gradient of $F$ is the horizontal lift of a vector field through a submersion of ${\mathrm{grad}}_{N}f$, sincefor horizontal $X$, and $dF(U)=0$ for vertical $U$. In particular no derivative of $f$ in a vertical direction occurs.

Here is the needed connection fact, which also follows directly from the Koszul formula: for horizontal lifts $X,Y$ of vector fields $\bar{X},\bar{Y}$ on $N$, the horizontal component of ${\nabla }_{X}Y$ projects to ${\nabla }_{\bar{X}}^N\bar{Y}$. To see this, pair the Koszul formula with a third horizontal lift $Z$. The horizontal inner products are pulled back from $N$, and the horizontal components of their Lie brackets of vector fields project to the brackets on $N$. Thus all six terms are the pullbacks of the corresponding terms on $N$.

Choose an adapted Riemannian orthonormal frame $X_1,\dots ,X_n,U_1,\dots ,U_r$, with the $X_a$ horizontal lifts. Using $\mathrm{Hess}F(A,A)=A(AF)-({\nabla }_{A}A)F$, the connection fact givesFor a vertical $U_b$, both $U_{b}F=0$ and $({\nabla }_{U_b}{U}_b)F=0$, the latter because the fibres are totally geodesic submanifolds. Taking the negative metric trace of the Riemannian Hessian therefore proves the basic-function Laplacian identityThis identity is local and does not require compactness. Vanishing mean curvature of the fibres would already suffice; total geodesicity makes each vertical summand vanish separately.

For the discrete eigenspace assertion, assume the two Riemannian manifolds are closed manifolds. Without a discrete spectral realization, an unrestricted noncompact version need not have an eigenbasis. The projections of the Riemannian product $M\times N$ are Riemannian submersions with totally geodesic submanifolds as fibres. Its Levi-Civita connection splits into the two factor connections. Consequently its positive Laplace-Beltrami operator isThe cross term in the product rule for the positive Laplace-Beltrami operator is zero because the two factor Riemannian gradients are orthogonal.

We use the standard compact elliptic compact elliptic spectral theorem: the positive Laplace-Beltrami operator on a closed manifold is self-adjoint, has compact resolvent, and has a complete orthonormal eigenbasis of smooth eigenfunctions, with finite-dimensional eigenspaces and eigenvalues tending to infinity. Let ${\Delta }_M{u}_i={\mu }_i{u}_i$ and ${\Delta }_N{v}_j={\nu }_j{v}_j$. Fubini's theorem and completeness on each factor show that $u_i(x)v_j(y)$ form a complete orthonormal basis of $L^2(M\times N)$. For example, a function orthogonal to all these products has, for each $i$, zero $u_i$-coefficient as an $L^2(N)$ function, hence is zero.

The displayed operator identity makes $u_i{v}_j$ an eigenfunction with eigenvalue ${\mu }_i+{\nu }_j$. Conversely, if $\Delta w=\lambda w$, self-adjointness of the positive Laplace-Beltrami operator givesAll other coefficients vanish. Only finitely many pairs can have ${\mu }_i+{\nu }_j=\lambda$, since both spectra are nonnegative and have finitely many eigenvalues below any fixed bound. Thus the product Laplacian eigenspace decomposition isThe tensor product summands are mutually orthogonal; their elements are actual smooth eigenfunctions, so this is an equality of eigenspaces, not just a formal expansion.

Let $M$ be a closed manifold with a Riemannian metric and use $\Delta =-\mathrm{div}\mathrm{grad}$, the positive Laplace-Beltrami operator. This is the setting in which the final discrete eigenfunction expansion applies without additional boundary or noncompact spectral hypotheses. A Riemannian heat kernel is a smooth $H(t,x,y)$ for $t>0$ such thatfor every smooth initial function $f$. Here $\nu$ solves the heat equation. The limit is the statement that the initial kernel is the Dirac delta distribution on the diagonal. In the closed setting this uniquely determines the heat kernel; it is symmetric, preserves constants, is nonnegative, and satisfies the semigroup property.

A heat parametrix $P(t,x,y)$ is an approximate version with the same delta initial limit and with errorregular enough at $t=0$ to correct by a convergent integral series. One may construct it by the local Gaussian ansatzwhere $\chi =1$ near the diagonal and is supported in a convex normal neighborhood, and the smooth coefficients solve the usual radial transport equations. The leading coefficient accounts for the Riemannian volume form. Taking $q$ sufficiently large makes the residual extend continuously, with any prescribed finite number of derivatives, to $t=0$. Alternatively a smooth asymptotic summation of all transport coefficients gives a residual vanishing to every order. These local construction, extension and differentiability facts are used here as subsidiary results, as permitted.

For the correction proof take a heat parametrix whose residual is smooth up to $t=0$ on a short interval $[0,T]$. We also use its uniform integral bound ${\sup }_{0<t\le T,x}{\int }_M|P(t,x,y)|d{V}_y<\infty$, its approximate-identity limit, and the following standard differentiability property: convolution with a smooth time-dependent kernel gives a smooth kernel for $t>0$, and differentiating it yields the identity below. These properties follow from the Gaussian estimates and the delta initial limit; no conclusion about the final exact heat kernel is assumed.

Define the Volterra convolution of kernels byIt is associative wherever these integrals converge. The delta term at the upper endpoint givesSeek $H=P+P\ast Q$. Then its error vanishes exactly when $Q+R+R\ast Q=0$. The solution is the Volterra parametrix correctionThe signs matter: the first correction is $-P\ast R$.

To prove convergence, let $|R|\le C$ and $V=\mathrm{vol}(M)$. The time variables in $R^{\ast k}$ range over a simplex of volume $t^{k-1}/(k-1)!$, soThe series for $Q$ therefore converges uniformly on $[0,T]\times M\times M$. Derivatives obey analogous bounds, with polynomial factors in $k$, by the quoted residual extension and differentiation properties. Thus the series can be convolved and differentiated as above. Associativity and absolute convergence give $Q+R+R\ast Q=0$, proving $({\partial }_t+{\Delta }_x)H=0$. The correction $P\ast Q$ has integral norm $O(t)$ by the integral bound on $P$ and boundedness of $Q$. It has zero initial limit, so $H$ has precisely the required delta initial data.

For uniqueness, a smooth solution $w$ with zero initial data satisfies the heat equation energy identityIts initial norm is zero, hence $w=0$. Applying uniqueness to consecutive evolutions proves the semigroup property. Choose $s_0<T$ and, for arbitrary $t>0$, compose enough kernels $H(t/r)$ that $t/r<s_0$. The resulting kernel is independent of the subdivision by uniqueness, is smooth for positive time by the quoted differentiability property, and extends the construction to every $t>0$. This proves a parametrix determines the global heat kernel in the closed setting. The heat equation maximum principle gives nonnegativity and applying uniqueness to the constant initial function gives conservation of total mass. Self-adjointness of the positive Laplace-Beltrami operator $\Delta$ gives symmetry.

Finally let $\{{\varphi }_j\}$ be a complete complex orthonormal eigenbasis of $\Delta$, with eigenvalues ${\lambda }_j\ge 0$, repeated by multiplicity. As subsidiary analytic facts we use the compact elliptic compact elliptic spectral theorem, elliptic regularity bounds making each fixed derivative of ${\varphi }_j$ grow at most polynomially in $1+{\lambda }_j$, and a polynomial eigenvalue-counting bound. The spectral expansion of the Riemannian heat kernel converges because exponential decay makes the following series converge with every derivative when $t\ge \epsilon >0$:Indeed, evolving initial data ${\varphi }_j$ gives $e^{-t{\lambda }_j}{\varphi }_j$ by direct substitution into the heat equation and uniqueness. For general $f\in L^2(M)$, completeness and the heat equation energy identity giveThe smoothly convergent kernel series represents this same operator, and equality for all smooth $f$ identifies it pointwise with the constructed heat kernel. The complex conjugate is required by the inner product; for a real eigenbasis it may be omitted. A noncompact Riemannian manifold may instead require a spectral integral, and the unqualified discrete formula should not be asserted there.

Use the given orientation throughout, and the usual convention that the manifold has no boundary. In an oriented manifold chart the metric volume form isIf $y$ is another oriented chart and $A=\partial x/\partial y$, then $g_y=A^t{g}_{x}A$, so $\sqrt{\det g_y}=(\det A)\sqrt{\det g_x}$ because $\det A>0$. This is exactly the transformation of the top exterior product. The formulas therefore agree on overlaps and define a smooth positive volume form. Equivalently it takes value one on any positively oriented orthonormal tangent frame.

The metric induces an inner product on $p$-forms by making the increasing exterior products of an orthonormal coframe orthonormal. The Hodge star operator is the unique pointwise linear map $\ast :{\Lambda }^p{T}^{\ast }M\to {\Lambda }^{n-p}{T}^{\ast }M$ satisfyingFor an increasing multi-index $I$, $\ast e^I$ is the complementary wedge with the sign making $e^I\wedge \ast e^I={\omega }_g$. Swapping the blocks of $p$ and $n-p$ factors introduces $(-1{)}^{p(n-p)}$. ThereforeIn particular $\ast 1={\omega }_g$ and $\ast {\omega }_g=1$.

For compactly supported smooth $f$, the $(n-1)$-form $f\ast \alpha$ has compact support. Stokes theorem and the graded Leibniz rule giveA top form $\tau$ equals $(\ast \tau ){\omega }_g$. Consequently this Hodge integration by parts for one-forms becomes exactlyCompact support of $f$ is enough; $\alpha$ need not itself have compact support.

Define the codifferential on $p$-forms by $\delta =(-1{)}^{n(p+1)+1}\ast d\ast$ and the Hodge Laplacian by $\Delta =d\delta +\delta d$. A harmonic differential form is a smooth form in $\ker \Delta$. On functions this is the positive Laplace-Beltrami operator, $\Delta f=\delta df=-{\mathrm{div}}_g{\mathrm{grad}}_{g}f$. This sign convention is required by the product identity; it is the negative of the $\mathrm{div}\mathrm{grad}$ convention also commonly used for the Laplace-Beltrami operator.

For a $p$-form $\beta$, the square of the Hodge star operator and the codifferential formula giveApplying the second formula to $d\beta$ and the first to $\delta \beta$ givesThus Hodge star commutes with the Hodge Laplacian. Since $\ast$ is invertible, $\beta$ is harmonic if and only if $\ast \beta$ is harmonic. This holds without compactness; we have not used the generally false noncompact implication that a harmonic form must be closed and coclosed.

For a smooth function $h$ and one-form $\alpha$, the same graded Leibniz rule gives $\delta (h\alpha )=h\delta \alpha -\ast (dh\wedge \ast \alpha )=h\delta \alpha -⟨dh,\alpha {⟩}_g$. Apply this to $d(f_1{f}_2)=f_{2}d{f}_1+f_{1}d{f}_2$ to obtain the product rule for the positive Laplace-Beltrami operator

The Hodge decomposition theorem states that on a compact oriented boundaryless Riemannian manifold, smooth forms have the $L^2$-orthogonal decompositionand each de Rham cohomology class has a unique harmonic differential form representative. To spell out the last conclusion, a harmonic form is closed and coclosed because $⟨\Delta \eta ,\eta ⟩=\Vert d\eta {\Vert }^2+\Vert \delta \eta {\Vert }^2$. If a closed form decomposes as $h+du+\delta v$, then $d\delta v=0$; integration by parts gives $\Vert \delta v{\Vert }^2=⟨v,d\delta v⟩=0$. Thus it represents $h$. A harmonic exact form has zero norm by adjointness, proving uniqueness. The analytic existence of the decomposition is the stated Hodge decomposition theorem.

Now let $M$ be compact, connected and oriented. A harmonic function satisfies $0=⟨f,\Delta f⟩=\Vert df{\Vert }^2$, so it is constant. The Hodge star operator identifies ${\mathcal{H}}^0(M)$ with ${\mathcal{H}}^n(M)$, hence ${\mathcal{H}}^n(M)=\mathbb{R}{\omega }_g$. A Riemannian metric exists by Question 3 even if none was initially chosen. Using the harmonic representative of each class, we obtain the top de Rham cohomology of a compact connected oriented manifoldThe class is nonzero also directly from Stokes theorem, since ${\int }_M{\omega }_g>0$ while every exact top form has zero integral. Boundarylessness matters: a compact interval has zero first de Rham cohomology.

For closed manifolds with the product Riemannian metric, the positive Laplace-Beltrami operator is the sum of the two factor operators. Products of their orthonormal eigenbases form a complete basis, and eigenvalues add. Hence the product eigenspace at $\lambda$ is the displayed finite orthogonal sum of tensor products. The discrete assertion requires a discrete spectral setting.

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.

The symmetric covariant two-tensor defined using the Levi-Civita connection by ${\mathrm{Hess}}_{g}f(X,Y)=X(Yf)-({\nabla }_{X}Y)f$. Its negative metric trace is the positive Laplace-Beltrami operator applied to $f$. Along a geodesic it gives the second derivative of $f$.

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.

For the positive Laplace-Beltrami operator $P$ on a closed manifold, the wave trace is the distribution ${\sum }_j{e}^{it\sqrt{{\lambda }_j}}$, counting multiplicities. The series is interpreted after testing against smooth compactly supported functions, rather than as an ordinarily convergent sum. Its singularities relate the spectrum to the length spectrum of closed geodesics, through the Duistermaat-Guillemin trace formula.

For a Laplace-type operator on a closed manifold, let $a_{t^0}$ be the constant coefficient of its full heat trace expansion. In the Mellin transform defining the spectral zeta function, the reduced trace has constant coefficient $a_{t^0}-\dim \ker L$. Its contribution is $(a_{t^0}-\dim \ker L)/s$, and $1/\Gamma (s)=s+O(s^2)$ proves the identity. For the scalar positive Laplace-Beltrami operator on a closed surface, $a_{t^0}=\int RdV/(24\pi )$ and $\dim \ker L$ is the number of connected components. In particular, a connected sphere has value $-2/3$ and a connected flat torus has value $-1$.