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.

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.

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.

Work with a finite-dimensional complex semisimple Lie algebra. The Dynkin diagram is obtained from its intrinsic root system, with several structure theorems entering the construction.

First choose a Cartan subalgebra $\mathfrak{h}$. The semisimple structure theory gives the root-space decompositionThe nonzero roots span ${\mathfrak{h}}^{\ast }$, root spaces have dimension one, and every root has its negative. The rank is $r=\dim \mathfrak{h}$.

Next use the Killing form. Its nondegeneracy on $\mathfrak{h}$ and its positive-definite restriction to the real coroot span identify the real root span $E=\mathbb{R}\Phi$ with a Euclidean space. Explicitly, if $t_{\lambda }={\mathcal{B}}^{-1}(\lambda )$, put $(\lambda ,\mu )=B(t_{\lambda },t_{\mu })$. The sl2 subalgebra associated with a root supplies its normalized coroot $H_{\alpha }$ and the integral pairingThe root-string theorem makes the reflection $s_{\alpha }(\beta )=\beta -\beta (H_{\alpha })\alpha$ permute $\Phi$. Together with the absence of multiples other than $\pm \alpha$, this shows that $\Phi$ is a finite reduced root system and a crystallographic root system. These properties come from semisimple Lie algebra theory, not from the diagram's definition.

Choose a vector in $E$ that is not perpendicular to any root. Define the positive roots by positive inner product with this vector. The simple roots are those positive roots that cannot be expressed as sums of two positive roots. Root-system theory shows that they form a basis $\Delta =\{{\alpha }_1,\dots ,{\alpha }_r\}$, and each root is an integer combination of them with coefficients all of one sign. Distinct simple roots have nonpositive inner product.

Use the standard Cartan matrix conventionThere is one diagram vertex for each simple root. For $i\ne j$, the integral off-diagonal entries are nonpositive andConnect the two vertices with respectively zero, one, two or three edges. The corresponding angles are $90^{\circ },120^{\circ },135^{\circ },150^{\circ }$. A multiple edge carries an arrow toward the shorter root. For joined vertices, the length ratio is determined by $A_{ij}/A_{ji}=({\alpha }_i,{\alpha }_i)/({\alpha }_j,{\alpha }_j)$; hence the diagram records both angle and relative length.

For example, ${\mathfrak{sl}}_3$ has simple roots $L_1-L_2$ and $L_2-L_3$, giving two equal-length vertices joined by one edge: type $A_2$. For ${\mathfrak{sp}}_4$, take ${\alpha }_1=L_1-L_2$ and ${\alpha }_2=2L_2$. Thenso the two vertices have a double edge directed toward ${\alpha }_1$, the shorter root: type $C_2$.

Finally, orthogonal irreducible components of the root system correspond to the simple ideals of a Lie algebra in the semisimple Lie algebra, so its diagram is the disjoint union of their connected diagrams. Conjugacy of Cartan subalgebras and the Weyl group action on choices of positive roots make the resulting diagram independent of these choices up to isomorphism. Thus the vertices, edge multiplicities, arrows and connected components encode the simple-root geometry and the simple-ideal decomposition.

Work over $\mathbb{C}$ with finite-dimensional Lie algebras and Lie algebra representation spaces. The adjoint form of Engel theorem isHere a nilpotent Lie algebra has terminating lower central series, ${\gamma }_1\mathfrak{g}=\mathfrak{g}$, ${\gamma }_{j+1}\mathfrak{g}=[\mathfrak{g},{\gamma }_j\mathfrak{g}]$. We prove the stronger linear statement: if $V\ne 0$, $L\subseteq \mathfrak{gl}(V)$ and every member of $L$ is a nilpotent endomorphism, then there is a nonzero $v\in V$ killed by all of $L$, and $L$ can be simultaneously represented by strictly upper triangular matrices.

For the proof of Engel theorem by induction and normalizers, first establish a useful nilpotence fact. If $x^m=0$, left and right multiplication by $x$ on $\mathrm{End}(V)$ commute. Expanding their difference givesFor $N=2m-1$, every summand contains a power of $x$ at least $m$, so $(\mathrm{ad}x{)}^{2m-1}=0$. Restrictions and induced quotient maps remain nilpotent.

Now prove the common-kernel statement, the Engel lemma, by induction on $\dim L$. The zero algebra is immediate. Choose a maximal proper Lie subalgebra $M$. Its action on $L/M$ by commutators consists of nilpotent maps by the preceding fact. Applying the lower-dimensional induction hypothesis to its image gives a nonzero coset $y+M$ with $[M,y]\subseteq M$. HenceThe normalizer of a Lie subalgebra is a Lie subalgebra by the Jacobi identity. Maximality implies $N_L(M)=L$, so $M$ is an ideal of a Lie algebra. Moreover $L/M$ has no proper nonzero Lie subalgebra, and every one-dimensional subspace is a Lie subalgebra. Thus $\dim L/M=1$.

Induction also provides a nonzero common annihilatorIt is stable under $L$: for $m\in M$, $x\in L$, $v\in W$,because $[m,x]\in M$. Choose $x$ representing a basis of $L/M$. The restriction of the nilpotent endomorphism $x$ to the nonzero space $W$ has nonzero kernel. A vector in this kernel is killed by both $M$ and $x$, hence by all of $L$. This finishes the induction. When an action is not faithful, induction is applied to its image, whose dimension is no larger than $\dim M$; no faithfulness assumption is hidden.

Starting with the common-kernel line, apply the same statement to each successive quotient Lie algebra representation space. It constructs a complete flag of invariant subspacesAn adapted basis represents every member of $L$ by a strictly upper triangular matrix. Every product of $d$ such matrices is zero, so every iterated commutator with $d$ entries is zero as well.

Finally, if every $\mathrm{ad}x$ on $\mathfrak{g}$ is nilpotent, apply the linear statement to the Adjoint representation. With $d=\dim \mathfrak{g}$, every product of $d$ adjoint operators vanishes, so ${\gamma }_{d+1}\mathfrak{g}=0$. Conversely, if ${\gamma }_{c+1}\mathfrak{g}=0$, then $(\mathrm{ad}x{)}^{c}y=0$ for all $x,y$, proving the other implication. This proves Engel's theorem and its simultaneous-triangularization form. The hypothesis concerns every element of the algebra, not merely a chosen set of generators.

A smooth manifold alone has no distinguished notion of an isometry; one must first choose a Riemannian metric $g$. We therefore interpret the essay as concerning a smooth connected Riemannian manifold $(M,g)$. Its isometry group consists of diffeomorphisms $f$ with $f^{\ast }g=g$. In coordinates this isAn isometry preserves lengths, distances, the Levi-Civita connection, geodesics and all curvature tensors. The Myers-Steenrod theorem ensures that the group is a finite-dimensional Lie group acting smoothly, and that even distance-preserving bijections have the required smoothness.

The crucial rigidity is determination by one first jet. An isometry fixing $p$ satisfieswherever the exponential chart is defined. If also $d{f}_p=I$, it is the identity on a neighborhood. Continuing this equality through overlapping exponential neighborhoods proves it on connected $M$. Therefore the isotropy group at $p$ embeds in $O(T_{p}M,g_p)$. Its dimension is at most $n(n-1)/2$, while an orbit has dimension at most $n$. ConsequentlyOn a compact manifold the group is compact. Indeed distance preservation makes a family of isometries equicontinuous; compactness and the same argument for their inverses give subsequential limits which are inverse distance-preserving bijections. Smoothness and the Lie-group structure then follow from Myers-Steenrod theorem.

Infinitesimal isometries are Killing vector fields. Differentiating $f_t^{\ast }g=g$ at zero givesOn a complete Riemannian manifold these fields have complete flows and constitute the Lie algebra of the isometry group; on an incomplete manifold one must retain only fields whose flows are globally defined. A Killing vector field is fixed by its value and its skew derivative at one point, consistent with the same dimension bound. Along a geodesic, its restriction is a Jacobi field, so local infinitesimal data propagate rigidly.

In dimension three there is no independent Weyl curvature tensor: the Riemann curvature tensor is determined by the Ricci tensor. Large isotropy therefore strongly constrains the geometry. If the full dimension six is attained, the group is transitive with three-dimensional isotropy. That isotropy rotates all tangent two-planes transitively, forcing equal sectional curvatures at each point; transitivity makes the value constant. The simply connected complete models are the round $S^3$, Euclidean ${\mathbb{R}}^3$ and hyperbolic ${\mathbb{H}}^3$. Their groups are respectively $O(4)$, ${\mathbb{R}}^3⋊O(3)$ and $O^{+}(3,1)$, all of dimension six. The orientation-preserving hyperbolic group is $PS{L}_2(\mathbb{C})$.

The standard Thurston geometries show how symmetry drops when a preferred direction is present:For the product models, the surface factor has a three-dimensional isometry group, and translations along the line add one parameter. For Nil and the universal covering of $S{L}_2(\mathbb{R})$ with their standard geometric metrics, a three-dimensional transitive group is enlarged by a circle of isotropy. The Sol metricis preserved by $(u,v,w)\mapsto (u_0+e^{-t}u,v_0+e^{t}v,w+t)$, giving three continuous parameters; its continuous isotropy is trivial because the distinguished expanding and contracting directions cannot be rotated into each other. Additional discrete symmetries may still occur. These dimensions refer to the model metrics, not arbitrary left-invariant perturbations.

A quotient generally has fewer global isometries than its local model. If $X$ is the complete simply connected universal cover and $M=X/\Gamma$, the isometry group of a Riemannian quotient isTo prove this, lift an isometry to the universal cover. Conjugating a deck transformation by the lift again gives a deck transformation, so the lift normalizes $\Gamma$. Conversely a normalizing isometry descends, and two lifts differ by a deck transformation. A flat three-torus retains its three-dimensional translation group, with additional discrete lattice symmetries. A round spherical quotient retains only the normalizer of its finite deck group; the full six-dimensional group does not automatically descend.

For a closed hyperbolic three-manifold, there are no nonzero Killing vector fields. The integrated Killing identity isBoth nonnegative terms must vanish. The compact isometry group is therefore zero-dimensional and hence finite. The Mostow rigidity theorem additionally says that isomorphisms between the fundamental groups of complete finite-volume hyperbolic three-manifolds are induced by unique hyperbolic isometries; this fixes the complete hyperbolic metric up to isometry and contrasts with the moduli of flat tori. An arbitrary smooth three-manifold need not admit a single homogeneous metric: prime decomposition and torus splitting may lead to several geometric pieces. Thus local homogeneous symmetry, global quotient symmetry and the chosen metric must be distinguished when discussing isometries in dimension three.