In this paper, we show that the diffraction of the primes is absolutely continuous, showing no bright spots (Bragg peaks). We introduce the notion of counting diffraction, extending the classical notion of (density) diffraction to sets of density zero. We develop the counting diffraction theory and give many examples of sets of zero density of all possible spectral types.

Consider the extended hull of a weak model set together with its natural shift action. Equip the extended hull with the Mirsky measure, which is a certain natural pattern frequency measure. It is known that the extended hull is a measure-theoretic factor of some group rotation, which is called the underlying torus. Among other results, in the article Periods and factors of weak model sets, we showed that the extended hull is isomorphic to a factor group of the torus, where certain periods of the window of the weak model set have been factored out. This was proved for weak model sets having a compact window. In this note, we argue that the same results hold for arbitrary measurable and relatively compact windows. Our arguments crucially rely on Moody’s work on uniform distribution in model sets. We also discuss implications for the diffraction of such weak model sets and discuss a new class of examples which are generic for the Mirsky measure.

This paper gives a brief overview of some configurations in which high-frequency wave propagation modelled by Helmholtz equation gives rise to solutions that vary rapidly across thin layers. The configurations are grouped according to their mathematical structure and tractability and one of them concerns a famous open problem of mathematical physics.

We investigate the Fano resonance in grating structures using coupled resonators. The grating consists of a perfectly conducting slab with periodically arranged subwavelength slit holes, where inside each period, a pair of slits sit very close to each other. The slit holes act as resonators and are strongly coupled. It is shown rigorously that there exist two groups of resonances corresponding to poles of the scattering problem. One sequence of resonances has imaginary part in the order of ε, where ε is the size of the slit aperture, while the other sequence has imaginary part in the order of ε2. When coupled with the incident wave at resonant frequencies, the narrow-band resonant scattering induced by the latter will interfere with the broader background resonant radiation induced by the former. The interference of these two resonances generates the Fano-type transmission anomaly, which persists in the whole radiation continuum of the grating structure as long as the slit aperture size is small compared to the incident wavelength.

Applications of a WKBJ-type ‘ray ansatz’ to obtain asymptotic solutions of the Helmholtz equation in the high-frequency limit are now standard and underpin the construction of ‘geometrical optics’ ray diagrams in many electromagnetic, acoustic and elastic reflection, transmission and other scattering problems. These applications were subsequently extended by Keller to include other types of rays – called ‘diffracted’ rays – to provide an accessible and impressively accurate theory which is relevant in wide-ranging sets of circumstances. Friedlander and Keller then introduced a modified ray ansatz to extend yet further the scope of ray theory and its applicability to certain other classes of diffraction problems (tangential ray incidence upon an obstructing boundary, for instance) and did so by the inclusion of an extra term proportional to a power of the wave number within the exponent of the initial ansatz. Our purpose here is to generalise this further still by the inclusion of several such terms, ordered in a natural sequence in terms of strategically chosen fractional powers of the large wave number, and to derive a systematic sequence of boundary value problems for the coefficient phase functions that arise within this generalised exponent, as well as one for the leading-order amplitude occurring as a pre-exponential factor. One particular choice of fractional power is considered in detail, and waves with specified radially symmetric or planar wavefronts are then analysed, along with a boundary value problem typifying two-dimensional radiation whereby arbitrary phase and amplitude variations are specified on a prescribed boundary curve. This theory is then applied to the scattering of plane and cylindrical waves at curved boundaries with small-scale perturbations to their underlying profile.

Periodic structures involving crossed arrays of cylinders appear as special three-dimensional photonic crystals and cross-stacked gratings. Such a structure consists of a number of layers where each layer is periodic in one spatial direction and invariant in another direction. They are relatively simple to fabricate and have found valuable applications. For analyzing scattering properties of such structures, general computational electromagnetics methods can certainly be used, but special methods that take advantage of the geometric features are often much more efficient. In this paper, an efficient method based on operators mapping electromagnetic field components between two spatial directions is developed to analyze structures with crossed arrays of circular cylinders. The method is much simpler than an earlier method based on similar ideas, and it does not require evaluating slowly converging lattice sums.

An algebraic multilevel method is proposed for efficiently simulating linear wave propagation using higher-order numerical schemes. This method is used in conjunction with the Finite Volume Time Domain (FVTD) technique for the numerical solution of the time-domain Maxwell’s equations in electromagnetic scattering problems. In the multilevel method the solution is cycled through spatial operators of varying orders of accuracy, while maintaining highest-order accuracy at coarser approximation levels through the use of the relative truncation error as a forcing function. Higher-order spatial accuracy can be enforced using the multilevel method at a fraction of the computational cost incurred in a conventional higher-order implementation. The multilevel method is targeted towards electromagnetic scattering problems at large electrical sizes which usually require long simulation times due to the use of very fine meshes dictated by point-per-wavelength requirements to accurately model wave propagation over long distances.

A spectral-element method is developed to solve the scattering problem for time-harmonic sound waves due to an obstacle in an homogeneous compressible fluid. The method is based on a boundary perturbation technique coupled with an efficient spectral-element solver. Extensive numerical results are presented, in order to show the accuracy and stability of the method.

Consider the acoustic wave scattering by an impenetrable obstacle in two dimensions, where the wave propagation is governed by the Helmholtz equation. The scattering problem is modeled as a boundary value problem over a bounded domain. Based on the Dirichlet-to-Neumann (DtN) operator, a transparent boundary condition is introduced on an artificial circular boundary enclosing the obstacle. An adaptive finite element based on a posterior error estimate is presented to solve the boundary value problem with a nonlocal DtN boundary condition. Numerical experiments are included to compare with the perfectly matched layer (PML) method to illustrate the competitive behavior of the proposed adaptive method.

In this paper, hp-adaptive finite element methods are studied for time-harmonic Maxwell’s equations. We propose the parallel hp-adaptive algorithms on conforming unstructured tetrahedral meshes based on residual-based a posteriori error estimates. Extensive numerical experiments are reported to investigate the efficiency of the hp-adaptive methods for point singularities, edge singularities, and an engineering benchmark problem of Maxwell’s equations. The hp-adaptive methods show much better performance than the h-adaptive method.

We extend classic Sommerfeld and Silver-Müller radiation conditions for bounded scatterers to acoustic and electromagnetic fields propagating over three isotropic homogeneous layers in three dimensions. If X= (x1,x2,x3)ϵℝ3, with x3 denoting the direction orthogonal to the layers, standard conditions only hold for the outer layers in the region ∣x3∣ > ∣∣x∣γ, for γϵ(1/4,1/2) and x large. For ∣x3∣ < ∣∣x∣∣γ and inside the slab, asymptotic behavior depends on the presence of surface or guided modes given by the discrete spectrum of the associated operator.

We show how to apply convolution quadrature (CQ) to approximate the time domain electric field integral equation (EFIE) for electromagnetic scattering. By a suitable choice of CQ, we prove that the method is unconditionally stable and has the optimal order of convergence. Surprisingly, the resulting semi discrete EFIE is dispersive and dissipative, and we analyze this phenomena. Finally, we present numerical results supporting and extending our convergence analysis.

In a two-dimensional (2D) photonic crystal (PhC) composed of circular cylinders (dielectric rods or air holes) on a square or triangular lattice, various PhC devices can be created by removing or modifying some cylinders. Most existing numerical methods for PhC devices give rise to large sparse or smaller but dense linear systems, all of which are expensive to solve if the device is large. In a previous work [Z. Hu et al., Optics Express, 16 (2008), 17383-17399], an efficient Dirichlet-to-Neumann (DtN) map method was developed for general 2D PhC devices with an infinite background PhC to take full advantage of the underlying lattice structure. The DtN map of a unit cell is an operator that maps the wave field to its normal derivative on the cell boundary and it allows one to avoid computing the wave field in the interior of the unit cell. In this paper, we extend the DtN map method to PhC devices with a finite background PhC. Since there is no bandgap effect to confine the light in a finite PhC, a different technique for truncating the domain is needed. We enclose the finite structure with a layer of empty boundary and corner unit cells, and approximate the DtN maps of these cells based on expanding the scattered wave in outgoing plane waves. Our method gives rise to a relatively small and sparse linear systems that are particularly easy to solve.