🔬 applied physics

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

This paper presents a reduced-order model based on a quadratic eigenvalue problem to characterize surface Bloch waves in truncated periodic media, enabling the analysis of their dispersion, waveforms, and skin depth to facilitate the optimal design of surface cuts for manipulating these boundary layers.

Original authors: Bojan B. Guzina, Shixu Meng, Prasanna Salasiya, Long Nguyen

Published 2026-02-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Bojan B. Guzina, Shixu Meng, Prasanna Salasiya, Long Nguyen

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a giant, perfectly repeating wall made of tiny, identical bricks. In physics, this is called a periodic medium. When waves (like sound or vibrations) travel through this infinite wall, they behave in a very predictable, rhythmic way. Scientists have long known how to predict these "body waves" that travel deep inside the wall.

However, this paper asks a different question: What happens when you cut a slice off the side of this wall?

When you create a free surface (a cut edge), a special kind of wave can get "stuck" right at the edge, skimming along the surface and dying out as it goes deeper into the wall. The authors call these Surface Bloch (SB) waves. Think of them like a surfer riding the very top edge of a wave, rather than swimming in the deep water.

Here is the breakdown of their discovery, using simple analogies:

1. The "Rational Slope" Rule

The wall isn't just a smooth block; it's made of a specific grid of bricks (a lattice). If you cut the wall at a random, messy angle, the math gets incredibly messy and hard to solve.

The authors discovered a trick: If you cut the wall at a "rational slope" (meaning the cut line hits the grid points in a neat, repeating pattern, like going 3 bricks right and 1 brick up), the problem becomes solvable.

Analogy: Imagine walking on a tiled floor. If you walk straight along the tiles, it's easy. If you walk at a weird angle that never lines up with the tile edges, it's chaotic. The authors only study the "neat" angles where the cut lines up perfectly with the repeating pattern of the material.

2. The "Magic Recipe" (The Quadratic Eigenvalue Problem)

To figure out how these surface waves behave, the authors developed a mathematical "recipe" called a Quadratic Eigenvalue Problem (QEP).

The Ingredients: You give the recipe two things: the frequency (how fast the wave vibrates) and the speed along the surface.
The Output: The recipe spits out a list of "ghost waves." These are waves that could exist but are designed to fade away (evanesce) as they go deeper into the material.
The Analogy: Imagine you have a box of different-sized tuning forks. You shake the box (the QEP), and it tells you exactly which forks will vibrate in a way that matches the edge of your wall.

3. Building the Wave (The Superposition)

Once the recipe gives them the list of "ghost waves," the authors show how to mix them together to create the actual surface wave.

The Analogy: Think of the surface wave as a smoothie. The "ghost waves" are the individual fruits (strawberries, bananas, blueberries). You can't eat just one fruit and get the smoothie; you have to blend them in the right proportions. The authors figured out the exact "recipe" (the amplitudes) to blend these ghost waves so that they cancel out perfectly at the very edge of the cut, satisfying the rule that the edge must be free of stress (like a loose piece of paper flapping in the wind).

4. The "Skin Depth" (How Deep Does It Go?)

One of the most important findings is about skin depth. This is a measure of how deep the wave penetrates into the material before it becomes too weak to matter.

The Analogy: Imagine dipping a sponge in water. The "skin depth" is how wet the sponge gets before it feels dry in the middle.
The authors found that by changing the shape of the cut (making it wavy or undulating instead of straight), you can control this depth. A wavy cut can make the wave stay very close to the surface or dive deeper, depending on the design.

5. Why This Matters (The "What-If" Machine)

The biggest advantage of their method is speed and simplicity.

Old Way: To study these waves, scientists usually had to simulate the entire huge wall on a computer, which is slow and computationally heavy.
New Way: The authors' method reduces the problem to a tiny, low-dimensional puzzle.
The Analogy: Instead of simulating a whole ocean to see how a wave hits a specific rock, they built a small, portable "wave tank" that predicts the result instantly. This allows engineers to quickly test thousands of different surface shapes (undulations) to find the perfect one for a specific job, like trapping energy or blocking vibrations.

Summary

In short, this paper provides a blueprint for understanding and designing waves that travel along the edges of repeating materials. By focusing on "neat" cuts and using a clever mathematical shortcut (the QEP), they turned a complex, messy problem into a simple, fast calculation. This allows us to "tune" the edge of a material to control how waves move along it, much like tuning a guitar string to get the perfect note.

Technical Summary: On the Surface Bloch Waves in Truncated Periodic Media

Problem Statement
While the dispersive characteristics of waves in infinite periodic media (Floquet-Bloch body waves) are well-established, the analysis of boundary layers in truncated periodic media remains a significant challenge. Existing approaches often treat evanescent surface or interface waves as implicit computational byproducts or "waves to avoid," lacking explicit analytical frameworks to describe their dispersion, waveforms, and energetic relevance (skin depth). Furthermore, while topological mechanics offers tools for identifying protected interfacial states, these methods often fail to provide quantitative information regarding the penetration depth or power flow of such boundary layers. The paper addresses the need for a comprehensive framework to understand, simulate, and control surface waves propagating along the free surface of a periodic medium, specifically those characterized by homogeneous Neumann (traction-free) boundary conditions.

Methodology
The authors adopt a two-dimensional scalar wave equation with periodic coefficients as a test bed, modeling antiplane shear waves in phononic crystals. The core of the methodology involves a geometric interplay between the "mother" Bravais lattice of the periodic medium and the orientation of the surface cut.

Rational Surface Cuts: The analysis is restricted to surface cuts with a "rational slope" relative to the Bravais lattice. This condition allows the definition of a new, effective unit cell (a "multi-cell") and a companion Brillouin zone tailored to the surface wave direction.
Quadratic Eigenvalue Problem (QEP): The authors formulate a QEP over this effective unit cell. Unlike standard linear eigenvalue problems that relate frequency and wavenumber for propagating waves, this QEP seeks complex wavenumbers ( $\kappa$ ) for a given excitation frequency ( $\omega$ ) and a prescribed real wavenumber ( $k_\beta$ ) along the surface cut. The imaginary part of the complex eigenvalue $\kappa$ quantifies the exponential decay (evanescence) of the wave away from the cut plane.
Reduced-Order Model (ROM): The surface Bloch (SB) wave is constructed as a superposition of the evanescent eigenstates derived from the QEP. The relative amplitudes of these modes are determined by imposing the homogeneous Neumann boundary condition on the cut surface.
Surface-Wave Eigenproblem: To solve for the specific surface wavenumber and modal amplitudes, the traction-free condition is projected onto a finite set of QEP eigenmodes. This results in a low-dimensional matrix eigenvalue problem (or a singular value decomposition search) that identifies the surface-wave eigenvalues ( $k_e^*$ ) where the condition number of the system matrix peaks.
Handling Undulations: The framework extends to undulated (periodic) surface cuts by mapping the boundary condition onto the undulated surface using Fourier series expansions or, for non-bijective undulations, discrete prolate spheroidal wave (DPSW) functions to ensure numerical stability.

Key Contributions

Explicit QEP Formulation: The paper derives a specific QEP for the effective unit cell that yields complex wavenumbers, serving as the fundamental building blocks for evanescent boundary layers.
Low-Dimensional ROM: By leveraging the QEP eigenspectrum, the evaluation of SB waves (dispersion and waveforms) is reduced to a low-dimensional eigenvalue problem, avoiding the computational cost of full boundary value problems or supercell numerical evaluations.
Quantitative Energetic Metrics: The framework explicitly calculates the power flow density and "skin depth" of surface Bloch waves. The skin depth is defined energetically as the depth at which the contribution to the averaged power flow drops below a specific threshold, providing a physically meaningful measure of wave penetration.
Design Capability: The reduced-order nature of the model enables rapid exploration of the design space, specifically the effect of periodic surface undulations on the dispersion and penetration depth of surface waves.

Results
Numerical examples on a periodic medium with a square unit cell demonstrate the framework's capabilities:

Dispersion Characteristics: The method successfully generates dispersion diagrams for SB waves, identifying pass bands and band gaps. For specific cut depths (e.g., at unit cell interfaces), the SB dispersion coincides with that of body waves; however, for generic cut depths, the dispersion deviates significantly, exhibiting unique surface-wave branches.
Effect of Cut Depth: The skin depth and surface wavenumber are shown to be periodic functions of the cut depth. Changing the cut depth can significantly alter the band gap width and the shape of dispersion branches.
Effect of Undulation: Introducing sinusoidal undulations to the surface cut was shown to nearly double the width of the first band gap and significantly alter the skin depth, particularly in low-wavenumber regimes where the first QEP eigenvalue has a small imaginary part.
Rational Slopes: The framework successfully handles rational slopes (e.g., 1:1) by constructing appropriate multi-cells, demonstrating that the dispersion characteristics depend heavily on the geometric interplay between the lattice and the cut orientation.
Convergence: The reduced-order model was shown to converge rapidly, with the first few evanescent QEP eigenmodes dominating the solution.

Significance
The paper claims to provide the first comprehensive analytical and computational framework for elucidating boundary layers in periodic media. By moving beyond the "waves to avoid" paradigm, the authors establish a method to explicitly compute and manipulate surface Bloch waves. The significance lies in the ability to:

Predict and Control: Explicitly determine the dispersion and skin depth of surface waves, which are critical for the energetic relevance of boundary layers.
Optimize Design: Facilitate the optimal design of surface cuts (including undulations) to manipulate wave propagation, such as creating surface-wave band gaps or controlling penetration depth.
Generalize: While demonstrated on scalar waves, the framework is presented as readily extendable to other governing equations (e.g., Navier, Maxwell) and boundary conditions, offering a unified approach to wave motion in bounded periodic structures.

The authors emphasize that their approach bridges the gap between the theoretical analysis of infinite periodic media and the practical engineering of bounded metamaterials, providing a tool to manipulate boundary layers via the geometric design of the surface cut.