Waveguiding in systems of high contrast resonators: Theory and fast computations
Abstract
In this work, we study guided modes in systems of high-contrast resonators near nonzero interior Neumann frequencies, beyond the subwavelength regime. In the regular exterior regime, where the exterior Dirichlet problem is well-posed at the reference wavenumber, we introduce an infinite-dimensional frequency-dependent capacitance operator obtained by compressing the exterior Helmholtz Dirichlet-to-Neumann map to the traces of the interior resonant Neumann eigenspaces. We prove the norm-resolvent convergence of the continuous problem to this discrete effective operator as the contrast δ→0, and derive first-order asymptotic formulas for compact-defect frequencies and line-defect band functions. We then establish exponential off-diagonal decay of the capacitance coefficients by a Combes--Thomas argument, yielding an exponentially accurate truncation of the discrete operator, and show that its retained coefficients can be computed from local Helmholtz problems. At the physical frequency, this local approximation converges exponentially under a uniform stability assumption for the growing finite-cluster problems. The stability assumption can be removed by introducing a vanishing complex absorption together with a Hermitian symmetrization. In particular, an absorption parameter of order √δ, together with interaction truncation and patch radii of order |logδ|, suffices to preserve the O(δ2) accuracy of the first-order high-contrast expansion of the defect eigenfrequencies, yielding a fast computational method. Numerical experiments for dipole and quadrupole resonances illustrate the accuracy, exponential locality, and applicability of the discrete model to straight and bent waveguides generated by material or geometric detuning.
Related
- Solver Agent: an Agentic AI Framework for Theoretical Physics Computations Applied to F-theory Uplifts of O3-planes and S-folds
- Dissipative Systems Analysis and Control. Theory and Applications: Erratum/Addendum. Springer Nature Switzerland AG, Communications and Control Engineering, Third uncorrected Edition, 2020
- Spectral computations on lamplighter groups and Diestel-Leader graphs
- Delta Redistribution: a variation of State Redistribution for steady-state embedded boundary computations
- Linear Response Theory for Jump-Driven Stochastic Systems: Transient Statistics and Escape Dynamics
BibTeX
@misc{ammari2026waveguiding,
title = {{Waveguiding in systems of high contrast resonators: Theory and fast computations}},
author = {Habib Ammari and Bowen Li and Borui Miao and Jiayu Qiu and Lara Vrabac},
howpublished = {arXiv},
year = {2026},
month = aug,
doi = {10.48550/arxiv.2608.26906},
eprint = {2608.26906},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.26906},
}
RIS
TY - UNPB
TI - Waveguiding in systems of high contrast resonators: Theory and fast computations
AU - Habib Ammari
AU - Bowen Li
AU - Borui Miao
AU - Jiayu Qiu
AU - Lara Vrabac
PY - 2026
DA - 2026/08/31
PB - arXiv
DO - 10.48550/arxiv.2608.26906
UR - https://arxiv.org/abs/2608.26906
AB - In this work, we study guided modes in systems of high-contrast resonators near nonzero interior Neumann frequencies, beyond the subwavelength regime. In the regular exterior regime, where the exterior Dirichlet problem is well-posed at the reference wavenumber, we introduce an infinite-dimensional frequency-dependent capacitance operator obtained by compressing the exterior Helmholtz Dirichlet-to-Neumann map to the traces of the interior resonant Neumann eigenspaces. We prove the norm-resolvent convergence of the continuous problem to this discrete effective operator as the contrast δ→0, and derive first-order asymptotic formulas for compact-defect frequencies and line-defect band functions. We then establish exponential off-diagonal decay of the capacitance coefficients by a Combes--Thomas argument, yielding an exponentially accurate truncation of the discrete operator, and show that its retained coefficients can be computed from local Helmholtz problems. At the physical frequency, this local approximation converges exponentially under a uniform stability assumption for the growing finite-cluster problems. The stability assumption can be removed by introducing a vanishing complex absorption together with a Hermitian symmetrization. In particular, an absorption parameter of order √δ, together with interaction truncation and patch radii of order |logδ|, suffices to preserve the O(δ^2) accuracy of the first-order high-contrast expansion of the defect eigenfrequencies, yielding a fast computational method. Numerical experiments for dipole and quadrupole resonances illustrate the accuracy, exponential locality, and applicability of the discrete model to straight and bent waveguides generated by material or geometric detuning.
ER -