Codex @codex  0

For a compact smooth Riemannian surface with smooth boundary and a pure Dirichlet boundary condition or pure Neumann boundary condition, the constant coefficient in the heat trace isThe equality is the Gauss-Bonnet theorem. Thus such isospectral manifolds have the same Euler characteristic. A simply connected compact connected bordered surface is a disk with $\chi =1$, whereas a nonorientable surface with $c\ge 1$ crosscaps and $b\ge 1$ boundary components has $\chi =2-c-b\le 0$. These cannot be isospectral in this setting. Corners and mixed boundary conditions contribute additional heat coefficients, so the smooth pure-boundary-condition hypothesis is essential.

Two isospectral manifolds remain isospectral after taking a Riemannian product with the same closed factor: eigenvalues add, with multiplicities given by the product Laplacian eigenspace decomposition. For nonisometric closed connected factors $M_1,M_2$, choose $\epsilon <\delta <{\min }_i\mathrm{inj}(M_i)$. Every closed product geodesic of length less than $\delta$ has constant projection to $M_i$, since a nonconstant closed geodesic shorter than the injectivity radius is impossible. The short coordinate geodesics of the small flat torus span exactly its tangent factor at every point. Any product Riemannian isometry must preserve that intrinsically defined smooth distribution and its orthogonal complement, hence restrict to an isometry of the original factors. Thus these products are still nonisometric. Taking $M_i$ to be four-dimensional tori proves the higher-dimensional stabilization statement without assuming a cancellation theorem for arbitrary products.

Embed $C_4\times C_4$ and $Q_8\times C_2$ in $S_{16}$ by their regular representations. Both have one element of order one, three of order two, and twelve of order four. An order-$r$ element in a regular action has cycle type $r^{16/r}$. Their intersections with every conjugacy class of $S_{16}$ consequently have the same size, so they are Gassmann equivalent. The first group is abelian and the second is not, so their underlying groups are not isomorphic. Finite-cover constructions from the Sunada theorem therefore give isospectral manifolds with different fundamental groups.

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 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.

Use the following skew-symmetric matrix, whose columns reduce modulo three to four of the specified nonzero tetracode words:Set $B_{+}=3I+S$, $B_{-}=-3I+S$, and define the full Euclidean lattices $L_{\pm }^0=B_{\pm }{\mathbb{Z}}^4$. Their determinants have absolute value $144$, sinceFor arbitrary positive real numbers $a_1,\dots ,a_4$, put $D=\mathrm{diag}(\sqrt{a_1},\dots ,\sqrt{a_4})$ and $L_{\pm }=D{L}_{\pm }^0$. This is the tetracode length-isospectral lattice construction. We prove equality of the length multisets by an explicit bijection, not just equality of the shortest lengths.

Modulo three, both lattices have image $\mathrm{im}S=T$. The matrix $S$ has rank two over ${\mathbb{F}}_3$; its four column lines are precisely the four lines in $T$. Also $S^2=0$ modulo three, so $\ker S=\mathrm{im}S=T$. Each nonzero tetracode word has exactly one zero coordinate. Every $v\in L_{\pm }^0$ therefore has some coordinate divisible by three. Moreover, all four coordinates of $v$ have the same parity, because every entry of either $B_{\pm }$ is odd.

For $v\in L_{+}^0$, choose the least index $i$ with $3\mid v_i$ and set $w=v-2v_i{e}_i$. The coordinate reflection bijection for tetracode lattices will show that $w\in L_{-}^0$. First, $w\equiv v(\mathrm{mod}3)$, so $Sw\equiv 0(\mathrm{mod}3)$. Second, positive-lattice membership implies $Sv\equiv 3v\equiv -v(\mathrm{mod}4)$. HenceThe expression in parentheses has $i$th coordinate zero and every other coordinate $v_j+v_i{S}_{ji}$ even, by the common parity and $S_{ji}=\pm 1$. Thus $(I-S)w\equiv 0(\mathrm{mod}4)$. Together the congruences modulo three and four say $(-3I-S)w\equiv 0(\mathrm{mod}12)$, which is precisely $B_{-}^{-1}w\in {\mathbb{Z}}^4$.

Conversely, if $v\in L_{-}^0$, the same reflection has unchanged reduction modulo three, and now $Sv\equiv v(\mathrm{mod}4)$. ThusIts $i$th coordinate is $4v_i$ and its other coordinates are even multiples of two, so it vanishes modulo four; it vanishes modulo three as well. The reflected vector is in $L_{+}^0$. The chosen index is unchanged by reflection, including for vectors reducing to zero, so this map is an involution between the two lattices. It preserves every $|v_j|$ and therefore preserves $\Vert Dv{\Vert }^2={\sum }_j{a}_j{v}_j^2$. The two lattices have the same length spectrum, including every multiplicity, for every positive choice of weights.

To prove a genuine failure of isometry, chooseThese four numbers are linearly independent over $\mathbb{Q}$: changing the signs of the independent square roots in the multiquadratic field $\mathbb{Q}(\sqrt{2},\sqrt{3},\sqrt{5})$ isolates each coefficient of a rational relation. If an orthogonal transformation $O$ maps $L_{+}$ to $L_{-}$, thenand $P^T\mathrm{diag}(a_j)P=\mathrm{diag}(a_j)$. Comparing rational coefficients of the independent $a_j$ in the diagonal entries gives $P_{ji}^2={\delta }_{ji}$. Consequently $P$ is diagonal with entries $\pm 1$.

Such a coordinate sign change must preserve $T$ modulo three. For each zero coordinate $i$, there are exactly two opposite nonzero tetracode words. Preserving their line requires the other three signs to be equal. Taking two different zero coordinates forces all four signs to be equal. Thus $P=I$ or $P=-I$, either of which would imply $L_{+}^0=L_{-}^0$. ButThis contradiction proves $L_{+}$ and $L_{-}$ are not isometric. Allowing affine isometries changes nothing: translating the image of zero back to zero reduces them to the orthogonal case.

For the flat tori, take period lattices equal to the dual lattices:The standard spectrum of a flat torus lemma states that its characters of a real torus are indexed by the dual of its period lattice. Here they are $e^{2\pi i⟨v,x⟩}$ for $v\in L_{\pm }$, with Laplacian eigenvalue $4{\pi }^2\Vert v{\Vert }^2$. These characters give a complete orthogonal basis, so our length-preserving bijection proves that $X_{+}$ and $X_{-}$ are isospectral manifolds. The second standard lemma is that an isometry of flat tori lifts to an affine isometry of their Euclidean universal covers. Its orthogonal part maps the period lattices to one another, and therefore maps their dual lattices to one another. This would give the impossible isometry $L_{+}\cong L_{-}$. Hence $X_{+}$ and $X_{-}$ are isospectral non-isometric flat four-dimensional tori.

The musical isomorphisms are the pointwise maps induced by the Riemannian metric:In local coordinates they are $(X^{♭}{)}_i=g_{ij}{X}^j$ and $({\alpha }^{♯}{)}^i=g^{ij}{\alpha }_j$. The inverse metric makes these inverse bundle maps, hence inverse maps between smooth vector fields and smooth differential one-forms.

For a covariant two-tensor $\beta$, its metric trace is ${\mathrm{tr}}_g\beta =g^{ij}{\beta }_{ij}={\sum }_i\beta (e_i,e_i)$ in a Riemannian orthonormal frame. The contraction is independent of the chosen frame. There is a terminology distinction here: if $\beta$ is literally an alternating differential two-form, this trace is zero. The object needed below is the covariant two-tensor $D(X^{♭})$, generally not alternating. Interpreting that tensor as an alternating form would make the stated divergence construction impossible.

Let $D$ be the Levi-Civita connection. The Riemannian gradient, divergence and positive Laplace-Beltrami operator areThus $g({\mathrm{grad}}_{g}f,Y)=Yf$, and in coordinatesThis choice makes the Laplacian eigenvalues nonnegative on a closed Riemannian manifold. The alternative convention $+\mathrm{div}\mathrm{grad}$ reverses all eigenvalue signs, but none of the comparisons below.

For the Riemannian submersion, let $F=f\circ \pi$ and choose orthonormal horizontal lifts $E_i$ of a base frame and an orthonormal vertical frame $V_{\alpha }$. The horizontal part of $D_{E_i}{E}_i$ projects to the base Levi-Civita connection; this follows directly from its Koszul formula applied to basic horizontal fields, because their inner products and projected brackets are pulled back from the base. Therefore the horizontal trace of the Riemannian Hessian of $F$ is the pulled-back base Hessian trace. Since $V_{\alpha }F=0$, its vertical terms areWriting $H={\sum }_{\alpha }(D_{V_{\alpha }}{V}_{\alpha }{)}^{\mathrm{hor}}$ for the unnormalised fibre mean curvature vector givesThe totally geodesic submanifold condition on the fibres makes $H=0$. Consequently eigenfunction pullback by a Riemannian submersion givesHere the Riemannian submersion is onto, so the map is injective. Its image is exactly the members of the upstairs eigenspace constant on each whole fibre: such a function descends smoothly by local sections of the submersion, and the displayed identity forces the descended function to have the same eigenvalue. In particular, whenever these eigenspaces are finite dimensional, the downstairs multiplicity is no larger than the upstairs multiplicity. The upstairs eigenspace may also contain functions varying along fibres.

A flat torus is ${\mathbb{R}}^d/\Lambda$ with the metric induced by the Euclidean metric, where $\Lambda$ is a full-rank Euclidean lattice. Its dual lattice is ${\Lambda }^{\ast }=\{w:⟨w,v⟩\in \mathbb{Z}\text{ for all }v\in \Lambda \}$. The characters $e_w(x)=e^{2\pi i⟨w,x⟩}$ descend to the quotient and form an orthogonal Fourier basis of $L^2$ on it. Direct differentiation gives the spectrum of a flat torus:with multiplicity equal to the number of dual lattice vectors of each length. In particular zero has multiplicity one.

Let the given four-dimensional isospectral manifolds be ${\mathbb{R}}^4/{\Lambda }_1$ and ${\mathbb{R}}^4/{\Lambda }_2$. For $d=4+k$, take their products with the same $k$ circles of length $\epsilon$, equivalently the lattices ${\Lambda }_i\oplus \epsilon {\mathbb{Z}}^k$. The dual lattices are ${\Lambda }_i^{\ast }\oplus {\epsilon }^{-1}{\mathbb{Z}}^k$, so their Laplacian eigenvalues areThe given equality of base spectra with multiplicities makes these product spectra identical.

Choose $\epsilon$ smaller than the shortest nonzero vector length in either original Euclidean lattice. By short-vector cancellation for flat tori, the shortest nonzero vectors in the product lattice are exactly $\pm \epsilon e_j$ in its new factor. To spell out the nonisometry argument, any isometry of the product flat tori lifts to an isometry of their universal Euclidean covers. This is affine, and its orthogonal linear part maps one product lattice onto the other. It must map the span of their shortest vectors, namely the new $k$-dimensional factor, onto itself; orthogonality then preserves the original four-dimensional factor and maps ${\Lambda }_1$ onto ${\Lambda }_2$. That would give an isometry of the original flat tori, contrary to the assumption. Isospectral nonisometric flat tori therefore exist in every dimension $d\ge 4$.

Use the positive Laplace-Beltrami operator on a compact Riemannian manifold without boundary. Its heat kernel is the smooth integral kernel $K(t,x,y)$, $t>0$, of the solution operators for the heat equation:The last condition holds, in particular, in $L^2$ for smooth initial functions. The heat trace isThe heat invariants are the coefficients in the small-time heat kernel expansion of this trace:Equivalently, the local diagonal coefficients $u_j(x,x)$ integrate to $a_j$; $R_g$ is the scalar curvature. Some normalizations absorb $(4\pi {)}^{-d/2}$ into the coefficients. These definitions here use the displayed normalization consistently.

To prove uniqueness, let $K_1,K_2$ satisfy the heat kernel conditions and let $w$ be the difference of their solutions with the same smooth initial function. It solves $w_t=-\Delta w$ and tends to zero in $L^2$ as $t\downarrow 0$. The heat equation energy identity, using integration by parts on the closed manifold, gives for $t>0$Integrating from $s>0$ to $t$ gives $\Vert w(t){\Vert }_{L^2}\le \Vert w(s){\Vert }_{L^2}$; sending $s\downarrow 0$ proves $w(t)=0$. Thus for every smooth test function $f$,For fixed positive $t$ and $x$, smoothness makes this difference a smooth function of $y$. Taking its conjugate as test function shows its squared norm is zero. The heat kernel is unique.

The compact elliptic spectral theorem for the self-adjoint elliptic Laplace-Beltrami operator gives a complete orthonormal basis of smooth eigenfunctions ${\varphi }_j$, with discrete eigenvalues $0\le {\lambda }_0\le {\lambda }_1\le \cdots$ counted with multiplicity. Solving the heat equation for initial data ${\varphi }_j$ gives $e^{-t{\lambda }_j}{\varphi }_j$, so the unique heat kernel has its eigenfunction expansion and its heat trace isThe series converges for every $t>0$. This proves that the spectrum with multiplicities determines the heat trace.

Conversely, its large-time behaviour recovers the distinct eigenvalues and their multiplicities successively. The multiplicity of zero is $m_0={\lim }_{t\to \infty }{Z}_M(t)$. After subtracting the already determined terms, letbe the positive remaining sum, and let $\mu$ be its smallest remaining eigenvalue. ThenIndeed, after factoring $e^{-t\mu }$, the next distinct eigenvalue is separated from $\mu$ by a positive gap. The remaining tail tends to zero: for fixed $t_0>0$ its absolute sum is bounded by a constant times $e^{-({\mu }_{\mathrm{next}}-\mu )(t-t_0)}$, using finiteness of $Z_M(t_0)$. Induction recovers every spectral term. Thus the heat trace and the spectrum with multiplicities determine one another.

Finally the leading coefficient among the heat invariants is positive, soThe common heat trace of isospectral manifolds therefore forces them to have the same dimension. Once the dimension is recovered, its leading coefficient also recovers their common volume.

Let $p=|a|$, $q=|b|$ and $r=|ab|$ denote element orders. Remove three points from an oriented sphere. Its fundamental group has peripheral generators ${\gamma }_a,{\gamma }_b,{\gamma }_c$ with ${\gamma }_a{\gamma }_b{\gamma }_c=1$. Send them to $a,b,(ab{)}^{-1}$. Since $a,b$ generate $T$, the kernel defines a connected regular covering of the punctured sphere with degree $|T|$ and deck transformation group $T$.

Complete the lifted punctures by disks. Above the three removed points the monodromy cycles have lengths $p,q,r$, respectively, and the disk projection has the local form $z\mapsto z^p$, $z^q$ or $z^r$. This gives a closed connected oriented topological surface $M$ and an orientation-preserving $T$-action. At the three kinds of filled-in points the stabiliser subgroups are conjugates of $⟨a⟩$, $⟨b⟩$ and $⟨ab⟩$, respectively. Every other point has trivial stabiliser subgroup. Orders equal to one simply correspond to unbranched filled points.

The punctured sphere has Euler characteristic $-1$, so its cover has Euler characteristic $-|T|$. Filling the ends adds $|T|/p+|T|/q+|T|/r$ disks. This proves the formula for the branched surface action from two generators:

For a subgroup $U$, its quotient map is always regular as a branched quotient. It is an ordinary unbranched normal covering map exactly when $U$ acts freely, that is,Necessity follows from the listed stabiliser subgroups. Conversely, these conditions eliminate every possible nontrivial point stabiliser. A free finite action gives a covering space, and its deck transformations are $U$, acting transitively on each fibre; hence the cover is normal. In this caseThe condition $U⊴T$ concerns regularity of the different map $U\backslash M\to T\backslash M$; it is not needed for $M\to U\backslash M$ when the action of $U$ is free.

For the specified order-$96$ group, $p=q=3$. Set $c=ab$ and $z=c^3$. Since $z$ has order two, $c^6=1$ and $c^3\ne 1$, so $|c|$ is either two or six. The order-two alternative would make the displayed formula give $\chi (M)=96(1/3+1/3+1/2-1)=16$. That is impossible for a connected closed oriented topological surface, whose Euler characteristic is $2-2g\le 2$. Thus $|c|=6$, and

Each $U_i$ has order eight. By Lagrange's theorem, it contains no element of order three or six, so it meets all conjugates of $⟨a⟩$ and $⟨b⟩$ trivially. In $⟨c⟩$, the only nonidentity element whose order divides eight is $z=c^3$. Since $z$ is central, it is also the unique such element in every conjugate of $⟨c⟩$, and it is excluded from both $U_i$ by hypothesis. Both subgroup actions are therefore free. The quotients are closed oriented surfaces with

The branched-cover construction supplies a smooth structure, with the stabilisers acting by rotations in the disk charts. Average a smooth Riemannian metric over $T$ to make the whole action isometric. The Gassmann equivalence of $U_1,U_2$ then gives isospectral manifolds: for every upstairs eigenspace with character ${\chi }_{\lambda }$, its downstairs multiplicities are $|U_i{|}^{-1}{\sum }_{u\in U_i}{\chi }_{\lambda }(u)$, and these are equal because the two subgroups have equal intersection sizes with every conjugacy class. This proves the spectral assertion for any such invariant metric.

To guarantee nonisometry, choose the invariant metric with curvature markers distinguishing finite-cover quotients. Here are the geometric details, since nonconjugacy alone would not guarantee nonisometry for an arbitrary metric. Choose a small disk in the free locus of the $T$-action, with all its translates disjoint. In that disk choose three noncollinear points in a geodesically convex neighbourhood. Supported conformal smooth bump functions can give them three distinct high values as strict Gaussian curvature maxima, occurring nowhere else in the unmodified metric. Arrange these three maxima as one tight cluster, much closer to each other than to any translated cluster, and copy the same modification to every translate. Curvature and its first few derivatives can be adjusted by the second and higher derivatives of the conformal factor; narrow bumps give high distinguished curvature with small change in metric distances. Further small supported perturbations remove accidental maxima of these exact heights. Thus the resulting metric remains smooth and $T$-invariant, and its marked clusters are intrinsic and uniquely labelled by their three curvature values.

Any proposed quotient isometry $F:U_1\backslash M\to U_2\backslash M$ must carry a whole labelled cluster to another. Lift its restriction near that cluster to the corresponding two marked disks in $M$, and identify these disks using the appropriate group element $t\in T$. The resulting local isometry fixes all three marked points. Its differential at the first point fixes the two independent directions of the unique geodesics to the other two; hence it is the identity, and the local isometry itself is the identity near that point. If $p_i:M\to U_i\backslash M$ are the quotient maps, we therefore haveon an open set. Both sides are local isometries on the connected complete Riemannian surface $M$. A local isometry is determined along each geodesic by its value and differential at the initial point, so the equality propagates to all of $M$.

For $u\in U_1$, this equality implies $p_2\circ tu=p_2\circ t$. Evaluating at a free point of the $T$-action shows $tu{t}^{-1}\in U_2$; equivalently these maps differ by a deck transformation of $p_2$. Thus $t{U}_1{t}^{-1}\subseteq U_2$, and equality of subgroup orders gives equality. This contradicts their assumed nonconjugacy. For this marked invariant metric the two quotients are nonisometric, isospectral Riemannian surfaces of genus two. No constant-curvature requirement is imposed or used.

Fix a bounded tile in Euclidean space $\Omega$ with labelled flat boundary faces. For each face label $s$, encode an assembly of $n$ copies by an involutive real symmetric signed permutation matrix $M_s$: a paired interface contributes a two-by-two swap, an unglued face with a Neumann boundary condition contributes $+1$ on the diagonal, and an unglued face with a Dirichlet boundary condition contributes $-1$. Use the same face-identification convention on each copy. If a second assembly has matrices $N_s$ and there is an invertible constant matrix $C$ satisfyingthen the map $u=(u_1,\dots ,u_n)\mapsto Cu$ is an isomorphism between corresponding Laplace-Beltrami operator eigenspaces. This is the transplantation theorem.

To prove it, restrict an eigenfunction to the tile copies and identify each copy with $\Omega$. Along a face, let $a$ be its vector of boundary values and $b$ its vector of outward normal derivatives. At a glued interface the values agree and the outward derivatives are opposite. At an unglued face with a Neumann boundary condition the derivative vanishes, and at a face with a Dirichlet boundary condition the value vanishes. All three cases are encoded byThe intertwining identities give $N_s(Ca)=Ca$ and $N_s(Cb)=-Cb$, so the transplanted function satisfies the second assembly's interface and boundary conditions. On each tile its Laplace-Beltrami operator is the same constant linear combination of the original Laplacians, hence it has the same eigenvalue. Across an interface, matching values and opposite normal derivatives make it a weak solution; elliptic regularity, or the allowed reflection principle, makes it a genuine smooth solution there. At corners the operator is understood by its finite-energy, Friedrichs extension; the same matching argument applies to its quadratic form domain. Applying $C^{-1}$ reverses the construction. Thus dimensions of all eigenspaces agree, including multiplicities, and the two specified boundary-value spectra coincide.

The requested pair exists with explicitly specified mixed boundary conditions and polygonal boundary. Take four copies $Q_0,Q_1,Q_2,Q_3$ of an octagonal regular polygon in Euclidean space. Label four alternate sides $a,b,c,d$ in cyclic order; the four intervening sides always have a Neumann boundary condition. The labelled sides are disjoint even at their endpoints. Parameterize every side in the positive boundary direction of the reference octagon. When pairing same-labelled sides, identify equal parameters. This produces a locally flat surface with boundary: a neighbourhood of an interior point of an interface consists of two half-disks glued along their diameters. All tile vertices remain boundary points, with boundary angles $3\pi /4$ or $3\pi /2$; there are no interior cone points.

Here are the complete gluing and boundary conditions. An entry $ij$ pairs the corresponding sides of $Q_i,Q_j$. The indicated unpaired sides have a Dirichlet boundary condition; every other unpaired side has a Neumann boundary condition.Both assemblies are compact and connected. The Laplace-Beltrami operator is the nonnegative operator $P=-\Delta$ associated with the finite-energy quadratic form $\int |\nabla u{|}^2$, with zero Sobolev trace on the specified Dirichlet arcs. Thus the mixed boundary condition is defined without requiring classical differentiability at a corner or at a Dirichlet–Neumann junction.

The four signed permutation matrices for $D$ areFor $B$ they areUse the Hadamard matrixThe intertwining relations can be checked on all columns at once. Writing the columns of $C$ as $h_0,h_1,h_2,h_3$, direct application of the displayed $N_s$ givesFor the four intervening faces both assemblies have the identity face matrix, so these intertwining relations hold as well. Also $C^{T}C=C^2=4I$. This four-tile mixed-boundary transplantation between a disk and a nonorientable surface gives a bijection of every pair of eigenspaces, preserving eigenvalues and multiplicities.

There is also a direct argument on the quadratic form domains that handles all boundary corners. The orthogonal matrix $U=C/2$ acts on tile restrictions, preserves their L2 inner product and the sum of their Dirichlet energies, and commutes with taking a Sobolev trace. The condition $M_{s}u=u$ on each face encodes equality of the two traces at a paired interface and vanishing trace at an unpaired Dirichlet face. Intertwining takes these conditions exactly to $N_s(Uu)=Uu$. Consequently $U$ maps the whole form domain onto the other form domain and preserves its quadratic form. The associated self-adjoint operators are unitarily equivalent, proving isospectrality for the stated mixed boundary conditions without any assumption of smoothness at a Dirichlet–Neumann junction.

We now verify the actual topology using the topology of disk tiles glued along disjoint boundary arcs. Because the gluing arcs in each tile are disjoint, the resulting surface with boundary is a thickening of its gluing graph, with one vertex for each closed disc tile and one edge for each paired arc. One can construct this homotopy equivalence by taking a small central closed disc in each tile and narrow strips to its paired arcs. What remains is a collar neighbourhood of the boundary, which retracts onto this connected system of closed discs and strips. Shrinking each central closed disc and each strip then gives the gluing graph. Equivalently, begin with a spanning tree of the gluing graph: its closed disc tiles joined along boundary arcs form one closed disc, and every remaining pair of boundary arcs adds one band. Thus for $n$ tiles and $e$ paired arcs the fundamental group is a free group of rank $e-n+1$, and the Euler characteristic is $n-e$.

For $D$, the gluing graph is the pathGluing a new closed disc to a closed disc along one boundary arc again gives a closed disc; applying this three times proves that $D$ is homeomorphic to a closed disc. In particularFor $B$, the five paired arcs give the graph edges $13,12,02,01,23$. It is connected, and henceHere $F_2$ is the free group on two generators, so $B$ is not simply connected. To decide orientation of a surface, give every reference tile its counterclockwise orientation. Our same-parameter gluing preserves boundary directions, so a global orientation requires opposite orientation signs on paired tiles: the two induced boundary orientations at an interior interface must be opposite. Such signs exist exactly when the gluing graph is bipartite. In $B$ the cyclehas length three. Following it forces the sign to reverse three times and return to its starting tile with the opposite sign. This proves that $B$ is a nonorientable surface. It has nonempty boundary because every intervening octagon side remains unpaired. Thus $D$ is a simply-connected orientable surface with boundary, whereas its isospectral partner $B$ is neither simply connected nor orientable, exactly as required, for the specified mixed boundary problem.

The boundary qualification is essential. With smooth boundary and the same pure Dirichlet boundary condition or pure Neumann boundary condition on both surfaces, the following existing obstruction applies. Polygonal corners and the mixed boundary conditions in the construction above place it outside that obstruction.

With the interpretation of compact smooth Riemannian surfaces with the same pure Dirichlet boundary condition or pure Neumann boundary condition, the printed simply-connected/nonorientable pair cannot exist. The two-dimensional heat trace has the following expansion (the smooth-boundary coefficients are given in the heat-coefficient formulas, Theorems 5.1–5.2):where the boundary sign is minus for Dirichlet and plus for Neumann. The Gauss-Bonnet theorem identifies the constant term with $\chi (M)/6$. Isospectrality therefore forces equal Euler characteristics. A compact connected simply-connected surface with boundary is a disk and has $\chi =1$. A connected nonorientable surface with $c\ge 1$ crosscaps and $b\ge 1$ boundary components has $\chi =2-c-b\le 0$. Thus no such smooth pure-Dirichlet or pure-Neumann pair exists. This is the Euler characteristic from the heat trace of a bordered surface.

Piecewise-smooth boundaries introduce corner terms, and mixed boundary conditions change the heat coefficients, as the explicit construction above demonstrates. Moreover, in the usual pure-Neumann polygonal reflection-transplantation setup, the pure Neumann reflection transplantation preserves Euler characteristic argument applies even with corners. The numbers of tiles agree. For each face label the number of edge classes is $(n+\mathrm{tr}{M}_s)/2$, so these agree too. For each tile-vertex label, vertex classes are orbits of the subgroup generated by incident-face involutions. Their number is the dimension of invariant vectors in that permutation representation, which is unchanged by the intertwiner. Hence the numbers of vertices agree, proving equal $V-E+F$. Thus this standard tiled construction also cannot supply a disk on one side and a nonorientable surface on the other.

For comparison, the following existing construction gives an orientable and a nonorientable pair with the same spectrum for the pure Neumann boundary condition. Both have nontrivial fundamental groups; this is a separate example, not the construction used to meet the simply-connected requirement above. This eight-tile Neumann transplantation across orientability makes the comparison explicit. The finite seed is the affine modulo-eight example in the original orientability construction; the intertwiner and the topology checks are derived below. Take eight copies of a regular Euclidean hexagon and label three pairwise nonadjacent sides $\sigma ,t,u$. Every remaining side is a side with a Neumann boundary condition. Index the copies by $j\in \mathbb{Z}/8\mathbb{Z}$. For the first assembly pair same-labelled faces according tofor the second, usewith all arithmetic modulo eight. A fixed point means an unglued side with a Neumann boundary condition. Pair faces by the reflection convention, matching the same endpoint parameter in the reference tiles. The distinguished sides are disjoint, so no interior vertex singularities are created; the resulting surfaces have polygonal boundary corners.

Let $R{e}_j=e_{j+1}$ on the eight-dimensional tile-index space and setIt commutes with the permutation matrices for $t$ and $\sigma$, because its coefficients are invariant under reversing cyclic order. For $u_1,u_2$, write $C_{ij}=c_{i-j}$, where $c_0=c_4=1$, $c_1=c_7=2$, and all other $c_r$ vanish. The required relation $C{M}_{u_1}=M_{u_2}C$ is exactly $c_r=c_{3r+4}$, verified from these four nonzero entries. Thus all three intertwining relations hold.

The cyclic Fourier vectors diagonalize $C$. Its eigenvalues arenamely $6,2\sqrt{2},2,-2\sqrt{2},-2,-2\sqrt{2},2,2\sqrt{2}$. None is zero, so $C$ is invertible and the two assemblies are isospectral manifolds for the Neumann boundary condition.

We prove their topology as well. Their gluing graphs are connected: the first contains the spanning path $0,1,3,6,2,7,5,4$, and the second contains the spanning cycle $0,1,7,2,6,3,5,4,0$. Each has eight tile vertices and ten paired-face edges. Since tiles are disks and gluing arcs are disjoint, each surface retracts up to homotopy to this multigraph; thus its Euler characteristic is $8-10=-2$ and its fundamental group is a free group of rank three.

An orientation assigns a sign to each reference tile, and a reflection gluing requires opposite signs at the ends of every graph edge. The second graph admits this assignment: take sign $+$ on $\{0,5,6,7\}$ and sign $-$ on $\{1,2,3,4\}$. Every paired face joins opposite signs, proving orientability. In the first graph the cyclehas odd length. Transporting orientation around it reverses the initial sign, proving nonorientability. Hence the first surface is also not simply connected. The second surface is orientable but likewise not simply connected; the disk required in the question is supplied instead by the four-tile mixed-boundary construction above.

The four-tile construction resolves the printed simply-connected requirement under explicit mixed boundary conditions. The eight-tile comparison only concerns orientation of a surface under a pure Neumann boundary condition.

Let ${\Lambda }_1,{\Lambda }_2\subset {\mathbb{R}}^d$ be full-rank Euclidean lattices, and choose $\epsilon$ smaller than every nonzero vector length in either. The shortest vectors of ${\Lambda }_i\oplus \epsilon {\mathbb{Z}}^k$ are exactly the new coordinate vectors and their negatives. Any isometry between the associated flat tori lifts to an affine Euclidean isometry, whose orthogonal part maps the lattices. It must preserve the span of the shortest vectors, hence also its orthogonal complement, and therefore maps ${\Lambda }_1$ to ${\Lambda }_2$. Consequently a nonisometric pair remains nonisometric after adding sufficiently short equal circle factors. The spectrum of a flat torus shows that adding the same factors preserves isospectral manifolds.

Spectral geometry studies which properties of a Riemannian manifold can be recovered from the spectrum of geometric differential operators, especially the Laplace-Beltrami operator. Explicit spherical harmonics, Fourier series on a flat torus, and isospectral manifolds illustrate three different aspects of the problem.

Let a finite group $T$ be generated by elements of exact orders $p,q$ whose product has exact order $r$, with $p^{-1}+q^{-1}+r^{-1}<1$. The epimorphism from the hyperbolic triangle group sending its generators to these elements has torsion-free kernel. Its quotient of the hyperbolic plane is a closed surface with a $T$ action. A subgroup acts freely precisely when it meets no conjugate of a nonidentity power of the three elliptic generators. Free almost conjugate subgroups give isospectral manifolds by the Sunada theorem.