## Biharmonic Split Ring Resonator Metamaterial: Artificially dispersive effective density in thin periodically perforated plates

Mohamed Farhat, 1, 2 Stefan Enoch, 3 and Sebastien Guenneau 3

1 Institute of Condensed Matter Theory and Solid State Optics, Abbe Center of Photonics, Friedrich-Schiller-Universit¨t Jena, D-07743 Jena, Germany a 2 Division of Computer, Electrical, and Mathematical Sciences and Engineering, 4700 King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia 3 Aix-Marseille Universit´, CNRS, Centrale Marseille, Institut Fresnel, e Campus universitaire de Saint-J´rˆme, 13013 Marseille, France e o

(Dated: August 4, 2014)

## Abstract

We present in this paper a theoretical and numerical analysis of bending waves localized on the boundary of a platonic crystal whose building blocks are split ring resonators (SRR). We first derive the homogenized parameters of the structured plate using a three-scale asymptotic expansion in the linearized biharmonic equation. In the limit when the wavelength of the bending wave is much larger than the typical heterogeneity size of the platonic crystal, we show that it behaves as an artificial plate with an anisotropic effective Young modulus and a dispersive effective mass density. We then analyze dispersion diagrams associated with bending waves propagating within an infinite array of SRR, for which eigen-solutions are sought in the form of Floquet-Bloch waves. We finally demonstrate that this structure displays the hallmarks of All-Angle-Negative-Refraction(AANR) and it leads to superlensing and ultrarefraction effects, interpreted thanks to our homogenization model as a consequence of negative and vanishing effective density, respectively.

PACS numbers: 43.40.+s, 46.40.Cd, 62.30.+d, 81.05.Xj

## I. INTRODUCTION

Left Handed Materials (LHM) are a new kind of materials which were theoretically envisionned by Veselago 1 as early as in 1967. Such materials have simultaneously negative relative permittivity ( ε r ) and negative relative permeability ( µ r ). This theoretical curiosity became a real field of research in 2000 after Pendry showed the potential of LHM to overcome the diffraction limit 2 and Smith et al. 3 proposed a first realization for such extra-ordinary materials based on periodic lattices combining Split Ring Resonators (SRR: Concentric annular rings with splits) and wires. The latter work can be considered as the experimental foundation of LHM (as the first experimental evidence of negative refraction) and it is based on a theoretical study by Pendry et al. which shows that negative permittivity could be obtained by a periodic arrangement of parallel wires and that a periodic lattice of SRR had a negative magnetic response around its resonance frequency . 4

In the recent years, there has been a keen interest in wave propagation in periodically structured media . 5 Investigation of photonic crystals has paved the way to the theoretical prediction and experimental realization of photonic band gaps 6-12 i.e. ranges of frequencies for which light, or a light polarization, is disallowed to propagate. Soon after, the focus has been extended to the study of acoustic waves in periodic media, and the existence of phononic band gaps has been verified both theoretically and experimentally 13-18 . Recently, the interest was even extended to the study of different types of waves, e.g. liquid surface waves 19-23 or biharmonic waves 24-28 in perforated thin plates. It has been shown that complete bandgaps also exist for these waves when propagating through a periodic lattice of vertically standing rods or over a periodically perforated thin plate 22,26 . In addition, many interesting phenomena have been reported, including negative refraction 1,29-33 , the superlensing effect and cloaking 34-38 . The essential condition for the AANR effect is that the constant frequency surfaces (EFS: equifrequency surfaces) should become convex everywhere about some point in the reciprocal space, and the size of these EFS should shrink with increasing frequency 9-11 .

In this paper, we focus on the application of split ring resonator (SRR) structures 4,39,40 to the domain of elastic waves. We first derive the homogenized governing equations of bending waves propagating within a thin-plate with a doubly periodic square array of freely

vibrating holes shaped as SRR, from the generalized biharmonic equation, and an asymptotic analysis involving three scales (one for the thickness of the thin-cut of each SRR, one for the array-pitch, and one for the wave wavelength). We then present an analysis of dispersion curves. To do this, we set the spectral problem for the biharmonic operator within a doubly periodic square array of SRR, homogeneous stress-free boundary conditions are prescribed on the contour of each resonator and the standard Floquet-Bloch conditions are set on the boundary of an elementary cell of the periodic structure. Such a structure presents an elastic bandgap at low frequencies. It turns out that the asymptotic analysis of our structure allows us to get analytically the frequency of the first localized mode and then the frequency of the first band gap. The aim of our work is actually to demonstrate the AANR effect at low frequencies for elastic thin perforated plates as well as their superlensing properties. Ultrarefraction is also considered and shows the versatility and power of using such structured media to realize new functionalities for surface elastic waves.

## II. HOMOGENIZATION OF A THIN-PLATE WITH AN ARRAY OF STRESSFREE SRR INCLUSIONS NEAR RESONANCE

The equations for bending of plates can be found in many textbooks 41,42 . The wavelength λ is supposed to be large enough compared to the thickness of the plate H and small compared to its in-plane dimension L , i.e. H /lessmuch λ /lessmuch L . In this case we can adopt the hypothesis of the theory of Von-Karman 41,42 . In this way, the mathematical setup is essentially two-dimensional, the thickness H of the plate appearing simply as a parameter in the governing equation.

We would like to homogenize a periodically structured thin-plate involving resonant elements. The resonances are associated with fast-oscillating displacement fields in thin-bridges of perforations shaped as split ring resonators (SRR), and we filter these oscillations by introducing a third scale in the usual two-scale expansion. We start with the Kirchhoff-Love equation and we consider an open bounded region Ω f ∈ R 2 . This region is e.g. a slab lens consisting of a square array of SRR shaped as the letter C .

When the bending wave penetrates the structured area Ω f of the plate whose geometry is shown in Figs. 1(c)-(d), it undergoes fast periodic oscillations. To filter these oscillations, we consider an asymptotic expansion of the associated vertical displacement U η solution of the

FIG. 1: (a) Geometry of a split ring resonator C consisting of a disc Σ connected to a thin ligament Π η of length l and thickness ηh where 0 &lt; η /lessmuch 1 in a unit cell Y; (b) Helmholtz resonator consisting of a mass connected to a wall via a spring which models resonances of SRR in (a); (c) Doubly periodic square lattice of SRR with the first Brillouin zone ΓXM in reciprocal space; (d) Geometry of the thin plate of thickness H , with a source on the left side of a platonic crystal (PC) slab occupying the region Ω f .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-150aa54d-949e-4bd5-b188-d6ccf2d162a4.jpg)

biharmonic equation given in (1) in terms of a macroscopic (or slow) variable x = ( x , x 1 2 ) and a microscopic (or fast) variable x η = x /η , where η is a small positive real parameter. With all the above assumptions, the out-of-plane displacement u η = (0 0 , , U η ( x , x 1 2 )) in the x 3 -direction (along the vertical axis) is solution of (assuming a time-harmonic dependence

exp( -iωt ) with ω the angular wave frequency):

$$\begin{array} { c } \nu _ { \eta } \, \max \, \nu \, \max \, \nu \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \, \nu _ { \eta } \, \max \\ + \frac { \partial ^ { 2 } } { \partial x _ { 2 } \partial x _ { 2 } } \left ( D _ { \eta } \left ( \frac { \partial ^ { 2 } } { \partial x _ { 2 } \partial x _ { 2 } } + \nu _ { \eta } \frac { \partial ^ { 2 } } { \partial x _ { 1 } \partial x _ { 1 } } \right ) \right ) U _ { \eta } \\ + 2 \frac { \partial ^ { 2 } } { \partial x _ { 1 } \partial x _ { 2 } } \left ( D _ { \eta } \left ( 1 - \nu _ { \eta } \right ) \frac { \partial ^ { 2 } } { \partial x _ { 1 } \partial x _ { 2 } } \right ) U _ { \eta } - \, \beta _ { \eta } ^ { 4 } U _ { \eta } = 0 \,, \\ \text{he heterogeneous isotropic region } \Omega _ { f } \left ( \text{the platonic crystal, PC, in Fig. } 1 ( d ) \right ), \text{where} \end{array}$$

inside the heterogeneous isotropic region Ω f (the platonic crystal, PC, in Fig. 1(d)), where

$$D _ { \eta } = D ( \frac { \mathbf x } { \eta } ) \, \ \nu _ { \eta } = \nu ( \frac { \mathbf x } { \eta } ) \ \text{ and } \ \beta _ { \eta } ^ { 4 } = \beta _ { 0 } ^ { 4 } \rho ( \frac { \mathbf x } { \eta } )$$

are nondimensionalized spatially varying parameters related to the flexural rigidity of the plate, its Poisson ratio and the wave frequency, respectively. In most cases, D and ν take piecewise constant values, with D &gt; 0 and -1 / 2 &lt; ν &lt; 1 / 2. Note that β 2 0 = ω √ ρ H/D 0 0 , where D 0 is the flexural rigidity of the plate outside the platonic crystal, ρ 0 its density and H its thickness.

Remark that (1) is written in weak form and we notably retrieve the classical boundary conditions for a homogeneous plate with stress-free inclusions (vanishing of bending moments and shearing stress for vanishing D η and ν η in the soft phase) 42 . Since there is only one phase in the problem which we consider (homogeneous medium outside freely-vibrating inclusions), it is also possible to recast (1) as

$$( \sqrt { D _ { \eta } } \nabla ^ { 2 } U _ { \eta } + \beta _ { \eta } ^ { 2 } ) ( \sqrt { D _ { \eta } } \nabla ^ { 2 } U _ { \eta } - \beta _ { \eta } ^ { 2 } ) U _ { \eta } = 0 \, \ \text{in} \ \Omega _ { f } \ \, \overline { \Theta } _ { \eta } \, \\ \intertext { e D _ { \ } v a n i s h e s i n s i d e t h e i n c i s i n s o n s o n } \Theta _ { \ } = 1 1.. \in \{ n ( i + C \), \, \text{and} \, i t i s a \, \text{constant in the} }$$

$$& \text{at} \left ( ( 1 + \nu _ { \eta } ) \left ( \frac { \partial ^ { 2 } } { \partial x _ { 1 } ^ { 2 } } + \frac { \partial ^ { 2 } } { \partial x _ { 2 } ^ { 2 } } \right ) + 2 ( 1 - \nu _ { \eta } ) \left ( \frac { \partial ^ { 2 } } { \partial x _ { 1 } \partial x _ { 2 } } \right ) \right ) U _ { \eta } = 0 \, \\ & \left ( ( 3 - \nu _ { \eta } ) \left ( \frac { \partial ^ { 3 } } { \partial x _ { 1 } ^ { 3 } } + \frac { \partial ^ { 3 } } { \partial x _ { 2 } ^ { 3 } } \right ) + ( 1 + \nu _ { \eta } ) \left ( \frac { \partial ^ { 3 } } { \partial x _ { 1 } ^ { 2 } \partial x _ { 2 } } + \frac { \partial ^ { 3 } } { \partial x _ { 1 } \partial x _ { 2 } ^ { 2 } } \right ) \right ) U _ { \eta } = 0 \, \\ ; & \text{the boundary } \partial \Theta _ { \eta } \text{ of } \Theta _ { \eta }, \text{ which is consistent with our former work on thin perforated}$$

since D η vanishes inside the inclusions Θ η = ⋃ i ∈ Z 2 { η i ( + C ) } and it is a constant in the matrix. Bear in mind that the number of SRR in Ω f is an integer which scales as η -2 . Note also that the vanishing of bending moment and shearing stress deduced from (1) requires that

At the boundary ∂ Θ η of Θ , which is consistent with our former work on thin perforated η plates 27 .

In the present case, perforations are shaped as split ring resonators, and each SRR C can be modeled as

$$C = \{ a < \sqrt { x _ { 1 } ^ { 2 } + x _ { 2 } ^ { 2 } } < b \} \, \langle \, \overline { \Pi _ { \eta } } \, \\ \, \rangle$$

where a and b are functions of variables x , x 1 2 , unless the ring is circular and

$$\Pi _ { \eta } = \left \{ ( x _ { 1 }, x _ { 2 } ) \, \colon 0 < x _ { 1 } < l \,, \, | \, x _ { 2 } \, | < \eta h / 2 \right \}, \, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \Pi _ { \eta } = \left \} & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & + & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & = & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & 0 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ^ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & b & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & | & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

is a thin ligament of length l = b -a between the ends of the letter C , see Fig. 1(a).

Our aim is to show that the homogenized multi-structured platonic structure within Ω f is characterized by an effective density which can take negative values near the fundamental resonance of the SRR. To do this, we need to perform a homogenization of a periodic structure involving resonant elements. The resonances are associated with fast-oscillating fields in thin-bridges of SRR perforations of the plate, and we filter these oscillations by introducing a third scale in the usual two-scale expansion:

$$\forall \mathbf x \in \Omega _ { f }, \ U _ { \eta } ( \mathbf x ) \, = \, & U _ { 0 } ( \mathbf x, \frac { \mathbf x } { \eta }, \frac { x _ { 2 } } { \eta ^ { 2 } } ) + \eta U _ { 1 } ( \mathbf x, \frac { \mathbf x } { \eta }, \frac { x _ { 2 } } { \eta ^ { 2 } } ) \\ & + \eta ^ { 2 } U _ { 2 } ( \mathbf x, \frac { \mathbf x } { \eta }, \frac { x _ { 2 } } { \eta ^ { 2 } } ) + \dots \\ \wedge \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd�$$

↦

where U i : Ω f × Y × -[ h/ , h/ 2 2] -→ C is a smooth function of variables ( x y , , ξ ) = ( x , x , y 1 2 1 , y 2 , ξ ), independent of η , such that ∀ x ∈ Ω f , φ i ( x , · , ξ ) is Y -periodic, and h denotes the thickness of the thin-cut of the SRR.

The differential operator is rescaled accordingly as ∇ = ∇ x + 1 η ∇ y + 1 η 2 ∇ ξ , so that Eq. (2) can be reexpressed as

$$\ a n \ v o c \colon \ e q p e s c u a s & \quad \\ & \left \{ \sqrt { D _ { \eta } } \left ( \nabla _ { x } + \frac { 1 } { \eta } \nabla _ { y } + \frac { 1 } { \eta ^ { 2 } } \nabla _ { \xi } \right ) \cdot \left ( \nabla _ { x } + \frac { 1 } { \eta } \nabla _ { y } + \frac { 1 } { \eta ^ { 2 } } \nabla _ { \xi } \right ) + \beta _ { \eta } ^ { 2 } \right \} \\ & \times \left \{ \sqrt { D _ { \eta } } \left ( \nabla _ { x } + \frac { 1 } { \eta } \nabla _ { y } + \frac { 1 } { \eta ^ { 2 } } \nabla _ { \xi } \right ) \cdot \left ( \nabla _ { x } + \frac { 1 } { \eta } \nabla _ { y } + \frac { 1 } { \eta ^ { 2 } } \nabla _ { \xi } \right ) - \beta _ { \eta } ^ { 2 } \right \} \\ & \left \{ \sum _ { i = 0 } ^ { \infty } \eta ^ { i U _ { i } ( x, \frac { x } { \eta }, \frac { x _ { 2 } } { \eta } ) } \right \} = 0 \, \\ & \text{ith the notations } \nabla _ { x } = ( \frac { \partial } { \partial x _ { i } }, \frac { \partial } { \partial x _ { i } } ), \nabla _ { y } = ( \frac { \partial } { \partial u }, \frac { \partial } { \partial u _ { n } } ) \text{ and } \nabla _ { \xi } = ( 0, \frac { \partial } { \partial \xi } ).$$

with the notations ∇ x = ( ∂ ∂x 1 , ∂ ∂x 2 ), ∇ y = ( ∂ ∂y 1 , ∂ ∂y 2 ) and ∇ ξ = (0 , ∂ ∂ξ ). Collecting terms of same power of η in Eq. (7), we obtain the homogenized problem 43 ,

$$( \sqrt { D _ { h o m } } \nabla ^ { 2 } U _ { \eta } + \beta _ { 0 } ^ { 2 } \sqrt { \rho _ { h o m } } ) ( \sqrt { D _ { h o m } } \nabla ^ { 2 } U _ { \eta } - \beta _ { 0 } ^ { 2 } \sqrt { \rho _ { h o m } } ) U _ { \eta } = 0 \, \ \text{in} \ \Omega _ { f } \ \langle \overline { \Theta _ { \eta } } \, \quad \quad ( 8 ) \\ \text{in the limit when $n$ tends to zero. in the natonic crystal $\Omega_{\cdot}$ This is the homoserized}$$

in the limit when η tends to zero, in the platonic crystal Ω . f This is the homogenized biharmonic equation with the homogenized plate rigidity D hom = ( ∫ Y ∗ D -1 ( y , y 1 2 ) dy dy 1 2 ) -1 where Y ∗ = Y \ ¯ . C We point out that the circular geometry of SRR is essential, if one would take elliptical SRR, then D hom would be a rank-2 anisotropic tensor.

The homogenized density ρ hom is given in the form

$$\sqrt { \rho _ { h o m } } ( \beta ) = 1 - \sum _ { m = 1 } ^ { \infty } \frac { \beta ^ { 2 } } { \beta ^ { 2 } - \beta _ { m } ^ { 2 } } \| V _ { m } \| _ { L ^ { 2 } ( 0, l ) } ^ { 2 } \,, \\ \intertext { m a t i n d } \lambda \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{matrix}$$

where the eigen-solutions V m correspond to longitudinal vibrations of the thin cut Π η within the split ring resonator.

It is clear from Eq. (9) that √ ρ hom ( β ) takes negative values near resonances β 2 = β 2 m . For our purpose, it is enough to look at the first few resonant frequencies β 2 m (the higher the frequency the worse the asymptotic approximation). These frequencies are associated with vibrations V m of the thin domain Π η 39 :

$$V _ { m } ^ { \prime \prime } ( x _ { 1 } ) + \beta _ { m } ^ { \ 2 } V _ { m } ( x _ { 1 } ) \, = \, 0 \,, \, 0 < x _ { 1 } < l \,,$$

$$V _ { m } ( 0 ) \, = \, 0$$

$$\eta h V _ { m } ^ { \prime } ( l ) \, = \, \text{area} ( \Xi ) \beta _ { m } ^ { \, 2 } V _ { m } ( l ) \,,$$

where ηh and l are the thickness and the length of the thin ligament Π , and Ξ is the central η disc within the SRR. The ligament Π η is connected to Ξ, hence V ( ) = l V , where V is the vibration of the stress-free air cavity Ξ. Note that the derivation of Eq. (12) required a boundary layer analysis, and we refer to 44 for more details.

The solution of the problem (10)-(12) has the form

$$V _ { m } ( x _ { 1 } ) = A \sin ( \beta _ { m } x _ { 1 } )$$

where β m (square root of frequency) is given as the solution of the following equation

$$\eta h \cot ( \beta _ { m } l ) = \text{area} ( \Xi ) \beta _ { m } \.$$

Note also that a similar problem to Eqs. (10)-(12) is deduced from Eq. (7) with a minus sign in (10), which leads to a sinh function in (14), but this does not give any additional resonant frequencies. The resonant frequencies of √ ρ hom also lead to negative values of ρ hom , in the same way that a product of two simultaneously negative square roots give a negative refractive index in the metamaterial's literature 1,33 .

In the sequel, we focus our analysis on the Floquet-Bloch bending wave problem and give some numerical results and comparisons with the homogenization theory.

## III. NUMERICAL ANALYSIS OF SRR PLATONIC CRYSTALS

## A. Dispersive properties

FIG. 2: (a) Band diagram for the doubly periodic square lattice representing the normalized frequency β 2 versus the Bloch vector k in the first Brillouin zone ΓXM. (b) Map of the eigen-mode corresponding to the eigen-frequency β 2 = 2 95 in a unit cell for a SRR with an outer radius 0 48 . . d and inner radius 0 1 . d and thickness of the thin ligament equal to 0 05 . d .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-49b43365-8139-4939-9ad3-b9ed3344b01d.jpg)

To investigate numerically the stop-band properties for out-of-plane bending waves propagating within the array of SRR, we use the Finite Element Method (by implementing the weak form of (2) in the commercial software COMSOL and its corresponding Floquet-Bloch boundary conditions as well as the PMLs: perfectly matched layers).

In Fig. 2(a), we give the band diagram for normalized eigen-frequencies β 2 as a function of the projection of the Bloch vector k on the first Brillouin zone ΓXM. We consider a square array of normalized pitch d =1 (take for example d = 1cm for comparisons with feasible experiments) with embedded SRR of inner radius 0 1 . d and outer radius 0 48 . d with a thin cut (ligament) of thickness 0 05 . d .

This figure displays two full platonic band gaps for the range of normalized frequencies [2 39 . , 2 92] and [4 23 . . , 24 32]. . As depicted in a previous study 39 , the low frequency band gap (the first one) is associated with the resonant modes of a single SRR (microscopic structure),

FIG. 3: (a) Band diagram of surface mode along the ΓX direction. (b) 1D profile along the x 2 direction of the mode. (c) Supercell showing the configuration where the plasmon mode was excited from the top layer with PMLs in the bottom layer to avoid any reflection and Floquet-Bloch boundary conditions on the left and right sides to account for the infinite array in x 1 direction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-2f3a8ba4-a5ad-4ed4-956f-fb0eed090c88.jpg)

while the second one is due to a Bragg scattering phenomenon (macroscopic structure). We also plot in Fig. 2(b) the corresponding localized eigen-function (the first eigen-mode) which corresponds in the context of continuum mechanics, to oscillations of the central region of the SRR as a rigid solid connected by the thin cut Π ε to the fixed rigid region around it 44 . In our context of bending waves, this can be interpreted as a collective vibration of the plate elements together, up and down in each SRR.

The transcendental equation (14) can be further simplified if we look at the first low

frequency, for which we deduce the explicit asymptotic approximation

$$\beta _ { 1 } ^ { 2 } \sim \frac { \eta h } { \text{area} ( \Xi ) l } \.$$

Our numerical estimate is

$$\psi _ { 1 } ^ { 2 } \sim \frac { 0. 0 5 } { 0. 3 8 } \frac { 1 } { \pi 0. 1 ^ { 2 } } = 4. 1 9$$

which is in excellent agreement with the finite element value for the plasmon frequency β 2 = 4 23 occurring at the . M point when k = ( π, π ).

It is interesting to also analyze the dispersive properties of surface elastic modes (waves exponentially localized at the interface of a platonic crystal). These modes originate from interference effect in phononic or platonic crystals and are important, for example, for superlensing effect. The dispersion diagram of these modes is shown in Fig. 3(a), where the frequency is plotted versus Bloch vector in the ΓM direction. The red dotted-dashed line reprsents the surface mode's dispersion and shows that's standing in the elastic bandgap (between the first and the second propagating modes). This means, for instance, that this mode is of evanescent nature. Figure 3(b) shows the flexural displacement at the interface between the platonic crystal and free-space [Fig. 3(c)], validating the decaying nature of the mode.

## B. (Super-)Lensing effect via negative effective density of a platonic crystal

We then consider a 2D finite platonic crystal perforating a thin elastic plate. The crystal consists of 230 holes with a C-shaped cross-section with identical parameters as those of Fig. 2(b). The crystal has a four-fold (square) symmetry with pitch d = 1. A flexural waves source of wavenumber β = 1 98 is located at a distance 0 1 . . d above the top of the crystal. The field-maps in Fig. 4 (a) show the existence of an image of the source, located twice the width of the crystal away in the lower side. As we can see in Fig. 4 (a), the amplitude of the field is the highest near the thin ligaments Π ε of the SRR. Formula (9) shows that the resonance of the field with the microscopic structure of SRR is responsible for the appearance of negative refraction through negative effective density ρ hom . Finally, the resolution of the image δ is enhanced compared to that obtained using an array of perforations with a circular or square cross-section, as demonstrated in Fig. 4(b), (c), whereby the full width at half maximum of the image point δ ≈ λ/ 3.

FIG. 4: (a) Pattern of the energy of the bending wave at the frequency β = 1 98 where AANR . occurs. (b) and (c) give respectively the profile of the energy along the dotted white line at the source (red curve) and the image (green curve) axis, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-72fad5ad-f77e-425e-8b70-1026d703a8be.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-7939fd0b-9665-4e02-93b6-73b16c83ed2c.jpg)

We note that at the plasmon frequency, the vibration of the thin ligaments is enhanced, which is in accordance with the behavior of the field as observed on Fig. 4. It is therefore important to incorporate the resonances within the thin-bridges as we did, to model the frequency at which negative refraction occurs.

## C. Gaussian beam at oblique incidence on a platonic crystal

The interaction of a Gaussian beam (GB) with the platonic crystal around the AANR is considered in this section. Refraction of GB (whose size is in the order of the wavelength) obliquely incident on the SRR platonic crystal is demonstrated to be negative.

As known, considering the interaction of a GB with any optical or mechanical media, reveals much of their propagating and scattering properties 45,46 . These are paraxial solutions of the biharmonic scalar wave equation (∆ 2 -β 4 ) W = 0, and thus have the form

$$U ( x, y ) = u ( x, y ) e ^ { - j \beta x } + n. p.$$

FIG. 5: Bending displacement when the platonic crystal is excited by a gaussian beam at an incidence angle θ inc = 25 degrees. (b) Snapshot of the bending wave in the presence of an effective medium with negative elastic index n e = -1 5 + 10 . -3 j , with n e = √ 12 ρ (1 -ν 2 ) / EH ( 2 ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-88aa71f2-fcfe-4bb5-91e1-eb0a18d287e6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-f4fb9efe-9482-4f52-9a02-cb5e496b2989.jpg)

9 with the wavenumber β , the coefficient u x ( = 0 , y ) has the well know Gaussian distribution in the x -y plane. n.p. denotes a non propagating contribution due to the fourth order nature of the biharmonic equation. It could be neglected in the real implementation. The expression of the obtained elastic beam is thus

$$u ( x, y ) = u _ { 0 } \frac { \omega _ { 0 } } { \omega ( x ) } e ^ { \left ( - \frac { y ^ { 2 } } { \omega ^ { 2 } ( x ) } - i k \frac { y ^ { 2 } } { 2 R ( x ) } + i \zeta ( x ) \right ) } \,,$$

with ω x ( ) = ω 0 √ 1 + ( x x R ) 2 , x R = πω 2 0 λ and ω 0 is called the waist of the beam. The radius of the curvature is given by R x ( ) = x [ 1 + ( x R x ) 2 ] and the Gouy phase shift is ζ ( x ) = atan x x R . Figure 5 gives the snapshot of the displacement field U in the presence of the platonic crystal at oblique incidence with angle θ inc = 25 degrees; a reflected beam could be distinguished as well as a refracted one which turns out to be in the negative refraction regime. To further verify this claim, another simulation using, now, an homogeneous effective slab of the same size, with negative elastic index n e = ω/β 2 = √ 12 ρ (1 -ν 2 ) / EH ( 2 ), is shown in Fig. 5(b). The patterns of the wave are identical in 5(a) and 5(b) and show that the elastic index of

the SRR crystal could be described as negative: n e = -1 5 + 10 . -3 j . The robustness of the AANR was further verified by changing the angle of incidence θ inc and the wavelength of operation around their main values, and the effect still could be observed.

## D. Ultrarefraction effect via vanishing effective density of a platonic crystal

FIG. 6: Pattern of the elastic displacement (a) and energy of the bending wave (b) at the frequency where ultrarefraction occurs when a point source is located at the center of the platonic crystal.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-84e91dd1-9302-49cc-8a41-c3af0bab1466.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0004v1-figure-33cdbd46-ddc3-44d7-8f64-55a2301f94aa.jpg)

To conclude this section, let us consider another anomalous refractive effect, namely, the ultra-refraction of bending waves. This phenomenon generally occurs when the effective elastic index n e defined above, tends to zero or equally speaking the group velocity v g ≈ 0. In this case, independently of the angle of incidence of an incident ray, propagating towards the crystal, it will be refracted normally to it (parallel to the normal). And the operation point could be easily deduced from the dispersion diagram by picking a point where the mode's dispersion curve becomes flat ( ∇ k β 2 ≈ 0). This allows us to build an elastic antenna (for instance by placing a cylindrical source of flexural waves in the center of the crystal). The wave will be then refracted along the normal to the crystal, permitting thus highly directive beam as could be seen from Fig. 6(a) and (b) where the displacement snapshot and its energy are respectively plotted around a wavelength λ = 2 π/β ≈ 3. This corresponds

as previously stated to the M point in the diagram of Fig. 2(a) at the maximum. Using formula (9), one can indeed achieve ultrarefraction when the effective density of the platonic crystal ρ hom is close to zero, which happens for β 2 2 ∼ 4 37. . To get this estimate, one need replace the boundary condition (11) by V m (0) = -area(Σ) / area(Y -C) which amounts to assuming that the field takes the same non zero values on the opposite edges of the basic cell. This leads to the new frequency estimate

$$\beta _ { 2 } ^ { 2 } & \sim \frac { \eta h } { \text{area} ( \Xi ) l } \left ( 1 + \frac { \text{area} ( \Sigma ) } { \text{area} ( Y - C ) } \right ) \,, \\. \quad \tilde { \lambda } \, \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \tilde { \lambda } \, \quad \ t \, \end{$$

which is a refinement of (15). The numerical result β 2 2 ∼ 4 37 is in . good agreement with finite element computations (4.385).

## IV. CONCLUDING REMARKS

In conclusion, we have proposed an original route towards elastic negative refraction (superlensing) and ultrarefraction based on excitation of flexural surface modes. To analyze these effects, we have derived the homogenized biharmonic equation using a multi-scale asymptotic approach. We found that the homogenized elastic parameters are described by a scalar flexural plate rigidity (since the perforations have a circular geometry owing to the fact that thin-bridges come unseen in the second scale asymptotics) and a scalar density, the latter taking negative values near thin-bridges resonances (unveiled by the third scale). We then performed numerical computations based on the finite element method which confirmed that the asymptotic model gives accurate predictions on dispersive properties of the elastic Lamb modes.

We believe that such a micro-structured plate could be manufactured easily, having in mind some potential applications in superlensing or directive elastic antennas. The range of industrial applications is vast, and our proof of concept should foster research efforts in this emerging area of acoustic and seismic metamaterials.

## Acknowledgments

S.G. would like to acknowledge a funding of the European Research Council through ERC grant ANAMORPHISM.

- 1 V. G. Veselago, Usp. Fiz. Nauk 92 , 517 (1967); V. G. Veselago, Sov. Phys. Usp. 10 , 509 (1968).
- 2 J.B. Pendry, Phys. Rev. Lett. 86 , 3966-3969 (2000).
- 3 D.R. Smith and N. Kroll, Phys. Rev. Lett. 85 , 2933 (2000).
- 4 J.B. Pendry, A.J. Holden, D.J. Robbins, W.J. Stewart, IEEE Trans. Microwave Theory Tech. 47 , 2075 (1999).
- 5 S. Datta, Classical Wave Propagation in Periodic and Random Media (Iowa State University, 1994).
- 6 E. Yablonovitch, Phys. Rev. Lett. 58 , 2059 (1987).
- 7 S. John, Phys. Rev. Lett. 58 , 2486 (1987).
- 8 T. F. Krauss, R. M. DeLaRue, and S. Brand, Nature 383 , 6602 (1996).
- 9 R. Zengerle, J. Mod. Opt. 34 , 1589 (1987).
- 10 B. Gralak, S. Enoch, and G. Tayeb, J. Opt. Soc. Am. A 17 , 1012 (2000).
- 11 M. Notomi, Opt. Quantum Electron. 34 , 133 (2002).
- 12 C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, Phys. Rev. B 65 , 201104(R) (2002).
- 13 C.G. Poulton, A.B. Movchan, R.C. McPhedran, N.A. Nicorovici, and Y.A. Antipov, Proc. R. Soc. London, Ser. A 456 , 2543-2559 (2000).
- 14 C.G. Poulton, R.C. McPhedran,N.A. Nicorovici, L.C. Botten, and A.B. Movchan, IUTAM Symposium on Mechanical and Eelectromagnetic Waves in Structured Media, (2001).
- 15 L. Feng, X. Liu, Y. Chen, Z. Huang, Y. Mao, Y. Chen, J. Z, and Y. Zhu, Phys. Rev. B 72 , 033108 (2005).
- 16 X. Zhang and Z. Liu, Appl. Phys. Lett. 85 , 341 (2004).
- 17 X. Hu, Y. Shu, X. Liu, R. Fu, and J. Zi, Phys. Rev. E 69 , 030201(R) (2004).
- 18 L. Feng, X. P. Liu, M. H. Lu, Y. B. Chen, Y. F. Chen, Y. W. Mao, J. Zi, Y. Y. Zhu, S. N. Zhu, and N. B. Ming, Phys. Rev. B 73 , 193101 (2006).

- 19 T. Chou, Phys. Rev. Lett. 79 , 4802 (1997).
- 20 T. Chou, J. Fluid Mech. 369 , 333 (1998).
- 21 P. McIver, J. Fluid Mech. 424 , 101 (2000).
- 22 M. Torres, J. P. Adrados, F. R. Montero de Espinosa, D. Garcia-Pablos, and J. Fayos, Phys. Rev. E 63 , 011204 (2000).
- 23 X. Hu, Y. Shen, X. Liu, R. Fu, J. Zi, X. Jiang, and S. Feng, Phys. Rev. E 68 , 037301 (2003).
- 24 N.V. Movchan, R.C. McPhedran, A.B. Movchan, G.G. Poulton, Proc. R. Soc. A-Math. Phys. Eng. Sci. 465 , 3383-3400 (2009).
- 25 A.B. Movchan, N.V. Movchan, R.C. McPhedran, Proc. R. Soc. A-Math. Phys. Eng. Sci. 463 , 2505-2518 (2007).
- 26 M. Farhat, S. Guenneau, S. Enoch, A.B. Movchan, and G. Petursson, Appl. Phys. Lett. 96 , 081909 (2010).
- 27 M. Farhat, S. Guenneau, S. Enoch, and A.B. Movchan, EPL-Europhys. Lett. 91 , 54003 (2010).
- 28 M. Farhat, S.Guenneau, and S. Enoch, Phys. Rev. Lett. 103 , 024301 (2009).
- 29 J. B. Pendry, Phys. Rev. Lett. 85 , 3966 (2000).
- 30 D. Maystre and S. Enoch, J. Opt. Soc. Am. 21 , 122 (2004).
- 31 S. Guenneau, A. C. Vutha, and S. A. Ramakrishna, New J. Phys. 7 , 164 (2005).
- 32 D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, Phys. Rev. Lett. 84 , 4184 (2000).
- 33 S. A. Ramakrishna, Rep. Prog. Phys. 68 , 449 (2005).
- 34 S.A. Cummer, D. Schurig, New J. Phys. 9 , 45 (2007).
- 35 A.N. Norris, Proc. R. Soc. Lond. A 464 2411 (2008).
- 36 M. Farhat, S. Enoch, S. Guenneau, and A. B. Movchan, Phys. Rev. Lett. 101 , 134501 (2008).
- 37 G. Dupont, M. Farhat, A. Diatta, S. Guenneau, and S. Enoch, Wave Motion 48 , 483-496 (2011).
- 38 M. Farhat, S. Guenneau, S. Enoch, G. Tayeb, A. B. Movchan, and N. V. Movchan, Phys. Rev. E 77 , 046308 (2008).
- 39 A.B. Movchan, S. Guenneau, Phys. Rev. B 70 , 125116 (2004).
- 40 S. Guenneau, A.B. Movchan, N.V. Movchan, Physica B 394 , 301648 (2006).
- 41 S. Timoshenko, Theory of plates and shells (McGraw-Hill, New York, 1940).
- 42 K.F. Graff, Wave motion in elastic solids (Dover, New York, 1975).
- 43 V.V. Jikhov, S.M. Kozlov, and O. A. Oleinik, Homogenization of Differential Operators and

Integral Functionals (Springer Verlag, New York, 1994).

- 44 V. Kozlov, V. Mazya, and A. B. Movchan, Fields in Multistructures. Asymptotic Analysis (Oxford University Press, Oxford, 1999).
- 45 R.W. Ziolkowski, Opt. Express 11 , 662-681 (2003).
- 46 L.W. Cai and J. Sanchez-Dehesa, New J. Phys. 9 , 450 (2007).
M. Farhat S. Enoch S. Guenneau
Biharmonic Split Ring Resonator Metamaterial: Artificially dispersive effective density in thin periodically perforated plates
We present in this paper a theoretical and numerical analysis of bending
waves localized on the boundary of a platonic crystal whose building blocks are
split ring resonators (SRR). We first derive the homogenized parameters of the
structured plate using a three-scale asymptotic expansion in the linearized
biharmonic equation. In the limit when the wavelength of the bending wave is
much larger than the typical heterogeneity size of the platonic crystal, we
show that it behaves as an artificial plate with an anisotropic effective Young
modulus and a dispersive effective mass density. We then analyze dispersion
diagrams associated with bending waves propagating within an infinite array of
SRR, for which eigen-solutions are sought in the form of Floquet-Bloch waves.
We finally demonstrate that this structure displays the hallmarks of
All-Angle-Negative-Refraction(AANR) and it leads to superlensing and
ultrarefraction effects, interpreted thanks to our homogenization model as a
consequence of negative and vanishing effective density, respectively.