Abstract:This paper presents a novel algorithm, based on use of rational approximants of a randomly scalarized boundary integral resolvent in conjunction with an adaptive search strategy and an exponentially convergent secant-method termination stage, for the evaluation of acoustic and electromagnetic resonances in open and closed cavities. The desired cavity resonances are obtained as the poles of associated rational approximants; both the approximants and their poles are obtained by means of the recently introduced AAA rational-approximation algorithm. In fact, the proposed resonance-search method applies to any nonlinear eigenvalue problem associated with a given function $F: U \to \mathbb{C}^{d\times d}$, wherein, denoting $F(k) = F_k$, a complex value $k$ is sought for which $F_kw = 0$ for some nonzero $w\in \mathbb{C}^d$. For the scattering problems considered in this paper, $F_k$ is taken to equal a spectrally discretized version of a Green function-based boundary integral operator at spatial frequency $k$. In all cases, the scalarized resolvent is given by an expression of the form $u^* F_k^{-1} v$, where $u,v \in \mathbb{C}^d$ are fixed random vectors. The proposed adaptive search strategy relies on use of a rectangular subdivision of the resonance search domain which is locally refined to ensure that all resonances in the domain are captured. The approach works equally well in the case in which the search domain is an interval of the real line, in which case the rectangles used degenerate into subintervals of the search domain. A variety of numerical results are presented, including comparisons with well-known methods based on complex contour integration, and a discussion of the asymptotics that result as open cavities approach closed cavities -- in all, demonstrating the accuracy provided by the method, for low- and high-frequency states alike.

Efficient rational approximation of optical response functions with the AAA algorithm

Representation of Complex Optical Surfaces with Adaptive Radial Basis Functions

Description Of Complex Optical Surface Based On Sparse Radial Basis Function Approximation With Spatially Variable Shape Parameters

A framework to compute resonances arising from multiple scattering

Modal-based Nonlinear Optimization Algorithm for Wavefront Measurement with Under-Sampled Data.

Evaluation of resonances: adaptivity and AAA rational approximation of randomly scalarized boundary integral resolvents

Adaptive Finite Element Method for Simulation of Optical Nano Structures

Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy

Invited Article: Acousto-Optic Finite-Difference Frequency- Domain Algorithm for First-Principles Simulations of On-Chip Acousto-Optic Devices

Accurate near-field calculation in the rigorous coupled-wave analysis method

Highly Efficient Rigorous Coupled-Wave Analysis Implementation Without Eigensystem Calculation

Efficient numerical method for analyzing optical bistability in photonic crystal microcavities

Efficient Nonlinear Metasurfaces using Multiresonant High-Q Plasmonic Arrays

Data driven design of optical resonators

Computing resonant modes of accelerator cavities by solving nonlinear eigenvalue problems via rational approximation

Rational Design of On-Chip Refractive Index Sensors Based on Lattice Plasmon Resonances (presentation Recording)

A simulation method for determining the optical response of highly complex photonic structures of biological origin

Efficient and Parallel Implementation of Real and Complex Response Functions Employing the Second-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator

A Multi-Scale Approach to Simulate the Nonlinear Optical Response of Molecular Nanomaterials

Three-Dimensional Modeling of the Optical Switch Based on Guided-Mode Resonances in Photonic Crystals

Finite-Difference Time-Domain Algorithm For Modeling Sagnac Effect In Rotating Optical Elements