pith. sign in

arxiv: 1907.09746 · v1 · pith:CJ5IGOHMnew · submitted 2019-07-23 · 🧮 math.NA · cs.NA

Complex scaled infinite elements for exterior Helmholtz problems

Lothar Nannen , Markus Wess This is my paper

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

classification 🧮 math.NA cs.NA
keywords complex scalinginfinite elementsHelmholtz equationGalerkin methodexterior problemsresonance computationsuper-algebraic convergence
0
0 comments X

The pith

Galerkin discretization with infinite-support functions on complex-scaled domains yields super-algebraic error decay for exterior Helmholtz resonances.

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

The paper introduces a discretization for exterior Helmholtz resonance problems that applies complex coordinate stretching and then uses a Galerkin method whose basis functions extend to infinity. It proves that the error measured against the number of radial unknowns falls super-algebraically. Numerical tests also indicate that the same accuracy is reached with fewer degrees of freedom than with a standard perfectly matched layer. A reader would care because the approach directly addresses the cost of truncating unbounded domains while preserving high accuracy for resonance computations.

Core claim

The authors discretize complex scaled Helmholtz resonance problems using a Galerkin method with ansatz functions that have infinite support. They establish that the approximation error decays super algebraically with respect to the number of unknowns in radial direction.

What carries the argument

Complex coordinate stretching that produces exponentially decaying solutions, combined with a Galerkin method employing infinite-support ansatz functions.

If this is right

  • The discretization error in the radial direction decreases faster than any algebraic rate once the complex scaling is applied.
  • Numerical examples demonstrate higher efficiency than a standard perfectly matched layer method for the same accuracy level.
  • The technique applies to time-harmonic wave equations whose solutions become exponentially localized after complex scaling.

Where Pith is reading between the lines

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

  • The same infinite-element construction could be tested on three-dimensional exterior problems where radial degrees of freedom dominate the cost.
  • If the super-algebraic rate persists under mesh refinement in the angular directions, the method would offer a practical alternative for resonance search in scattering configurations.
  • The approach implicitly suggests that other unbounded-domain problems admitting exponential decay after scaling might benefit from analogous infinite-support bases.

Load-bearing premise

The complex stretching must generate solutions whose exponential decay is captured accurately by the chosen infinite-support basis functions inside the Galerkin formulation.

What would settle it

A sequence of computations that records the resonance error for successively more radial basis functions and finds only algebraic rather than super-algebraic decay would falsify the central claim.

Figures

Figures reproduced from arXiv: 1907.09746 by Lothar Nannen, Markus Wess.

Figure 1
Figure 1. Figure 1: Two dimensional example domains and exterior coordinates [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Approximation error of h0(ω(R + σ·)) for varying parameters. The dashed lines mark the predicted convergence rates from Section 5 [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Approximation error of hν(ω(R+iσ·)). The dashed lines mark the predicted convergence rates from Section 5. 23 [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Condition numbers of the discretization matrices of ˜s [PITH_FULL_IMAGE:figures/full_fig_p024_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Errors of the eigenvalue ω ≈ 2.903916 − 1.201866i obtained by solving the separated problem for ν = 3. The dashed lines mark the squared exponential and super algebraic convergence rates from Section 5 respectively. Note the different scalings on the horizontal axes. 6.2 Computational costs In this subsection we compare the computational costs of our infinite elements and a conventional PML by approximatin… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of errors against factorization times for eigenvalues [PITH_FULL_IMAGE:figures/full_fig_p026_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Potential functions [PITH_FULL_IMAGE:figures/full_fig_p027_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Inhomogeneous exterior problem with radial inhomogeneity. The lines mark the locations [PITH_FULL_IMAGE:figures/full_fig_p028_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: shows resonances of the same problem with an additional potential pˆ(ˆx) := ˆp(x, y, z) := z, (cf. Figures 7b and 7c) and varying values of ˆ. This problem is not separable any more, thus only a three dimensional simulation is possible. Due to the disturbed symmetry, the multiple eigenvalues fan out. 1.2 1.4 1.6 1.8 2 2.2 2.4 −0.85 −0.8 −0.75 −0.7 −0.65 −0.6 ν = 2 ν = 3 <(ω) ˆ = 0 ˆ = 0.5 1.2 1.4 1.6 1.… view at source ↗
Figure 10
Figure 10. Figure 10: Resonance functions corresponding to eigenvalues from Figure 8 with ˜ε [PITH_FULL_IMAGE:figures/full_fig_p029_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Resonance functions corresponding to eigenvalues from Figure 9 with ˆε [PITH_FULL_IMAGE:figures/full_fig_p029_11.png] view at source ↗
read the original abstract

The technique of complex scaling for time harmonic wave type equations relies on a complex coordinate stretching to generate exponentially decaying solutions. In this work, we use a Galerkin method with ansatz functions with infinite support to discretize complex scaled Helmholtz resonance problems. We show that the approximation error of the method decays super algebraically with respect to the number of unknowns in radial direction. Numerical examples underline the theoretical findings and show the superior efficiency of our method compared to a standard perfectly matched layer method.

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

Summary. The paper introduces a Galerkin discretization of complex-scaled exterior Helmholtz resonance problems that employs ansatz functions with infinite radial support. It establishes a theoretical result that the approximation error decays superalgebraically in the number of radial unknowns once the complex scaling has produced exponentially decaying solutions, and it presents numerical examples demonstrating superior efficiency relative to a standard perfectly matched layer discretization.

Significance. If the superalgebraic convergence holds under the stated hypotheses on the complex scaling, the approach supplies a theoretically justified and computationally efficient alternative to truncated PML methods for exterior resonance problems. The combination of an established transformation with an infinite-element ansatz that exploits the resulting exponential decay is a natural and potentially advantageous discretization strategy; the explicit rate result and the numerical comparisons are the primary contributions.

minor comments (2)
  1. The abstract and introduction would benefit from a brief statement of the precise function spaces and the precise hypotheses on the stretching function that are used in the convergence theorem.
  2. Figure captions should explicitly indicate the radial polynomial degree and the number of radial degrees of freedom employed in each run so that the superalgebraic decay can be directly compared with the theory.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript, the recognition of the theoretical result on superalgebraic convergence, and the recommendation to accept. No major comments were raised that require a point-by-point reply.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's central result is a mathematical proof that the Galerkin error for infinite-support ansatz functions on complex-scaled exterior Helmholtz problems decays superalgebraically in the number of radial degrees of freedom. This follows from the established complex-scaling transformation (which produces exponentially decaying solutions) combined with standard approximation theory for the chosen basis functions; the derivation does not reduce any claimed prediction or uniqueness statement to a fitted quantity, self-citation chain, or definitional tautology. The abstract and method description treat the stretching as an external input and the error bound as a consequence of that input under the stated hypotheses, with no load-bearing step that collapses to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review performed on abstract only; no explicit free parameters, additional axioms, or invented entities are identifiable from the given text.

axioms (1)
  • domain assumption Complex coordinate stretching generates exponentially decaying solutions for time-harmonic Helmholtz problems.
    Invoked in the first sentence of the abstract as the foundation of the technique.

pith-pipeline@v0.9.0 · 5593 in / 1123 out tokens · 26839 ms · 2026-05-24T17:28:20.252508+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

29 extracted references · 29 canonical work pages

  1. [1]

    Milton Abramowitz and Irene A. Stegun. Handbook of mathematical functions with formulas, graphs, and mathematical tables , volume 55 of National Bureau of Standards Applied Mathe- matics Series. For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, D.C., 1964

    work page 1964

  2. [2]

    Babuˇ ska and J

    I. Babuˇ ska and J. Osborn. Eigenvalue problems. In Handbook of numerical analysis, Vol. II , Handb. Numer. Anal., II, pages 641–787. North-Holland, Amsterdam, 1991

    work page 1991

  3. [3]

    A perfectly matched layer for the absorption of electromagnetic waves

    Jean-Pierre Berenger. A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. , 114(2):185–200, 1994

    work page 1994

  4. [4]

    Berm´ udez, L

    A. Berm´ udez, L. Hervella-Nieto, A. Prieto, and R. Rodr´ ı guez. An exact bounded perfectly matched layer for time-harmonic scattering problems. SIAM J. Sci. Comput. , 30(1):312–338, 2007/08

    work page 2007

  5. [5]

    Bramble and Joseph E

    James H. Bramble and Joseph E. Pasciak. Analysis of a Cartesian PML approximation to acoustic scattering problems in R2 and R3. J. Comput. Appl. Math. , 247:209–230, 2013

    work page 2013

  6. [6]

    The perfectly matched layer in curvilinear coordinates

    Francis Collino and Peter Monk. The perfectly matched layer in curvilinear coordinates. SIAM J. Sci. Comput. , 19(6):2061–2090 (electronic), 1998

    work page 2061

  7. [7]

    Inverse acoustic and electromagnetic scattering theory , vol- ume 93 of Applied Mathematical Sciences

    David Colton and Rainer Kress. Inverse acoustic and electromagnetic scattering theory , vol- ume 93 of Applied Mathematical Sciences. Springer-Verlag, Berlin, second edition, 1998

    work page 1998

  8. [8]

    Convergence of hardy space infinite elements for helmholtz scattering and resonance problems

    Martin Halla. Convergence of hardy space infinite elements for helmholtz scattering and resonance problems. SIAM Journal on Numerical Analysis , 54(3):1385–1400, 2016

    work page 2016

  9. [9]

    Hardy space infinite elements for time harmonic wave equations with phase and group velocities of different signs

    Martin Halla, Thorsten Hohage, Lothar Nannen, and Joachim Sch¨ oberl. Hardy space infinite elements for time harmonic wave equations with phase and group velocities of different signs. Numer. Math., 133(1):103–139, 2016

    work page 2016

  10. [10]

    Hardy space infinite elements for time-harmonic two- dimensional elastic waveguide problems

    Martin Halla and Lothar Nannen. Hardy space infinite elements for time-harmonic two- dimensional elastic waveguide problems. Wave Motion, 59:94–110, 2015

    work page 2015

  11. [11]

    Two scale Hardy space infinite elements for scalar waveguide problems

    Martin Halla and Lothar Nannen. Two scale Hardy space infinite elements for scalar waveguide problems. Adv. Comput. Math. , 44(3):611–643, 2018

    work page 2018

  12. [12]

    P. D. Hislop and I. M. Sigal. Introduction to spectral theory, volume 113 of Applied Mathemat- ical Sciences. Springer-Verlag, New York, 1996. With applications to Schr¨ odinger operators

    work page 1996

  13. [13]

    Hardy space infinite elements for scattering and reso- nance problems

    Thorsten Hohage and Lothar Nannen. Hardy space infinite elements for scattering and reso- nance problems. SIAM J. Numer. Anal. , 47(2):972–996, 2009

    work page 2009

  14. [14]

    Convergence of infinite element methods for scalar waveguide problems

    Thorsten Hohage and Lothar Nannen. Convergence of infinite element methods for scalar waveguide problems. BIT Numerical Mathematics , 55(1):215–254, 2015

    work page 2015

  15. [15]

    Solving time-harmonic scattering problems based on the pole condition

    Thorsten Hohage, Frank Schmidt, and Lin Zschiedrich. Solving time-harmonic scattering problems based on the pole condition. I. Theory. SIAM J. Math. Anal. , 35(1):183–210, 2003. 31

    work page 2003

  16. [16]

    Solving time-harmonic scattering problems based on the pole condition

    Thorsten Hohage, Frank Schmidt, and Lin Zschiedrich. Solving time-harmonic scattering problems based on the pole condition. II. Convergence of the PML method. SIAM J. Math. Anal., 35(3):547–560, 2003

    work page 2003

  17. [17]

    Seungil Kim and Joseph E. Pasciak. The computation of resonances in open systems using a perfectly matched layer. Math. Comp., 78(267):1375–1398, 2009

    work page 2009

  18. [18]

    Lassas and E

    M. Lassas and E. Somersalo. On the existence and the convergence of the solution of the pml equations. Computing, 60:229–241, 1998

    work page 1998

  19. [19]

    Quantum theory of resonances: Calculating energies, width and cross-sections by complex scaling

    N Moiseyev. Quantum theory of resonances: Calculating energies, width and cross-sections by complex scaling. Physics reports, 302:211–293, 1998

    work page 1998

  20. [20]

    Hardy space infinite elements for Helmholtz-type problems with unbounded inhomogeneities

    Lothar Nannen and Achim Sch¨ adle. Hardy space infinite elements for Helmholtz-type problems with unbounded inhomogeneities. Wave Motion, 48(2):116–129, 2010

    work page 2010

  21. [21]

    Computing scattering resonances using perfectly matched layers with frequency dependent scaling functions

    Lothar Nannen and Markus Wess. Computing scattering resonances using perfectly matched layers with frequency dependent scaling functions. BIT, 58(2):373–395, 2018

    work page 2018

  22. [22]

    Numerical methods for large eigenvalue problems , volume 66 of Classics in Applied Mathematics

    Yousef Saad. Numerical methods for large eigenvalue problems , volume 66 of Classics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2011. Revised edition of the 1992 original [ 1177405]

    work page 2011

  23. [23]

    A new approach to coupled interior-exterior Helmholtz-type problems: Theory and algorithms

    Frank Schmidt. A new approach to coupled interior-exterior Helmholtz-type problems: Theory and algorithms. Habilitation, Freie Universitt Berlin, 2002

    work page 2002

  24. [24]

    Discrete transparent boundary conditions for the numer- ical solution of Fresnel’s equation

    Frank Schmidt and Peter Deuflhard. Discrete transparent boundary conditions for the numer- ical solution of Fresnel’s equation. Computers Math. Appl. , 29:53–76, 1995

    work page 1995

  25. [25]

    Netgen - an advancing front 2d/3d-mesh generator based on abstract rules

    Joachim Sch¨ oberl. Netgen - an advancing front 2d/3d-mesh generator based on abstract rules. Comput.Visual.Sci, 1:41–52, 1997

    work page 1997

  26. [26]

    C++11 implementation of finite elements in ngsolve

    Joachim Sch¨ oberl. C++11 implementation of finite elements in ngsolve. Preprint 30/2014, Institute for Analysis and Scientific Computing, TU Wien, 2014

    work page 2014

  27. [27]

    Spectral methods, volume 41 of Springer Series in Computational Mathematics

    Jie Shen, Tao Tang, and Li-Lian Wang. Spectral methods, volume 41 of Springer Series in Computational Mathematics . Springer, Heidelberg, 2011. Algorithms, analysis and applica- tions

    work page 2011

  28. [28]

    Steinbach and G

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

    work page 2012

  29. [29]

    Nico M. Temme. Asymptotic methods for integrals , volume 6 of Series in Analysis . World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2015. 32

    work page 2015