Let denote the Half-range Fourier transforms of , and let denote those of the upward derivatives . Also set
Integrate the tangential derivatives in the three spectral functions by parts. For example, , and . All corner values cancel in . The result is the semistrip Laplace spectral global relation
For , use the sine transforms
Reality of makes and real, so taking imaginary parts removes the unknown left-side normal derivative. At and the two consequences of the global relation are
Subtracting these equations and then substituting back gives the semistrip Dirichlet-to-Neumann sine transforms entirely in terms of the prescribed Dirichlet boundary data:
These are the requested transforms of itself. The bottom outward Neumann boundary data have the opposite sign. The apparent division at is handled by a continuous limit whenever the corresponding moments exist; there is no singularity for any .
We construct the two domains explicitly from seven congruent triangles. Let a scalene reference triangle have sides opposite its three labelled vertices. Choose its three side lengths in a small open neighbourhood of , with a fixed strict ordering. Thus the three angles are distinct and lie between and . Each copy keeps these side labels.
For the first assembly glue the following tile pairs, reflecting across the indicated side:
Each coloured adjacency graph is a tree with central tile and three two-tile arms. Place tile in the plane and obtain every other tile by the prescribed reflections. At the equilateral reference shape the copies have disjoint interiors and form a polygonal disk; the same holds throughout a sufficiently small neighbourhood, since the only contacts between tiles are the prescribed edges and their incident vertices. In particular each assembly is a bounded connected planar domain.
Figure 1.
Seven-triangle point and line assemblies with matching Dirichlet spectra
. Details.
To prove equality of spectra, index tiles by the seven nonzero vectors , using with label . Let
These involutions act on points by and on nonzero covectors by . Their swaps are exactly the two columns of the table. Write for the corresponding permutation matrices. The Fano plane incidence matrix
satisfies , since incidence is preserved by a linear transformation and its inverse transpose. Each point lies on three lines, while two different points have one common line. Therefore
where is the all-ones matrix, proving that is invertible.
For the Dirichlet boundary condition, a glued side has a positive swap and an exposed side has diagonal entry . Both trees have the same bipartition, with sign matrix
Their signed reflection matrices are and . On a glued edge its endpoint signs are opposite, so the off-diagonal entries are ; on an exposed edge the entry is . Consequently the invertible matrix satisfies
We spell out why this gives isospectrality. Pull each function back from the seven tiles to the reference triangle. Its seven boundary values along side form a vector satisfying : this imposes equal values at glued edges and zero values at exposed edges. Outward normal derivatives of a smooth eigenfunction satisfy , giving opposite normal derivatives on glued edges. The identities above preserve both conditions under . Since is constant, it also commutes with the scalar Laplacian on each triangle.
One can include all weak eigenfunctions without any corner regularity issue by orthogonalization of a transplantation matrix. The matrix commutes with the symmetric , so is orthogonal and still intertwines . On the disjoint tiles, preserves the norm and the sum of the Dirichlet energies; it maps the trace-matching form domain bijectively to . Thus it is a unitary equivalence of the Dirichlet Laplacians, proving equality of every eigenvalue with its multiplicity.
Finally, the only reentrant vertices of either polygon are the three vertices of tile , with interior angles four times the corresponding reference-triangle angles. All other boundary angles are less than . These three distinct angles force any hypothetical isometry to map each labelled central vertex to its counterpart. An isometry between connected open planar domains is a restriction of a rigid motion: its derivative is orthogonal, and preservation of the flat Levi-Civita connection makes that derivative constant. After aligning the two central triangles, this rigid motion fixes three noncollinear points and is the identity. But the two domains do not coincide: the tile beyond side of the central triangle has its next attachment along in the point assembly and along in the line assembly. These occupy different open regions. This contradiction proves nonisometry.
Therefore the three freely varying side lengths give a three-parameter family of nonisometric Dirichlet-isospectral bounded planar domains. The parameters include scale; no congruence normalization removes it.
Use the data norm. The coefficients of the measurement error are , so Parseval's identity gives
An unregularized continuation of these measurements to replaces by in the outgoing spherical harmonic expansion. Whenever this continuation has a finite trace, its error is
Choose a measurement error supported on a single normalized spherical harmonic, with and all other coefficients zero. The data error has norm exactly , but
This perturbation theory is a finite expansion, so it defines an actual outgoing solution of the Helmholtz equation outside any sphere of positive radius. It is not merely an example of nonconvergent formal coefficients. Taking instead gives data errors tending to zero and trace errors of norm one, which directly proves discontinuity of inversion. General square-summable measurement noise can even violate the spherical far-field range condition, leaving no continuation to the chosen sphere.
For obstacle reconstruction one would continue the field toward the unknown boundary and look for a surface on which , the Dirichlet boundary condition. The unbounded spherical far-to-near-field continuation shows why small far-field errors can cause arbitrarily large errors in the field used to estimate that surface. Locally, a differentiable normal displacement of a candidate boundary must obey the linearized condition . Thus a large continued-field error, or a small normal derivative, defeats an uncontrolled boundary estimate.
There is a necessary qualification to the wording about errors in the scatterer itself. The calculation proves an unbounded field-continuation error. An arbitrary perturbed outgoing field need not correspond to any Dirichlet obstacle, and no admissible class or metric on obstacle shapes is specified. It therefore does not by itself prove an arbitrarily large geometric error in every constrained reconstruction method. Such a claim requires a shape model and an analysis of the boundary-extraction step. The instability displayed here is exactly the mechanism that an unregularized reconstruction must control.
Fix a bounded tile in Euclidean space with labelled flat boundary faces. For each face label , encode an assembly of copies by an involutive real symmetric signed permutation matrix : a paired interface contributes a two-by-two swap, an unglued face with a Neumann boundary condition contributes on the diagonal, and an unglued face with a Dirichlet boundary condition contributes . Use the same face-identification convention on each copy. If a second assembly has matrices and there is an invertible constant matrix satisfying
then the map 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 . Along a face, let be its vector of boundary values and 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 by
The intertwining identities give and , 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 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 of an octagonal regular polygon in Euclidean space. Label four alternate sides 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 or ; there are no interior cone points.
Here are the complete gluing and boundary conditions. An entry pairs the corresponding sides of . 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 associated with the finite-energy quadratic form , 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 are
For they are
Use the Hadamard matrix
The intertwining relations can be checked on all columns at once. Writing the columns of as , direct application of the displayed gives
For the four intervening faces both assemblies have the identity face matrix, so these intertwining relations hold as well. Also . 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 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 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 . Consequently 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 tiles and paired arcs the fundamental group is a free group of rank , and the Euler characteristic is .
For , the gluing graph is the path
Gluing a new closed disc to a closed disc along one boundary arc again gives a closed disc; applying this three times proves that is homeomorphic to a closed disc. In particular
For , the five paired arcs give the graph edges . It is connected, and hence
Here is the free group on two generators, so 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 the cycle
has length three. Following it forces the sign to reverse three times and return to its starting tile with the opposite sign. This proves that is a nonorientable surface. It has nonempty boundary because every intervening octagon side remains unpaired. Thus is a simply-connected orientable surface with boundary, whereas its isospectral partner 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 . Isospectrality therefore forces equal Euler characteristics. A compact connected simply-connected surface with boundary is a disk and has . A connected nonorientable surface with crosscaps and boundary components has . 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 , 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 . 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 . Every remaining side is a side with a Neumann boundary condition. Index the copies by . For the first assembly pair same-labelled faces according to
for the second, use
with 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 on the eight-dimensional tile-index space and set
It commutes with the permutation matrices for and , because its coefficients are invariant under reversing cyclic order. For , write , where , , and all other vanish. The required relation is exactly , verified from these four nonzero entries. Thus all three intertwining relations hold.
The cyclic Fourier vectors diagonalize . Its eigenvalues are
namely . None is zero, so 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 , and the second contains the spanning cycle . 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 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 and sign on . Every paired face joins opposite signs, proving orientability. In the first graph the cycle
has 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.
Work above the plane, with time dependence and away from grazing incidence. Put and . The incident field and its chosen flat Dirichlet reflection cancel in value at the plane, while their normal derivative is . Thus the diffracted correction has zero trace on the negative half-line and prescribed normal derivative on the positive half-line.
There is an incompatibility in the stated edge data: a normal derivative that is identically zero on the Neumann half-line cannot diverge as that half-line approaches the junction. The inverse-square-root normal-derivative singularity belongs to the Dirichlet side, ; on the total normal derivative remains zero. The following solution satisfies the mixed conditions and bounded-field edge behavior, with that corrected side.
Take the transform and its inverse as
Let be the transform of the diffracted boundary value, and the transform of its unknown normal derivative on the negative half-line. The outgoing transformed field is , where
Continue into the upper half-plane. Choose branches so that , outgoing propagating waves have positive vertical wavenumber, and evanescent waves decay. The branch point lies below and belongs to ; lies above and belongs to . Since , the boundary transform equation is
The forcing splits by subtracting its numerator at the incident pole:
The pole in is removable; is analytic below , while is analytic above it. Hence is entire. The bounded mixed-boundary edge solution has and , making both sides decay at infinity. These orders can also be obtained locally from the leading wedge solution , which is Neumann at and Dirichlet at . Liouville's theorem sets the entire remainder to zero. Therefore
The solution of mixed Dirichlet-Neumann half-plane diffraction is consequently
For , take from left to right inside and above the incident pole . Cuts from run into the lower half-plane and cuts from into the upper half-plane, without crossing . For example choose with . Then let , keeping the contour above and and below . Any subsequent deformation for asymptotic evaluation must preserve these bypasses or include the residues of crossed poles. This prescription fixes the reflected-wave contributions and the radiation condition; a principal-value integral alone would not do so.
As a boundary check, . On the positive half-line the lower-contour residue gives , cancelling the incident-plus-reflected derivative. The boundary trace vanishes on the negative half-line by upper-half-plane analyticity. Finally gives a square-root boundary value and hence a bounded field at the junction; the singular derivative occurs on its Dirichlet side.
Consider two finite assemblies of the same number of congruent bounded Euclidean tiles. Label their matching faces so the same prescribed face identifications are used in both assemblies. For each face type , encode the first assembly by a symmetric involution : a paired face exchanges its two tile indices, an exterior Dirichlet boundary condition contributes diagonal , and an exterior Neumann boundary condition contributes diagonal . Let encode the second assembly. Face coordinates are pulled back to the same reference face; if identifications have additional face isometries, these pullbacks must be included in the intertwining operators.
The transplantation theorem says that a constant matrix with for every face type sends the vector of tile restrictions of a Laplacian eigenfunction on the first assembly to one on the second, by . If is invertible, it is a bijection of each eigenspace and the assemblies are isospectral with multiplicities.
To prove it, each component solves the same interior equation , and a constant linear combination solves that equation too. On a face, let be the vector of boundary values and the vector of outward normal derivatives. All matching and exterior conditions are exactly
For an internal face these equations say that values agree and the two outward derivatives sum to zero. For diagonal they impose zero value, and for diagonal zero normal derivative. The intertwining identity gives and , so the transplanted functions match and satisfy the correct boundary conditions. Applying proves the eigenspace bijection. At corners the same reasoning is interpreted in the finite-energy weak domain: no value jump and cancellation of normal flux prevent an extra distributional source. It gives the usual self-adjoint Dirichlet/mixed Laplacians on polygonal assemblies.
Here is an explicit transplantation by reflection parity. Take a rectangle of width and height , with the Dirichlet boundary condition on every exterior side. Cut it at its vertical midline into two congruent half-rectangles, pulling the right-hand half back by reflection. The midline gluing matrix is
Replace the connected assembly by the disjoint union of two width-, height- rectangles: the first has Dirichlet conditions on all sides; the second has Dirichlet on three sides and Neumann on the side representing the cut. Its corresponding matrix is . The invertible, indeed orthogonal, matrix
intertwines the midline conditions. Every other face has matrix on both assemblies, which automatically intertwines. The transplantation theorem proves equality of the complete spectra, although one surface is connected and the other has two components with the required uniform and mixed conditions.
Equivalently, odd reflection eigenfunctions vanish on the cut, while even reflection eigenfunctions have zero normal derivative there. A direct separation-of-variables check gives the connected rectangle's eigenvalues
Even gives the fully Dirichlet half-rectangle spectrum, and odd , , gives the mixed half-rectangle spectrum. The split preserves every multiplicity, including coincidences between different pairs of indices.
Figure 1.
Reflection splits a connected Dirichlet rectangle into a disjoint Dirichlet rectangle and a mixed Dirichlet-Neumann rectangle
. Details.
The Cauchy problem for a partial differential equation here prescribes both the value and the normal derivative , together with . Only one of these traces would be boundary data for a usual elliptic boundary problem, rather than full Cauchy data.
Every real hypersurface is a non-characteristic hypersurface for the Laplace equation. In local real analytic coordinates flattening the real analytic hypersurface , the coefficient of the second transverse derivative is nonzero: its principal coefficient is the squared length of the conormal. The equation can therefore be solved for that second derivative. The normal derivative data determine the transverse first derivative, because the coefficient relating them is nonzero and the tangential first derivatives are already determined by .
The coefficients, flattened Cauchy data, and coordinate change are all real analytic. The Cauchy-Kovalevskaya theorem applies, giving a unique local real analytic solution around each point of . This is a local existence assertion, not a claim of stable dependence in arbitrary Sobolev space norms.
The Cauchy-Kovalevskaya theorem cannot be applied to merely , non- Cauchy data. It requires real analytic data.
There is also no solution of the Laplace equation on a neighborhood of a point where one of these prescribed traces fails to be . By interior elliptic regularity, any such solution would be smooth and real analytic. On the real analytic hypersurface retained from the preceding part, both its restriction and its normal derivative would then be real analytic, hence . This contradicts the prescribed trace. There is no solution on a neighborhood of all of with the stated non- data. This does not exclude solutions near other points where the data happen to be real analytic.
For this Neumann Poisson problem, interpret the forcing in , as is automatic if it is smooth up to the boundary. Literal interior smoothness alone does not ensure the integrals or bounded functionals required in this question: for example, on is interior smooth but even diverges. Classical regularity in the converse is likewise understood up to the boundary.
Suppose the weak solution is smooth on . Testing against compactly supported test functions gives in distributions and hence pointwise. Now the weak identity and Green's first identity imply
for every smooth on . Every smooth boundary function has such an extension, so on . This proves both the interior equation and the boundary condition.
Conversely, for satisfying the equation and zero normal derivative, Green's first identity gives the weak identity for all smooth on . The density of smooth functions in a Sobolev space and the Cauchy-Schwarz inequality extend it continuously to every . Thus the classical solution is a weak solution, and the smooth weak solution is classical.
The clamped second-order Sobolev space is . On a smooth bounded domain the Sobolev trace theorem characterizes it by zero value and zero normal derivative on the boundary. In particular, its whole first-order boundary jet is zero, since tangential derivatives of the zero trace also vanish. As above, assume and classical regularity up to the boundary.
For a smooth weak solution, compactly supported test functions and two integrations by parts give
so pointwise. Membership in supplies on . Thus it is a classical solution of the clamped biharmonic problem.
Conversely, a classical solution with these traces belongs to . For every compactly supported test function, two integrations by parts give . Both sides are continuous for the norm, so the defining density of in extends this equality to every required test. The two notions agree under the stated smoothness.
The dispersion symmetry elimination of a boundary trace uses the symmetry of the dispersion relation
For , , so the global relation is valid at . Since , it gives
Substitution into the contour integral representation produces an unwanted integral . Its integrand is analytic in and decays on closing the contour upwards; Jordan lemma makes this integral zero for . Hence the Fokas method eliminates the unknown normal derivative:
All quantities here are determined by the prescribed initial and Dirichlet boundary data up to time .
For verification and for numerical evaluation it is useful to evaluate the spectral contour integrals, giving a half-line drift reflection kernel. Put
and define the half-line drift boundary kernel
Fubini's theorem, the Gaussian Fourier transform and contour deformation give the equivalent causal formula
For the reflected initial term, on the contour; the evaluated Gaussian supplies the necessary large- decay. If is not integrable, truncate the initial conditions first, evaluate, and pass to the limit using Gaussian bounds. No extra exponential-decay assumption on the original data is needed for this kernel formula.
The boundary kernel follows particularly simply from
This calculation independently checks both the sign and the coefficient of the boundary forcing.