Search results

In this paper we consider scattering resonance computations in optics when the resonators consist of frequency dependent and lossy materials, such as metals at optical frequencies. The proposed computational approach combines a novel hp-FEM strategy, based on dispersion analysis for complex frequencies, with a fast implementation of the nonlinear eigenvalue solver NLEIGS. Numerical computations il

This article reviews both recent progress on the mathematics of dispersive and absorptive photonic crystals and well-established results on conservative photonic crystals. The focus is on properties of the photonic band structures and we also provide results that are of importance for the understanding of lossy metal-dielectric photonic crystals.

For analytic operator functions, we prove accumulation of branches of complex eigenvalues to the essential spectrum. Moreover, we show minimality and completeness of the corresponding system of eigenvectors and associated vectors. These results are used to prove sufficient conditions for eigenvalue accumulation to the poles and to infinity of rational operator functions. Finally, an application of

We present an efficient procedure for computing resonances and resonant modes of Helmholtz problems posed in exterior domains. The problem is formulated as a nonlinear eigenvalue problem (NEP), where the nonlinearity arises from the use of a Dirichlet-to-Neumann map, which accounts for modeling unbounded domains. We consider a variational formulation and show that the spectrum consists of isolated

In this paper we present equivalence results for several types of unbounded operator functions. A generalization of the concept equivalence after extension is introduced and used to prove equivalence and linearization for classes of unbounded operator functions. Further, we deduce methods of finding equivalences to operator matrix functions that utilizes equivalences of the entries. Finally, a met

In this paper, we discuss problems arising when computing resonances with a finite element method. In the pre-asymptotic regime, we detect for the one dimensional case, spurious solutions in finite element computations of resonances when the computational domain is truncated with a perfectly matched layer (PML) as well as with a Dirichlet-to-Neumann map (DtN). The new test is based on the Lippmann

In this paper we introduce an enclosure of the numerical range of a class of rational operator functions. In contrast to the numerical range the presented enclosure can be computed exactly in the infinite dimensional case as well as in the finite dimensional case. Moreover, the new enclosure is minimal given only the numerical ranges of the operator coefficients and many characteristics of the num

We establish new analytic results for a general class of rational spectral problems. They arise e.g. in modelling photonic crystals whose capability to control the flow of light depends on specific features of the eigenvalues. Our results comprise a complete spectral analysis including variational principles and two-sided bounds for all eigenvalues, as well as numerical implementations. They apply

We present a-posteriori analysis of higher order finite element approximations (hp-FEM) for quadratic Fredholm-valued operator functions. Residual estimates for approximations of the algebraic eigenspaces are derived and we reduce the analysis of the estimator to the analysis of an associated boundary value problem. For the reasons of robustness we also consider approximations of the associated in

We present an algorithm for approximating an eigensubspace of a spectral component of an analytic Fredholm valued function. Our approach is based on numerical contour integration and the analytic Fredholm theorem. The presented method can be seen as a variant of the FEAST algorithm for infinite dimensional nonlinear eigenvalue problems. Numerical experiments illustrate the performance of the algor

Metallic nano-antennas are devices used to concentrate the energy in light into regions that are much smaller than the wavelength. These structures are currently used to develop new measurement and printing techniques, such as optical microscopy with sub-wavelength resolution, and high-resolution lithography. Here, we analyze and design a nano-antenna in a two-dimensional setting with the source b

Galerkin spectral approximation theory for non-self-adjoint quadratic operator polynomials with periodic coefficients is considered. The main applications are complex band structure calculations in metallic photonic crystals, periodic waveguides, and metamaterials. We show that the spectrum of the considered operator polynomials consists of isolated eigenvalues of finite multiplicity with a nonzer

The meshless radial point interpolation method (RPIM) in frequency domain for electromagnetic scattering problems is presented. This method promises high accuracy in a simple collocation approach using radial basis functions. The treatment of high-order non-reflecting boundary conditions for open waveguides is discussed and implemented up to fourth-order. RPIM allows the direct calculation of high

Meshless methods are a promising new field in computational electromagnetics. Instead of relying on an explicit mesh topology, a numerical solution is computed on an unstructured set of collocation nodes. This allows to model fine geometrical details with high accuracy and facilitates the adaptation of node distributions for optimization or refinement purposes. The radial point interpolation metho

Band structure calculations for photonic crystals require the numerical solution of eigenvalue problems. In this paper, we consider crystals composed of lossy materials with frequency-dependent permittivities. Often, these frequency dependencies are modeled by rational functions, such as the Lorentz model, in which case the eigenvalue problems are rational in the eigenvalue parameter. After spatia

In this paper a high-order finite element method with curvilinear elements is proposed for the simulation of plasmonic structures. Most finite element packages use low order basis functions and non-curved elements, which is very costly for demanding problems such as the simulation of nanoantennas. To enhance the performance of finite elements, we use curvilinear quadrilateral elements to calculate

Meshless methods are numerical methods that have the advantage of high accuracy without the need of an explicitly described mesh topology. In this class of methods, the Radial Point Interpolation Method (RPIM) is a promising collocation method where the application of radial basis functions yields high interpolation accuracy for even strongly unstructured node distributions. For electromagnetic si

Three different meshless methods based on radial basis functions are investigated for the numerical solution of electromagnetic eigenvalue problems. The three algorithms, the non-symmetric Kansa approach, the symmetric Kansa method and the radial point interpolation method, are first described putting emphasis on the influence of their formalism on practical implementation. The convergence rate of

A meshless method based on a radial basis collocation approach is presented to calculate eigenvalues for the second-order wave equation. Instead of an explicit mesh topology only a node distribution is required to calculate electric fields, thus facilitating dynamic alteration of the discretization of an electromagnetic problem. An algorithm is presented that automatically adapts an initially very

The concept of an adaptive meshless eigenvalue solver is presented and implemented for two-dimensional structures. Based on radial basis functions, eigenmodes are calculated in a collocation approach for the second-order wave equation. This type of meshless method promises highly accurate results with the simplicity of a node-based collocation approach. Thus, when changing the discrete representat