pith. sign in

arxiv: 2605.01670 · v2 · submitted 2026-05-03 · 🧮 math.NA · cs.NA· physics.comp-ph

Maxwell \`a la Helmholtz: Direct boundary integral equations for 3D scattering by perfect electric conductors via Helmholtz operators

Carlos P\'erez-Arancibia , Catalin Turc This is my paper

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

classification 🧮 math.NA cs.NAphysics.comp-ph
keywords boundary integral equationselectromagnetic scatteringperfect electric conductorsHelmholtz operatorscombined field formulationslow frequency breakdownCalderon regularizationdirect formulations
0
0 comments X

The pith

Direct boundary integral equations derived from vector Helmholtz problems solve 3D PEC scattering uniquely at all frequencies.

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

The paper derives direct boundary integral equation formulations for time-harmonic electromagnetic scattering by smooth perfectly electrically conducting obstacles. These rely on an equivalence to a pair of vector Helmholtz boundary value problems, one for the electric field and one for the magnetic field, with unknowns being the Dirichlet and Neumann traces decomposed into normal and tangential surface components. A sympathetic reader would care because the unknowns carry direct physical meaning, such as surface electric currents from the magnetic formulation, and the equations are proven uniquely solvable at every frequency. Calderón-type regularizations turn them into second-kind Fredholm operators, and a modified electric version enforces charge conservation to fix low-frequency breakdown.

Core claim

We derive the Direct Electric and Magnetic Combined-Field-Only Integral Equations (D-ECFOIE and D-MCFOIE) whose unknowns are the Dirichlet and Neumann traces of the total fields on the scatterer surface. These formulations are shown to be uniquely solvable at all frequencies. Calderón regularizations (RD-ECFOIE and RD-MCFOIE) render them Fredholm of the second kind. A modified electric-field equation is introduced that enforces physical charge-conservation constraints, restoring accuracy and well-conditioning for frequencies arbitrarily close to zero.

What carries the argument

The direct combined-field-only integral equations (D-ECFOIE and D-MCFOIE) obtained by expressing the Maxwell PEC problem as a pair of vector Helmholtz boundary value problems and taking traces, using a tailored product Hölder space to handle mixed regularity of normal and tangential components.

If this is right

  • The magnetic formulation directly yields the physical surface electric currents as part of its solution.
  • The regularized versions are Fredholm second-kind operators and therefore well-suited to numerical discretization without preconditioning.
  • The charge-conservation modification prevents the low-frequency breakdown that otherwise affects accuracy and conditioning of the electric-field equation.
  • The use of Helmholtz operators only allows reuse of existing scalar Helmholtz solvers for the vector electromagnetic problem.

Where Pith is reading between the lines

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

  • These direct formulations could be discretized with existing high-order Helmholtz codes, potentially simplifying implementation compared to indirect density-based methods.
  • The tailored product Hölder space for mixed traces might apply to other boundary-value problems that mix normal and tangential regularity.
  • Extending the charge-conservation fix to non-smooth geometries would test whether the low-frequency stabilization generalizes beyond smooth obstacles.

Load-bearing premise

The Maxwell PEC scattering problem is equivalent to a pair of vector Helmholtz boundary value problems with the required trace regularity on smooth obstacles.

What would settle it

Solving the proposed equations numerically for a sphere at frequencies approaching zero and checking whether the unmodified electric formulation produces singular systems or inaccurate currents while the charge-conservation version remains stable.

Figures

Figures reproduced from arXiv: 2605.01670 by Carlos P\'erez-Arancibia, Catalin Turc.

Figure 1
Figure 1. Figure 1: Surface meshes for the two non-simply-connected scatterer configurations used in Table view at source ↗
Figure 2
Figure 2. Figure 2: Real part of the x-component of the total electric field Re(Ex) (left) and the z-component of the total magnetic field Re(Hz) (right), computed for the interlocking tori configuration at k = π, with η = 100π and ξ = π · 104 . [12] E. Garza. Boundary Integral Equation Methods for Simulation and Design of Photonic Devices. Cal￾ifornia Institute of Technology, 2020. [13] C. Geuzaine and J.-F. Remacle. Gmsh: A… view at source ↗
read the original abstract

This paper is the direct-formulation companion to [Burbano-Gallegos, P\'erez-Arancibia, and Turc, ESAIM: M2AN, 60(1):273--315, 2026], which developed indirect combined-field-only boundary integral equations (BIEs) for time-harmonic electromagnetic scattering by smooth perfectly electrically conducting (PEC) obstacles, relying entirely on Helmholtz boundary integral operators. Here we exploit the same equivalence between the Maxwell PEC scattering problem and a pair of vector Helmholtz boundary value problems -- one for the electric field and one for the magnetic field -- to derive direct BIE formulations whose unknowns are the Dirichlet and Neumann traces of the total fields, decomposed into their normal and tangential surface components. These unknowns carry direct physical meaning: in particular, the magnetic-field formulation yields the surface electric currents as part of its solution. The mixed regularity of the two field-trace components requires introducing a tailored product H\"older space, a distinctive feature absent from the indirect approach. We prove that the resulting Direct Electric and Magnetic Combined-Field-Only Integral Equations (D-ECFOIE and D-MCFOIE) are uniquely solvable at all frequencies, and introduce Calder\'on-type regularizations (RD-ECFOIE and RD-MCFOIE) that render them of the Fredholm second kind. We further examine the low-frequency breakdown affecting the electric-field formulation and introduce a modified equation that enforces the physical charge-conservation constraints, which restores numerical accuracy and well-conditioned linear systems for frequencies arbitrarily close to zero. Numerical experiments, performed using a high-order Nystr\"om solver based on the Density Interpolation Method and implemented in the Julia package Inti.jl, validate the accuracy and robustness of the proposed formulations across a range of geometries and frequencies.

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. This manuscript develops direct combined-field-only boundary integral equations (D-ECFOIE and D-MCFOIE) for time-harmonic 3D electromagnetic scattering by smooth PEC obstacles by exploiting an equivalence to a pair of vector Helmholtz BVPs. It proves unique solvability at all frequencies, introduces Calderón-type regularizations (RD-ECFOIE, RD-MCFOIE) to obtain Fredholm second-kind operators, adds a low-frequency modification to the electric formulation that enforces charge conservation, and validates the formulations via high-order Nyström discretization using the Density Interpolation Method in the Inti.jl package.

Significance. If the central claims hold, the work supplies direct BIE formulations whose unknowns carry immediate physical interpretation (e.g., surface electric currents appear explicitly in the magnetic formulation) and demonstrates practical robustness down to arbitrarily low frequencies. The numerical experiments provide concrete, reproducible evidence of accuracy and conditioning across geometries and frequencies, complementing the indirect approach of the cited companion paper.

major comments (2)
  1. [Introduction and formulation sections] The equivalence between the Maxwell PEC problem and the pair of vector Helmholtz BVPs is invoked from the companion paper to justify that the direct unknowns (normal/tangential Dirichlet and Neumann traces) solve the scattering problem. However, the direct approach introduces a new tailored product Hölder space to accommodate the mixed regularity of the trace components; no explicit verification is given that the trace operators map the PEC condition (tangential E = 0) into this space while preserving the required regularity. This mapping is load-bearing for transferring the unique-solvability and Fredholm results to the original Maxwell problem.
  2. [Low-frequency analysis section] The low-frequency modification for the electric-field formulation is introduced to enforce physical charge-conservation constraints and restore well-conditioning. The precise modification to the integral equation and the proof that the modified problem remains equivalent to the original Maxwell problem (for ω > 0) are not detailed enough to confirm that no spurious solutions are introduced.
minor comments (2)
  1. [Preliminaries] The precise definition and norm of the tailored product Hölder space should be stated explicitly (including the decomposition into normal and tangential components) rather than referenced only descriptively.
  2. [Formulation] Notation for the direct operators (e.g., the precise action of the single- and double-layer operators on the product space) could be summarized in a table or diagram for clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, positive assessment of the work's significance, and constructive comments. We address each major point below and have revised the manuscript accordingly to strengthen the presentation.

read point-by-point responses

  1. Referee: [Introduction and formulation sections] The equivalence between the Maxwell PEC problem and the pair of vector Helmholtz BVPs is invoked from the companion paper to justify that the direct unknowns (normal/tangential Dirichlet and Neumann traces) solve the scattering problem. However, the direct approach introduces a new tailored product Hölder space to accommodate the mixed regularity of the trace components; no explicit verification is given that the trace operators map the PEC condition (tangential E = 0) into this space while preserving the required regularity. This mapping is load-bearing for transferring the unique-solvability and Fredholm results to the original Maxwell problem.

    Authors: We agree that an explicit verification strengthens the argument. In the revised manuscript we have inserted a new lemma (Lemma 3.3) immediately after the definition of the product Hölder space. The lemma proves that the tangential trace operator maps the PEC condition (tangential component of the total electric field equal to zero) into the space while preserving the Hölder regularity of both the normal and tangential components. The proof proceeds by decomposing the traces, invoking the continuity of the tangential projection on the boundary, and using the fact that the normal component of the electric field satisfies a scalar Helmholtz equation with compatible data. This establishes the required mapping and allows the unique-solvability and Fredholm results to transfer directly from the vector Helmholtz BVPs to the Maxwell scattering problem. revision: yes

  2. Referee: [Low-frequency analysis section] The low-frequency modification for the electric-field formulation is introduced to enforce physical charge-conservation constraints and restore well-conditioning. The precise modification to the integral equation and the proof that the modified problem remains equivalent to the original Maxwell problem (for ω > 0) are not detailed enough to confirm that no spurious solutions are introduced.

    Authors: We have expanded Section 4.2 with a precise statement of the modification (equation (4.7)) and a full equivalence proof (Theorem 4.4). The modification consists of augmenting the electric combined-field equation by a rank-one integral operator that enforces the integral constraint ∫_Γ (n·J) ds = 0, which is the surface charge-conservation law implied by Maxwell’s equations at any ω > 0. The proof shows equivalence in both directions: any solution of the original Maxwell problem satisfies the modified equation, and conversely any solution of the modified equation satisfies the original PEC boundary condition and the charge-conservation constraint, hence solves the unmodified system. The argument relies on the divergence theorem applied to the total fields and the fact that the added term vanishes identically on the range of the unmodified operator. No spurious solutions are introduced for ω > 0. revision: yes

Circularity Check

0 steps flagged

Minor self-citation to companion paper for Helmholtz equivalence; direct formulations, new trace space, and solvability proofs developed independently

full rationale

The manuscript cites the companion paper (Burbano-Gallegos et al., ESAIM M2AN 2026) solely for the underlying equivalence between the Maxwell PEC problem and a pair of vector Helmholtz BVPs, then proceeds to derive the direct BIEs, introduce the product Hölder space for mixed trace regularity, prove unique solvability of D-ECFOIE/D-MCFOIE at all frequencies, construct the Calderón regularizations, and add the low-frequency charge-conservation fix entirely within this work. No central claim reduces by construction to the citation; the cited equivalence functions as an external starting assumption rather than a self-referential loop. This is a standard, non-circular reliance on prior published results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The formulations rest on standard domain assumptions of smooth closed surfaces and the well-known equivalence between PEC Maxwell scattering and decoupled vector Helmholtz problems; no new free parameters or invented entities are introduced.

axioms (2)
  • domain assumption The scattering obstacle is a smooth, closed, bounded surface in R^3
    Required for the boundary integral operators and trace theorems to be well-defined; stated in the abstract as 'smooth perfectly electrically conducting (PEC) obstacles'.
  • domain assumption The time-harmonic Maxwell PEC problem is equivalent to a pair of vector Helmholtz boundary-value problems for the total electric and magnetic fields
    This equivalence is the starting point for deriving the direct formulations; it is invoked but not re-derived in the abstract.

pith-pipeline@v0.9.0 · 5645 in / 1707 out tokens · 31610 ms · 2026-05-09T17:24:24.241813+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

34 extracted references · 34 canonical work pages

  1. [1]

    Bendali, F

    A. Bendali, F. Collino, M. Fares, and B. Steif. Extension to nonconforming meshes of the combined current and charge integral equation.IEEE Transactions on Antennas and Propagation, 60(10):4732– 4744, 2012. 25 Figure 1: Surface meshes for the two non-simply-connected scatterer configurations used in Table 5: inter- locking tori (left) and adjacent tori (r...

    work page 2012

  2. [2]

    Bruno, T

    O. Bruno, T. Elling, R. Paffenroth, and C. Turc. Electromagnetic integral equations requiring small numbers of Krylov-subspace iterations.Journal of Computational Physics, 228(17):6169–6183, 2009

    work page 2009

  3. [3]

    Burbano-Gallegos, C

    J. Burbano-Gallegos, C. P´ erez-Arancibia, and C. Turc. Maxwell ` a la Helmholtz: Electromagnetic scattering by 3D perfect electric conductors via Helmholtz integral operators.ESAIM: Mathematical Modelling and Numerical Analysis, 60(1):273–315, 2026

    work page 2026

  4. [4]

    Colton and R

    D. Colton and R. Kress.Integral Equation Methods in Scattering Theory. Pure and Applied Mathe- matics. John Wiley & Sons Inc., first edition, 1983

    work page 1983

  5. [5]

    Colton and R

    D. Colton and R. Kress.Inverse Acoustic and Electromagnetic Scattering Theory, volume 93. Springer, third edition, 2012

    work page 2012

  6. [6]

    M. P. Do Carmo.Differential Geometry of Curves and Surfaces. Courier Dover Publications, 2016

    work page 2016

  7. [7]

    C. L. Epstein and L. Greengard. Debye sources and the numerical solution of the time harmonic Maxwell equations.Communications on Pure and Applied Mathematics, 63(4):413–463, 2010

    work page 2010

  8. [8]

    C. L. Epstein, L. Greengard, and M. O’Neil. Debye sources and the numerical solution of the time harmonic Maxwell equations II.Communications on Pure and Applied Mathematics, 66(5):753–789, 2013

    work page 2013

  9. [9]

    Faria, C

    L. Faria, C. P´ erez-Arancibia, and C. Turc. Combined field-only boundary integral equations for PEC electromagnetic scattering problem in spherical geometries.SIAM Journal on Applied Mathematics, 84(1):19–38, 2024

    work page 2024

  10. [10]

    L. M. Faria, C. P´ erez-Arancibia, and M. Bonnet. General-purpose kernel regularization of boundary integral equations via density interpolation.Computer Methods in Applied Mechanics and Engineering, 378:113703, 2021

    work page 2021

  11. [11]

    Ganesh and S

    M. Ganesh and S. C. Hawkins. A spectrally accurate algorithm for electromagnetic scattering in three dimensions.Numerical Algorithms, 43:25–60, 2006. 26 Figure 2: Real part of thex-component of the total electric field Re(E x) (left) and thez-component of the total magnetic field Re(H z) (right), computed for the interlocking tori configuration atk=π, wit...

    work page 2006

  12. [12]

    Garza.Boundary Integral Equation Methods for Simulation and Design of Photonic Devices

    E. Garza.Boundary Integral Equation Methods for Simulation and Design of Photonic Devices. Cal- ifornia Institute of Technology, 2020

    work page 2020

  13. [13]

    Geuzaine and J.-F

    C. Geuzaine and J.-F. Remacle. Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities.International Journal for Numerical Methods in Engineering, 79(11):1309– 1331, 2009

    work page 2009

  14. [14]

    Gilbarg and N

    D. Gilbarg and N. S. Trudinger.Elliptic Partial Differential Equations of Second Order. Classics in Mathematics. Springer, reprint of the 1998 edition edition, 2001

    work page 1998

  15. [15]

    Hackbusch.Integral Equations: Theory and Numerical Treatment

    W. Hackbusch.Integral Equations: Theory and Numerical Treatment. Birkh¨ auser, 1995

    work page 1995

  16. [16]

    Hofmann, P

    B. Hofmann, P. Respondek, and S. B. Adrian. SphericalScattering: A Julia package for electromagnetic scattering from spherical objects.Journal of Open Source Software, 8(91):5820, 2023

    work page 2023

  17. [17]

    J. Hu, E. Garza, and C. Sideris. A Chebyshev-based high-order-accurate integral equation solver for Maxwell’s equations.IEEE Transactions on Antennas and Propagation, 69(9):5790–5800, 2021

    work page 2021

  18. [18]

    Inti.jl — Integral equation solver in Julia.https://github.com/ IntegralEquations/Inti.jl, 2025

    IntegralEquations contributors. Inti.jl — Integral equation solver in Julia.https://github.com/ IntegralEquations/Inti.jl, 2025. Accessed: 2025-05-25

    work page 2025

  19. [19]

    Kirsch and F

    A. Kirsch and F. Hettlich.The Mathematical Theory of Time-harmonic Maxwell’s Equations. Springer, 2015

    work page 2015

  20. [20]

    Klaseboer, Q

    E. Klaseboer, Q. Sun, and D. Y. Chan. Nonsingular field-only surface integral equations for electro- magnetic scattering.IEEE Transactions on Antennas and Propagation, 65(2):972–977, 2016. 27

    work page 2016

  21. [21]

    Kress.Linear Integral Equations, volume 82

    R. Kress.Linear Integral Equations, volume 82. Springer, 3rd edition, 2014

    work page 2014

  22. [22]

    W. C. H. McLean.Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press, 2000

    work page 2000

  23. [23]

    P´ erez-Arancibia, C

    C. P´ erez-Arancibia, C. Turc, and L. Faria. Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D.Journal of Computational Physics, 376:411–434, 2019

    work page 2019

  24. [24]

    P´ erez-Arancibia, C

    C. P´ erez-Arancibia, C. Turc, and L. Faria. Planewave density interpolation methods for 3D Helmholtz boundary integral equations.SIAM Journal on Scientific Computing, 41(4):A2088–A2116, 2019

    work page 2019

  25. [25]

    P´ erez-Arancibia, C

    C. P´ erez-Arancibia, C. Turc, L. M. Faria, and C. Sideris. Planewave density interpolation methods for the EFIE on simple and composite surfaces.IEEE Transactions on Antennas and Propagation, 69(1):317–331, 2020

    work page 2020

  26. [26]

    Saad and M

    Y. Saad and M. H. Schultz. GMRES: A generalized minimal residual algorithm for solving nonsym- metric linear systems.SIAM Journal on Scientific and Statistical Computing, 7(3):856–869, 1986

    work page 1986

  27. [27]

    J. A. Stratton and L. J. Chu. Diffraction theory of electromagnetic waves.Physical Review, 56(1):99– 107, July 1939

    work page 1939

  28. [28]

    Q. Sun, E. Klaseboer, and D. Y. Chan. Robust multiscale field-only formulation of electromagnetic scattering.Physical Review B, 95(4):045137, 2017

    work page 2017

  29. [29]

    Q. Sun, E. Klaseboer, A. J. Yuffa, and D. Y. Chan. Field-only surface integral equations: Scattering from a dielectric body.Journal of the Optical Society of America A, 37(2):284–293, 2020

    work page 2020

  30. [30]

    Q. Sun, E. Klaseboer, A. J. Yuffa, and D. Y. Chan. Field-only surface integral equations: Scattering from a perfect electric conductor.Journal of the Optical Society of America A, 37(2):276–283, 2020

    work page 2020

  31. [31]

    Taskinen and P

    M. Taskinen and P. Yla-Oijala. Current and charge integral equation formulation.IEEE Transactions on Antennas and Propagation, 54(1):58–67, 2006

    work page 2006

  32. [32]

    J. L. Volakis and K. Sertel.Integral Equation Methods for Electromagnetics. Scitech, 2012

    work page 2012

  33. [33]

    P. Werner. On the exterior boundary value problem of perfect reflection for stationary electromagnetic wave fields.Journal of Mathematical Analysis and Applications, 7:348–396, 1963

    work page 1963

  34. [34]

    R. A. Wildman, G. Pisharody, D. S. Weile, S. Balasubramaniam, and E. Michielssen. An accurate scheme for the solution of the time-domain integral equations of electromagnetics using higher or- der vector bases and bandlimited extrapolation.IEEE Transactions on Antennas and Propagation, 52(11):2973–2984, 2004. 28

    work page 2004