Plug-and-play analytical paradigm for the scattering of plane waves by 'layer-cake' periodic systems

We investigate the scattering of scalar plane waves in two dimensions by a heterogeneous layer that is periodic in the direction parallel to its boundary. On describing the layer as a union of periodic laminae, we develop a solution of the scattering problem by merging the concept of propagator matrices and that of Bloch eigenstates featured by the unit cell of each lamina. The featured Bloch eigenstates are obtained by solving the quadratic eigenvalue problem (QEVP) that seeks a complex-valued wavenumber normal to the layer boundary given (i) the excitation frequency, and (ii) real-valued wavenumber parallel to the boundary - that is preserved throughout the system. Spectral analysis of the QEVP reveals sufficient conditions for discreteness of the eigenvalue spectrum and the fact that all eigenvalues come in complex-conjugate pairs. By deploying the factorization afforded by the propagator matrix approach, we demonstrate that the contribution of individual eigenvalues (and so eigenmodes) to the solution diminishes exponentially with absolute value of their imaginary part, which then forms a rational basis for truncation of the factorized Bloch (FB)-wave solution. The proposed methodology caters for the optimal design of rainbow traps, energy harvesters and metasurfaces, whose potency to manipulate waves is decided not only by the individual dispersion characteristics of the component laminae, but also by ordering and generally fitting of the latter into a composite layer. Using the FB approach, evaluation of trial configurations - as generated by the permutation and window translation/stretching of the component laminae - can be accelerated by decades.