pith. sign in

arxiv: 1907.04282 · v1 · pith:NAH4XHTRnew · submitted 2019-07-09 · 🧮 math.SP · cs.NA· math-ph· math.AP· math.MP· math.NA

Boundary integral formulations of eigenvalue problems for elliptic differential operators with singular interactions and their numerical approximation by boundary element methods

Markus Holzmann , Gerhard Unger This is my paper

Pith reviewed 2026-05-24 23:56 UTC · model grok-4.3

classification 🧮 math.SP cs.NAmath-phmath.APmath.MPmath.NA
keywords eigenvalue problemssingular delta interactionsboundary integral operatorsboundary element methodselliptic differential operatorsself-adjoint operatorsnumerical approximationconvergence analysis
0
0 comments X

The pith

Eigenvalue problems for elliptic operators with delta and delta-prime interactions admit equivalent boundary integral formulations that support boundary element method computations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines discrete eigenvalues of second-order elliptic differential operators on the whole space that include singular delta and delta-prime interactions supported on lower-dimensional sets. It first establishes that the associated operators are self-adjoint. It then derives equivalent formulations of the eigenvalue problems that are expressed entirely in terms of boundary integral operators. These integral formulations are designed so that the discrete eigenvalues and eigenfunctions can be approximated by boundary element methods, and the paper supplies convergence results together with numerical illustrations.

Core claim

The eigenvalue problems for elliptic second-order differential operators in L2(R^n) with singular delta- and delta'-interactions are equivalent to formulations involving boundary integral operators; these equivalent formulations are suitable for numerical computation of the discrete eigenvalues and eigenfunctions by boundary element methods, and the numerical approximations converge.

What carries the argument

Boundary integral operators obtained from the fundamental solution of the elliptic differential operator, which reduce the eigenvalue problem on the interaction support to an integral equation on that support.

If this is right

  • Self-adjointness of the operators guarantees that all eigenvalues are real.
  • The boundary integral formulations allow computation of both eigenvalues and eigenfunctions without volume meshing.
  • Convergence of the boundary element approximations to the true discrete eigenvalues is guaranteed under the paper's assumptions.
  • Numerical examples confirm that the method produces accurate approximations for model problems in low dimensions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same reduction might apply to other classes of singular potentials whose Green's functions are known explicitly.
  • The approach could be combined with fast multipole or hierarchical matrix techniques to treat larger interaction supports.
  • Eigenvalue tracking with respect to interaction strength parameters becomes feasible once the integral formulation is in hand.

Load-bearing premise

The original differential eigenvalue problem remains equivalent to the boundary integral formulation when the elliptic operator and the supports of the singular interactions satisfy the stated function-space conditions.

What would settle it

A concrete choice of elliptic operator and interaction support for which the eigenvalues computed from the boundary integral equation differ from those obtained by direct discretization of the original differential operator.

Figures

Figures reproduced from arXiv: 1907.04282 by Gerhard Unger, Markus Holzmann.

Figure 1
Figure 1. Figure 1: Computed eigenfunctions of Aα, α “ ´6, in the xy￾plane for the unit ball. 4.3.2. Screen. For the second numerical example we have chosen a δ-potential sup￾ported on the non-closed surface Γ :“ r0, 1s ˆ r0, 1s ˆ t0u Ă R 3 , which is referred to as screen. The interaction strength α is defined by α “ ´15χΓ, where χΓ is the characteristic function on Γ given as χΓpxq :“ # 1, for x P Γ, 0, else . Such a proble… view at source ↗
Figure 2
Figure 2. Figure 2: Computed eigenfunctions of Aα in the xy-plane for α “ ´15χr0,1sˆr0,1sˆt0u . similar as for Aα in Section 4. First, in Section 5.1 we show the self-adjointness of Bβ in L 2 pR nq and obtain the Birman-Schwinger principle to characterize the dis￾crete eigenvalues of Bβ via boundary integral operators in Proposition 5.2. Then, in Section 5.2 we discuss how these boundary integral equations can be solved nu￾me… view at source ↗
Figure 3
Figure 3. Figure 3: plots of computed eigenfunctions of Bβ in the xy-plane are given where for each exact eigenvalue one approximated eigenfunction is selected. h ˇ ˇ ˇλ p0q h ´λ p0q ˇ ˇ ˇ |λp0q | eoc ˇ ˇ ˇλpp1q h ´λ p1q ˇ ˇ ˇ |λp1q | eoc ˇ ˇ ˇλpp2q h ´λ p2q ˇ ˇ ˇ |λp2q | eoc 0.2 3.232e-3 - 1.885e-3 - 6.745e-3 - 0.1 7.099e-4 2.19 3.926e-4 2.26 1.406e-3 2.26 0.05 1.635e-4 2.11 8.958e-5 2.13 3.054e-4 2.20 [PITH_FULL_IMAGE:figu… view at source ↗
Figure 4
Figure 4. Figure 4: Computed eigenfunctions of Bβ, β ´1 “ ´0.75, for the L-shape domain Ωi “ p´1, 1q 3 zpr0, 1s 2 ˆ r´1, 1sq in the xy-plane. eigenvalue problem inside this contour, namely λ p0q h “ ´5.54, λ p1q h “ ´4.41 and λ p2q h “ ´2.94 for the mesh-size h “ 0.1. Plots of the numerical approximations of the eigenfunctions in the xy-plane are given in [PITH_FULL_IMAGE:figures/full_fig_p037_4.png] view at source ↗
read the original abstract

In this paper the discrete eigenvalues of elliptic second order differential operators in $L^2(\mathbb{R}^n)$, $n \in \mathbb{N}$, with singular $\delta$- and $\delta'$-interactions are studied. We show the self-adjointness of these operators and derive equivalent formulations for the eigenvalue problems involving boundary integral operators. These formulations are suitable for the numerical computations of the discrete eigenvalues and the corresponding eigenfunctions by boundary element methods. We provide convergence results and show numerical examples.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper studies discrete eigenvalues of elliptic second-order differential operators in L²(ℝⁿ) with singular δ- and δ'-interactions supported on hypersurfaces. It establishes self-adjointness of the associated operators, derives equivalent boundary-integral formulations of the eigenvalue problems, shows that these formulations are amenable to discretization by boundary element methods, proves convergence of the discrete eigenvalues and eigenfunctions, and presents numerical examples.

Significance. If the claimed equivalence and convergence hold, the work supplies a mathematically justified BEM framework for computing eigenvalues of Schrödinger operators with singular surface interactions, a setting that arises in quantum mechanics and scattering theory. The combination of functional-analytic self-adjointness proofs with numerical analysis is a concrete contribution to computational spectral theory.

major comments (2)
  1. [Abstract / derivation of formulations] The central equivalence claim (abstract) between the differential eigenvalue problem and the boundary-integral formulation must be verified in both directions for the δ' case. The reconstruction step that recovers the precise jump relations from a solution of the BIE relies on trace theorems and jump formulas for the normal derivative of the single-layer potential; the manuscript should state the precise Sobolev regularity assumed on the support and confirm that these formulas remain valid under the stated assumptions.
  2. [Convergence results] Convergence results for the discrete eigenvalues obtained by BEM discretization are stated in the abstract. The proof should explicitly identify the compactness or approximation properties of the boundary-integral operators that are used to pass to the limit, and should indicate whether the rate depends on the smoothness of the interaction support.
minor comments (2)
  1. Notation for the single- and double-layer operators and the associated boundary integral operators should be introduced with explicit reference to the underlying trace spaces.
  2. The numerical examples section would benefit from a brief statement of the mesh refinement strategy and the observed convergence order for at least one test case.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment below and indicate the revisions we intend to incorporate.

read point-by-point responses

  1. Referee: [Abstract / derivation of formulations] The central equivalence claim (abstract) between the differential eigenvalue problem and the boundary-integral formulation must be verified in both directions for the δ' case. The reconstruction step that recovers the precise jump relations from a solution of the BIE relies on trace theorems and jump formulas for the normal derivative of the single-layer potential; the manuscript should state the precise Sobolev regularity assumed on the support and confirm that these formulas remain valid under the stated assumptions.

    Authors: We agree that the equivalence must be established in both directions for the δ' interaction. Section 3 derives the boundary-integral formulation from the differential eigenvalue problem by applying the jump relations of the single-layer potential. The converse direction reconstructs the solution in the distributional sense and verifies that the resulting function satisfies the original jump conditions via the trace theorems. To strengthen the presentation we will add an explicit statement of the Sobolev regularity assumed on the hypersurface Γ (C^{1,1} or smoother) and a short remark confirming that the required jump formulas remain valid in the trace spaces used throughout the paper. revision: yes

  2. Referee: [Convergence results] Convergence results for the discrete eigenvalues obtained by BEM discretization are stated in the abstract. The proof should explicitly identify the compactness or approximation properties of the boundary-integral operators that are used to pass to the limit, and should indicate whether the rate depends on the smoothness of the interaction support.

    Authors: The convergence analysis in Section 4 proceeds from the approximation properties of the boundary-element spaces together with the compactness of the boundary-integral operators that appear after reformulation. We will revise the proof to name these compactness and collective-compactness arguments explicitly when passing to the limit. We will also add a remark stating that the convergence rate depends on the smoothness of the interaction support Γ, with higher regularity permitting higher-order rates in the BEM error estimates. revision: yes

Circularity Check

0 steps flagged

Derivations rely on standard potential theory and trace theorems; no reductions to self-inputs or self-citations

full rationale

The paper establishes self-adjointness of the operators with δ/δ'-interactions and derives equivalent boundary integral eigenvalue formulations via layer potential representations. These steps invoke classical jump relations, Green's identities, and Sobolev trace theorems under the stated regularity assumptions on the interaction supports. No equation reduces by construction to a fitted parameter renamed as a prediction, nor does any load-bearing premise collapse to a self-citation chain. The numerical BEM convergence analysis is independent of the analytic equivalence claim. The derivation chain is therefore self-contained against external functional-analytic benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard background results from elliptic operator theory and boundary integral operators; no free parameters, invented entities, or ad-hoc axioms are introduced in the abstract.

axioms (2)
  • domain assumption Elliptic second-order differential operators with delta and delta-prime interactions on hypersurfaces are self-adjoint in L2(R^n) under appropriate conditions on the interaction supports.
    Invoked to justify the eigenvalue problem setup before deriving the integral formulations.
  • standard math Boundary integral operators arising from the Green's function of the elliptic operator are well-defined and compact or Fredholm in the appropriate trace spaces.
    Required for the equivalence between the differential and integral formulations.

pith-pipeline@v0.9.0 · 5622 in / 1235 out tokens · 17134 ms · 2026-05-24T23:56:45.907783+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Integral formulation of Dirac singular waveguides

    math-ph 2023-12 unverdicted novelty 6.0

    Derives and analyzes a boundary integral equation for the massive Dirac equation with discontinuous mass, proves uniqueness, extends to two interfaces, and provides numerical examples of surface wave scattering.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Albeverio, F

    S. Albeverio, F. Gesztesy, R. Høegh-Krohn, and H. Holden, Solvable Models in Quantum Mechanics. With an Appendix by Pavel Exner, 2nd ed., Amer. Math. Soc. Chelsea Publish- ing, Providence, RI (2005)

    work page 2005

  2. [2]

    Antoine, F

    J. Antoine, F. Gesztesy and J. Shabani, Exactly solvable models of sphere interactions in quantum mechanics, J. Phys. A 20 (1987), 3687–3712

    work page 1987

  3. [3]

    Behrndt, P

    J. Behrndt, P. Exner, M. Holzmann, and V. Lotoreichik, Approximation of Schr¨ odinger operators with δ-interactions supported on hypersurfaces , Math. Nachr. 290(8-9) (2017), 1215–1248

    work page 2017

  4. [4]

    The Landau Hamiltonian with $\delta$-potentials supported on curves

    J. Behrndt, P. Exner, M. Holzmann, and V. Lotoreichik, The Landau Hamiltonian with δ-potentials supported on curves , arxiv:1812.09145 (2018)

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  5. [5]

    Behrndt, P

    J. Behrndt, P. Exner, M. Holzmann, V. Lotoreichik, On Dirac operators in R3 with electro- static and Lorentz scalar δ-shell interactions, To appear in Quantum Stud.: Math. Found., DOI: https://doi.org/10.1007/s40509-019-00186-6

    work page doi:10.1007/s40509-019-00186-6

  6. [6]

    Behrndt, G

    J. Behrndt, G. Grubb, M. Langer, and V. Lotoreichik, Spectral asymptotics for resolvent dif- ferences of elliptic operators with δ and δ1-interactions on hypersurfaces, J. Spectr. Theory 5(4) (2015), 697–729

    work page 2015

  7. [7]

    Behrndt, M

    J. Behrndt, M. Langer, and V. Lotoreichik, Schr¨ odinger operators withδ and δ1-potentials supported on hypersurfaces, Ann. Henri Poincar´ e14 (2013), 385–423

    work page 2013

  8. [8]

    Behrndt and J

    J. Behrndt and J. Rohleder, An inverse problem of Calder´ on type with partial data, Comm. Partial Differential Equations 37(6) (2012), 1141–1159

    work page 2012

  9. [9]

    Behrndt and J

    J. Behrndt and J. Rohleder, Spectral analysis of selfadjoint elliptic differential operators, Dirichlet-to-Neumann maps, and abstract Weyl functions , Adv. Math. 285 (2015), 1301– 1338

    work page 2015

  10. [10]

    W. J. Beyn, An integral method for solving nonlinear eigenvalue problems , Linear Algebra Appl. 432 (2012), 3839–3863

    work page 2012

  11. [11]

    Brasche, P

    J. Brasche, P. Exner, Y. Kuperin, and P. ˇSeba, Schr¨ odinger operators with singular inter- actions, J. Math. Anal. Appl. 184 (1994), 112–139

    work page 1994

  12. [12]

    Brasche, R

    J. Brasche, R. Figari, and A. Teta, Singular Schr¨ odinger operators as limits of point inter- action Hamiltonians, Potential Anal. 8 (1998), 163–178

    work page 1998

  13. [13]

    Brown, M

    M. Brown, M. Marletta, S. Naboko, and I. Wood, Boundary triplets and M-functions for non-selfadjoint operators, with applications to elliptic PDEs and block operator matrices , J. Lond. Math. Soc. (2) 77(3) (2008), 700–718

    work page 2008

  14. [14]

    Brasche and K

    J. Brasche and K. Oˇ zanov´ a,Convergence of Schr¨ odinger operators, SIAM J. Math. Anal. 39 (2007), 281–297. 38 M. HOLZMANN AND G. UNGER

    work page 2007

  15. [15]

    Exner, Leaky quantum graphs: a review, in: Analysis on graphs and its applications

    P. Exner, Leaky quantum graphs: a review, in: Analysis on graphs and its applications . Selected papers based on the Isaac Newton Institute for Mathematical Sciences programme, Cambridge, UK, 2007. Proc. Symp. Pure Math. 77 (2008), 523–564

    work page 2007

  16. [16]

    Exner and H

    P. Exner and H. Kovaˇ r´ ık, Quantum Waveguides. Theoretical and Mathematical Physics, Springer (2015)

    work page 2015

  17. [17]

    Exner and K

    P. Exner and K. Nˇ emcov´ a,Leaky quantum graphs: approximations by point-interaction Hamiltonians, J. Phys. A 36 (2003), 10173–10193

    work page 2003

  18. [18]

    Figotin and P

    A. Figotin and P. Kuchment, Band-gap structure of spectra of periodic dielectric and acous- tic media. II. Two-dimensional photonic crystals , SIAM J. Appl. Math. 56(6) (1996), 1561– 1620

    work page 1996

  19. [19]

    I. C. Gohberg, S. Goldberg, M. A. Kaashoek, Classes of Linear Operators. Vol. I. Birkh¨ auser Verlag, Basel (1990)

    work page 1990

  20. [20]

    I. C. Gohberg and E. I. Sigal, An operator generalization of the logarithmic residue theorem and Rouch´ e’s theorem, Math. USSR-Sb. 13 (1971), 603–625

    work page 1971

  21. [21]

    Karma, Approximation in eigenvalue problems for holomorphic Fredholm operator func- tions

    O. Karma, Approximation in eigenvalue problems for holomorphic Fredholm operator func- tions. I, Numer. Funct. Anal. Optim. 17 (1996), 365–387

    work page 1996

  22. [22]

    Karma, Approximation in eigenvalue problems for holomorphic Fredholm operator func- tions

    O. Karma, Approximation in eigenvalue problems for holomorphic Fredholm operator func- tions. II. (Convergence rate) , Numer. Funct. Anal. Optim. 17 (1996), 389–408

    work page 1996

  23. [23]

    Kato, Perturbation theory for linear operators

    T. Kato, Perturbation theory for linear operators. Springer-Verlag, Berlin, reprint of the 1980 edition (1995)

    work page 1980

  24. [24]

    Kimeswenger, O

    A. Kimeswenger, O. Steinbach, and G. Unger, Coupled finite and boundary element methods for fluid-solid interaction eigenvalue problems , SIAM J. Numer. Anal. 52(5) (2014), 2400– 2414

    work page 2014

  25. [25]

    Kleefeld, A numerical method to compute interior transmission eigenvalues , Inverse Problems 29(10) (2013), 104012, 20

    A. Kleefeld, A numerical method to compute interior transmission eigenvalues , Inverse Problems 29(10) (2013), 104012, 20

    work page 2013

  26. [26]

    Kozlov and V

    V. Kozlov and V. Maz 1ya, Differential Equations with Operator Coefficients with Applica- tions to Boundary Value Problems for Partial Differential Equations. Springer-Verlag, Berlin (1999)

    work page 1999

  27. [27]

    R. de L. Kronig and W. Penney, Quantum mechanics of electrons in crystal lattices , Proc. Roy. Soc. Lond. 130 (1931), 499–513

    work page 1931

  28. [28]

    McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge Uni- versity Press, Cambridge (2000)

    W. McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge Uni- versity Press, Cambridge (2000)

    work page 2000

  29. [29]

    Mantile, A

    A. Mantile, A. Posilicano, and M. Sini, Self-adjoint elliptic operators with boundary condi- tions on not closed hypersurfaces , J. Differential Equations 261(1) (2016), 1–55

    work page 2016

  30. [30]

    Oˇ zanov´ a,Approximation by point potentials in a magnetic field , J

    K. Oˇ zanov´ a,Approximation by point potentials in a magnetic field , J. Phys. A 39 (2006), 3071–3083

    work page 2006

  31. [31]

    S. A. Sauter and C. Schwab, Boundary element methods. Springer-Verlag, Berlin (2011)

    work page 2011

  32. [32]

    ´Smigaj, S

    W. ´Smigaj, S. Arridge, T. Betcke , J. Phillips, and J. Schweiger, Solving Boundary Integral Problems with BEM++ , ACM Trans. Math. Software 41(6) (2015), 1–40

    work page 2015

  33. [33]

    Steinbach, Numerical approximation methods for elliptic boundary value problems

    O. Steinbach, Numerical approximation methods for elliptic boundary value problems. Finite and boundary elements. Springer, New York (2008)

    work page 2008

  34. [34]

    Steinbach and G

    0. Steinbach and G. Unger, Convergence analysis of a Galerkin boundary element method for the Dirichlet Laplacian eigenvalue problem , SIAM J. Numer. Anal. 50 (2012), 710–728

    work page 2012

  35. [35]

    L. H. Thomas, The interaction between a neutron and a proton and the structure of H 3, Phys. Rev., II. Ser. 47 (1935), 903–909

    work page 1935

  36. [36]

    Unger, Analysis of Boundary Element Methods for Laplacian Eigenvalue Problems

    G. Unger, Analysis of Boundary Element Methods for Laplacian Eigenvalue Problems. PhD thesis, Graz University of Technology, Graz, 2009. Institut f¨ur Angewandte Mathematik, Technische Universit¨at Graz, Steyrergasse 30, 8010 Graz, Austria, E-mail: holzmann@math.tugraz.at Institut f¨ur Angewandte Mathematik, Technische Universit¨at Graz, Steyrergasse 30, ...

    work page 2009