Lorenz-Mie theory for 2D scattering and resonance calculations

This PhD tutorial is concerned with a description of the two-dimensional generalized Lorenz-Mie theory (2D-GLMT), a well-established numerical method used to compute the interaction of light with arrays of cylindrical scatterers. This theory is based on the method of separation of variables and the application of an addition theorem for cylindrical functions. The purpose of this tutorial is to assemble the practical tools necessary to implement the 2D-GLMT method for the computation of scattering by passive scatterers or of resonances in optically active media. The first part contains a derivation of the vector and scalar Helmholtz equations for 2D geometries, starting from Maxwell's equations. Optically active media are included in 2D-GLMT using a recent stationary formulation of the Maxwell-Bloch equations called steady-state ab initio laser theory (SALT), which introduces new classes of solutions useful for resonance computations. Following these preliminaries, a detailed description of 2D-GLMT is presented. The emphasis is placed on the derivation of beam-shape coefficients for scattering computations, as well as the computation of resonant modes using a combination of 2D-GLMT and SALT. The final section contains several numerical examples illustrating the full potential of 2D-GLMT for scattering and resonance computations. These examples, drawn from the literature, include the design of integrated polarization filters and the computation of optical modes of photonic crystal cavities and random lasers.Comment: This is an author-created, un-copyedited version of an article published in Journal of Optics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from i

Numerical Analysis of Three-dimensional Acoustic Cloaks and Carpets

We start by a review of the chronology of mathematical results on the Dirichlet-to-Neumann map which paved the way towards the physics of transformational acoustics. We then rederive the expression for the (anisotropic) density and bulk modulus appearing in the pressure wave equation written in the transformed coordinates. A spherical acoustic cloak consisting of an alternation of homogeneous isotropic concentric layers is further proposed based on the effective medium theory. This cloak is characterised by a low reflection and good efficiency over a large bandwidth for both near and far fields, which approximates the ideal cloak with a inhomogeneous and anisotropic distribution of material parameters. The latter suffers from singular material parameters on its inner surface. This singularity depends upon the sharpness of corners, if the cloak has an irregular boundary, e.g. a polyhedron cloak becomes more and more singular when the number of vertices increases if it is star shaped. We thus analyse the acoustic response of a non-singular spherical cloak designed by blowing up a small ball instead of a point, as proposed in [Kohn, Shen, Vogelius, Weinstein, Inverse Problems 24, 015016, 2008]. The multilayered approximation of this cloak requires less extreme densities (especially for the lowest bound). Finally, we investigate another type of non-singular cloaks, known as invisibility carpets [Li and Pendry, Phys. Rev. Lett. 101, 203901, 2008], which mimic the reflection by a flat ground.Comment: Latex, 21 pages, 7 Figures, last version submitted to Wave Motion. OCIS Codes: (000.3860) Mathematical methods in physics; (260.2110) Electromagnetic theory; (160.3918) Metamaterials; (160.1190) Anisotropic optical materials; (350.7420) Waves; (230.1040) Acousto-optical devices; (160.1050) Acousto-optical materials; (290.5839) Scattering,invisibility; (230.3205) Invisibility cloak

Direct and Inverse Computational Methods for Electromagnetic Scattering in Biological Diagnostics

Scattering theory has had a major roll in twentieth century mathematical physics. Mathematical modeling and algorithms of direct,- and inverse electromagnetic scattering formulation due to biological tissues are investigated. The algorithms are used for a model based illustration technique within the microwave range. A number of methods is given to solve the inverse electromagnetic scattering problem in which the nonlinear and ill-posed nature of the problem are acknowledged.Comment: 61 pages, 5 figure

Embedding approach to modeling electromagnetic fields in a complex two-dimensional environment

An approach is presented to combine the response of a two-dimensionally inhomogeneous dielectric object in a homogeneous environment with that of an empty inhomogeneous environment. This allows an efficient computation of the scattering behavior of the dielectric cylinder with the aid of the CGFFT method and a dedicated extrapolation procedure. Since a circular observation contour is adopted, an angular spectral representation can be employed for the embedding. Implementation details are discussed for the case of a closed 434 MHz microwave scanner, and the accuracy and efficiency of all steps in the numerical procedure are investigated. Guidelines are proposed for choosing computational parameters such as truncation limits and tolerances. We show that the embedding approach does not increase the CPU time with respect to the forward problem solution in a homogeneous environment, if only the fields on the observation contour are computed, and that it leads to a relatively small increase when the fields on the mesh are computed as well

A finite element-boundary integral method for scattering and radiation by two- and three-dimensional structures

A review of a hybrid finite element-boundary integral formulation for scattering and radiation by two- and three-dimensional composite structures is presented. In contrast to other hybrid techniques involving the finite element method, the proposed one is in principle exact and can be implemented using a low O(N) storage. This is of particular importance for large scale applications and is a characteristic of the boundary chosen to terminate the finite element mesh, usually as close to the structure as possible. A certain class of these boundaries lead to convolutional boundary integrals which can be evaluated via the fast Fourier transform (FFT) without a need to generate a matrix; thus, retaining the O(N) storage requirement. The paper begins with a general description of the method. A number of two- and three-dimensional applications are then given, including numerical computations which demonstrate the method's accuracy, efficiency, and capability

Numerical methods for computing Casimir interactions

We review several different approaches for computing Casimir forces and related fluctuation-induced interactions between bodies of arbitrary shapes and materials. The relationships between this problem and well known computational techniques from classical electromagnetism are emphasized. We also review the basic principles of standard computational methods, categorizing them according to three criteria---choice of problem, basis, and solution technique---that can be used to classify proposals for the Casimir problem as well. In this way, mature classical methods can be exploited to model Casimir physics, with a few important modifications.Comment: 46 pages, 142 references, 5 figures. To appear in upcoming Lecture Notes in Physics book on Casimir Physic