Ray Lattice: Acoustic & Elastic Wave Control

• Ray lattices are periodic structures comprising Rayleigh beams that integrate both transverse and rotational inertia to uniquely control wave dispersion.
• They exhibit distinctive phenomena such as Dirac cones, flat bands, and directional anisotropy, which facilitate applications in waveguides, filters, and negative refraction devices.
• The framework uses Floquet–Bloch methods and tailored beam geometries to design elastic and acoustic metamaterials for controlled wave localization and transmission.

A Ray lattice is a periodic structural or acoustic network whose wave physics and dispersion properties are governed or substantially influenced by Rayleigh modes—vibrational modes exhibiting nontrivial effects of rotational inertia. The concept encompasses elastic flexural beam systems represented by Rayleigh beam models and acoustic arrays supporting generalised Rayleigh–Bloch waves. Ray lattices enable phenomena such as Dirac cones, flat bands, directional anisotropy, negative refraction, wave localization, and radiative lattice resonances, with broad implications for wave manipulation in metamaterial design.

1. Elastic Ray Lattice Structures: Geometric and Physical Foundations

Ray lattices in elasticity are formalized as planar networks of Rayleigh beams. Two canonical geometries appear in the literature:

• Honeycomb lattice: The arrangement consists of beam ligaments of length $h$, joined in a hexagonal tiling with primitive vectors $a_1 = (\sqrt{3}/2\, h,\ 3/2\, h)$ and $a_2 = (-\sqrt{3}/2\, h,\ 3/2\, h)$. The elementary cell is a parallelogram spanned by these vectors, containing two inequivalent junctions (E and F), each connecting three beams.
• Square lattice: Beams meet at right angles with vectors $a_1 = (h, 0)$, $a_2 = (0, h)$, and each node symmetrically connects four beams.

The structural unit (Rayleigh beam) incorporates both transverse inertia ($\rho A$) and additional rotational inertia ($\rho I$) per unit length. The governing equation for time-harmonic flexural waves $w(x, t) = w(x) e^{i \omega t}$ in a beam (with $P=0$, $\beta=0$) is

$$EI\, w''''(x) - \rho I\, \omega^2\, w''(x) - \rho A\, \omega^2\, w(x) = 0$$

This equation differentiates Rayleigh beams from classical Euler–Bernoulli beams (where the $\rho I$ term is absent) and underpins all Ray lattice wave phenomena (Cabras et al., 2017, Piccolroaz et al., 2017).

2. Floquet–Bloch Wave Methods and Dispersion Relations

Periodic Ray lattices enable Floquet–Bloch analysis for the determination of wave spectra. In honeycomb and square geometries, the general solution for displacement along each beam is expanded over the four characteristic roots $\kappa_{1,2,3,4}$ of the Rayleigh beam equation:

$$w_q(x_q) = \sum_{p=1}^4 C_{pq}\, e^{i\kappa_{pq} x_q}$$

Applying Bloch boundary conditions and continuity/equilibrium at lattice junctions leads to a homogeneous linear system

$A(\omega, k) \cdot C = 0$

The dispersion surfaces $\omega = \omega(k)$ are implicitly defined by $\det A(\omega, k) = 0$.

In honeycomb Ray lattices with normalized parameters, an explicit factorized dispersion relation is constructed leveraging $D_6$ symmetry (Cabras et al., 2017). The equation incorporates branches corresponding to both propagating and standing wave (flat-band) solutions and encodes the interplay of geometry and inertia.

3. Dirac Cones, Flat Bands, and Directional Anisotropy

Ray lattices support distinctive features in their band structures:

• Dirac cones: Conical intersections of dispersion sheets occur at K-points of the Brillouin zone, giving rise to massless-like wave behavior. In the Rayleigh-beam honeycomb lattice, rotational inertia lowers the frequency $\omega_D$ of the Dirac points compared to the Euler–Bernoulli case.
• Flat bands: The explicit factorization of the dispersion equation yields conditions ($\Omega_1 h = n\pi$) for flat bands, representing standing waves localized on individual beams. The Rayleigh and Euler–Bernoulli frequencies for such modes differ by inclusion of $I$ in the denominator.
• Directional anisotropy: Saddle points in the dispersion surfaces produce non-convex isofrequency contours, manifesting as directional wave propagation and preferential energy flux along beam axes, especially pronounced in square geometries.

These phenomena enable advanced control over elastic wave propagation, including the realization of edge modes and localized states via geometric and inertial design (Cabras et al., 2017, Piccolroaz et al., 2017).

4. Structured Interfaces, Transmission Phenomena, and Wave Localization

Embedding a slab of Rayleigh beams within an Euler–Bernoulli lattice leads to frequency-dependent transmission, reflection, and waveguiding behavior:

• Low-frequency transparency: For $\omega \rightarrow 0$, Rayleigh and Euler–Bernoulli lattices are dynamically indistinguishable, resulting in perfect transmission across interfaces.
• Negative refraction and focusing: Near saddle-point frequencies, the Rayleigh slab exhibits hyperbolic slowness contours, inducing negative refraction, flat-lens focusing, and wave mirroring.
• Waveguide and edge modes: At frequencies where the Rayleigh slab admits propagating or resonant modes (but the host does not), flexural energy localizes within the slab or along its interface, giving rise to trapped waveguides and edge-localized waves, especially at Dirac-point frequencies.
• Resonant transmission and band gaps: High-frequency regions can yield resonance-enhanced transmission, wide band gaps, and near-zero group velocities due to the downward shift and flattening of Rayleigh lattice bands.

These effects underpin the design of elastic metamaterials, waveguides, filters, and phononic devices capable of tailored wave control (Cabras et al., 2017, Piccolroaz et al., 2017).

5. Acoustic Ray Lattices and Generalised Rayleigh–Bloch Waves

In airborne acoustic systems, the Ray lattice paradigm centers on 1D periodic arrays of Neumann (sound-hard) cylinders in air described by the Helmholtz equation

$(\Delta + k^2)\, \phi = 0$

with Floquet–Bloch periodicity and outgoing-radiation conditions. Rayleigh–Bloch (RB) waves are eigenmodes exponentially localized to the array, existing both below and—crucially, for sufficiently small cylinder spacing—above the first cut-off frequency (the radiative regime). Generalised RB waves above cut-off exhibit complex Bloch phase $\beta$ and radiate into the ambient medium, a phenomenon confirmed both theoretically (via T-matrix multiple-scattering and lattice-sum techniques) and experimentally (Chaplain et al., 2024).

• Dispersion relations: RB branches in $k$–$\beta$ space are shaped by array aspect ratio $r/a$. For $r/a<0.33$, a radiative branch peels off above cut-off; for $r/a=0.35,$ it vanishes, as the inter-cylinder dipole moment is suppressed.
• Finite array resonances: Short arrays exhibit sharp lattice resonances below and (for $r/a=0.15$) above cut-off; long arrays manifest propagating generalised RB waves corresponding to these resonances.
• Nomenclature unification: Infinite Ray lattices correspond to RB eigenmodes (Floquet–Bloch periodic, exponentially localized); special modes at zone edges are termed Neumann and Dirichlet; finite arrays admit broadened, lossy resonances due to end reflection.

This acoustic Ray lattice theory generalizes Rayleigh–Bloch physics to non-resonant scatterer arrays and connects the eigenmode picture to observable macroscopic wave control (Chaplain et al., 2024).

6. Device and Metamaterial Implications, Design Principles

The wave-manipulation capabilities of Ray lattices—elastic and acoustic—suggest device concepts encompassing:

• Flexible waveguides: The existence of edge and waveguide modes in structured interfaces enables low-loss, tunable guiding of waves.
• Transparent coatings or filters: Low-frequency perfect transmission and high-frequency band gaps can be harnessed in phononic filtering and “invisible” coatings.
• Flat acoustic and elastic lenses: Negative refraction and hyperbolic dispersion allow subwavelength focusing via simple lattice engineering.
• Metasurface antennas and delay lines: Radiative generalised RB waves and lattice resonances form the basis for highly directive radiators and efficient signal delay.

Ray lattice design proceeds via controlled selection of geometric parameters (e.g., lattice spacing, node connectivity), mechanical properties (mass, rotational inertia), and array aspect ratio, in line with the criteria laid out for existence of radiative RB branches (e.g., $r/a<0.33$ in acoustics). All functionalities derive from wave phenomena inherent in the Floquet–Bloch structure and the nontrivial role of rotational inertia.

A plausible implication is that Ray lattice concepts unify a broad class of wave phenomena, guiding rational design of next-generation metamaterials across acoustics and elasticity (Cabras et al., 2017, Piccolroaz et al., 2017, Chaplain et al., 2024).