Modal analysis of a domain decomposition method for Maxwell's equations in a waveguide
Pith reviewed 2026-05-23 02:33 UTC · model grok-4.3
The pith
Limiting spectra of Toeplitz matrices govern convergence of one-level Schwarz methods for Maxwell equations in general waveguides
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The convergence of the one-level Schwarz iteration for Maxwell's equations in strip-wise decomposed waveguides is governed by the limiting spectrum of associated Toeplitz matrices combined with the modal decomposition of the electromagnetic solutions, allowing extension to arbitrary cross-sections and impedance or PML transmission conditions.
What carries the argument
Limiting spectrum of Toeplitz matrices arising from strip-wise decomposition, analyzed together with modal decomposition of Maxwell solutions
Load-bearing premise
The limiting spectrum of the Toeplitz matrices from strip-wise decomposition together with modal decomposition continues to govern convergence of the one-level Schwarz iteration for arbitrary cross-sections and impedance or PML conditions.
What would settle it
A numerical computation of the actual iteration operator spectrum for a waveguide with non-rectangular cross-section that deviates substantially from the predicted limiting spectrum.
read the original abstract
Time-harmonic wave propagation problems, especially those governed by Maxwell's equations, pose significant computational challenges due to the non-self-adjoint nature of the operators and the large, non-Hermitian linear systems resulting from discretization. Domain decomposition methods, particularly one-level Schwarz methods, offer a promising framework to tackle these challenges, with recent advancements showing the potential for weak scalability under certain conditions. In this paper, we analyze the weak scalability of one-level Schwarz methods for Maxwell's equations in strip-wise domain decompositions, focusing on waveguides with general cross sections and different types of transmission conditions such as impedance or perfectly matched layers (PMLs). By combining techniques from the limiting spectrum analysis of Toeplitz matrices and the modal decomposition of Maxwell's solutions, we provide a novel theoretical framework that extends previous work to more complex geometries and transmission conditions. Numerical experiments confirm that the limiting spectrum effectively predicts practical behavior even with a modest number of subdomains. Furthermore, we demonstrate that the one-level Schwarz method can achieve robustness with respect to the wave number under specific domain decomposition parameters, offering new insights into its applicability for large-scale electromagnetic wave problems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to develop a novel theoretical framework for the weak scalability of one-level Schwarz domain decomposition methods applied to time-harmonic Maxwell's equations in waveguides. The framework combines limiting-spectrum analysis of Toeplitz matrices with modal decomposition of the solutions, extending prior results to waveguides with arbitrary cross-sections and to impedance or PML transmission conditions. Numerical experiments are said to confirm that the limiting spectrum predicts practical convergence behavior even for modest numbers of subdomains, and that robustness with respect to the wave number is achievable under suitable domain-decomposition parameters.
Significance. If the structural preservation of the limiting spectrum under arbitrary cross-sections and PML stretching can be rigorously established, the work would supply a useful analytic tool for predicting convergence of non-overlapping Schwarz methods in electromagnetic waveguide problems. The explicit combination of Toeplitz limiting spectra with modal analysis, together with the reported numerical confirmation for small subdomain counts, would constitute a concrete advance over purely numerical studies of the same methods.
major comments (2)
- [Abstract / theoretical framework] Abstract (paragraph on the novel theoretical framework): the central claim requires that the limiting spectrum of the Toeplitz matrices arising from strip-wise decomposition, together with modal decomposition, continues to control one-level Schwarz convergence when the waveguide cross-section is arbitrary. For non-rectangular cross-sections the transverse eigenmodes generally lack the separability that produces a clean block-Toeplitz structure along the guide axis; the manuscript must show explicitly (e.g., in the derivation of the interface symbol) that the limiting-spectrum formulas survive this loss of separability without additional assumptions on the transverse operator.
- [Abstract / theoretical framework] Abstract (paragraph on PML transmission conditions): PML stretching perturbs the interface symbols that enter the Toeplitz analysis. The manuscript must demonstrate that the limiting-spectrum formulas remain valid under this perturbation; otherwise the extension to PML conditions rests on an unverified structural assumption that is load-bearing for the claimed generality.
minor comments (1)
- [Numerical experiments] The abstract states that numerical experiments confirm the predictions but supplies no information on mesh sizes, number of subdomains tested, or quantitative error measures; these details should be added to the numerical section for reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address the two major comments point by point below. In both cases we agree that additional explicit derivations are needed to fully substantiate the claims and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Abstract / theoretical framework] Abstract (paragraph on the novel theoretical framework): the central claim requires that the limiting spectrum of the Toeplitz matrices arising from strip-wise decomposition, together with modal decomposition, continues to control one-level Schwarz convergence when the waveguide cross-section is arbitrary. For non-rectangular cross-sections the transverse eigenmodes generally lack the separability that produces a clean block-Toeplitz structure along the guide axis; the manuscript must show explicitly (e.g., in the derivation of the interface symbol) that the limiting-spectrum formulas survive this loss of separability without additional assumptions on the transverse operator.
Authors: The modal decomposition is performed with respect to the eigenfunctions of the transverse Maxwell operator defined on the (possibly non-rectangular) cross-section. These eigenfunctions form a complete basis that diagonalizes the transverse operator, thereby reducing the original problem to a countable set of independent one-dimensional modal problems along the waveguide axis. For each fixed mode the strip-wise decomposition produces a block-Toeplitz operator whose symbol depends only on the modal wave number and the chosen transmission condition; the transverse geometry enters solely through the modal eigenvalue and does not affect the longitudinal Toeplitz structure. Consequently the limiting-spectrum formulas derived for the one-dimensional modal problems apply directly. We will add an explicit derivation of the interface symbol in modal coordinates to the revised manuscript, confirming that the argument requires only the standard spectral properties of the transverse operator and no further separability assumptions. revision: yes
-
Referee: [Abstract / theoretical framework] Abstract (paragraph on PML transmission conditions): PML stretching perturbs the interface symbols that enter the Toeplitz analysis. The manuscript must demonstrate that the limiting-spectrum formulas remain valid under this perturbation; otherwise the extension to PML conditions rests on an unverified structural assumption that is load-bearing for the claimed generality.
Authors: We agree that the perturbation induced by PML stretching must be treated explicitly. In the modal setting the PML appears as a complex stretching of the longitudinal coordinate near the artificial interfaces, which modifies the transmission coefficients but leaves the block-Toeplitz character of the discrete operator intact. The limiting spectrum is obtained from the adjusted symbol that incorporates the stretched transmission parameters. We will insert a dedicated subsection deriving this modified symbol and verifying that the spectral bounds (and the resulting wave-number robustness under suitable parameter choices) continue to hold. This derivation will be placed immediately after the impedance-case analysis so that the extension to PML is fully rigorous. revision: yes
Circularity Check
No circularity; derivation extends external Toeplitz limiting-spectrum and modal-analysis techniques
full rationale
The provided abstract and claims describe a framework obtained by combining standard limiting-spectrum analysis of Toeplitz matrices with modal decomposition of Maxwell solutions, then extending the combination to arbitrary cross-sections and impedance/PML transmission conditions. No equation, definition, or claim is shown that reduces a derived quantity to a fitted parameter taken from the same data, nor does any load-bearing step rest on a self-citation whose content is itself unverified within the paper. Numerical experiments are presented only as confirmation of the already-derived limiting spectrum, not as the source of the spectrum formulas. The derivation is therefore self-contained against the cited external mathematical techniques.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Limiting spectrum analysis applies to the Toeplitz matrices generated by the strip-wise discretization of the waveguide operator
- domain assumption Modal decomposition of Maxwell solutions remains valid for general cross sections under the chosen transmission conditions
Reference graph
Works this paper leans on
-
[1]
I. M. Babu s ka and S. A. Sauter , Is the pollution effect of the FEM avoidable for the H elmholtz equation considering high wave numbers? , SIAM J. Numer. Anal., 34 (1997), pp. 2392--2423
work page 1997
-
[2]
M. Bonazzoli, V. Dolean, I. G. Graham, E. A. Spence, and P.-H. Tournier , Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption , Math. Comp., 88 (2019), pp. 2559--2604
work page 2019
-
[3]
N. Bootland, S. Borzooei, V. Dolean, and P.-H. Tournier , Numerical assessment of pml transmission conditions in a domain decomposition method for the helmholtz equation , in International Conference on Domain Decomposition Methods, Springer, 2022, pp. 445--453
work page 2022
-
[4]
N. Bootland, V. Dolean, A. Kyriakis, and J. Pestana , Analysis of parallel schwarz algorithms for time-harmonic problems using block toeplitz matrices , Electronic Transactions on Numerical Analysis, 55 (2022), pp. 112--141, https://doi.org/10.1553/etna_vol55s112
-
[5]
J. Bramble and J. Pasciak , Analysis of a finite pml approximation for the three dimensional time-harmonic maxwell and acoustic scattering problems , Mathematics of Computation, 76 (2007), pp. 597--614
work page 2007
- [6]
-
[7]
E. Cances, Y. Maday, and B. Stamm , Domain decomposition for implicit solvation models , J. Chem. Phys., 139(5) (2013), https://doi.org/10.1063/1.4816767
-
[8]
G. Ciaramella and M. J. Gander , Analysis of the parallel S chwarz method for growing chains of fixed-sized subdomains: P art I , SIAM J. Numer. Anal., 55 (2017), pp. 1330--1356, https://doi.org/10.1137/16M1065215, https://doi.org/10.1137/16M1065215
-
[9]
G. Ciaramella, M. Hassan, and B. Stamm , On the scalability of the S chwarz method , SMAI J. Comput. Math., 6 (2020), pp. 33--68, https://doi.org/10.5802/smai-jcm.61, https://doi.org/10.5802/smai-jcm.61
-
[10]
A.-S. B.-B. Dhia and \'E . Lun \'e ville , Propagation dans les guides d'ondes , (2021)
work page 2021
- [11]
- [12]
-
[13]
V. Dolean, M. J. Gander, S. Lanteri, J.-F. Lee, and Z. Peng , Effective transmission conditions for domain decomposition methods applied to the time-harmonic curl--curl M axwell's equations , J. Comput. Phys., 280 (2015), pp. 232--247, https://doi.org/10.1016/j.jcp.2014.09.024
-
[14]
V. Dolean, P. Jolivet, and F. Nataf , An Introduction to Domain Decomposition Methods: Algorithms, Theory, and Parallel Implementation , SIAM, Philadelphia, 2015, https://doi.org/10.1137/1.9781611974065.ch1
- [15]
-
[16]
M. El Bouajaji , V. Dolean, M. J. Gander, and S. Lanteri , Optimized S chwarz methods for the time-harmonic M axwell equations with damping , SIAM J. Sci. Comput., 34 (2012), pp. A2048--A2071
work page 2012
-
[17]
M. J. Gander, L. Halpern, and F. Magoul \`e s , An optimized S chwarz method with two-sided R obin transmission conditions for the H elmholtz equation , Internat. J. Numer. Methods Fluids, 55 (2007), pp. 163--175
work page 2007
-
[18]
M. J. Gander and A. Tonnoir , Analysis of schwarz methods for convected helmholtz-like equations , SIAM Journal on Scientific Computing, 46 (2024), pp. A1--A22
work page 2024
-
[19]
S. Gong, I. G. Graham, and E. A. Spence , Domain decomposition preconditioners for high-order discretizations of the heterogeneous H elmholtz equation , IMA J. Numer. Anal., (2020). draa080
work page 2020
-
[20]
I. G. Graham, E. A. Spence, and E. Vainikko , Recent results on domain decomposition preconditioning for the high-frequency helmholtz equation using absorption , in Modern Solvers for Helmholtz Problems, D. Lahaye, J. Tang, and K. Vuik, eds., Geosystems Mathematics, Birkh \"a user, Cham, 2017, pp. 3--26
work page 2017
-
[21]
I. G. Graham, E. A. Spence, and J. Zou , Domain decomposition with local impedance conditions for the H elmholtz equation , SIAM J. Numer. Anal., 58 (2020), pp. 2515--2543
work page 2020
- [22]
-
[23]
W. F. Hall and A. V. Kabakian , A sequence of absorbing boundary conditions for maxwell’s equations , Journal of Computational Physics, 194 (2004), pp. 140--155
work page 2004
-
[24]
Hecht , New development in freefem++ , J
F. Hecht , New development in freefem++ , J. Numer. Math., 20 (2012), pp. 251--265, https://freefem.org/
work page 2012
-
[25]
V. Mattesi, M. Darbas, and C. Geuzaine , A high-order absorbing boundary condition for 2D time-harmonic elastodynamic scattering problems , Comput. Math. Appl., 77 (2019), pp. 1703--1721
work page 2019
- [26]
-
[27]
P.-H. Tournier, P. Jolivet, and F. Nataf , FFDDM : Freefem domain decomposition method . https://doc.freefem.org/documentation/ffddm/index.html , 2019
work page 2019
-
[28]
" write newline "" before.all 'output.state := FUNCTION fin.entry add.period write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1 'skip if FUNCTION new.block.checka empty 'skip 'new.block if FUNCTION field.or.null duplicate empty pop "" 'skip ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.