On the surface Bloch waves in truncated periodic media: scalar-wave primer

Much like their counterparts in homogenous elastic solids, waves in periodic media can be broadly classified into Floquet-Bloch body waves, and evanescent surface waves. Our goal is to elucidate the latter boundary layers, termed surface Bloch (SB) waves, affiliated with rational surface cuts and homogeneous Neumann data. To this end we adopt a two-dimensional (2D) scalar wave equation with periodic coefficients (describing anticline shear waves in phonoic crystals) as a test bed and develop a unit cell-of-periodicity-based, reduced order model of the SB waves that is capable of describing both their dispersion, waveforms, and skin depth''. The centerpiece of our analysis is a quadratic eigenvalue problem (QEP) for the effective unit cell of periodicity -- deriving from a geometric interplay between the mother Bravais lattice and orientation of the surface cut -- that seeks the complex wavenumber normal to the cut plane given (i) the excitation frequency and (ii) wavenumber in the direction of the cut plane. In this way the sought boundary layer is derived via superposition of the evanescent QEP eigenstates, whose relative amplitudes are obtained by imposing the homogeneous boundary condition. With the QEP eigenspectrum at hand, evaluation of an SB wave -- in terms of both dispersion characteristics and evanescent waveforms -- entails only a low-dimensional eigenvalue problem. This feature caters for rapid exploration of the effect of (periodic) surface undulations, and so enables manipulation of the SB waves via optimal design of the surface cut. Our analysis also includes an account for the power flow and skin depth’’ of a surface Bloch wave, both of which are critical for the energetic relevance of boundary layers.

💡 Research Summary

The paper presents a comprehensive theoretical and computational framework for surface Bloch (SB) waves that arise when a periodic medium is truncated, focusing on a two‑dimensional scalar wave model with periodic coefficients. The authors begin by motivating the study: while most metamaterial research concentrates on bulk Bloch dispersion in infinite periodic structures, practical applications inevitably involve boundaries, and the associated evanescent surface layers have received comparatively little analytical attention. Recent advances in handling complex wave vectors through quadratic eigenvalue problems (QEP) in band‑gap analysis suggest a promising route for surface‑layer modeling.

The governing equation is a scalar anti‑plane shear wave equation, ∇·(G(x)∇u)+ω²ρ(x)u=0, where G(x) and ρ(x) are Y‑periodic shear modulus and density, respectively. By applying the Floquet‑Bloch theorem, the field is written as u(x)=ϕ(x) e^{i k·x} with a Y‑periodic envelope ϕ. The wave vector k is split into two components: k_β, parallel to the cut plane, and k_α, normal to it. The latter controls exponential decay away from the surface.

The core of the analysis is the derivation of a QEP for the normal component κ≡k_α. Starting from a weak variational form, the authors define three self‑adjoint operators A, B, and C acting on the periodic function space H¹_p(Y). The QEP reads A ϕ + κ B ϕ + κ² C ϕ = 0. Lemmas establish that B is compact and that the combined operator T(λ)=A+λB+λ²C is Fredholm of index zero, guaranteeing a discrete spectrum consisting of complex‑conjugate pairs. This mathematical foundation ensures that only a finite number of evanescent modes need to be considered for any given frequency and parallel wave number.

To obtain physically admissible surface waves, the authors impose a “rational slope” condition on the cut plane S and the propagation direction e: the direction must be a rational linear combination of the Bravais lattice vectors (e ∥ q₁a₁+q₂a₂ with integer q₁,q₂). This condition preserves Bloch periodicity along the surface and allows the surface wave to be expressed as a superposition of the QEP eigenstates with Im(κ_i)>0:

u_S(x)=∑{i=1}^N α_i ϕ_i(x) e^{i(k_β·x∥+κ_i x_⊥)}.

The homogeneous Neumann boundary condition (traction‑free surface) supplies a low‑dimensional linear system for the amplitudes α_i. Consequently, the full surface‑wave field is obtained by solving only a small eigenvalue problem on the unit cell, rather than a large global boundary‑value problem.

Beyond dispersion, the paper defines and computes the power flow associated with each evanescent mode and introduces a “skin depth” metric that quantifies how far the surface wave penetrates into the bulk. By varying the surface geometry—e.g., adding periodic undulations or changing the cut angle—the QEP spectrum changes, enabling systematic control of phase velocity, attenuation rate, and energy transmission. The authors demonstrate that the reduced‑order QEP approach dramatically speeds up parametric sweeps and design optimization compared with conventional finite‑element simulations.

In summary, the work establishes a rigorous, low‑dimensional QEP‑based methodology for predicting and engineering surface Bloch waves in truncated periodic media. It bridges the gap between bulk Bloch analysis and boundary‑layer physics, offering a practical tool for metamaterial designers who wish to tailor surface‑wave characteristics through rational cut orientations and engineered surface patterns.