pith. sign in

arxiv: 2605.01230 · v1 · submitted 2026-05-02 · 🧮 math.NA · cs.NA

Inverse Scattering for Dirac Equations Arising in Waveguide Arrays

John C. Schotland , Shenwen Yu This is my paper

Pith reviewed 2026-05-09 18:48 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords inverse scatteringDirac equationswaveguide arraysinverse Born serieswell-posednessnumerical reconstructionerror estimates
0
0 comments X

The pith

Inverse Born series and its reduced form solve inverse scattering for Dirac equations in waveguide arrays with convergence guarantees.

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

The paper examines inverse scattering problems for Dirac equations that model waveguide arrays. It first shows that the forward scattering problem is well-posed. For the inverse problem, it introduces the inverse Born series and a reduced version to recover the potential from scattering data. The authors prove that these series converge under appropriate conditions and derive rigorous bounds on the reconstruction error. Numerical experiments illustrate that the algorithms work effectively in practice.

Core claim

For the inverse problems, the inverse Born series and the reduced inverse Born series converge, providing analysis of convergence and rigorous error estimates for recovering the coefficients in the Dirac equations arising from waveguide arrays.

What carries the argument

The inverse Born series, an iterative method for reconstructing the unknown potential in the Dirac operator from the scattering data.

If this is right

  • The forward scattering model for Dirac equations in waveguide arrays is well-posed.
  • The inverse Born series converges when the scattering data meet the required conditions.
  • Rigorous error estimates quantify the accuracy of the reconstructed potential.
  • Numerical tests confirm the practical utility of the method for continuum waveguide models.

Where Pith is reading between the lines

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

  • This framework might be adapted to inverse problems for other first-order hyperbolic systems.
  • Direct reconstruction via series could complement optimization methods in designing photonic devices based on waveguide arrays.
  • Error bounds may help assess stability in noisy measurement scenarios for similar inverse problems.

Load-bearing premise

The scattering data must satisfy conditions that ensure the convergence of the inverse Born series.

What would settle it

A case where the inverse Born series diverges or produces reconstructions with errors exceeding the derived bounds for a known Dirac potential in a waveguide array.

Figures

Figures reproduced from arXiv: 2605.01230 by John C. Schotland, Shenwen Yu.

Figure 1
Figure 1. Figure 1: Reconstructions of low contrast disk (chiral model) 10 view at source ↗
Figure 2
Figure 2. Figure 2: Reconstructions of medium contrast disk (chiral model) 11 view at source ↗
Figure 3
Figure 3. Figure 3: Reconstructions of high contrast disk (chiral model) 12 view at source ↗
Figure 4
Figure 4. Figure 4: Reconstructions of two low contrast Gaussian scatterers (chiral model) 13 view at source ↗
Figure 5
Figure 5. Figure 5: Reconstructions of two medium contrast Gaussian scatterers (chi￾ral model) 14 view at source ↗
Figure 6
Figure 6. Figure 6: Reconstructions of two high contrast Gaussian scatterers (chiral model) 15 view at source ↗
Figure 7
Figure 7. Figure 7: Reconstructions for the chiral model using the RIBS The numerical experiments for the chiral model made use of the computational domain Ω = [0, 2.0] × [−8.0, 8.0], with wave number k = 2.0. A total of 64 incident waves were generated, whose initial conditions are prescribed by ψ (j) (0, y) =   exp − (y − y (j) c ) 2 0.08 ! 0   , y ∈ [−8, 8], where the centers y (j) c are uniformly distributed in the … view at source ↗
Figure 8
Figure 8. Figure 8: Reconstructions of a low contrast disk (anti-chiral model) points. At the measurement location x = 2, 1600 detectors are uniformly placed along the transverse direction. To avoid the inverse crime, the synthetic data are generated on a finer grid with 800 time steps in the x-direction and 3200 spatial grid points in the y-direction. The Tikhonov regularization parameter is fixed as 0.001 throughout all exp… view at source ↗
Figure 9
Figure 9. Figure 9: Reconstructions of a medium contrast disk (anti-chiral model) We investigate the reconstruction of two types of scattering potentials with low, medium, and high contrast. The first is a disk defined by V (x, y) = ( σ, (x − 1.0)2 + y 2 ≤ (0.4)2 , 0, otherwise, where the contrast parameter σ is chosen as 0.1, 0.5, and 2.0. The corresponding numerical reconstructions are presented in Figures 1, 2, and 3. 18 view at source ↗
Figure 10
Figure 10. Figure 10: Reconstructions of a high contrast disk (anti-chiral model) The second potential consists of two Gaussian components, V (x, y) = σ  exp − (x − 0.5)2 + y 2 0.08  + exp − (x − 1.5)2 + y 2 0.08  , where the parameter σ is again selected as 0.1, 0.5, and 2.0. The corresponding results are shown in Figures 4, 5, and 6. In addition, we apply the RIBS with 10 terms with σ = 0.5. The results are displayed in view at source ↗
Figure 11
Figure 11. Figure 11: Reconstructions of two low contrast Gaussians (anti-chiral model) Two classes of scattering potentials are considered. The first is a circular disk given by V (x, y) = ( σ, (x − 12.8)2 + (y − 12.8)2 ≤ (6.4)2 , 0, otherwise, where the contrast σ is set to 0.1, 0.2, and 0.4. The corresponding reconstructions are shown in Figures 8, 9, and 10. 20 view at source ↗
Figure 12
Figure 12. Figure 12: Reconstructions of two medium contrast Gaussian scatterers (anti-chiral model) The second consists of two Gaussian-shaped scatterers, V (x, y) = σ " exp − (x − 16.64)2 + (y − 8.96)2 20.48  + exp − (x − 8.96)2 + (y − 16.64)2 10.83  # , where the intensity σ is chosen as 0.1, 0.3, and 0.6. The corresponding numerical results are presented in Figures 11, 12, and 13. Finally, we test the RIBS algorithm fo… view at source ↗
Figure 13
Figure 13. Figure 13: Reconstructions of two high contrast Gaussian scatterers (anti￾chiral model) For both models, the numerical results exhibit the same qualitative behavior. In the low￾contrast regime, one or two IBS terms are sufficient to accurately recover the scattering potential. As the contrast increases to a moderate level, a larger number of IBS terms is required for convergence, leading to improved reconstruction q… view at source ↗
Figure 14
Figure 14. Figure 14: Reconstructions of two medium contrast Gaussian scatterers for the anti-chiral model using RIBS 6. Discussion In conclusion, we have investigated the IBS and the RIBS for two Dirac equations arising from the quantum optics of chiral and anti-chiral quantum waveguide arrays. We established the convergence properties of the proposed series and conducted extensive numerical exper￾iments to demonstrate the ef… view at source ↗
read the original abstract

We investigate inverse scattering problems for Dirac equations that arise as continuum models of waveguide arrays. We first establish the well-posedness of the forward models. For the associated inverse problems, we develop the inverse Born series and the reduced inverse Born series, providing analysis of convergence and rigorous error estimates. Numerical experiments are presented to validate the proposed algorithms and demonstrate their effectiveness.

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 / 1 minor

Summary. The manuscript investigates inverse scattering problems for Dirac equations that arise as continuum models of waveguide arrays. It first establishes well-posedness of the forward scattering models. For the inverse problems it develops both the inverse Born series and a reduced inverse Born series, supplies convergence analysis together with rigorous error estimates, and validates the methods through numerical experiments.

Significance. If the stated well-posedness, convergence, and error-estimate results hold, the work supplies a mathematically grounded reconstruction framework for potentials or coefficients in Dirac waveguide models. The combination of the standard and reduced inverse Born series, explicit error bounds, and numerical validation constitutes a concrete contribution to inverse scattering theory in this setting.

minor comments (1)
  1. The abstract states that convergence analysis and rigorous error estimates are provided, yet the precise hypotheses on the scattering data (e.g., smallness or regularity requirements) are not summarized; adding a single sentence listing the main assumptions would improve readability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript on inverse scattering for Dirac equations in waveguide arrays, including the recognition of the well-posedness results, convergence analysis, error estimates, and numerical validation. The recommendation for minor revision is noted.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper establishes well-posedness of the forward Dirac waveguide models, then constructs the inverse Born series and reduced inverse Born series with convergence analysis and error estimates, followed by numerical validation. No load-bearing step reduces by the paper's own equations to a fitted input, self-definition, or self-citation chain. The inverse series are derived from the forward operator under standard conditions on scattering data, with independent mathematical content and external numerical checks. This is the normal case of a self-contained analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the work rests on standard well-posedness assumptions for Dirac scattering operators that are not detailed here.

pith-pipeline@v0.9.0 · 5342 in / 1014 out tokens · 25619 ms · 2026-05-09T18:48:21.647822+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. Reconstruction methods for inverse scattering problems with phaseless data

    math-ph 2026-05 unverdicted novelty 6.0

    Extensions of the inverse Born series with Fourier and polarization techniques are proposed to reconstruct scattering potentials from phaseless data.

Reference graph

Works this paper leans on

15 extracted references · 15 canonical work pages · cited by 1 Pith paper

  1. [1]

    J. G. Hoskins, M. Rachh, J. C. Schotland, Quantum electrodynamics of chiral and antichiral waveguide arrays, Optics Letters 48 (5) (2023) 1232–1235.doi:10.1364/OL.477807

    work page doi:10.1364/ol.477807 2023

  2. [2]

    Moskow, J

    S. Moskow, J. C. Schotland, Convergence and stability of the inverse scattering series for diffuse waves, Inverse Problems 24 (6) (2008) 065005.doi:10.1088/0266-5611/24/6/065005

    work page doi:10.1088/0266-5611/24/6/065005 2008

  3. [3]

    J. G. Hoskins, J. C. Schotland, Analysis of the inverse Born series: An approach through geometric function theory, Inverse Problems 38 (7) (2022) 074001.doi:10.1088/1361-6420/ac661f

    work page doi:10.1088/1361-6420/ac661f 2022

  4. [4]

    V. A. Markel, J. C. Schotland, Reduced inverse Born series: A computational study, JOSA A 39 (12) (2022) C179–C189.doi:10.1364/JOSAA.473683

    work page doi:10.1364/josaa.473683 2022

  5. [5]

    B. A. Watson, Inverse spectral problems for weighted Dirac systems, Inverse Problems 15 (3) (1999) 793.doi:10.1088/0266-5611/15/3/311

    work page doi:10.1088/0266-5611/15/3/311 1999

  6. [6]

    Ya. V. Mykytyuk, D. V. Puyda, Inverse spectral problems for Dirac operators on a finite interval, Journal of Mathematical Analysis and Applications 386 (1) (2012) 177–194.doi:10.1016/j.jmaa. 2011.07.061

    work page doi:10.1016/j.jmaa 2012

  7. [7]

    M. Salo, L. Tzou, Inverse problems with partial data for a Dirac system: A Carleman estimate approach, Advances in Mathematics 225 (1) (2010) 487–513.doi:10.1016/j.aim.2010.03.003. 25

    work page doi:10.1016/j.aim.2010.03.003 2010

  8. [8]

    Kurylev, M

    Y. Kurylev, M. Lassas, Inverse problems and index formulae for Dirac operators, Advances in Mathe- matics 221 (1) (2009) 170–216.doi:10.1016/j.aim.2008.12.001

    work page doi:10.1016/j.aim.2008.12.001 2009

  9. [9]

    Hachem, The partial-bar approach to inverse scattering for Dirac operators, Inverse Problems 11 (1) (1995) 123.doi:10.1088/0266-5611/11/1/007

    G. Hachem, The partial-bar approach to inverse scattering for Dirac operators, Inverse Problems 11 (1) (1995) 123.doi:10.1088/0266-5611/11/1/007

    work page doi:10.1088/0266-5611/11/1/007 1995

  10. [10]

    Isozaki, Inverse scattering theory for Dirac operators, Annales de l’I.H.P

    H. Isozaki, Inverse scattering theory for Dirac operators, Annales de l’I.H.P. Physique th´ eorique 66 (2) (1997) 237–270

    work page 1997

  11. [11]

    Kawamoto, M

    A. Kawamoto, M. Yamamoto, Determination of an electromagnetic potential for the Dirac equation, Inverse Problems 28 (11) (2012) 115012.doi:10.1088/0266-5611/28/11/115012

    work page doi:10.1088/0266-5611/28/11/115012 2012

  12. [12]

    Grebert, Inverse scattering for the Dirac operator on the real line, Inverse Problems 8 (5) (1992) 787

    B. Grebert, Inverse scattering for the Dirac operator on the real line, Inverse Problems 8 (5) (1992) 787. doi:10.1088/0266-5611/8/5/007

    work page doi:10.1088/0266-5611/8/5/007 1992

  13. [13]

    M. I. Belishev, V. S. Mikhailov, Inverse problem for a one-dimensional dynamical Dirac system (BC- method), Inverse Problems 30 (12) (2014) 125013.doi:10.1088/0266-5611/30/12/125013

    work page doi:10.1088/0266-5611/30/12/125013 2014

  14. [14]

    Yi, Bounded time inverse scattering for semilinear Dirac equation, Inverse Problems 41 (6) (2025) 065024.doi:10.1088/1361-6420/ade471

    Y. Yi, Bounded time inverse scattering for semilinear Dirac equation, Inverse Problems 41 (6) (2025) 065024.doi:10.1088/1361-6420/ade471

    work page doi:10.1088/1361-6420/ade471 2025

  15. [15]

    Sharp Bounds for the Modulus and Phase of

    P. Freitas, Sharp bounds for the modulus and phase of Hankel functions with applications to Jaeger integrals, Mathematics of Computation 87 (309) (2018) 289–308.doi:10.1090/mcom/3267. Department of Mathematics and Department of Physics, Yale University, New Haven, CT, USA Email address:john.schotland@yale.edu Department of Mathematical Sciences, Tsinghua ...

    work page doi:10.1090/mcom/3267 2018