pith. sign in

arxiv: 2511.11016 · v2 · pith:ENTZLZLEnew · submitted 2025-11-14 · 🧮 math.NA · cs.NA· math-ph· math.MP· math.SP

Bifurcations in Interior Transmission Eigenvalues: Theory and Computation

Davide Pradovera , Alessandro Borghi , Lukas Pieronek , Andreas Kleefeld This is my paper

Pith reviewed 2026-05-21 18:30 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath-phmath.MPmath.SP
keywords interior transmission eigenvaluesbifurcationsnonlinear eigenproblemscontour integrationinverse scatteringspectral analysisradial symmetryparameter continuation
0
0 comments X

The pith

The interior transmission eigenvalue spectrum can bifurcate and become non-smooth even when the refractive index varies smoothly in the underlying PDE.

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

The interior transmission eigenvalue problem is central to inverse scattering and the spectral analysis of inhomogeneous media. The paper shows that the spectral map from material parameters to eigenpairs can nevertheless display non-smooth behavior and bifurcations. It supplies sufficient conditions for this phenomenon on general domains and sharper characterizations when the geometry is radially symmetric. The analysis is paired with a numerical method that treats the problem as a parametric nonlinear eigenproblem and tracks eigenvalue trajectories with an adaptive contour solver. Experiments confirm the predicted non-smooth effects and uncover additional ones.

Core claim

The authors establish a theoretical framework that identifies sufficient conditions for non-smooth spectral behavior and bifurcations in the interior transmission eigenvalue problem on general domains. For radially symmetric geometries they obtain a more precise description of the bifurcation points. They reformulate the problem as a parametric discrete nonlinear eigenproblem and apply a match-based adaptive contour eigensolver to follow eigenvalue paths under parameter variation; numerical tests verify the theory and reveal further non-smooth spectral phenomena.

What carries the argument

The match-based adaptive contour eigensolver applied to the parametric discrete nonlinear eigenproblem formulation of the interior transmission problem, together with the theoretical conditions that guarantee non-smooth spectral behavior.

If this is right

  • Bifurcation points in the spectrum can be located theoretically before any numerical computation is performed.
  • Eigenvalue trajectories remain continuous except at explicitly characterizable parameter values in radial geometries.
  • The adaptive contour method can be used to monitor and detect these non-smooth transitions in practice.
  • Non-smooth spectral effects must be accounted for when interpreting scattering data from inhomogeneous media.

Where Pith is reading between the lines

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

  • Bifurcations may introduce sudden changes in the number or location of eigenvalues that affect the conditioning of inverse scattering reconstructions.
  • The same theoretical conditions could be checked in related transmission or scattering eigenvalue problems outside the interior transmission setting.
  • Numerical tests on non-radial domains with moderate asymmetry would provide a direct check on how far the general-domain conditions extend.

Load-bearing premise

The match-based adaptive contour eigensolver accurately and completely tracks eigenvalue trajectories under continuous parameter changes without missing or mischaracterizing bifurcations.

What would settle it

A computation on a radially symmetric domain in which the theory predicts a bifurcation point but the numerical solver shows a smooth crossing of eigenvalues, or conversely a smooth crossing where the theory predicts a bifurcation.

Figures

Figures reproduced from arXiv: 2511.11016 by Alessandro Borghi, Andreas Kleefeld, Davide Pradovera, Lukas Pieronek.

Figure 1
Figure 1. Figure 1: Examples of non-smooth eigenpair behavior. We set [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Seven ITE trajectories (five of which are purely real) for the unit disk in the parameter [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Real (left) and imaginary (center) parts of the ITE trajectories for the unit disk. [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Indicator I near three bifurcations for the ITE trajectories for the unit disk. Line colors and markers are the same as in [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Ten ITE trajectories (eight of which are purely real) with [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Real (left) and imaginary (center) parts of the ITE trajectories for the annulus with [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Six ITE trajectories with m = 1 for the annulus with r = 0.1 in the parameter range p ∈ [4, 25]. Arrows are used to indicate the direction of travel along the dotted magenta curve with triangular markers. When the dotted curve is on the real axis, the travel direction is always towards the right. The symbol “⊗” marks the position of the Laplacian eigenvalue κ ⋆ . In [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Real (left) and imaginary (center) parts of the ITE trajectories for the annulus with [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Indicator I near eight bifurcations for the ITE trajectories for the annulus with m = 1 and r = 0.1. which prescribes convergence of non-real trajectories to eigenvalues of the negative Dirichlet Laplacian on D as p → ∞. We show real and imaginary parts of the eigenvalues in [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Indicator I near eight bifurcations for the ITE trajectories for the annulus with m = 1 and r = 0.1. Line colors and markers are the same as in [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Some ITE trajectories with m = 1 for the unit disk. functions, focusing on their smoothness with respect to variations in the refractive index. Our general results advance the understanding of ITP spectra and provide simple, practical criteria to verify the local smoothness of real ITP eigenpair trajectories. These criteria enable a quick diagnostic for identifying bifurcations and other exceptional point… view at source ↗
read the original abstract

The interior transmission eigenvalue problem (ITP) plays a central role in inverse scattering theory and in the spectral analysis of inhomogeneous media. Despite its smooth dependence on the refractive index at the PDE level, the corresponding spectral map from material parameters to eigenpairs may exhibit non-smooth or bifurcating behavior. In this work, we develop a theoretical framework identifying sufficient conditions for such non-smooth spectral behavior in the ITP on general domains. We further specialize our analysis to some radially symmetric geometries, enabling a more precise characterization of bifurcations in the spectrum. Computationally, we formulate the ITP as a parametric, discrete, nonlinear eigenproblem and use a match-based adaptive contour eigensolver to accurately and efficiently track eigenvalue trajectories under parameter variation. Numerical experiments confirm the theoretical predictions and reveal novel non-smooth spectral 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

1 major / 1 minor

Summary. The paper claims to develop a theoretical framework identifying sufficient conditions for non-smooth spectral behavior in the ITP on general domains, with specialization to radially symmetric geometries for precise bifurcation characterization. It formulates the ITP as a parametric discrete nonlinear eigenproblem solved via a match-based adaptive contour eigensolver to track eigenvalue trajectories, and uses numerical experiments to confirm predictions and reveal novel effects.

Significance. If validated, this could advance the spectral analysis in inverse scattering by explaining non-smooth behaviors in transmission eigenvalues. The theoretical specialization to radial cases and the computational tracking approach are notable strengths, provided the solver's accuracy at bifurcation points is established.

major comments (1)
  1. The formulation of the parametric nonlinear eigenproblem and the match-based adaptive contour eigensolver lacks discussion of error bounds or robustness near non-smooth parameter values where bifurcations occur. Without this, the numerical confirmation of the theoretical predictions on non-smooth effects cannot be considered conclusive, as the solver might miss or mischaracterize bifurcations.
minor comments (1)
  1. The abstract could briefly mention one example of the novel non-smooth spectral effect to better engage readers.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed review and constructive feedback on our work concerning bifurcations in interior transmission eigenvalues. We address the major comment below and will revise the manuscript to strengthen the numerical validation.

read point-by-point responses

  1. Referee: The formulation of the parametric nonlinear eigenproblem and the match-based adaptive contour eigensolver lacks discussion of error bounds or robustness near non-smooth parameter values where bifurcations occur. Without this, the numerical confirmation of the theoretical predictions on non-smooth effects cannot be considered conclusive, as the solver might miss or mischaracterize bifurcations.

    Authors: We agree that a dedicated discussion of error control and robustness near bifurcation points would enhance the rigor of the numerical results. The match-based adaptive contour eigensolver adjusts the integration contour dynamically according to eigenvalue estimates and employs a matching criterion to detect spectral crossings, which in practice allows reliable tracking through non-smooth transitions. However, explicit a-priori error bounds tailored to these points were not provided in the original manuscript. In the revised version we will add a subsection on the solver's numerical analysis, including bounds derived from the contour radius, the adaptive step-size control, and the tolerance of the matching procedure. We will also include supplementary numerical experiments that quantify the solver's accuracy at the predicted bifurcation values, thereby confirming that the observed non-smooth effects are not artifacts of the discretization or contour integration. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation builds new conditions and solver on established ITP theory

full rationale

The paper introduces sufficient conditions for non-smooth spectral behavior in the ITP on general domains, specializes the analysis to radially symmetric cases for bifurcation characterization, reformulates the ITP as a parametric discrete nonlinear eigenproblem, and applies a match-based adaptive contour eigensolver to track trajectories. Numerical experiments then confirm the predictions. No load-bearing step reduces by construction to a fitted parameter, self-definition, or self-citation chain; the central claims add independent theoretical content and computational formulation without the outputs being equivalent to the inputs by definition. The work is self-contained against external ITP benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Relies on standard domain assumptions from inverse scattering theory; no free parameters or invented entities explicitly introduced in abstract.

axioms (1)
  • domain assumption The interior transmission eigenvalue problem exhibits smooth dependence on the refractive index at the PDE level.
    Explicitly stated as background in the abstract.

pith-pipeline@v0.9.0 · 5683 in / 1244 out tokens · 49977 ms · 2026-05-21T18:30:54.556585+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

43 extracted references · 43 canonical work pages

  1. [1]

    Abramowitz and I

    M. Abramowitz and I. A. Stegun , Handbook of mathematical functions with formulas, graphs, and mathematical tables, vol. 55, US Government printing office, 1964

    work page 1964

  2. [2]

    Amsallem and C

    D. Amsallem and C. F arhat , An online method for interpolating linear parametric reduced-order models, SIAM Journal on Scientific Computing, 33 (2011), pp. 2169–2198, https://doi.org/10.1137/ 100813051. ISBN: 9550101053

    work page 2011

  3. [3]

    A. L. Andrew, K.-W. E. Chu, and P. Lancaster , Derivatives of eigenvalues and eigenvectors of matrix functions, SIAM Journal on Matrix Analysis and Applications, 14 (1993), pp. 903–926, https://doi.org/10.1137/0614061

    work page doi:10.1137/0614061 1993

  4. [4]

    Beyn, An integral method for solving nonlinear eigenvalue problems, Linear Algebra and its Ap- plications, 436 (2012), pp

    W.-J. Beyn, An integral method for solving nonlinear eigenvalue problems, Linear Algebra and its Ap- plications, 436 (2012), pp. 3839–3863, https://doi.org/10.1016/j.laa.2011.03.030. 26 D. PRADOVERA, A. BORGHI, L. PIERONEK, AND A. KLEEFELD

    work page doi:10.1016/j.laa.2011.03.030 2012

  5. [5]

    Cakoni, M

    F. Cakoni, M. Çayören, and D. Colton , Transmission eigenvalues and the nondestructive testing of dielectrics, Inverse Problems, 24 (2008), 065016 (15 pages), https://doi.org/10.1088/0266-5611/24/ 6/065016

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

  6. [6]

    Cakoni and D

    F. Cakoni and D. Colton , Qualitative Methods in Inverse Scattering Theory - An Introduction, Springer, Berlin, 2006

    work page 2006

  7. [7]

    Cakoni, D

    F. Cakoni, D. Colton, and D. Gintides , The interior transmission eigenvalue problem, SIAM Journal on Mathematical Analysis, 42 (2010), pp. 2912–2921, https://doi.org/10.1137/100793542

    work page doi:10.1137/100793542 2010

  8. [8]

    Cakoni, D

    F. Cakoni, D. Colton, and H. Haddar , Inverse Scattering Theory and Transmission Eigenvalues, Second Edition, Society for Industrial and Applied Mathematics, Philadelphia, PA, 2022, https://doi. org/10.1137/1.9781611977424

    work page doi:10.1137/1.9781611977424 2022

  9. [9]

    Cakoni and H

    F. Cakoni and H. Haddar , On the existence of transmission eigenvalues in an inhomogeneous medium, Applicable Analysis, 88 (2009), pp. 475–493, https://doi.org/10.1080/00036810802713966

    work page doi:10.1080/00036810802713966 2009

  10. [10]

    Cakoni, H

    F. Cakoni, H. Haddar, and D. Gintides , The existence of an infinite discrete set of transmission eigenvalues, SIAM Journal on Mathematical Analysis, 42 (2010), pp. 237–255, https://doi.org/10. 1137/090769338

    work page 2010

  11. [11]

    Colton and Y.-J

    F. Cakoni and R. Kress , A boundary integral equation method for the transmission eigenvalue problem, Applicable Analysis, 96 (2017), pp. 23–38, https://doi.org/10.1080/00036811.2016.1189537

    work page doi:10.1080/00036811.2016.1189537 2017

  12. [12]

    Colton, A

    D. Colton, A. Kirsch, and L. Päivärinta , Far-field patterns for acoustic waves in an inhomogeneous medium, SIAM Journal on Mathematical Analysis, 20 (1989), pp. 1472–1483, https://doi.org/10. 1137/0520096

    work page 1989

  13. [13]

    Colton and Y.-J

    D. Colton and Y.-J. Leung , The existence of complex transmission eigenvalues for spherically stratified media, Applicable Analysis, 96 (2017), pp. 39–47, https://doi.org/10.1080/00036811.2016.1210788

    work page doi:10.1080/00036811.2016.1210788 2017

  14. [14]

    Colton, Y.-J

    D. Colton, Y.-J. Leung, and S. Meng , Distribution of complex transmission eigenvalues for spher- ically stratified media, Inverse problems, 31 (2015), 035006 (19 pages), https://doi.org/10.1088/ 0266-5611/31/3/035006

    work page 2015

  15. [15]

    Colton and P

    D. Colton and P. Monk , The inverse scattering problem for time-harmonic acoustic waves in an inhomogeneous medium, The Quarterly Journal of Mechanics and Applied Mathematics, 41 (1988), pp. 97–125, https://doi.org/10.1093/qjmam/41.1.97

    work page doi:10.1093/qjmam/41.1.97 1988

  16. [16]

    Colton, P

    D. Colton, P. Monk, and J. Sun , Analytical and computational methods for transmission eigenvalues, Inverse Problems, 26 (2010), 045011 (16 pages), https://doi.org/10.1088/0266-5611/26/4/045011

    work page doi:10.1088/0266-5611/26/4/045011 2010

  17. [17]

    Colton, L

    D. Colton, L. Päiv ärinta, and J. Syl vester , The interior transmission problem, Inverse Problems and Imaging, 1 (2007), pp. 13–28, https://doi.org/10.3934/ipi.2007.1.13

    work page doi:10.3934/ipi.2007.1.13 2007

  18. [18]

    Cossonnière, Valeurs propres de transmission et leur utilisation dans l’identification d’inclusions à partir de mesures électromagnétiques, PhD thesis, Université de Toulouse, 2011

    A. Cossonnière, Valeurs propres de transmission et leur utilisation dans l’identification d’inclusions à partir de mesures électromagnétiques, PhD thesis, Université de Toulouse, 2011

    work page 2011

  19. [19]

    Cossonnière and H

    A. Cossonnière and H. Haddar , Surface integral formulation of the interior transmission problem, Journal of Integral Equations and Applications, 25 (2013), pp. 341–376, https://doi.org/10.1216/ JIE-2013-25-3-341

    work page 2013

  20. [20]

    Hickmann , Interior transmission eigenvalue problem with refractive index havingC 2-transition to the background medium, Applicable Analysis, 91 (2012), pp

    K. Hickmann , Interior transmission eigenvalue problem with refractive index havingC 2-transition to the background medium, Applicable Analysis, 91 (2012), pp. 1675–1690, https://doi.org/10.1080/ 00036811.2011.577741

    work page arXiv 2012

  21. [21]

    Hryniv and P

    R. Hryniv and P. Lancaster , On the perturbation of analytic matrix functions, Integral Equations and Operator Theory, 34 (1999), pp. 325–338, https://doi.org/10.1007/BF01300582, https://doi.org/ 10.1007/BF01300582

    work page doi:10.1007/bf01300582 1999

  22. [22]

    Ji and J

    X. Ji and J. Sun , A multi-level method for transmission eigenvalues of anisotropic media, Journal of Computational Physics, 255 (2013), pp. 422–435, https://doi.org/10.1016/j.jcp.2013.08.030

    work page doi:10.1016/j.jcp.2013.08.030 2013

  23. [23]

    Perturbation Theory for Linear Operators

    T. Kato, Perturbation Theory for Linear Operators, Springer Berlin, Heidelberg, 1995, https://doi.org/ 10.1007/978-3-642-66282-9

    work page doi:10.1007/978-3-642-66282-9 1995

  24. [24]

    Kirsch, The denseness of the far field patterns for the transmission problem, IMA Journal of Applied Mathematics, 37 (1986), pp

    A. Kirsch, The denseness of the far field patterns for the transmission problem, IMA Journal of Applied Mathematics, 37 (1986), pp. 213–225, https://doi.org/10.1093/imamat/37.3.213

    work page doi:10.1093/imamat/37.3.213 1986

  25. [25]

    Kirsch and F

    A. Kirsch and F. Hettlich , The Mathematical Theory of Time-Harmonic Maxwell’s Equations: Expansion-, Integral-, and Variational Methods, vol. 190 of Applied Mathematical Sciences, Springer International Publishing, Cham, 2015, https://doi.org/10.1007/978-3-319-11086-8

    work page doi:10.1007/978-3-319-11086-8 2015

  26. [26]

    Kleefeld, A numerical method to compute interior transmission eigenvalues, Inverse Problems, 29 (2013), 104012 (20 pages), https://doi.org/10.1088/0266-5611/29/10/104012

    A. Kleefeld, A numerical method to compute interior transmission eigenvalues, Inverse Problems, 29 (2013), 104012 (20 pages), https://doi.org/10.1088/0266-5611/29/10/104012. BIFURCATIONS IN INTERIOR TRANSMISSION EIGENVALUES 27

    work page doi:10.1088/0266-5611/29/10/104012 2013

  27. [27]

    Kleefeld and L

    A. Kleefeld and L. Pieronek , Computing interior transmission eigenvalues for homogeneous and an- isotropic media, Inverse Problems, 34 (2018), p. 105007, https://doi.org/10.1088/1361-6420/aad7c4

    work page doi:10.1088/1361-6420/aad7c4 2018

  28. [28]

    Kleefeld and L

    A. Kleefeld and L. Pieronek , The method of fundamental solutions for computing acoustic interior transmission eigenvalues, Inverse Problems, 34 (2018), 035007 (28 pages), https://doi.org/10.1088/ 1361-6420/aaa72d

    work page 2018

  29. [29]

    Kubrusly, Spectral theory of operators on Hilbert spaces, Birkhauser, New York, 2012

    C. Kubrusly, Spectral theory of operators on Hilbert spaces, Birkhauser, New York, 2012

    work page 2012

  30. [30]

    N. N. Lebedev and R. A. Sil verman , Special functions and their applications, Courier Corporation, 1972

    work page 1972

  31. [31]

    Leung and D

    Y.-J. Leung and D. Colton , Complex transmission eigenvalues for spherically stratified media, Inverse Problems, 28 (2012), 075005 (9 pages), https://doi.org/10.1088/0266-5611/28/7/075005

    work page doi:10.1088/0266-5611/28/7/075005 2012

  32. [32]

    Päiv ärinta and J

    L. Päiv ärinta and J. Syl vester, Transmission eigenvalues, SIAM Journal on Mathematical Analysis, 40 (2008), pp. 738–753, https://doi.org/10.1137/070697525

    work page doi:10.1137/070697525 2008

  33. [33]

    Peters and A

    S. Peters and A. Kleefeld , Numerical computations of interior transmission eigenvalues for scatter- ing objects with cavities, Inverse Problems, 32 (2016), 045001 (28 pages), https://doi.org/10.1088/ 0266-5611/32/4/045001

    work page 2016

  34. [34]

    Pieronek, The method of fundamental solutions for computing interior transmission eigenvalues, PhD thesis, BTU Cottbus-Senftenberg, 2020

    L. Pieronek, The method of fundamental solutions for computing interior transmission eigenvalues, PhD thesis, BTU Cottbus-Senftenberg, 2020

    work page 2020

  35. [35]

    Pieronek and A

    L. Pieronek and A. Kleefeld , On trajectories of complex-valued interior transmission eigenvalues, Inverse Problems and Imaging, 18 (2024), pp. 480–516, https://doi.org/10.3934/ipi.2023041

    work page doi:10.3934/ipi.2023041 2024

  36. [36]

    Pradovera and A

    D. Pradovera and A. Borghi ,Match-based solution of general parametric eigenvalue problems, Journal of Computational Physics, 519 (2024), 113384 (19 pages), https://doi.org/10.1016/j.jcp.2024.113384

    work page doi:10.1016/j.jcp.2024.113384 2024

  37. [37]

    Bifurcations in Interior Transmission Eigenvalues: Theory and Computation

    D. Pradovera, A. Borghi, L. Pieronek, and A. Kleefeld , Code and results for numerical experi- ments in “Bifurcations in Interior Transmission Eigenvalues: Theory and Computation” (version v1), November 2025, https://doi.org/10.5281/zenodo.17597807

    work page doi:10.5281/zenodo.17597807 2025

  38. [38]

    Rellich, Perturbation theory of eigenvalue problems, Gordon and Breach, 1969

    F. Rellich, Perturbation theory of eigenvalue problems, Gordon and Breach, 1969. Series Title: Notes on mathematics and its applications

    work page 1969

  39. [39]

    B. P. Rynne and B. D. Sleeman , The interior transmission problem and inverse scattering from inhomogeneous media, SIAM Journal on Mathematical Analysis, 22 (1991), pp. 1755–1762, https: //doi.org/10.1137/0522109

    work page doi:10.1137/0522109 1991

  40. [40]

    N. G. Stephen , On veering of eigenvalue loci, Journal of Vibration and Acoustics, 131 (2009), https: //doi.org/10.1115/1.3147130

    work page doi:10.1115/1.3147130 2009

  41. [41]

    Sun and A

    J. Sun and A. Zhou , Finite Element Methods for Eigenvalue Problems, CRC Press, Boca Raton, 2017

    work page 2017

  42. [42]

    Xu, C.-F

    X.-C. Xu, C.-F. Yang, S. A. Buterin, and V. A. Yurko , Estimates of complex eigenvalues and an inverse spectral problem for the transmission eigenvalue problem, Electronic Journal of Qualitative Theory of Differential Equations, 38 (2019), pp. 1–15, https://doi.org/10.14232/ejqtde.2019.1.38

    work page doi:10.14232/ejqtde.2019.1.38 2019

  43. [43]

    F. Zeng, J. Sun, and L. Xu , A spectral projection method for transmission eigenvalues, Science China Mathematics, 59 (2016), pp. 1613–1622, https://doi.org/10.1007/s11425-016-0289-8

    work page doi:10.1007/s11425-016-0289-8 2016