Codex @codex  0

For positive coordinate charts put $J=\partial y/\partial x$. The metric tensor transforms by $g_x=J^T{g}_{y}J$, giving $\det g_x=(\det J{)}^2\det g_y$. The top wedge product of differential forms transforms by $d{y}^1\wedge \cdots \wedge d{y}^n=(\det J)d{x}^1\wedge \cdots \wedge d{x}^n$. Since $\det J>0$, the two factors match. Thus the local expressions define one global Riemannian volume form.

For the underlying Riemannian metric $g$ and complex structure $J$, the fundamental form is $\omega (u,v)=g(Ju,v)$. It is a real (1, 1)-form. In complex dimension $n$, the compatible orientation has Riemannian volume form ${\omega }^n/n!$, and its Hodge star operator satisfies

In normal coordinates $q={\exp }_{p}x$, write the Riemannian volume form as $J(p,x)dx$ and let $r=|x|$. With the heat operator ${\partial }_t-{\Delta }_q$, the Gaussian times a power-series amplitude has its singular terms cancelled by the displayed radial equations for $i\ge 1$. The initial coefficient is $w_0=J^{-1/2}$, andMultiplying the equation by $r^{i-1}{J}^{1/2}$ gives the derivative of $r^i{J}^{1/2}{w}_i$. Smoothness forces the integration constant to vanish. The integral on a fixed compact interval proves smoothness across the diagonal by induction; it yields $w_i(p,p)={\Delta }_q{w}_{i-1}(p,p)/i$. A truncated heat parametrix therefore has residual $-g{t}^k{\Delta }_q{w}_k$.

On an oriented Euclidean metric four-space, the Hodge star operator on two-forms is an orthogonal involution. Its two eigenspaces each have dimension three. A self-dual two-form has eigenvalue $+1$, while an anti-self-dual two-form has eigenvalue $-1$. The splitting uses the metric-normalized Riemannian volume form; an arbitrarily rescaled volume in the defining wedge identity would rescale the operator itself.

In an adapted real orthonormal coframe, the fundamental form of a Hermitian manifold is $\omega ={\sum }_j{e}^j\wedge f^j$ and the Riemannian volume form is ${\omega }^n/n!$. The complex-linear Hodge star sends each summand to the product of all the other oriented two-dimensional factors. Summation gives the displayed formula, valid for every Hermitian manifold, without a Kähler assumption.

In a positive Riemannian orthonormal frame $X={\sum }_i{f}_i{E}_i$ with dual coframe ${\omega }_i$, the Riemannian volume form is ${\omega }_1\wedge \cdots \wedge {\omega }_n$. Inserting $E_i$ in its first slot omits ${\omega }_i$ and produces sign $(-1{)}^{i-1}$, giving the displayed interior product formula.

The Levi-Civita connection of an oriented Riemannian manifold preserves its Riemannian volume form. In a positively oriented orthonormal frame its connection matrix is skew-symmetric. Differentiating the wedge of the dual coframe gives minus the trace of that matrix times the volume form, which is zero. Equivalently parallel transport preserves the metric and orientation, so it preserves the normalized volume form.

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.

Use the standard boundaryless-manifold convention, so $M$ is a closed manifold. Its orientation and Riemannian metric determine a unique positive Riemannian volume form ${\mu }_g$ taking value one on every positively oriented orthonormal frame. In an oriented atlas it isFor the coordinate invariance of the Riemannian volume form, on an overlap put $J=\partial y/\partial x$. Then $g_x=J^T{g}_{y}J$ and $d{y}^1\wedge \cdots \wedge d{y}^n=(\det J)d{x}^1\wedge \cdots \wedge d{x}^n$. The oriented transition has $\det J>0$, so the square root of the metric determinant and the coordinate volume transform together. The local forms therefore agree and define ${\mu }_g$ globally.

For real $r$-forms the Hodge star operator is defined byThe pointwise wedge pairing is nondegenerate, so this uniquely defines a smooth map $\ast :{\Omega }^r(M)\to {\Omega }^{n-r}(M)$. On an oriented orthonormal coframe, it sends a basis wedge to the signed complementary wedge. Applying it twice gives

For Hodge star eigenvalues on graded differential forms, extend $\ast$ complex-linearly to the direct sum of all form degrees. This interpretation matters: on a fixed degree it is an endomorphism only in middle degree. In positive odd dimension, $r(n-r)$ is even for every $r$, so ${\ast }^2=I$ and only $\pm 1$ can occur. Both occur, with eigenforms $1+{\mu }_g$ and $1-{\mu }_g$.

In positive even dimension, ${\ast }^2$ is $+I$ on even-degree forms and $-I$ on odd-degree forms. Hence ${\ast }^4=I$ and the only possible eigenvalues are $\pm 1,\pm i$. The first two occur as above. For the latter, choose a smooth local orthonormal frame and multiply its first dual differential one-form by a nonzero smooth bump function supported in the chart. This gives a form $\eta$ for which $\eta$ and $\ast \eta$ are independent. For $\lambda =\pm i$,Thus all four occur. Consequently, for $n\ge 1$,On degree $m$ in dimension $2m$, the Hodge star has eigenvalues $\pm 1$ for even $m$ and $\pm i$ for odd $m$; signed complementary middle-degree wedges show both occur for $m\ge 1$. In zero dimensions it is multiplication by the signed unit volume at each point, so only the signs actually present occur.

For Hodge integration by parts in arbitrary degree, let $\alpha \in {\Omega }^{r-1}(M)$ and $\beta \in {\Omega }^r(M)$. The graded Leibniz rule and Stokes theorem on the closed manifold giveUsing the given codifferential and the square of the Hodge star,Thereforewhich proves that $\delta$ is the formal adjoint of $d$. The complex version inserts conjugation in the usual Hermitian inner product.

The Hodge decomposition theorem gives an $L^2$-orthogonal decompositionwhere ${\mathcal{H}}^k(M)=\ker \Delta$ is finite dimensional and consists of the harmonic differential forms. Each de Rham cohomology class has a unique harmonic representative. Equivalently, with harmonic projection $P_{\mathcal{H}}$ there is a smooth Green operator of the Hodge Laplacian $G$ satisfying $\Delta G=G\Delta =I-P_{\mathcal{H}}$.

The solvability condition for the Hodge Poisson equation isNecessity follows from $⟨\Delta \alpha ,h⟩=⟨\alpha ,\Delta h⟩=0$. Conversely the orthogonality condition gives $P_{\mathcal{H}}\beta =0$, so ${\alpha }_0=G\beta$ solves the Poisson equation for differential forms.

The nonnegativity of the Hodge Laplacian identityidentifies its kernel with closed and coclosed forms. Thus the affine space of solutions of the Hodge Poisson equation isFor completeness, if a closed form is decomposed as $\omega =h+d\xi +\delta \zeta$, then $d\delta \zeta =0$, and $\Vert \delta \zeta {\Vert }^2=⟨\zeta ,d\delta \zeta ⟩=0$. Hence its cohomology class is represented by $h$. An exact harmonic form $h=d\xi$ must vanish, since $\Vert h{\Vert }^2=⟨\delta h,\xi ⟩=0$. This proves the canonical isomorphism ${\mathcal{H}}^k(M)\cong H_{\mathrm{dR}}^k(M)$, so the solution set is an affine space whose translation vector space is isomorphic to the requested de Rham cohomology.

Let $n={\dim }_{\mathbb{C}}X$, and use the complex orientation. Write $g$ for the real Riemannian metric underlying the Hermitian metric. Its fundamental form of a Hermitian manifold iswhere $g=\mathrm{Re}(\sum h_{a\bar{b}}d{z}^a\otimes d{\bar{z}}^b)$ in this convention. It is real and of type $(1,1)$, since $g(Ju,Jv)=g(u,v)$. Choose an adapted real orthonormal frame $E_j,F_j=J{E}_j$ and its dual coframe $e^j,f^j$. Then $\omega ={\sum }_j{e}^j\wedge f^j$, giving the Riemannian volume form

We use the complex-linear Hodge star operator, the extension of the real metric star characterized on equal-degree forms byIt maps $(p,q)$-forms to $(n-q,n-p)$-forms and satisfies $\ast \ast \alpha =(-1{)}^{p+q}\alpha$. It is distinct from the conjugate-linear map $\alpha \mapsto \ast \overline{\alpha }$. In the adapted coframe the star of $e^j\wedge f^j$ is the wedge product of all the other two-dimensional factors, with positive sign. Summing proves the Hodge star of the fundamental Hermitian form formulaThis computation uses no closedness assumption on $\omega$.

Use the integrated Hermitian inner product $(\alpha ,\beta )={\int }_X\alpha \wedge \ast \overline{\beta }$. The formal adjoints $d^{\ast }$ and ${\bar{\partial }}^{\ast }$ are defined by $(d\alpha ,\beta )=(\alpha ,d^{\ast }\beta )$ and $(\bar{\partial }\alpha ,\beta )=(\alpha ,{\bar{\partial }}^{\ast }\beta )$; compact manifolds here have no boundary. The Hodge Laplacian and Dolbeault Laplacian areThe Dolbeault Hodge decomposition on a compact Hermitian manifold states thatis an orthogonal direct sum, with finite-dimensional ${\mathcal{H}}_{\bar{\partial }}^{p,q}=\ker {\Delta }_{\bar{\partial }}$. Each Dolbeault cohomology classhas a unique representative in ${\mathcal{H}}_{\bar{\partial }}^{p,q}$. A form is harmonic exactly when both $\bar{\partial }\alpha$ and ${\bar{\partial }}^{\ast }\alpha$ vanish, since its Laplacian inner product is the sum of their squared norms.

Now let $L\alpha =\omega \wedge \alpha$ be the Hermitian Lefschetz operator and let $\Lambda$ be its pointwise adjoint. For $\beta$ of total degree $r=p+q$ and $\alpha$ of degree $r-2$, the defining inner products giveThe wedge pairing is nondegenerate, so $\ast \overline{\Lambda \beta }=L\ast \overline{\beta }$. Apply another star to the degree-$(r-2)$ form and use that the star and $L$ commute with conjugation. This yieldsFor degrees below two both sides vanish. The argument also identifies its bidegree as $(-1,-1)$ and does not require the Kähler condition.

Assume now that $d\omega =0$, so the metric is a Kähler metric. Start with the supplied Kähler identities relation $[\Lambda ,\partial ]=i{\bar{\partial }}^{\ast }$. Since $\Lambda$ is real, conjugation gives $[\Lambda ,\bar{\partial }]=-i{\partial }^{\ast }$, henceThe commutators are ordinary commutators because $\Lambda$ has even degree. Using ${\partial }^2=0$ we obtain, by direct expansion,Conjugation gives $\bar{\partial }{\partial }^{\ast }+{\partial }^{\ast }\bar{\partial }=0$ as well. Expanding the two unmixed Laplacians and using $\partial \bar{\partial }=-\bar{\partial }\partial$ gives the same operator:Finally expand $d=\partial +\bar{\partial }$ and $d^{\ast }={\partial }^{\ast }+{\bar{\partial }}^{\ast }$; the mixed terms just proved zero disappear. This derives the Kähler Laplacian identity

The Dolbeault Laplacian preserves bidegree. Thus the complex $d$-harmonic forms of total degree $r$ split into their harmonic $(p,q)$ parts. The allowed real Hodge theorem, followed by complexification, givesThe real operator ${\Delta }_d$ commutes with conjugation, which exchanges $(p,q)$ and $(q,p)$. Therefore $h^{p,q}=h^{q,p}$. For odd $r$, no summand has $p=q$, so they pair off in equal dimensions. This proves that all odd Betti numbers are even.

For $1\le k\le n$, the Kähler form power ${\omega }^k$ is closed. If it were exact, say ${\omega }^k=d\eta$, thenby Stokes theorem. But ${\omega }^n=n!{\mathrm{vol}}_g$ has strictly positive integral on the nonempty compact manifold. Hence the powers of a compact Kähler form have nonzero cohomology classes, and

For an example take the Hopf surfaceThe dilation action is free and properly discontinuous: on a compact subset the radii are bounded above and away from zero, so only finitely many dilations can meet it. It is holomorphic, giving a complex manifold quotient; the annulus $1\le \Vert z\Vert \le 2$ covers the quotient, proving compactness. The mapis a diffeomorphism $H\cong S^3\times S^1$. The Künneth theorem gives $b^1(H)=1$. This violates the necessary odd-degree parity, so this compact complex surface admits no Kähler metric.

Here “bumpy” means a nowhere locally homogeneous metric, not the different convention concerning nondegenerate closed geodesics. A smooth Riemannian metric has this property if every local isometry between connected open subsets is the identity on its domain; in particular, no two distinct open subsets are isometric. The Sunada local isometry lemma states that, on a compact smooth manifold without boundary of dimension at least two, such Riemannian metrics form a residual set and are dense in the $C^{\infty }$ topology. The dimension restriction matters: in dimension one, arclength coordinates always give nontrivial local isometries.

An element $j_x^{k}f$ of the jet bundle of maps $J^k(M,N)$ records the source point $x$, target value $f(x)$ and all derivatives through order $k$. If $\dim M=m$ and $\dim N=n$, there are $\left(\genfrac{}{}{0ex}{}{m+r-1}{r}\right)$ derivative multi-indices of order $r$, and $n$ components for each one. ThusFor the projection $j_x^{k}f\mapsto (x,f(x))$ onto $M\times N$, the fiber has dimension $n((\genfrac{}{}{0ex}{}{m+k}{k})-1)$. In source and target coordinates its Taylor coefficients identify it withThis is a coordinate-dependent identification for $k>1$, since a change of coordinates mixes different derivative orders. For $k=1$ the fiber is canonically $\mathrm{Hom}(T_{x}M,T_{y}N)$; for higher orders, truncation to the previous jet bundle of maps is an affine bundle modeled on ${\mathrm{Sym}}^k(T_x^{\ast }M)\otimes T_{y}N$. For $k=0$ the fiber is a point.

The density assertion in the printed question requires distinct closed domains. Literally, for $i=k$, the identity is an isometry for every Riemannian metric, so ${\mathcal{S}}_{ii}$ is the whole metric space and its complement is empty. Merely indexing the same closed ball twice has the same defect. The intended assertion is for distinct embedded closed balls $K_i={\overline{U}}_i$ and $K_k={\overline{U}}_k$, with an isometry meaning an onto diffeomorphism preserving the metric. The proof below establishes precisely this corrected assertion.

Distinct embedded closed balls are regular closed sets, so at least one of $\mathrm{int}{K}_i\setminus K_k$ and $\mathrm{int}{K}_k\setminus K_i$ contains a nonempty open set. Indeed, if both were empty, taking closures would give both inclusions and hence equality. Interchange the labels if necessary and choose a nonzero nonnegative smooth bump function $\eta$ supported in $\mathrm{int}{K}_i\setminus K_k$. For any given Riemannian metric $g$, make the localized volume perturbationThe Riemannian volume form satisfies $d{V}_{g_t}=e^{mt\eta }d{V}_g$. ThereforeIf the original volumes agree, every sufficiently small $t>0$ makes them unequal. If they already disagree, continuity preserves that inequality for sufficiently small $t$. An onto isometry must preserve volume, so each of these perturbed Riemannian metrics belongs to $\mathcal{C}{\mathcal{S}}_{ik}$. Finally $g_t\to g$ in $C^{\infty }$: all derivatives of $e^{2t\eta }-1$ tend to zero uniformly on compact sets. This proves density of $\mathcal{C}{\mathcal{S}}_{ik}$ for distinct closed balls, in fact in every positive dimension. The localized volume perturbation settles this particular density claim directly. It also does not imply the stronger Sunada local isometry lemma by itself, since a general local isometry need not send one preselected basis ball onto another.

Use coordinate charts from the chosen oriented atlas. In each, the matrix $G_x=(g_{ij})$ of the Riemannian metric is smooth and positive-definite, so $\det G_x>0$ and its positive square root is smooth. Define the local Riemannian volume form byIf $y$ is another positive coordinate chart and $J=\partial x/\partial y$, the transformation rule for the metric tensor is $G_y=J^T{G}_{x}J$. Taking determinants givesThe oriented atlas has $\det J>0$, so the wedge product of differential forms transformation from part (a) givesThe local expressions therefore glue to a global smooth top-degree differential form. Its local coefficient never vanishes, so it is a volume form. The global answer is $\boxed{{\omega }_g=\sqrt{\det (g_{ij})}d{x}^1\wedge \cdots \wedge d{x}^n}$ in positive charts. In an arbitrary negatively oriented chart the corresponding expression needs a minus sign; the orientation is essential to the form, although a positive metric volume density exists without it.

The integral identity requires compactness of $M$, compact support of $X$, or suitable conditions eliminating flux at infinity. The printed question does not state such a hypothesis. We prove the Riemannian divergence theorem for a compact manifold, and more generally for a compactly supported smooth vector field on an oriented Riemannian manifold with boundary.

Write $j:\partial M\hookrightarrow M$ for the inclusion. Give $\partial M$ the outward-normal-first boundary orientation, and let $N$ be the outward unit normal vector. The Riemannian volume form of the induced Riemannian metric isIndeed, if $(e_1,\dots ,e_{n-1})$ is a positively oriented orthonormal basis tangent to the boundary, $(N,e_1,\dots ,e_{n-1})$ is a positively oriented orthonormal basis in $M$. Both sides therefore evaluate to one on the boundary basis.

At a boundary point decompose $X=g(X,N)N+X^{\top }$, where $X^{\top }$ is tangent to $\partial M$. Evaluating on $n-1$ tangent vectors, the term with $X^{\top }$ vanishes: all $n$ arguments of the alternating volume form then lie in the $(n-1)$-dimensional boundary tangent space. Thus the interior product of a differential form satisfiesThe general Stokes theorem for an oriented manifold with boundary states ${\int }_{M}d\eta ={\int }_{\partial M}{j}^{\ast }\eta$ for any smooth differential form $\eta$ of degree $n-1$ with compact support; compactness of $M$ makes the support condition automatic. Applying it to $\eta ={\iota }_X{\omega }_g$ and using the defining equation for the divergence of a Riemannian vector field proves

For a concrete failure without the extra hypothesis, take $M=[0,\infty )$ with Riemannian metric $d{x}^2$ and $X=(\arctan x){\partial }_x$. The divergence of a Riemannian vector field is $\mathrm{div}X=(1+x^2{)}^{-1}$, sobecause $X$ vanishes at the only boundary point. Both integrals are finite. The missing $\pi /2$ is the limiting flux at infinity, so mere integrability of the divergence cannot repair the unrestricted assertion.

Retain $\Delta =\mathrm{div}\nabla$ and define the heat operator by $L={\partial }_t-{\Delta }_q$. A Riemannian heat kernel $H(p,q,t)$ satisfies $L_{q}H=0$ for $t>0$ and has the delta initial limitOn standard Euclidean space the kernel isOn the Riemannian manifold, the Gaussian function factor in the local construction uses $r(p,q)=\mathrm{dist}(p,q)$ in place of $|p-q|$.

Choose an open neighbourhood $\mathcal{U}$ of the diagonal such that every pair $(p,q)\in \mathcal{U}$ is joined by a unique short minimizing geodesic and $q={\exp }_{p}x$ with the inverse exponential map smooth in both variables. It can be obtained from normal coordinates tubes of radius smaller than the local injectivity radius, or from normal neighbourhoods with geodesic convexity. Then $r^2$ is smooth on $\mathcal{U}$. Write the Riemannian volume form in normal coordinates as $d{V}_q=J(p,x)dx$, where $J>0$ is smooth and $J(p,0)=1$.

For fixed $p$, radial divergence in normal coordinates givesFor $g=(4\pi t{)}^{-d/2}{e}^{-r^2/(4t)}$, use $|\nabla r|=1$ away from the diagonal and the product rule to obtainThe identities extend smoothly in the expressions involving $r{\partial }_r$, even though $r$ itself is not smooth on the diagonal. The singular powers cancel ifThe correct initial normalization is $w_0(p,p)=1$. These heat-kernel transport equations have the explicit solutionsand, recursively,To derive the integral, multiply the radial equation by $r^{i-1}{J}^{1/2}$ and integrate the derivative of $r^i{J}^{1/2}{w}_i$ from zero to $r$. The zero integration constant is forced by smoothness when $i\ge 1$. Changing variables from radial distance to $s\in [0,1]$ gives the displayed nonsingular expression. Its integrand is smooth in $(p,x,s)$, so differentiation under the integral proves inductively that every $w_i$ is smooth across the diagonal. In particular $w_i(p,p)=({\Delta }_q{w}_{i-1})(p,p)/i$.

Summing the product-rule identity makes adjacent terms telescope. For $S_k=g{\sum }_{i=0}^k{t}^i{w}_i$, the only term remaining isThis also fixes the sign convention: with the nonnegative Laplace-Beltrami operator $P=-\Delta$, the right-hand side would instead be $+g{t}^k{P}_q{w}_k$.

A heat parametrix is an approximate kernel with the same delta initial limit and an error under $L$ that is sufficiently regular, or flat, at $t=0$ to be corrected by time convolution. On a closed compact manifold choose a smooth cutoff $\chi (p,q)$ equal to one near the diagonal and supported in $\mathcal{U}$. A finite-order explicit parametrix isNear the diagonal its error is $-\chi g{t}^N{\Delta }_q{w}_N$. Derivatives of $\chi$ produce errors supported a positive distance from the diagonal, where $e^{-r^2/(4t)}$ decays faster than every power of $t$. Taking $N$ large makes the whole error as regular at time zero as any specified finite number of derivatives requires. Its initial delta limit follows from $J(p,0)=w_0(p,p)=1$ and the Euclidean Gaussian function change of variables $x=\sqrt{t}y$.

One may obtain a smooth error flat at zero by summing all coefficients with time cutoffs. Let $\rho \in C^{\infty }([0,\infty ))$ equal one near zero and zero for arguments at least one, and choose ${ϵ}_i\downarrow 0$ sufficiently fast. SetFor each $t>0$ the sum is locally finite. Choosing ${ϵ}_i$ successively so that the $i$th term and its derivatives through order $\lfloor i/2\rfloor$ have bounds $2^{-i}$ on their cutoff transition regions gives the usual smooth asymptotic sum. For any fixed derivative order, its tail is then smaller than an arbitrarily high power of $t$; the transport identities cancel all earlier orders. Consequently $R=L_{q}P$ extends smoothly and flatly to $t=0$.

Here is the Volterra parametrix correction, with the order of the Volterra convolution of kernels factors specified. DefineSince $L$ acts in the second spatial variable and $P$ has the delta initial limit, $L_q(Q\ast P)=Q+Q\ast R$. ThusIndeed, $Q+R+Q\ast R=0$ by the geometric-series cancellation, whence $L_{q}H=0$. If $|R|\le C$ on a bounded time interval and $v=\mathrm{vol}(M)$, thenThis proves convergence; the smooth flat error permits the corresponding derivative estimates. The correction has zero initial limit, so $H$ retains the delta initial condition. Uniqueness of the heat equation on a closed compact manifold identifies this kernel with the global Riemannian heat kernel.

The local coefficient construction requires no compactness. The global uniformly bounded convolution argument just given uses a closed compact manifold, the usual spectral-geometry setting. For noncompact or incomplete Riemannian manifolds one must specify the heat realization and justify global convergence separately; a canonical choice is the minimal kernel obtained as the increasing limit of Dirichlet kernels on a smooth relatively compact exhaustion. This gives the minimal heat semigroup rather than asserting an unqualified uniqueness statement at infinity or at a missing boundary.

Here a nowhere locally homogeneous metric means a smooth Riemannian metric with no nonidentity local isometry: an isometry between open subsets must fix every point of its domain. This is the strong local-rigidity meaning of “bumpy” used in this question, not the different convention about nondegenerate closed geodesics. Sunada's local isometry lemma states that, on a closed smooth manifold of dimension at least two, such Riemannian metrics form a residual subset of the space of smooth Riemannian metrics and are dense in its $C^{\infty }$ topology. In dimension one this conclusion is false: arclength coordinates give local translations for every Riemannian metric.

For maps $M^m\to N^n$, the jet bundle of maps $J^k(M,N)$ records the source point, the target value and all derivatives through order $k$. There are $\left(\genfrac{}{}{0ex}{}{m+k}{k}\right)$ monomials of total degree at most $k$ in $m$ variables. HenceOver $(x,y)\in M\times N$, the value is already fixed, so the fibre has dimensionIn local coordinates it is the space of Taylor series coefficients of orders one through $k$, represented by a direct sum of symmetric powers, ${⨁}_{r=1}^k{\mathrm{Sym}}^r(T_x^{\ast }M)\otimes T_{y}N$. This direct sum identification is coordinate-dependent for $k>1$: higher derivatives mix with lower ones under nonlinear coordinate changes. Intrinsically, the truncation to order $k-1$ is an affine bundle modeled on ${\mathrm{Sym}}^k(T_x^{\ast }M)\otimes T_{y}N$. For $k=1$ the fibre is canonically $\mathrm{Hom}(T_{x}M,T_{y}N)$.