Source-linked AI summary
Waveguiding in systems of high contrast resonators: Theory and fast computations
Habib Ammari, Bowen Li, Borui Miao, Jiayu Qiu, Lara Vrabac
TL;DR
Waveguiding near nonzero interior Neumann frequencies needs an effective description beyond the subwavelength regime. The paper constructs a frequency-dependent capacitance operator, proves its spectral and locality properties, and develops stabilized local computations that retain first-order accuracy. The framework applies to straight and bent waveguides, while exterior-spectrum and real-frequency stability questions remain open.
Problem
Waveguiding beyond the subwavelength regime involves multidimensional interior resonant eigenspaces, Helmholtz exterior fields, and leading-order radiation effects requiring a different effective model.
Method
The paper compresses the exterior Helmholtz DtN map onto interior resonant Neumann traces, proves norm-resolvent convergence and exponential coefficient decay, and approximates retained interactions with stabilized local Helmholtz problems.
Results
γ = √δ and interaction and patch radii N, R = O(|log δ|) preserve O(δ^2) accuracy of the first-order effective model for defect eigenfrequencies.
Takeaways & Limitations
The effective operator is exponentially compressible and locally computable, supporting fast modeling of nonperiodic configurations such as bent waveguides.
Takeaways & Limitations
The analysis excludes cases where the exterior wavenumber lies in the Dirichlet spectrum and leaves direct real-frequency stability under minimal assumptions open.
Abstract
from arXiv · showhide
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 $δ\to0$, 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 $\sqrtδ$, 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.
1. Introduction
The paper develops a frequency-dependent discrete framework for waveguiding near nonzero interior Neumann frequencies, where resonant eigenspaces may have multiple degrees of freedom and exterior fields are governed by Helmholtz equations. It proves effective reduction, exponential locality, and locally computable approximations with controlled accuracy.
- Motivation: Non-subwavelength resonances require a frequency-dependent model acting on possibly multidimensional interior Neumann eigenspaces.The exterior problem is Helmholtz rather than harmonic, and radiation effects may contribute at leading order.
- Effective model: The frequency-dependent capacitance operator is obtained by compressing the exterior Helmholtz DtN map to traces of interior resonant Neumann eigenspaces.For periodic structures, quasi-periodic capacitance matrices determine leading-order Bloch band asymptotics.
- Spectral reduction: Norm-resolvent reduction shows that the effective operator determines first-order asymptotics for compact-defect frequencies and line-defect guided bands.The active resonant spaces can produce a block operator, and line defects are treated fiberwise through quasi-periodic capacitance matrices.
- Exponential compression: The capacitance coefficients decay exponentially off diagonal, so interaction truncation with radius N = O(log ε^-1) achieves operator accuracy O(ε).For a finite line-like defect, only O(Nact log ε^-1) blocks are retained instead of O(N^2).
- Local computation: Retained coefficients can be computed from finite local Helmholtz patches, with exponential convergence under uniform stability of growing real-frequency clusters.Complex absorption and Hermitian symmetrization remove this stability requirement; γ = √δ and N, R = O(|log δ|) preserve O(δ^2) accuracy.
- Numerical illustrations: Numerical experiments for dipole and quadrupole resonances show accurate straight-waveguide bands, exponential interaction decay, finite-patch convergence, and applicability to sharply bent waveguides.The framework does not require global periodicity and supports material or geometric detuning.
2. High contrast crystal with defects
The paper formulates high-contrast resonator systems with compact or extended defects without requiring global periodicity, then eliminates the exterior field through a Dirichlet-to-Neumann map.
- Defect configurations: The framework allows compact, straight line, and completely nonperiodic defects, including shape and material-parameter perturbations.Defect sets are classified by whether they are finite, lower-dimensional periodic, or nonperiodic.
- Whole-space problem: The whole-space spectral problem is expressed using transmission conditions across resonator boundaries, including continuity of the field and scaled normal flux.The normal orientation and exterior/interior trace conventions are specified on each resonator boundary.
- Exterior reduction: Eliminating the exterior field produces an interior problem on δ-independent spaces and an operator pencil Aδ(z).This reformulation avoids the δ-dependent domain of the original whole-space operator.
- Exterior reduction: The exterior Dirichlet-to-Neumann map is defined for frequencies whose exterior wavenumber lies outside the Dirichlet spectrum.For real frequencies in this resolvent set, the map is self-adjoint under the boundary duality pairing.
3. Frequency-dependent capacitance operator and spectral asymptotics
The frequency-dependent capacitance operator reduces waveguiding near nonzero interior Neumann frequencies to active resonant degrees of freedom, including block-valued interactions and quasi-periodic line-defect matrices.
- Frequency-dependent capacitance operator: The discrete model acts on interior Neumann eigenspaces at a reference frequency, allowing several active degrees of freedom per resonator.Active sites are those with nontrivial eigenspaces, and their multiplicities determine the block structure.
- Frequency assumptions: The exterior gap condition places the reference exterior wavenumber in the resolvent set of the exterior Dirichlet operator.This condition supports analyticity and self-adjointness properties needed for the capacitance construction.
- Frequency-dependent capacitance operator: The capacitance operator is built by compressing exterior Helmholtz interactions to traces of the interior resonant Neumann eigenfunctions.Each resonant mode is extended through an exterior Helmholtz problem with prescribed boundary traces.
- Line defects: For a straight line defect, partial Floquet–Bloch transformation converts the full capacitance operator into Hermitian quasi-periodic matrices Cα(ω0).The matrices act on the resonant modes of the defect resonator and are related to the full operator by direct-integral decomposition.
The unperturbed bulk crystal.
The unperturbed bulk crystal is described through a quasi-periodic capacitance matrix whose spectral data determine leading bulk Bloch-frequency asymptotics.
- Bulk capacitance matrix: The bulk capacitance matrix is defined from an orthonormal basis of the relevant interior Neumann eigenspace of a reference resonator.Its quasi-periodic formulation is indexed by α in the Brillouin torus.
- Bulk spectral asymptotics: The bulk Bloch eigenfrequencies satisfy a first-order expansion with O(δ^2) remainder uniformly in α.The expansion is used to compute the bulk bands and to illustrate exponential decay of capacitance interactions.
- Bulk spectral asymptotics: The real-space bulk capacitance operator is recovered from its quasi-periodic symbol through the inverse Floquet–Bloch transform.This connects the bulk matrix representation with the operator’s spatial interaction coefficients.
3.3. Norm-resolvent convergence and spectral asymptotics
The continuous high-contrast problem converges, after reduction to interior resonant modes, to a frequency-dependent capacitance operator that governs leading spectral behavior.
- Spectral asymptotics: The eigenvalues of C(ω0) determine leading-order shifts of compact-defect eigenfrequencies, while the corresponding periodic formulation governs line-defect band asymptotics.The reduction applies under the interior-resonance and exterior-gap conditions.
- Norm-resolvent convergence: The rescaled resolvent δR(ω0^2+δζ,δ) converges in operator norm to the effective resolvent involving 2ω0C(ω0).The convergence holds uniformly for ζ in compact subsets of the resolvent set of the capacitance operator.
- Block reduction: The effective degrees of freedom are the interior-resonant Neumann modes, while the complementary block is generally indefinite at nonzero frequencies.Uniform invertibility of that complementary block is established using spectral calculus rather than coercivity.
- Spectral decomposition: The resonant-mode projection and its complement provide the decomposition used to formulate the reduced operator on ℓ2(I(ω0)) and Q.The complementary operator is boundedly invertible under the stated spectral condition.
- Block reduction: The proof decomposes the interior space into resonant and complementary components, estimates the four operator blocks, and eliminates the complementary block by a Schur complement.This block procedure yields the effective discrete resolvent after inversion.
3.4. Characterizations of defect eigenvalues and bands
The section characterizes compact-defect eigenfrequencies and line-defect bands through frequency-dependent capacitance operators, and places them in bulk spectral gaps for sufficiently small contrast.
- A bulk spectral gap opens around ω0^2 for sufficiently small contrasts under the exterior gap condition.
- Compact defects: For a compact defect, an eigenvalue λ of C(ω0) with multiplicity m produces exactly m eigenfrequencies with ωδ,k = ω0 + δλ + O(δ^2).
- Compact defects: These compact-defect eigenfrequencies lie in a spectral gap of the unperturbed bulk crystal.
- Line defects: For a straight line defect, quasi-periodic capacitance matrices yield m* real, continuous band functions bifurcating from ω0 uniformly in the quasi-momentum.
- Line defects: The line-defect bands lie in a projected bulk gap uniformly in α and j.
3.5. Exponential decay of the off-diagonal coefficients
The frequency-dependent capacitance operator is exponentially local: the exterior spectral gap causes interactions between distant resonators to decay exponentially.
- The capacitance coefficient measuring flux on one resonator from a source on another decays exponentially with their separation.This locality follows from a Combes–Thomas estimate for the exterior Dirichlet resolvent.
- The exponential bound holds uniformly across resonator indices, while same-resonator coefficients remain uniformly bounded.
4. Patch approximation and local computation
Exponential locality enables finite-range and local-patch approximations of the capacitance operator. Complex absorption and Hermitian symmetrization remove the need for real-frequency finite-cluster stability while preserving the first-order eigenfrequency accuracy.
- Complex-frequency stabilization: A vanishing complex absorption makes finite-cluster outgoing problems uniformly invertible, while Hermitian symmetrization removes the first-order regularization bias.
- Operator compression: For a straight line defect, truncation produces a block-banded operator of bandwidth N with exponentially small operator error.
- Real-frequency patches: At ω0, patch approximations converge exponentially in truncation radius N and patch radius R when growing local Helmholtz problems satisfy a uniform resolvent bound.
- Error and complexity: Interaction truncation and patch radii of order |log δ| preserve the O(δ^2) remainder in first-order defect-eigenfrequency expansions.
- Error and complexity: For finite straight or bent line defects, storage and matrix–vector costs are linear in Nact up to a polylogarithmic factor, excluding local PDE discretization costs.
5. Numerical results
The numerical results construct frequency-dependent capacitance matrices for dipole and quadrupole resonances, then compare their bands, locality, and patch approximations with full formulations. Across straight and bent waveguides, the discrete model shows close band agreement, exponential decay, and localized modal fields.
- Operator construction: Dipole and quadrupole resonances have multiplicity two, so a line waveguide uses a 2 × 2 Floquet symbol for the active modes.The symbol incorporates all longitudinal interactions through its quasimomentum dependence.
- Operator construction: The effective capacitance operator retains only active resonant modes, while inactive cladding and fringe resonators remain auxiliary unknowns in the boundary-integral system.For the illustrated finite patch, five active resonators produce a 5 × 5 block capacitance matrix.
- Operator construction: The quasi-periodic capacitance matrix is computed with a supercell containing finitely many cladding resonators, then transformed into the frequency-dependent operator.The numerical setup uses 64 quadrature points per resonator, 31 quasimomenta, 13 projected quasimomenta, and δ = 0.001.
- Tight-binding approximation: The capacitance operator is a block Laurent operator whose entries depend on longitudinal separation and are exponentially approximable by banded operators.This structure supports truncating long-range interactions in the effective model.
- Straight waveguides: The capacitance-based dipole and quadrupole defect bands agree excellently with the full characteristic-value results, producing two bands for each resonance type.The defect frequencies are selected deep inside bulk band gaps.
- Geometric detuning: For radius-detuned resonators, the bands span larger frequency intervals, while the quadrupole bands deviate at lower frequencies farthest from ω0 = 6.7143.This comparison identifies the range where the first-order frequency-dependent approximation is less accurate.
- Locality and bent waveguides: Patch and inverse Floquet–Bloch computations show exponential decay and nearly perfect overlap for both dipole and quadrupole interactions, alongside localized waveguide fields.The same behavior is reported for the bent-waveguide examples.
6. Concluding remarks
The paper develops a discrete framework for guided modes near nonzero interior Neumann frequencies and establishes its accuracy, locality, and computational applicability. Numerical results support the framework for dipole and quadrupole resonances, while several extensions remain open.
- Concluding framework: The continuous transmission problem reduces to a frequency-dependent capacitance operator acting on active interior Neumann modes.Norm-resolvent approximation connects compact-defect frequencies and line-defect band functions to the effective operator and its partial Floquet–Bloch transform.
- Computational structure: Exponential locality of the effective interaction enables exponentially accurate truncation and local patch approximation under a uniform stability assumption.The exterior spectral gap drives exponential decay of capacitance coefficients and supports local Helmholtz computations.
- Numerical evidence: Numerical experiments demonstrate accurate first-order approximations and rapid spatial decay for dipole and quadrupole resonances.Material and geometric detuning are treated within a common active-mode framework; finite outgoing computations produce small radiative imaginary parts distinct from the real infinite-operator spectrum.
- Open problems: A remaining scope boundary is the singular exterior case, where the exterior DtN map is singular and the leading splitting is expected to be of order δ1/2.Extending the framework to line defects in this regime is identified as an open problem.
- Open problems: Direct real-frequency convergence under minimal assumptions remains open because the uniform stability assumption for outgoing finite-cluster approximation has not been removed or verified.Complex-frequency stabilization provides an unconditional computational alternative.
A.1. Boundary integral operators
This section recalls boundary-integral ingredients for Helmholtz problems, including outgoing Green functions, layer potentials, boundary operators, jump relations, and quasi-periodic representations.
- Boundary integral operators: The boundary-integral formulation uses the outgoing Helmholtz Green function in two dimensions, expressed through the Hankel function of the first kind and order zero.The Green function is associated with the single- and double-layer constructions.
- Boundary integral operators: Single-layer and double-layer potentials, together with the Neumann–Poincaré operator, are introduced for boundary densities on the resonator boundary.The operators act between boundary function spaces including L2 and H1 spaces.
- Quasi-periodic formulation: For quasi-periodic problems, the frequency excludes reciprocal-lattice thresholds, and the quasi-periodic Green function admits an augmented-plane-wave representation via Poisson summation.The excluded frequencies satisfy ω ≠ |q + α| for reciprocal-lattice vectors q.
- Boundary relations: The formulation relies on jump relations for the double-layer potential and the normal derivative of the single-layer potential.These relations are stated on the boundary and support the associated boundary operators.
A.2. Characterizations of defect and band defect modes and their associated eigenmodes
This section characterizes defect and band-defect modes through operator-valued characteristic values and relates their eigenmodes to frequency-dependent capacitance matrices. It also distinguishes material and radius detuning cases.
- Detuning cases: Material and radius detuning are treated through distinct operator formulations, with radius detuning replacing the reference domain by a scaled domain and associated boundary maps.The appendix covers these two pure detuning cases, while simultaneous shape-and-parameter detuning is not used in the legacy formulas.
- Defect modes: Defect frequencies bifurcating from ω0 are characterized as characteristic values of an operator-valued function.The formulation applies within the band gap around ω0 and uses an operator Aα(ω,δ).
- Band-defect modes: For a line defect, α1-quasi-periodic Bloch eigenfrequencies bifurcating from ω0 are characterized similarly by an α1-quasi-periodic operator-valued function.The characterization is made separately for each α1.
- Associated eigenmodes: Eigenmodes are approximated by eigenvectors of the frequency-dependent capacitance matrix, and outgoing patch convergence yields finite-cluster reconstruction under Theorem 4.3 assumptions.The root functions associated with the characteristic values represent the defect and band-defect modes.