pith. sign in

arxiv: 2312.16701 · v1 · submitted 2023-12-27 · 🧮 math-ph · cs.NA· math.AP· math.MP· math.NA

Integral formulation of Dirac singular waveguides

Guillaume Bal , Jeremy Hoskins , Solomon Quinn , Manas Rachh This is my paper

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

classification 🧮 math-ph cs.NAmath.APmath.MPmath.NA
keywords Dirac equationboundary integral equationsurface wavesmass jumpholomorphic perturbation theorywaveguidesscattering
0
0 comments X

The pith

A boundary integral equation for the two-dimensional massive Dirac equation with a mass jump has a unique solution for almost all parameters.

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

The paper derives a boundary integral formulation for the massive Dirac equation in two dimensions when the mass term jumps across a one-dimensional interface. This setup models transitions between insulating materials that support surface waves traveling along the interface while decaying exponentially away from it. The authors prove that the resulting integral equation admits a unique solution for almost all choices of parameters by applying holomorphic perturbation theory. They extend the same approach to the case of two such interfaces and implement a fast numerical solver that computes solutions and scattering patterns. The reduction makes it possible to study these waveguide problems without discretizing the full two-dimensional domain.

Core claim

After deriving the boundary integral equation from the Dirac operator with a discontinuous mass, the paper establishes that this equation has a unique solution for almost all parameter values by means of holomorphic perturbation theory. The same uniqueness result holds when the formulation is extended to two parallel interfaces. A fast numerical method based on the integral equation is then used to compute explicit solutions and to illustrate scattering effects.

What carries the argument

Boundary integral equation obtained by reducing the Dirac equation across a mass jump, with uniqueness proved via holomorphic perturbation theory.

If this is right

  • Scattering and propagation along the interface can be computed by solving only a one-dimensional integral equation.
  • The same reduction and uniqueness proof apply when two separate interfaces are present.
  • Numerical examples demonstrate that the method captures both guided modes and scattering at defects or junctions.
  • Parameter studies become feasible because the integral formulation avoids repeated full-domain discretizations.

Where Pith is reading between the lines

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

  • The holomorphic perturbation technique used for uniqueness might extend to other first-order systems with singular coefficients.
  • The boundary reduction could be tested on related models such as the Schrödinger equation with discontinuous potentials.
  • Fast solvers of this type open the possibility of inverse problems that recover interface shapes from observed scattering data.

Load-bearing premise

The surface waves induced by the mass jump must decay exponentially in the transverse direction, which permits the reduction from the full PDE to a boundary integral equation.

What would settle it

An explicit set of parameters for which the integral equation either has no solution or has more than one solution would show that the uniqueness claim fails.

Figures

Figures reproduced from arXiv: 2312.16701 by Guillaume Bal, Jeremy Hoskins, Manas Rachh, Solomon Quinn.

Figure 1
Figure 1. Figure 1: Solution for the flat interface, Γ = {𝑥2 = 0}, where 𝑚 = 1, 𝐸 = 0.5 and the right-hand side is 𝑓1 = (0, 0), 𝑓2 = (𝛿𝑥0 , 0) with source location 𝑥0 = (0, 2). The full solution is illustrated by the second and third rows. The top left panel plots the real parts of the computed densities 𝜌 and 𝜏. Recall that 𝜏1 = 𝜌1. The top right panel contains the relative error of the solution 𝑢1 at the four specified poin… view at source ↗
Figure 2
Figure 2. Figure 2: A component of the solution 𝑢 of the Dirac equation (6) for various interfaces Γ. In each plot, the interface is illustrated by the solid curve. These examples all correspond to 𝑚 = 4, 𝐸 = 1 and a right-hand side given by 𝑓1 = (0, 0), 𝑓2 = (𝛿𝑥0 , 0), with source location of 𝑥0 = (−1, 1). Dirac equation (6) with energy 𝐸. But since solutions of the Dirac equation cannot get trapped (as discussed above), suc… view at source ↗
Figure 3
Figure 3. Figure 3: The solution 𝑢 is evaluated at the point (−10, 1), which lies just above the interface. As expected, the solution is smooth in 𝐸 with no resonance detected. The right panel shows the corresponding Klein-Gordon solution, which exhibits several sharp peaks. This direct comparison between the Dirac and Klein-Gordon solutions illustrates that there are resonances near the real axis for the Klein-Gordon equatio… view at source ↗
Figure 3
Figure 3. Figure 3: Top row: absolute value of (the first component of) the solution of the Dirac equation (6), with [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left: a component of the solution of the Dirac equation as a function of 0 [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Top row and bottom left: plots of ℜ(𝑢1) for the two-interface problem with various interface separations 𝑑. Bottom right panel: The transmission coefficient 𝑇𝑅 as a function of 𝑑. In these examples, 𝑚 = 2, 𝐸 = 0.8, and the right-hand side is given by 𝑓0 = (𝛿𝑥0 , 0), 𝑓1 = (0, 0), 𝑓2 = (0, 0) with source location 𝑥0 at the center of the red dot (at a fixed distance away from Γ1). where we assume that 0 < 𝐸2 … view at source ↗
Figure 6
Figure 6. Figure 6: Same set-up as in Figure 5, only with different interfaces. [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The 𝑑 → 0 limit of the transmission coefficient as a function of 𝜃, for interfaces Γ𝑗 parametrized by (38) with 𝑎 = tan 𝜃. and 𝐺 𝑗(𝑥, 𝑦) := 1 2𝜋 (−𝑖𝜕𝑥1𝜎1 − 𝑖𝜕𝑥2𝜎2 + (−1) 𝑗𝑚 𝑗𝜎3 + 𝐸)𝐾0 (𝜔𝑗 |𝑥 − 𝑦|), with 𝜔𝑗 := √︃ 𝑚2 𝑗 − 𝐸2 . We again wish to solve for 𝜇 such that 𝑢 is continuous across Γ. Recalling the shorthand ˆ𝑛(𝑡) · 𝜎 := 𝑛ˆ1 (𝑡)𝜎1 + 𝑛ˆ2 (𝑡)𝜎2, standard properties of layer potentials imply that [ [𝑢𝑠]] (… view at source ↗
read the original abstract

This paper concerns a boundary integral formulation for the two-dimensional massive Dirac equation. The mass term is assumed to jump across a one-dimensional interface, which models a transition between two insulating materials. This jump induces surface waves that propagate outward along the interface but decay exponentially in the transverse direction. After providing a derivation of our integral equation, we prove that it has a unique solution for almost all choices of parameters using holomorphic perturbation theory. We then extend these results to a Dirac equation with two interfaces. Finally, we implement a fast numerical method for solving our boundary integral equations and present several numerical examples of solutions and scattering effects.

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

0 major / 3 minor

Summary. The manuscript derives a boundary integral equation for the 2D massive Dirac operator with a mass term that jumps across a smooth curve (modeling an interface between insulating regions). It establishes that the resulting integral operator is Fredholm and applies holomorphic perturbation theory to prove that the equation admits a unique solution for all but a discrete (or measure-zero) set of parameters. The analysis is extended to the case of two parallel interfaces, and a Nyström-type discretization is implemented to compute surface-wave solutions and scattering examples.

Significance. If the derivation and uniqueness argument hold, the work supplies a mathematically justified boundary-integral framework for surface modes of gapped Dirac operators, a setting of current interest in mathematical physics and condensed-matter modeling. The explicit use of holomorphic perturbation theory on trace spaces to obtain uniqueness outside an exceptional set is a clear technical strength, as is the extension to multiple interfaces and the accompanying numerical illustrations.

minor comments (3)
  1. [§2.2] §2.2, after Eq. (2.7): the precise definition of the trace space H^{1/2}(Γ) and the associated single-layer operator should include the explicit norm and the mapping properties used in the subsequent Fredholm argument.
  2. [§4] §4, Theorem 4.1: the statement that uniqueness holds for 'almost all' parameters would be strengthened by a brief remark on whether the exceptional set is discrete or merely of measure zero, together with the dependence on the curve geometry.
  3. [§5.3] Figure 5 and §5.3: the convergence plot for the two-interface case reports relative L² errors but does not indicate the reference solution or the quadrature order used; adding this information would make the numerical validation self-contained.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the thorough summary and positive assessment of our work on the boundary integral formulation for the massive Dirac equation with discontinuous mass. The recommendation of minor revision is appreciated. No specific major comments were provided in the report, so we have no individual points requiring detailed rebuttal or revision at this stage.

Circularity Check

0 steps flagged

No significant circularity; derivation and uniqueness proof are independent

full rationale

The paper presents a derivation of a boundary integral equation for the 2D massive Dirac operator with a mass jump, based on the standard property of exponential decay of surface modes away from the essential spectrum. Uniqueness for almost all parameters is obtained by applying holomorphic perturbation theory to the Fredholm family of integral operators on trace spaces; this is a standard technique that does not reduce the result to any fitted quantity, self-definition, or prior self-citation chain. The extension to two interfaces and the numerical solver are likewise independent of the claimed results. No load-bearing self-citations, ansatzes, or renamings of known results are present in the central claims.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central modeling choice is the discontinuous mass across a 1D interface; no free parameters, invented entities, or additional ad-hoc axioms are stated in the abstract.

axioms (1)
  • domain assumption Mass term jumps across a one-dimensional interface modeling transition between two insulating materials
    Core physical modeling assumption that induces the surface waves and enables the boundary reduction.

pith-pipeline@v0.9.0 · 5633 in / 1173 out tokens · 27730 ms · 2026-05-24T05:17:36.359117+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages · 1 internal anchor

  1. [1]

    B. A. Bernevig, T. L. Hughes, Topological insulators and topological superconductors, Prince- ton university press, 2013

    work page 2013

  2. [2]

    Fefferman, M

    C. Fefferman, M. Weinstein, Honeycomb lattice potentials and Dirac points, Journal of the American Mathematical Society 25 (4) (2012) 1169–1220. 29

    work page 2012

  3. [3]

    Helsing, A

    J. Helsing, A. Ros ´en, Dirac integral equations for dielectric and plasmonic scattering, Integral Equations and Operator Theory 93 (5) (2021) 1–41

    work page 2021

  4. [4]

    Witten, Three lectures on topological phases of matter, La Rivista del Nuovo Cimento 39 (7) (2016) 313–370

    E. Witten, Three lectures on topological phases of matter, La Rivista del Nuovo Cimento 39 (7) (2016) 313–370

    work page 2016

  5. [5]

    Schulz-Baldes, J

    H. Schulz-Baldes, J. Kellendonk, T. Richter, Simultaneous quantization of edge and bulk hall conductivity, Journal of Physics A: Mathematical and General 33 (2) (2000) L27

    work page 2000

  6. [6]

    Elbau, G.-M

    P. Elbau, G.-M. Graf, Equality of bulk and edge hall conductance revisited, Communications in mathematical physics 229 (3) (2002) 415–432

    work page 2002

  7. [7]

    G. M. Graf, M. Porta, Bulk-edge correspondence for two-dimensional topological insulators, Communications in Mathematical Physics 324 (3) (2013) 851–895

    work page 2013

  8. [8]

    Prodan, H

    E. Prodan, H. Schulz-Baldes, Bulk and boundary invariants for complex topological insulators: from 𝑘-theory to physics, Mathematical Physics Studies (2016)

    work page 2016

  9. [9]

    Drouot, Microlocal analysis of the bulk-edge correspondence, Communications in Mathe- matical Physics 383 (3) (2021) 2069–2112

    A. Drouot, Microlocal analysis of the bulk-edge correspondence, Communications in Mathe- matical Physics 383 (3) (2021) 2069–2112

    work page 2021

  10. [10]

    Bal, Topological invariants for interface modes, Communications in Partial Differential Equations 47 (8) (2022) 1636–1679

    G. Bal, Topological invariants for interface modes, Communications in Partial Differential Equations 47 (8) (2022) 1636–1679

    work page 2022

  11. [11]

    Bal, Topological charge conservation for continuous insulators, Journal of Mathematical Physics 64 (3) (2023)

    G. Bal, Topological charge conservation for continuous insulators, Journal of Mathematical Physics 64 (3) (2023)

    work page 2023

  12. [12]

    Quinn, G

    S. Quinn, G. Bal, Approximations of interface topological invariants, arXiv preprint arXiv:2112.02686 (2021)

    work page arXiv 2021

  13. [13]

    Hatsugai, Edge states in the integer quantum hall effect and the riemann surface of the bloch function, Physical Review B 48 (16) (1993) 11851

    Y. Hatsugai, Edge states in the integer quantum hall effect and the riemann surface of the bloch function, Physical Review B 48 (16) (1993) 11851

    work page 1993

  14. [14]

    G. E. Volovik, The universe in a helium droplet, Vol. 117, Oxford University Press on Demand, 2003

    work page 2003

  15. [15]

    A. M. Essin, V. Gurarie, Bulk-boundary correspondence of topological insulators from their respective green’s functions, Physical Review B 84 (12) (2011) 125132

    work page 2011

  16. [16]

    Fukui, K

    T. Fukui, K. Shiozaki, T. Fujiwara, S. Fujimoto, Bulk-edge correspondence for chern topolog- ical phases: A viewpoint from a generalized index theorem, Journal of the Physical Society of Japan 81 (11) (2012) 114602

    work page 2012

  17. [17]

    Bal, Continuous bulk and interface description of topological insulators, Journal of Mathe- matical Physics 60 (8) (2019) 081506

    G. Bal, Continuous bulk and interface description of topological insulators, Journal of Mathe- matical Physics 60 (8) (2019) 081506

    work page 2019

  18. [18]

    Thaller, The Dirac equation, Springer Science & Business Media, 2013

    B. Thaller, The Dirac equation, Springer Science & Business Media, 2013

    work page 2013

  19. [19]

    G. Bal, J. Hoskins, S. Quinn, M. Rachh, Integral formulation of Klein-Gordon singular waveg- uides, arXiv preprint arXiv:2212.12619 (2022)

    work page arXiv 2022

  20. [20]

    Raether, Surface plasmons on smooth surfaces, Surface plasmons on smooth and rough surfaces and on gratings (1988) 4–39

    H. Raether, Surface plasmons on smooth surfaces, Surface plasmons on smooth and rough surfaces and on gratings (1988) 4–39. 30

    work page 1988

  21. [21]

    Tzarouchis, A

    D. Tzarouchis, A. Sihvola, Light scattering by a dielectric sphere: Perspectives on the mie resonances, Applied Sciences 8 (2) (2018) 184

    work page 2018

  22. [22]

    Helsing, K.-M

    J. Helsing, K.-M. Perfekt, The spectra of harmonic layer potential operators on domains with rotationally symmetric conical points, Journal de Math ´ematiques Pures et Appliqu ´ees 118 (2018) 235–287

    work page 2018

  23. [23]

    Helsing, A

    J. Helsing, A. Karlsson, An extended charge-current formulation of the electromagnetic trans- mission problem, SIAM Journal on Applied Mathematics 80 (2) (2020) 951–976

    work page 2020

  24. [24]

    Darbas, E

    M. Darbas, E. Darrigrand, Y. Lafranche, Combining analytic preconditioner and fast multipole method for the 3-d helmholtz equation, Journal of Computational Physics 236 (2013) 289–316

    work page 2013

  25. [25]

    Antoine, Advances in the on-surface radiation condition method: Theory, numerics and applications, Computational Methods for Acoustics Problems (2008) 169–194

    X. Antoine, Advances in the on-surface radiation condition method: Theory, numerics and applications, Computational Methods for Acoustics Problems (2008) 169–194

    work page 2008

  26. [26]

    Antoine, M

    X. Antoine, M. Darbas, Alternative integral equations for the iterative solution of acoustic scattering problems, The Quarterly Journal of Mechanics and Applied Mathematics 58 (1) (2005) 107–128

    work page 2005

  27. [27]

    Antoine, M

    X. Antoine, M. Darbas, Y. Y. Lu, An improved surface radiation condition for high-frequency acoustic scattering problems, Computer Methods in Applied Mechanics and Engineering 195 (33-36) (2006) 4060–4074

    work page 2006

  28. [28]

    Chaillat, M

    S. Chaillat, M. Darbas, F. Le Lou ¨er, Approximate local dirichlet-to-neumann map for three- dimensional time-harmonic elastic waves, Computer Methods in Applied Mechanics and En- gineering 297 (2015) 62–83

    work page 2015

  29. [29]

    Chaillat, M

    S. Chaillat, M. Darbas, F. Le Lou ¨er, A new analytic preconditioner for the iterative solution of dirichlet exterior scattering problems in 3d elasticity (2014)

    work page 2014

  30. [30]

    G. A. Kriegsmann, A. Taflove, K. R. Umashankar, A new formulation of electromagnetic wave scattering using an on-surface radiation boundary condition approach, IEEE Transactions on Antennas and Propagation 35 (2) (1987) 153–161

    work page 1987

  31. [31]

    G. Bal, J. G. Hoskins, Z. Wang, Asymmetric transport computations in Dirac models of topological insulators, Journal of Computational Physics 487 (2023) 112151

    work page 2023

  32. [32]

    G. Bal, S. Becker, A. Drouot, C. F. Kammerer, J. Lu, A. B. Watson, Edge state dynamics along curved interfaces, SIAM Journal on Mathematical Analysis 55 (5) (2023) 4219–4254

    work page 2023

  33. [33]

    G. Bal, S. Becker, A. Drouot, Magnetic slowdown of topological edge states, To appear in C.P.A.M.; arXiv:2201.07133 (2023)

    work page arXiv 2023

  34. [34]

    Drouot, Topological insulators in semiclassical regime, arXiv preprint arXiv:2206.08238 (2022)

    A. Drouot, Topological insulators in semiclassical regime, arXiv preprint arXiv:2206.08238 (2022)

    work page arXiv 2022

  35. [35]

    Bal, Semiclassical propagation along curved domain walls, To appear in SIAM M.M.S.; arXiv:2206.09439 (2022)

    G. Bal, Semiclassical propagation along curved domain walls, To appear in SIAM M.M.S.; arXiv:2206.09439 (2022)

    work page arXiv 2022

  36. [36]

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

    M. Holzmann, G. Unger, Boundary integral formulations of eigenvalue problems for elliptic differential operators with singular interactions and their numerical approximation by boundary element methods, arXiv preprint arXiv:1907.04282 (2019). 31

    work page internal anchor Pith review Pith/arXiv arXiv 1907

  37. [37]

    Abramowitz, I

    M. Abramowitz, I. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover, 1964

    work page 1964

  38. [38]

    chunkie, https://github.com/fastalgorithms/chunkie, GitHub Repository

    work page

  39. [39]

    Bremer, V

    J. Bremer, V. Rokhlin, I. Sammis, Universal quadratures for boundary integral equations on two-dimensional domains with corners, Journal of Computational Physics 229 (22) (2010) 8259–8280

    work page 2010

  40. [40]

    Bremer, Z

    J. Bremer, Z. Gimbutas, V. Rokhlin, A nonlinear optimization procedure for generalized gaussian quadratures, SIAM Journal on Scientific Computing 32 (4) (2010) 1761–1788

    work page 2010

  41. [41]

    Rokhlin, Rapid solution of integral equations of scattering theory in two dimensions, Journal of Computational physics 86 (2) (1990) 414–439

    V. Rokhlin, Rapid solution of integral equations of scattering theory in two dimensions, Journal of Computational physics 86 (2) (1990) 414–439

    work page 1990

  42. [42]

    Greengard, V

    L. Greengard, V. Rokhlin, A fast algorithm for particle simulations, Journal of computational physics 73 (2) (1987) 325–348

    work page 1987

  43. [43]

    fmm2d, https://github.com/flatironinstitute/fmm2d, GitHub Repository

    work page

  44. [44]

    Jiang, L

    S. Jiang, L. Greengard, Approximating the Gaussian as a Sum of Exponentials and its Appli- cations to the Fast Gauss Transform, Communications in Computational Physics 31 (1) (2022) 1–26

    work page 2022

  45. [45]

    Hammer, W

    R. Hammer, W. P ¨otz, Dynamics of domain-wall dirac fermions on a topological insulator: A chiral fermion beam splitter, Physical Review B 88 (23) (2013) 235119

    work page 2013

  46. [46]

    Z. Qiao, J. Jung, C. Lin, Y. Ren, A. H. MacDonald, Q. Niu, Current partition at topological channel intersections, Physical Review Letters 112 (20) (2014) 206601

    work page 2014

  47. [47]

    G. Bal, P. Cazeaux, D. Massatt, S. Quinn, Mathematical models of topologically protected transport in twisted bilayer graphene, Multiscale Modeling & Simulation 21 (3) (2023) 1081– 1121

    work page 2023

  48. [48]

    Kato, Perturbation theory for linear operators, Vol

    T. Kato, Perturbation theory for linear operators, Vol. 132, Springer Science & Business Media, 2013

    work page 2013

  49. [49]

    Colton, R

    D. Colton, R. Kress, Integral equation methods in scattering theory, SIAM, 2013. 32

    work page 2013