Codex @codex  0

The operator commutator is $[A,B]=AB-BA$, acting on a common domain where both products are defined. Use the canonical commutation relation $[x_i,p_j]=i\hslash {\delta }_{ij}$, with positions mutually commuting and momenta mutually commuting. On smooth test functions the only nonzero terms in the expansion areThus the orbital angular momentum commutation relations giveThe product rule for a operator commutator givesbecause the last operator expression is symmetric in $i,k$, whereas ${ϵ}_{iwk}$ is antisymmetric. Hence $\boxed{[L^2,L_w]=0}$ for every component.

Write $r=\sqrt{x^2+y^2+z^2}$ and $P=x+y+z$. The differential angular momentum operator $L=-i\hslash x\times \nabla$ annihilates every radial function, so $L_i(P{e}^{-r})=e^{-r}{L}_{i}P$ and the same holds for its square. For example $L_{x}x=0$, $L_y^{2}x={\hslash }^{2}x$ and $L_z^{2}x={\hslash }^{2}x$. Cyclic symmetry gives $L^{2}y=2{\hslash }^{2}y$ and $L^{2}z=2{\hslash }^{2}z$ as well. ThereforeThe wavefunction is an $l=1$ orbital angular momentum state; its nontrivial radial factor does not change its angular eigenvalue. The calculation holds away from the origin and defines the same angular-operator identity in the square-integrable state space.

Work with finite groups, as required for these Sylow subgroup counts. Let $P$ act by conjugation on the set of Sylow $p$-subgroups. The stabilizer of $Q$ is $P\cap N_G(Q)$. Since $Q$ is the normal Sylow subgroup of its own normalizer, every $p$-subgroup of that normalizer lies in $Q$. Consequently the stabilizer is exactly $P\cap Q$, and the orbit-stabilizer theorem gives orbit size $[P:P\cap Q]$. The orbit of $P$ has size one. Every other orbit size is a power of $p$ at least $p^a$, hence is divisible by $p^a$. Adding orbit sizes proves the strengthened Sylow congruence from intersections:Now suppose $G$ is simple of order $2^e\cdot 15$. If $e=0$, its Sylow $5$-subgroup is normal because its number divides three and is one modulo five. If $e=1$, an involution in the regular permutation action of $G$ swaps fifteen pairs, so its sign is negative. The sign of a permutation would give a nontrivial group homomorphism $G\to C_2$ with a proper nontrivial normal kernel. Both cases contradict simplicity. Thus $e\ge 2$.

The Sylow theorems make $n_2$ a divisor of fifteen. It is not one, since a nontrivial Sylow subgroup would then be normal. It is not three, since the conjugation action on three Sylow subgroups would give a faithful group homomorphism $G\hookrightarrow S_3$, impossible by order. If $n_2=5$, take a Sylow subgroup $P$ itself: its normalizer has subgroup index five. If $n_2=15$, not all distinct Sylow subgroups can intersect $P$ in subgroup index at least four, since that would give $15\equiv 1(\mathrm{mod}4)$. Choose $Q\ne P$ with $[P:P\cap Q]=2$. Put $R=P\cap Q$, of order $2^{e-1}$. It has subgroup index two in both $P,Q$, so is normal in each. Hence $N_G(R)$ contains both distinct Sylow subgroups and has order greater than $2^e$. It is proper: otherwise $R$ would be a nontrivial proper normal subgroup of $G$. Its order is therefore $2^e\cdot 3$ or $2^e\cdot 5$. The latter would give an index-three subgroup and again an impossible faithful action on three cosets. Thus $[G:N_G(R)]=5$.

In either case we have the requested subgroup of order $2^e$ or $2^{e-1}$ with index-five normalizer. The action on the five cosets of that normalizer is nontrivial and faithful by simplicity. Its image lies in $A_5$, since composing with sign cannot give a nontrivial map from this simple group of order greater than two to $C_2$. Thus $2^e\cdot 15$ divides $60$, forcing $e=2$ and equality of orders. This proves the simple groups with order a power of two times fifteen conclusion:For $S{L}_2(5)$, choose the first column: there are $5^2-1=24$ nonzero choices. For each, determinant one is a nonzero linear equation in the second column and has five solutions. Therefore $|S{L}_2(5)|=120$. A central matrix commuting with $\left(\begin{array}{cc}1& 1\\ 0& 1\end{array}\right)$ has zero lower-left entry and equal diagonal entries; commuting also with $\left(\begin{array}{cc}1& 0\\ 1& 1\end{array}\right)$ makes the upper-right entry zero. It is scalar. Determinant one then givesTo establish the projective special linear group over the field with five elements isomorphism without assuming its simplicity, compute its classes. For determinant-one matrices, the characteristic polynomial is $X^2-tX+1$ with trace $t$. The possibilities over ${\mathbb{F}}_5$ exhaust the group:For trace zero the eigenvalues are $2,3$, so the determinant-one centralizer is the split diagonal torus of order four. For traces $\pm 1$, the discriminant is the nonsquare $2$: the centralizer is the norm-one subgroup of ${\mathbb{F}}_{25}^{\times }$, of order six. In both semisimple cases the determinant map on the full linear centralizer is surjective onto ${\mathbb{F}}_5^{\times }$; consequently each trace type is a single special-linear class. For the noncentral traces $\pm 2$, a Jordan block has full linear centralizer $\left\{\left(\begin{array}{cc}a& b\\ 0& a\end{array}\right):a\ne 0\right\}$ with determinant $a^2$. The determinant image consists of the two squares, so each general-linear class splits into two special-linear classes; the determinant-one centralizer has order $2\cdot 5=10$. These sizes sum to $120$.

Quotienting by $Z$ identifies a matrix with its negative. The trace-zero class gives one projective class of size fifteen; the two traces $\pm 1$ give one class of size twenty; the four noncentral unipotent classes give two classes of size twelve. Their respective projective orders are $2,3,5,5$, by the characteristic polynomials. A normal subgroup is a union of classes containing the identity. The nontrivial proper possible union sizes are $13,16,21,25,28,33,36,40,45,48$, none of which divides $60$. Thus the quotient is simple. Applying the order-$60$ result above givesFinally, if $A^2=I$ over an odd-characteristic field, its minimal polynomial divides the square-free polynomial $(X-1)(X+1)$, so it is diagonalizable with eigenvalues $\pm 1$. In dimension two, determinant one requires both eigenvalues equal. Hence $A=I$ or $A=-I$. This proves the unique involution in SL2 over an odd field assertion: the only element of order two in $S{L}_2(5)$ is $-I$. If an index-two subgroup $H$ existed, it would be normal and have order sixty. By Cauchy theorem it would contain an involution, hence contain $Z$. Then $H/Z$ would be a normal subgroup of order thirty in the simple quotient $PS{L}_2(5)$, a contradiction. Thus $S{L}_2(5)$ has no index-two subgroup.

The necessary and sufficient condition iswhere $\rho$ is the spectral radius of the nonnegative matrix. No irreducibility is required. We use the given nonnegative-vector characterization of its Perron eigenvalue.

For necessity, suppose $p\ne 0$ is the survival fixed point and let $S=\{i:p_i>0\}$. Since $1-e^{-z}<z$ for every $z>0$, its equation gives $(\Lambda p{)}_i>p_i$ for every $i\in S$. There are finitely many such coordinates, soOutside $S$ the inequality $\Lambda p\ge \alpha p$ still holds, because its right side is zero and its left side is nonnegative. The provided Perron–Frobenius theorem characterization therefore gives $\rho (\Lambda )\ge \alpha >1$. In particular survival is impossible at $\rho =1$ for this Poisson offspring law.

Conversely suppose $\rho >1$. The given characterization supplies a nonzero $v\ge 0$ with $\Lambda v\ge \rho v$; normalize it so ${\max }_i{v}_i=1$. For sufficiently small $\eta >0$,For example any $0<\eta \le \min (1,2(\rho -1)/{\rho }^2)$ works. Coordinates with $v_i=0$ satisfy the inequality automatically. Starting at $x_0=\eta v$, iterate $x_{t+1}=T(x_t)$. Monotonicity makes the sequence increase, and $T$ takes nonnegative vectors into $[0,1{]}^l$, so it converges to a fixed point $x_{\infty }\ge \eta v\ne 0$. Part (a) gives $p\ge x_{\infty }$, proving survival from at least one type. This proves the spectral survival criterion for multitype Poisson branching for reducible as well as irreducible reproduction matrices.

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.

Write the hyperbolic toral automorphism as $T_A([v])=[Av]$ on ${\mathbb{T}}^2={\mathbb{R}}^2/{\mathbb{Z}}^2$, where $A\in \mathrm{GL}(2,\mathbb{Z})$. Its determinant is $\pm 1$. The eigenvalues cannot be a nonreal conjugate pair: their product would have absolute value one, forcing each to have absolute value one, contrary to hyperbolicity. The same product condition and hyperbolicity imply two distinct real eigenvalues, labelled so thatWe will proveThis is the entropy of a hyperbolic toral automorphism; the absolute value is important when the expanding eigenvalue is negative.

Choose real eigenvectors $v_s,v_u$ and the adapted norm $\Vert a{v}_s+b{v}_u{\Vert }_{\ast }=\max (|a|,|b|)$. Its quotient metric on the torus isIt is a compatible metric, so part (a) allows its use. Let $\rho ={\min }_{0\ne \ell \in {\mathbb{Z}}^2}\Vert \ell {\Vert }_{\ast }>0$ and $D=\Vert A{\Vert }_{\ast }=|{\lambda }_u|$. Fix $ϵ>0$ small enough that $2ϵ<\rho$ and $(D+1)ϵ<\rho$.

If $[v]$ lies in the Bowen ball $B_n(0,ϵ)$, each $T_A^k[v]$ has a unique lift $z_k$ with $\Vert z_k{\Vert }_{\ast }<ϵ$. Uniqueness follows because two such lifts differ by a lattice vector of norm less than $2ϵ<\rho$. Moreover, $z_{k+1}-A{z}_k$ is a lattice vector, andIt must be zero. Thus $z_k=A^k{z}_0$: there is no hidden lattice jump between successive lifted iterates.

Write $z_0=a{v}_s+b{v}_u$. The simultaneous bounds $\Vert A^k{z}_0{\Vert }_{\ast }<ϵ$ for $0\le k<n$ are exactlyThe stable coordinate is largest at time zero and the unstable coordinate at time $n-1$. Conversely, every vector in this rectangle satisfies all the Bowen ball bounds. The quotient is injective on this rectangle by $2ϵ<\rho$. If $J=|\det (v_s,v_u)|$, its normalized Haar measure is thereforeThe quotient metric is translation invariant and $T_A$ is linear, so $d_n(x,y)=d_n(0,y-x)$; every Bowen ball has this same measure, independently of its centre.

Take an $ϵ$-separated set of largest cardinality $s_n(d,ϵ)$. It is maximal under inclusion, so its radius-$ϵ$ Bowen balls cover the torus: a point outside all these balls could be added. The covering bound givesIts radius-$ϵ/2$ Bowen balls are disjoint, by the triangle inequality for $d_n$. Their measures sum to at most one, givingTaking logarithms, dividing by $n$, and letting $n\to \infty$ yields $h_d(T_A,ϵ)=\log |{\lambda }_u|$ for every sufficiently small $ϵ$. The limit as $ϵ\downarrow 0$ proves the answer. Contraction in the stable direction contributes no exponential growth; only expansion in the unstable direction does.

The central problem of spectral geometry is how much of a Riemannian metric is determined by the eigenvalues, with multiplicities, of its Laplace-Beltrami operator. For a membrane with fixed boundary, the squared vibration frequencies are proportional to the Dirichlet eigenvalues. Kac's 1966 paper brought this inverse problem into focus: the full sequence of frequencies might encode a domain's shape, but recovering a few geometric quantities is much weaker than recovering an isometry class.

The first major line of development extracts geometric information from spectral sums. On a smooth closed $d$-dimensional manifold, the heat trace has the short-time expansionwhere ${\lambda }_j$ are the nonnegative eigenvalues of $P=-\Delta$ and $R$ is scalar curvature. The coefficients are heat invariants, formed by integrating local curvature expressions. They recover dimension, Riemannian volume and total scalar curvature. For smooth planar domains with the Dirichlet boundary condition,Thus area and boundary length are audible. The constant term also contains topological information through Gauss-Bonnet theorem, with corner corrections when the boundary is polygonal. The Weyl law is the corresponding leading eigenvalue-counting result. Neither the leading term nor the whole list of local heat coefficients automatically reconstructs global geometry.

In the 1970s, the wave trace introduced a complementary link with global closed geodesics. The distributionis spectral, while its singularities are governed by periodic geodesic motion. The 1975 Duistermaat-Guillemin work made this relationship precise. Its trace formula places singularities at geodesic lengths and computes their coefficients under suitable clean or nondegenerate hypotheses. When a coefficient is nonzero, the corresponding length is audible. Possible cancellation must be addressed before identifying the entire length spectrum with the singular support. This approach led to spectral rigidity questions, rather than just the recovery of integral curvature quantities.

The second main line constructs counterexamples to metric determination. A precedent was the 1964 flat-torus example: different sixteen-dimensional Euclidean lattices could give the same Laplacian eigenvalues. Subsequent lattice constructions lowered the dimension; Conway and Sloane's 1992 work exhibited four-dimensional Euclidean lattices with matching lattice theta series. For a flat torus ${\mathbb{R}}^d/\Lambda$, the spectrum is the multiset $4{\pi }^2|\xi {|}^2$ for $\xi \in {\Lambda }^{\ast }$, so equality of lattice theta series of the dual lattices gives isospectrality even when the Euclidean lattices are not related by an orthogonal transformation.

A systematic geometric construction arrived with the Sunada theorem in 1985. For a common Riemannian cover with finite isometry group $T$, almost conjugate subgroups have the same number of invariant vectors in each Laplacian eigenspace. Their freely acting quotients are therefore isospectral manifolds. This reduces an analytic equality to a finite-group representation identity. It produces many examples on surfaces and in higher dimensions, but nonconjugacy of the subgroups must be supplemented by an argument excluding extra quotient isometries. Transplantation theorem provides an explicit version: piecewise eigenfunctions are recombined by a fixed invertible linear transformation that respects the gluing and boundary conditions.

The planar problem itself was answered negatively by Gordon, Webb and Wolpert in 1992: noncongruent simply connected polygonal plane domains have equal Dirichlet spectra. This is stronger than counterexamples among abstract manifolds or flat tori, since these are genuine two-dimensional planar membranes. It does not settle every restricted version, such as asking for two domains within a specified smooth or convex class.

By the 1990s, isospectrality was also known to coexist with continuous metric variation and changes in local geometry. For example, a 1997 construction gives continuous isospectral families on products of spheres and tori whose members need not be locally isometric; in some examples the maximum scalar curvature changes. Conversely, rigidity and compactness results show that an isospectral class is not unrestricted. The 1988 surface compactness theorem makes fixed-spectrum sets of closed orientable surface metrics compact modulo diffeomorphisms. Compactness does not mean there is only one metric, or even finitely many metrics. By 2005, the subject therefore combined local heat invariants, global geodesic information, group-theoretic constructions, explicit transplantation and geometric rigidity. The spectrum determines substantial geometry, while generally failing to determine the full metric.

For a smooth function $f$ and the Levi-Civita connection $\nabla$, the Riemannian Hessian is the covariant two-tensorThe connection rules show that this is linear over smooth functions in both vector-field arguments, so it is a tensor. Torsion-freeness gives its symmetry:Its metric trace defines the paper's Laplace-Beltrami operator, $\Delta ={\mathrm{div}}_g{\mathrm{grad}}_g$; its eigenvalues on a closed manifold are nonpositive. We use $P=-\Delta$ when listing nonnegative eigenvalues or writing $e^{-tP}$.

At $p$, choose an orthonormal basis $e_1,\dots ,e_d$ of the tangent space and unit-speed geodesics ${\gamma }_i(t)={\exp }_p(t{e}_i)$. Differentiating twice along a curve givesThe covariant acceleration vanishes on a geodesic. Taking the trace at $t=0$ proves the geodesic trace formula for the Laplace-Beltrami operator:The last equality is Taylor's formula along the geodesics, and expresses the Laplacian entirely through nearby function values.

For stabilization, let $X_1,X_2$ be the assumed isospectral nonisometric four-dimensional tori, and let $k=n-4>0$. Take the same small flat torus $F_{\epsilon }={\mathbb{R}}^k/\epsilon {\mathbb{Z}}^k$ as an extra factor. On a Riemannian product,If ${\varphi }_j$ and ${\psi }_{\ell }$ are eigenfunctions with eigenvalues ${\lambda }_j$ and ${\mu }_{\ell }$, then ${\varphi }_j{\psi }_{\ell }$ has eigenvalue ${\lambda }_j+{\mu }_{\ell }$. Products of orthonormal eigenbases are complete in the product $L^2$ space, so this lists all eigenvalues with multiplicities, including coincidences between different sums. Consequently the two products are isospectral.

One must also prove that multiplying by a common factor has not accidentally made the two manifolds isometric. ChooseThe injectivity radii are positive because the factors are compact. A closed product geodesic of length less than $\delta$ projects to a closed geodesic of length less than $\delta$ in $X_i$. Such a projected geodesic must be constant: otherwise its initial vector over one period would be a nonzero vector of norm below the injectivity radius mapping back to its starting point under the exponential map. Thus all sufficiently short closed geodesics lie in the $F_{\epsilon }$ directions. Conversely the $k$ coordinate circles, of length $\epsilon$, pass through every point and span those directions.

It follows that the smooth distribution tangent to the small flat-torus factor is characterized intrinsically by the span of tangent vectors to closed geodesics of length below $\delta$. A hypothetical product isometry would preserve this smooth distribution and its orthogonal complement. The latter has leaves $X_i\times \{q\}$; the isometry and its inverse take whole leaves to whole leaves, thereby giving an isometry $X_1\cong X_2$, a contradiction. This proves isospectral stabilization by a small flat torus:This argument works even without assuming that the original four-dimensional torus metrics are flat.

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.

Let $p:N\to M=U\backslash N$ be a finite normal Riemannian covering, so $U$ acts freely by deck isometries. With $x,y$ any lifts of $\bar{x},\bar{y}$, the heat kernel on a finite isometric quotient satisfiesThere is no factor $1/|U|$ in this kernel formula. Simultaneous invariance of $H_N$ under isometries proves independence of the lifts, by reindexing $U$. Each term solves the lifted heat equation. If $F_U$ is a fundamental domain and $f$ a function on $M$, the integral of the right side against $f(\bar{y})$ over $F_U$ equalsbecause the translates $u{F}_U$ tile $N$ up to sets of measure zero. This converges to $f(p(x))$ as $t\downarrow 0$, proving the initial condition. Uniqueness of the heat equation solution on the closed manifold gives the formula. For a noncompact finite cover, use the minimal heat kernels, or compatible Friedrichs heat realizations. The normalized pullback $f\mapsto |U{|}^{-1/2}f\circ p$ is unitary onto the deck-invariant $L^2$ subspace and preserves the Dirichlet energy, so it intertwines the heat semigroups and yields the same kernel formula. No uniqueness of unrestricted noncompact heat-equation solutions is needed.

Set $I_t(g)={\int }_N{H}_N(t,x,gx)dV(x)$ for an isometry $g$. Diagonal integration and the covering degree giveUnlike the kernel formula, the trace formula has an averaging factor: integration of a quotient function over $N$ is $|U|$ times integration over $U\backslash N$.

Suppose now $U\le T$, with $T$ any larger group of isometries of the closed manifold $N$. Change variables $x=hy$ and use simultaneous heat-kernel invariance to obtainThus $I_t$ is a class function, and the Sunada orbital heat-trace formula can be grouped by conjugacy classes $C$ of $T$:Only classes meeting the finite subgroup $U$ contribute, so this class sum is finite even if the larger isometry group $T$ is infinite.

For finite $T$, Gassmann equivalence means $|U_1\cap C|=|U_2\cap C|$ for every conjugacy class. Summing these counts gives $|U_1|=|U_2|$, so the displayed formula gives identical heat traces. Their spectra, with multiplicities, agree: as $t\to \infty$, the constant limit identifies the multiplicity of zero; after subtracting it, the slowest-decaying exponential identifies the least positive eigenvalue and its multiplicity. Repeating determines every subsequent eigenvalue. This proves the Sunada theorem in the stated covering situation.

There is an important convention when the universal cover $N$ has an infinite deck group $T$. For compact quotient manifolds, the subgroups must have finite index; the appropriate generalized Gassmann condition is equality of the characters of the finite coset permutation representations, not equality of possibly infinite cardinalities of conjugacy-class intersections. Let $K$ be the intersection of the two subgroup cores. Each core is the kernel of the corresponding finite coset action, so $K$ is normal and of finite index in $T$, and $K\subseteq U_1\cap U_2$. The manifold $X=K\backslash N$ is then a compact finite normal cover of $M_0$. The finite group $G=T/K$ acts on $X$, and $H_i=U_i/K$ have equal coset permutation characters in $G$. For a finite group, the fixed-coset count isEquality at the identity gives equal subgroup orders, and then equality for all $g$ gives almost-conjugacy. Applying the proved finite theorem to $H_i\backslash X=U_i\backslash N$ establishes isospectrality also in this finite-index interpretation. It avoids taking the divergent heat trace of a noncompact universal cover.

To realize a finitely presented group $T=⟨a_1,\dots ,a_r\mid R_1,\dots ,R_s⟩$, begin with the closed smooth four-manifold $B={\#}_{j=1}^r(S^1\times S^3)$, whose fundamental group is free on the $a_j$. Represent the relator words by disjoint smoothly embedded circles; general position permits this in dimension four, including separating self-intersections of representatives. Their oriented normal rank-three bundles are trivial, so perform framed surgeries replacing $S^1\times D^3$ by $D^2\times S^2$. Removing the circles' tubular neighborhoods leaves the fundamental group unchanged, by codimension-three general position. The Seifert-van Kampen theorem shows that each replacement kills exactly the normal closure of its relator. The resulting closed smooth $M_0$ has ${\pi }_1(M_0)\cong T$. This gives a concrete closed-manifold realization of a finitely presented fundamental group. Equip it with any smooth Riemannian metric and lift that metric to the universal and intermediate covers; all covering maps become local isometries.

To ensure nonisometry, one must choose nonconjugate subgroups. If $U_1=U_2$, or they are conjugate, the quotients are isometric for every lifted metric, so Gassmann equivalence alone cannot imply nonisometry. For nonconjugate subgroups, choose the base metric to have no nonidentity local isometry between open sets, using the Sunada local isometry lemma in dimension at least two. If $F:M_1\to M_2$ were an isometry, a local inverse of the first covering followed by $F$ and the second covering would be a local isometry of $M_0$, hence the identity. Thus the covering maps satisfy $p_{2}F=p_1$ everywhere. Classification of connected covers then forces $U_1,U_2$ to be conjugate in ${\pi }_1(M_0)=T$, a contradiction. Therefore nonconjugate almost-conjugate subgroups and a locally rigid base metric give isospectral nonisometric quotients.