Mechanical metamaterials.

Mechanical metamaterials, also known as architected materials, are rationally designed composites, aiming at elastic behaviors and effective mechanical properties beyond ('meta') those of their individual ingredients-qualitatively and/or quantitatively. Due to advances in computational science and manufacturing, this field has progressed considerably throughout the last decade. Here, we review its mathematical basis in the spirit of a tutorial, and summarize the conceptual as well as experimental state-of-the-art. This summary comprises disordered, periodic, quasi-periodic, and graded anisotropic functional architectures, in one, two, and three dimensions, covering length scales ranging from below one micrometer to tens of meters. Examples include extreme ordinary linear elastic behavior from artificial crystals, e.g. auxetics and pentamodes, 'negative' effective properties, behavior beyond classical linear elasticity, e.g. arising from local resonances, chirality, beyond-nearest-neighbor interactions, quasi-crystalline mechanical metamaterials, topological band gaps, cloaking based on coordinate transformations and on scattering cancelation, seismic protection, nonlinear and programmable metamaterials, as well as space-time-periodic architectures.

Active exterior cloaking for the two-dimensional Helmholtz equation with complex wavenumbers and application to thermal cloaking.

We design sources for the two-dimensional Helmholtz equation that can cloak an object by cancelling out the incident field in a region, without the sources completely surrounding the object to hide. As in previous work for real positive wavenumbers, the sources are also determined by the Green identities. The novelty is that we prove that the same approach works for complex wavenumbers which makes it applicable to a variety of media, including media with dispersion, loss and gain. Furthermore, by deriving bounds on Graf's addition formulas with complex arguments, we obtain new estimates that allow to quantify the quality of the cloaking effect. We illustrate our results by applying them to achieve active exterior cloaking for the heat equation. This article is part of the theme issue 'Wave generation and transmission in multi-scale complex media and structured metamaterials (part 2)'.

Active thermal cloaking and mimicking.

We present an active cloaking method for the parabolic heat (and mass or light diffusion) equation that can hide both objects and sources. By active, we mean that it relies on designing monopole and dipole heat source distributions on the boundary of the region to be cloaked. The same technique can be used to make a source or an object look like a different one to an observer outside the cloaked region, from the perspective of thermal measurements. Our results assume a homogeneous isotropic bulk medium and require knowledge of the source to cloak or mimic, but are in most cases independent of the object to cloak.

Hybrid topological guiding mechanisms for photonic crystal fibers.

We create hybrid topological-photonic localisation of light by introducing concepts from the field of topological matter to that of photonic crystal fiber arrays. S-polarized obliquely propagating electromagnetic waves are guided by hexagonal, and square, lattice topological systems along an array of infinitely conducting fibers. The theory utilises perfectly periodic arrays that, in frequency space, have gapped Dirac cones producing band gaps demarcated by pronounced valleys locally imbued with a nonzero local topological quantity. These broken symmetry-induced stop-bands allow for localised guidance of electromagnetic edge-waves along the crystal fiber axis. Finite element simulations, complemented by asymptotic techniques, demonstrate the effectiveness of the proposed designs for localising energy in finite arrays in a robust manner.

Isotropic Chiral Acoustic Phonons in 3D Quasicrystalline Metamaterials.

The elastic properties of three-dimensional (3D) crystalline mechanical metamaterials, unlike those of amorphous structures, are generally strongly anisotropic-even in the long-wavelength limit and for highly symmetric crystals. Aiming at isotropic linear elastic wave propagation, we therefore study 3D periodic approximants of 3D icosahedral quasicrystalline mechanical metamaterials consisting of uniaxial chiral metarods. Considering the increasing order of the approximants, we approach nearly isotropic effective speeds of sound and isotropic acoustical activity. The latter is directly connected to circularly polarized 3D metamaterial chiral acoustic phonons-for all propagation directions in three dimensions.

Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings.

This study follows from Maurel et al., Phys. Rev. B 98, 134311 (2018), where we reported on direct numerical observations of out-of-plane shear surface waves propagating along an array of plates atop a guiding layer, as a model for a forest of trees. We derived closed form dispersion relations using the homogenization procedure and investigated the effect of heterogeneities at the top of the plates (the foliage of trees). Here, we extend the study to the derivation of a homogenized model accounting for heterogeneities at both endings of the plates. The derivation is presented in the time domain, which allows for an energetic analysis of the effective problem. The effect of these heterogeneous endings on the properties of the surface waves is inspected for hard heterogeneities. It is shown that top heterogeneities affect the resonances of the plates, hence modifying the cut-off frequencies of a wave mathematically similar to the so-called Spoof Plasmon Polariton (SPP) wave, while the bottom heterogeneities affect the behavior of the layer, hence modifying the dispersion relation of the Love waves. The complete system simply mixes these two ingredients, resulting in hybrid surface waves accurately described by our model.

Cloaking In-Plane Elastic Waves with Swiss Rolls.

We propose a design of cylindrical cloak for coupled in-plane shear waves consisting of concentric layers of sub-wavelength resonant stress-free inclusions shaped as Swiss rolls. The scaling factor between inclusions' sizes is according to Pendry's transform. Unlike the hitherto known situations, the present geometric transform starts from a Willis medium and further assumes that displacement fields u in original medium and u ' in transformed medium remain unaffected ( u ' = u ). This breaks the minor symmetries of the rank-4 and rank-3 tensors in the Willis equation that describe the transformed effective medium. We achieve some cloaking for a shear polarized source at specific, resonant sub-wavelength, frequencies, when it is located in close proximity to a clamped obstacle surrounded by the structured cloak. The structured medium approximating the effective medium allows for strong Willis coupling, notwithstanding potential chiral elastic effects, and thus mitigates roles of Willis and Cosserat media in the achieved elastodynamic cloaking.

Fresnel drag in space-time-modulated metamaterials.

A moving medium drags light along with it as measured by Fizeau and explained by Einstein's theory of special relativity. Here we show that the same effect can be obtained in a situation where there is no physical motion of the medium. Modulations of both the permittivity and permeability, phased in space and time in the form of traveling waves, are the basis of our model. Space-time metamaterials are represented by effective bianisotropic parameters, which can in turn be mapped to a moving homogeneous medium. Hence these metamaterials mimic a relativistic effect without the need for any actual material motion. We discuss how both the permittivity and permeability need to be modulated to achieve these effects, and we present an equivalent transmission line model.

Topological beam-splitting in photonic crystals.

We create a passive wave splitter, created purely by geometry, to engineer three-way beam splitting in electromagnetism in transverse electric and magnetic polarisation. We do so by considering arrangements of Indium Phosphide dielectric pillars in air, in particular we place several inclusions within a cell that is then extended periodically upon a square lattice. Hexagonal lattice structures are more commonly used in topological valleytronics but, as we discuss, three-way splitting is only possible using a square, or rectangular, lattice. To achieve splitting and transport around a sharp bend we use accidental, and not symmetry-induced, Dirac cones. Within each cell pillars are either arranged around a triangle or square; we demonstrate the mechanism of splitting and why it does not occur for one of the cases. The theory is developed and full scattering simulations demonstrate the effectiveness of the proposed designs.

Controlling frequency dispersion in electromagnetic invisibility cloaks.

Electromagnetic cloaking, as challenging as it may be to the physicist and the engineer has become a topical subject over the past decade. Thanks to the transformations optics (TO) invisibility devices are in sight even though quite drastic limitations remain yet to be lifted. The extreme material properties which are deduced from TO can be achieved in practice using dispersive metamaterials. However, the bandwidth over which a metamaterial cloak is efficient is drastically limited. We design and simulate a spherical cloak which takes into account the dispersive nature of relative permittivity and permeability tensors realized by plasma-like metamaterials. This spherical cloak works over a broad frequency-band even though these materials are of a highly dispersive nature. We establish two equations of state that link the eigenvalues of the permittivity and permeability tensors in every spherical cloak regardless of the geometrical transformation. Frequency dispersive properties do not disrupt cloaking as long as the equations of states are satisfied in the metamaterial cloak.

Frequency domain transformation optics for diffusive photon density waves' cloaking.

We make use of transformation optics technique to realize cloaking operation in the light diffusive regime, for spherical objects. The cloak requires spatially heterogeneous anisotropic diffusivity, as well as spatially varying speed of light and absorption. Analytic calculations of Photon's fluence confirm minor role of absorption in reduction of far-field scattering, and a monopole fluence field converging to a constant in the static regime in the invisibility region. The latter is in contrast to acoustic and electromagnetic cloaks, for which the field vanishes inside the core. These results are finally discussed in the context of mass diffusion, where cloaking can be achieved with a heterogeneous anisotropic diffusivity.

Flat lens effect on seismic waves propagation in the subsoil.

We show that seismic energy simulated by an artificial source that mainly propagates Rayleigh surface waves, is focused in structured soil made of a grid of holes distributed in the ground. We carry out large-scale field tests with a structured soil made of a grid consisting of cylindrical and vertical holes in the ground and a low frequency artificial source (<10 Hz). This allows the identification of a distribution of energy inside the grid, which can be interpreted as the consequence of a dynamic anisotropy akin to an effective negative refraction index. Such a flat lens reminiscent of what Veselago and Pendry envisioned for light opens avenues in seismic metamaterials to counteract partially or totally the most devastating components of seismic signals.

Spanning the scales of mechanical metamaterials using time domain simulations in transformed crystals, graphene flakes and structured soils.

We begin with a brief historical survey of discoveries of quasi-crystals and graphene, and then introduce the concept of transformation crystallography, which consists of the application of geometric transforms to periodic structures. We consider motifs with three-fold, four-fold and six-fold symmetries according to the crystallographic restriction theorem. Furthermore, we define motifs with five-fold symmetry such as quasi-crystals generated by a cut-and-projection method from periodic structures in higher-dimensional space. We analyze elastic wave propagation in the transformed crystals and (Penrose-type) quasi-crystals with the finite difference time domain freeware SimSonic. We consider geometric transforms underpinning the design of seismic cloaks with square, circular, elliptical and peanut shapes in the context of honeycomb crystals that can be viewed as scaled-up versions of graphene. Interestingly, the use of morphing techniques leads to the design of cloaks with interpolated geometries reminiscent of Victor Vasarely's artwork. Employing the case of transformed graphene-like (honeycomb) structures allows one to draw useful analogies between large-scale seismic metamaterials such as soils structured with columns of concrete or grout with soil and nanoscale biochemical metamaterials. We further identify similarities in designs of cloaks for elastodynamic and hydrodynamic waves and cloaks for diffusion (heat or mass) processes, as these are underpinned by geometric transforms. Experimental data extracted from field test analysis of soil structured with boreholes demonstrates the application of crystallography to large scale phononic crystals, coined as seismic metamaterials, as they might exhibit low frequency stop bands. This brings us to the outlook of mechanical metamaterials, with control of phonon emission in graphene through extreme anisotropy, attenuation of vibrations of suspension bridges via low frequency stop bands and the concept of transformed meta-cities. We conclude that these novel materials hold strong applications spanning different disciplines or across different scales from biophysics to geophysics.

Enhanced sensing and conversion of ultrasonic Rayleigh waves by elastic metasurfaces.

Recent years have heralded the introduction of metasurfaces that advantageously combine the vision of sub-wavelength wave manipulation, with the design, fabrication and size advantages associated with surface excitation. An important topic within metasurfaces is the tailored rainbow trapping and selective spatial frequency separation of electromagnetic and acoustic waves using graded metasurfaces. This frequency dependent trapping and spatial frequency segregation has implications for energy concentrators and associated energy harvesting, sensing and wave filtering techniques. Different demonstrations of acoustic and electromagnetic rainbow devices have been performed, however not for deep elastic substrates that support both shear and compressional waves, together with surface Rayleigh waves; these allow not only for Rayleigh wave rainbow effects to exist but also for mode conversion from surface into shear waves. Here we demonstrate experimentally not only elastic Rayleigh wave rainbow trapping, by taking advantage of a stop-band for surface waves, but also selective mode conversion of surface Rayleigh waves to shear waves. These experiments performed at ultrasonic frequencies, in the range of 400-600 kHz, are complemented by time domain numerical simulations. The metasurfaces we design are not limited to guided ultrasonic waves and are a general phenomenon in elastic waves that can be translated across scales.

Finite element analysis of electromagnetic waves in two-dimensional transformed bianisotropic media.

We analyze the wave propagation in two-dimensional bianisotropic media with the Finite Element Method (FEM). Starting from the Maxwell-Tellegen's equations in bianisotropic media, we derive some system of coupled Partial Differential Equations (PDEs) for longitudinal electric and magnetic field components. These PDEs are implemented in FEM using a solid mechanics formulation. Perfectly Matched Layers (PMLs) are also discussed to model unbounded bianisotropic media. The PDEs and PMLs are then implemented in a finite element software, and transformation optics is further introduced to design some bianisotropic media with interesting functionalities, such as cloaks, concentrators and rotators. In addition, we propose a design of metamaterial with concentric layers made of homogeneous media with isotropic permittivity, permeability and magnetoelectric parameters that mimic the required effective anisotropic tensors of a bianisotropic cloak in the long wavelength limit (homogenization approach). Our numerical results show that transformation based electromagnetic metamaterials can be extended to bianisotropic media.

A seismic metamaterial: The resonant metawedge.

Critical concepts from three different fields, elasticity, plasmonics and metamaterials, are brought together to design a metasurface at the geophysical scale, the resonant metawedge, to control seismic Rayleigh waves. Made of spatially graded vertical subwavelength resonators on an elastic substrate, the metawedge can either mode convert incident surface Rayleigh waves into bulk elastic shear waves or reflect the Rayleigh waves creating a "seismic rainbow" effect analogous to the optical rainbow for electromagnetic metasurfaces. Time-domain spectral element simulations demonstrate the broadband efficacy of the metawedge in mode conversion while an analytical model is developed to accurately describe and predict the seismic rainbow effect; allowing the metawedge to be designed without the need for extensive parametric studies and simulations. The efficiency of the resonant metawedge shows that large-scale mechanical metamaterials are feasible, will have application, and that the time is ripe for considering many optical devices in the seismic and geophysical context.

Transformation seismology: composite soil lenses for steering surface elastic Rayleigh waves.

Metamaterials are artificially structured media that exibit properties beyond those usually encountered in nature. Typically they are developed for electromagnetic waves at millimetric down to nanometric scales, or for acoustics, at centimeter scales. By applying ideas from transformation optics we can steer Rayleigh-surface waves that are solutions of the vector Navier equations of elastodynamics. As a paradigm of the conformal geophysics that we are creating, we design a square arrangement of Luneburg lenses to reroute Rayleigh waves around a building with the dual aim of protection and minimizing the effect on the wavefront (cloaking). To show that this is practically realisable we deliberately choose to use material parameters readily available and this metalens consists of a composite soil structured with buried pillars made of softer material. The regular lattice of inclusions is homogenized to give an effective material with a radially varying velocity profile and hence varying the refractive index of the lens. We develop the theory and then use full 3D numerical simulations to conclusively demonstrate, at frequencies of seismological relevance 3-10 Hz, and for low-speed sedimentary soil (vs: 300-500 m/s), that the vibration of a structure is reduced by up to 6 dB at its resonance frequency.

Forests as a natural seismic metamaterial: Rayleigh wave bandgaps induced by local resonances.

We explore the thesis that resonances in trees result in forests acting as locally resonant metamaterials for Rayleigh surface waves in the geophysics context. A geophysical experiment demonstrates that a Rayleigh wave, propagating in soft sedimentary soil at frequencies lower than 150 Hz, experiences strong attenuation, when interacting with a forest, over two separate large frequency bands. This experiment is interpreted using finite element simulations that demonstrate the observed attenuation is due to bandgaps when the trees are arranged at the sub-wavelength scale with respect to the incident Rayleigh wave. The repetitive bandgaps are generated by the coupling of the successive longitudinal resonances of trees with the vertical component of the Rayleigh wave. For wavelengths down to 5 meters, the resulting bandgaps are remarkably large and strongly attenuating when the acoustic impedance of the trees matches the impedance of the soil. Since longitudinal resonances of a vertical resonator are inversely proportional to its length, a man-made engineered array of resonators that attenuates Rayleigh waves at frequency ≤10 Hz could be designed starting from vertical pillars coupled to the ground with longitudinal resonance ≤10 Hz.

Molding acoustic, electromagnetic and water waves with a single cloak.

We describe two experiments demonstrating that a cylindrical cloak formerly introduced for linear surface liquid waves works equally well for sound and electromagnetic waves. This structured cloak behaves like an acoustic cloak with an effective anisotropic density and an electromagnetic cloak with an effective anisotropic permittivity, respectively. Measured forward scattering for pressure and magnetic fields are in good agreement and provide first evidence of broadband cloaking. Microwave experiments and 3D electromagnetic wave simulations further confirm reduced forward and backscattering when a rectangular metallic obstacle is surrounded by the structured cloak for cloaking frequencies between 2.6 and 7.0 GHz. This suggests, as supported by 2D finite element simulations, sound waves are cloaked between 3 and 8 KHz and linear surface liquid waves between 5 and 16 Hz. Moreover, microwave experiments show the field is reduced by 10 to 30 dB inside the invisibility region, which suggests the multi-wave cloak could be used as a protection against water, sonic or microwaves.

Single frequency microwave cloaking and subwavelength imaging with curved wired media.

We consider the cloaking properties of electromagnetic wired media deduced from arbitrary coordinate transformations. We propose an interpretation of invisibility via sub-wavelength imaging features. The quality of cloaking is assessed by the level of deformation of the image of a P-shaped source through the stretched wired media: the lesser the image deformation, the more effective the cloaking. We numerically and experimentally demonstrate a tetrahedral wired cloak with longer edge length about 7cm at a frequency of 1GHz (the cloak is thus subwavelength). The wired cloak has two functionalities: it can serve as a high-resolution imaging system over long distances, and it can also perform space transformations such as, but not limited to, cloaking at a single operation frequency.