Author Correction: Octet lattice-based plate for elastic wave control

Motivated by the importance of lattice structures in multiple fields, we numerically investigate the propagation of flexural waves in a thin reticulated plate augmented with two classes of metastructures for wave mitigation and guiding, namely metabarriers and metalenses. The cellular architecture of this plate invokes the well-known octet topology, while the metadevices rely on novel customized octets either comprising spherical masses added to the midpoint of their struts or variable node thickness. We numerically determine the dispersion curves of a doubly-periodic array of octets, which produce a broad bandgap whose underlying physics is elucidated and leveraged as a design paradigm, allowing the construction of a metabarrier effective for inhibiting the transmission of waves. More sophisticated effects emerge upon parametric analyses of the added masses and node thickness, leading to graded designs that spatially filter waves through an enlarged bandgap via rainbow trapping. Additionally, Luneburg and Maxwell metalenses are realized using the spatial modulation of the tuning parameters and numerically tested. Wavefronts impinging on these structures are progressively curved within the inhomogeneous media and steered toward a focal point. Our results yield new perspectives for the use of octet-like lattices, paving the way for promising applications in vibration isolation and energy focusing.

www.nature.com/scientificreports/ the octet supports a bandgap whose existence depends upon the aspect ratio of the struts, but do not investigate the underlying physics of its generating mechanism; only briefly suggested by Gerard et al. 23 . In addition, most of the state-of-the-art research on such a FCC cell resort to frequency gaps produced by attached resonators 14 or use it as a constitutive matrix of multi-phased materials 24,25 , without delving into its inherent resonant modes. The absence of an overarching dynamic assessment of this architecture in the existing literature emphasizes the novelty of our work, which proposes the octet topology as a metamaterial lattice endowed with an energy gap stemming from the bending local resonance of its beam-like members. Moreover, due to the arrangement and behavior of wave propagation through the struts, we show that the octet cells form a medium in which flexural resonances naturally coincide with the Bragg frequency -as observed within the eigenmodes permitted by Floquet-Bloch waves.
In aerospace and aircraft engineering, lattice plates are frequently used as structural components or filling in sandwich panels 1 that, whilst providing lightness to the structure 2 , are liable to significant vibrations. The propagation of such vibrations in solid thin plates, occurring in the form of Lamb waves with different polarization, has been extensively documented 26,27 resulting in a variety of metamaterial-based solutions. Prominent among these are metabarriers (or metawedges 28 ) encompassing clusters of rods whose longitudinal resonant mode couples with the flexural waves propagating in the plate and triggers the opening of bandgaps 29,30 . To enlarge the frequency band of such barriers, that owing to local resonance offer a fairly limited range of applicability, graded designs have been explored, where the spatial variation of the tuning parameter, rod height, yields a wider spectrum of attenuation associated to the so-called Rainbow trapping effect 31,32 . The theory underpinning these metastructures can be easily applied to octet-based plates, with the longitudinal resonances replaced by the bending mode of the struts and the rod height by the parameters of the beam-like members. Based on these considerations, a number of parameters exist in the literature to tune the bandgap by tailoring the metamaterials at a micro-structural level 33,34 , two of which are of particular interest for the present study, namely auxiliary point masses and joint stiffness. In their pioneering work, Krödel et al. 12 were among the first to unfold the effect of added micro-inertial elements on the dispersive behavior of frame-based lattices. By placing added masses at different positions in the constitutive lattice cell, they showed that these lead to either Bragg-type or local resonance bandgaps and that their tuning strongly influences the frequency content of the attenuation zone. Junyi et al. 14 proposed equipping the end of cantilever beams, added to primary 3D Bravais lattices, with point masses whose size governs the gap frequency independently of the geometry of the base lattice. Shortly before this, the same spherical masses were attached to the joints of a Kagome lattice by Liu et al. 35 to stiffen its connections and open a bandgap tied to the rod-node vibration.
Beyond widening the bandgaps, graded designs lay the foundation for a multitude of compelling wave control possibilities. One of these dates back to the nineteenth century and the seminal work of Maxwell on so-called flat or gradient-index (GRIN) lenses 36 . These are effective, aberration-free devices that leverage composite structures to adjust the spatial distribution of the refractive index within a bounded region and ultimately shape the incident wavefronts. The Luneburg lens 37 is a type of spherical GRIN lens that diverts any collimated beam colliding at its boundaries to a diametrically opposite focal point. Similarly to this system, underpinned by a refractive index that increases radially from the outer toward the center of the graded area, the Maxwell fish-eye lens, bends the trajectories of rays originated by a point source at its surface and focuses them in a mirrored image point. While building such lenses with ordinary composite materials appears a daunting task, employing phononic crystals [38][39][40] or metamaterials 41,42 considerably relieves this complexity. From a metamaterial perspective, the cluster of locally resonant rods with graded height, previously discussed, has proven an effective candidate for the design of GRIN lenses [43][44][45] . Its rod-resonance bandgap, besides providing attenuation, serves as tuning variable for the phase velocity of waves traveling in the metastructure and leads to applications at different scales, from acoustics 43 to seismology 46 . In terms of frame structures, the research on wave guiding is still at a primordial state, mainly confined to two-dimensional lattices [47][48][49][50] and yielding, only in a handful of cases, to more sophisticated effects. Xie et al. 51 and Zhao et al. 52 were among the few to design 2D and 3D flattened Luneburg lenses consisting of a three-dimensional truss-based square lattice, to manipulate waves in the ultrasound regime. A recent study 53 , on the other hand, has explored the propagation of elastic waves in three-dimensional lattices, including specifically the octet cell. Although similar to our investigation with respect to the choice of topology, no bandgap is reported in this case and only the anisotropic character of the structure is emphasized, with no explicit reference to shaping wave trajectories via graded octet lattices. Therefore, the lack of existing literature on gradient-index lenses designed with octet lattices further motivates our work.
Here we take the octet lattice analysis one step further by using this geometry to numerically realize a reticulated plate and two types of integrated devices for flexural wave mitigation and control. Following a dispersion analysis of the proposed standard octet cell, the bandgap and its underlying generation mechanism are revealed, resulting in the identification of two tailoring parameters. Auxiliary masses, embedded in the midpoint of the struts, or variable joint thickness are used to tune the frequency range with forbidden propagation as well as the phase velocity spatially within the lattice. To avoid redundancy of results, only those stemming from the addition of masses to the primary cell are reported here; varying the node thickness produces similar outcomes. Two metabarrier designs are discussed, the first characterized by constant auxiliary masses and the second relying on the grading strategy to trap and hinder wave transmission over a wider frequency band. Finally, Luneburg and Maxwell lenses are built by smoothly increasing the supplementary spherical masses radially from the boundaries to the center of the devices. We find the octet cell and subsequent customizations possess desirable dispersive properties easily coaxed to required frequencies, these properties are readily exploited to create a variety of devices able to manipulate the propagation of waves over a variety of scales, in both a subwavelength (lenses) and Bragg (metabarriers) scattering regime.

Results
Octet cell design, dispersion and parametric analyses. Consider a thin reticulated plate, of infinite extension, constructed from repeating unit cells, each comprising a strut network generating the well-known octet topology 20 . We choose a fundamental cell whose edges coincide with the midpoint of the octet struts, as depicted in Fig. 1a, where we define physical basis vectors a 1 and a 2 , assumed collinear with the Cartesian axes x and y. We consider struts of circular cross-section, of diameter t strut = 1 mm; the circular cross-section provides axial symmetry and prevents stress concentration that would instead arise in sharp-edged profiles, e.g., square cross-section. We take the length of our lattice vectors to be l = 3 cm and hence have a strut length of l strut = l √ 2/2 . For our computations, we use versatile plastic, with Young's modulus E = 1.7 GPa, density ρ = 450 kg/m 3 and Poisson's ratio ν = 0.3, as constitutive material of our octet cell. This choice of geometry and material facilitates the 3D printing of the designed octet cell via selective laser sintering, allowing for more extensive prototyping and experimental validation to be carried out in future studies.
It is sufficient to consider waves propagating through the fundamental cell as the phase shift from cell-to-cell is completely described by the dispersion relation between eigenfrequencies and Bloch-wavevector, corresponding to the red dotted line reported in Fig. 2a (see "Methods" section). Despite the complexity of this folded band structure, three main frequency bands can be clearly identified.
In the low-frequency limit, from 0 to 2 kHz, the spectrum is dominated by zero-order modes similar to the solutions of Lamb's characteristic equations in solid, homogeneous, and thin plates for long wavelengths 26 . The underlying lattice structure has subwavelength behavior and wave propagation closely resembles an equivalent continuum supporting three propagating waves, respectively two with transverse and one with longitudinal polarization 54 (see Sect. 2 in Supplementary Material). Among these, the asymmetric flexural mode, simply denoted as A 0 mode in Fig. 2a, is the most attractive for its ease of activation and amount of carried energy hence, the cornerstone of this study.
At higher frequencies (2 to 5 kHz), the wavelength reduces rapidly and waves interacting with the detailed geometry of the cell open a bandgap between 3.85 and 5.08 kHz. Waves with frequency content matching this range are prevented from penetrating the octet lattice while being trapped by the bending vibrations of its struts. These deformations emerge as flat bands in the lower bound of the forbidden interval and serve as evidence of www.nature.com/scientificreports/ its local resonance-based generating mechanism 55 . They cluster around the first flexural frequency of a clampedclamped Euler-Bernoulli (EB) beam, geometrically equivalent to the struts in the octet and defined by: The analytical flexural frequency in Eq. (1), equal to 3.85 kHz, corresponds to the bandgap opening bound and poses as a design variable for tailoring the bandgap frequency range. By considering the octet struts as equivalent EB beams with clamped edges and bending resonance 13 , the bandgap opening frequency can be estimated via Eq. (1) and, consequently, also the correlated attenuation zone 56 . This simplification reduces the complexity of the octet cell to individual beams clamped at the nodes, facilitating the tuning of design parameters. This clamped-clamped nodal behavior relies on a network of highly connected struts peculiar of the octet topology. As outlined in Wang et al. 56 highly connected lattices (here intersection points of 12 members at a node) yield a clamped-clamped behavior of the struts at the nodes. Similarly to Nolde et al. 57 , our flat bands at f clamped correspond to standing waves induced by the fundamental resonance of vibrations along clamped-clamped beams. This behavior is highlighted in the BG opening eigenmode displayed in Fig. 2b, where the nodes stand still. Furthermore, such a fundamental resonance fits exactly within the length of the primitive lattice vectors (a half our conventional unit cell), hence the propagation regime is identical to Kaina et al. 8 in which the wavelength of the fundamental clamped-clamped resonance is twice that of the modulus of the primitive lattice vectors; subsequently the Bragg frequency and f clamped coalesce, resulting in destructive interferences producing a broad bandgap. Wang et al. 56 remarks that bandgaps, in this setting, are generated by local resonances and not Bragg scattering; however, in the case where beams are parallel to and of the same length as the primitive lattice vectors, one expects a medium in which the fundamental resonances naturally match the Bragg frequencies, all of which are properties the underlying octet cell possesses. Notice that, a complete bandgap is obtained when computing the band structure over the full two-dimensional Irreducible Brillouin zone, whose lower and upper bounds are still enforced by the two modes in Fig www.nature.com/scientificreports/ bound of the attenuation zone (see top inset in Fig. 2b), in contrast to the two modes converging at 4.7 kHz, which exhibit, instead, in-plane deformations. One known limitation of the Bloch analysis is that only a minority of the computed modes are triggered owing to their strong source dependence. Hence, a finite octet-based plate, consisting of 120 × 40 cells equivalent to the one of the infinite study (Fig. 1b), is adopted to corroborate the predicted band structure via time-domain simulations (see "Methods" section). The outcome of this analysis is shown as a background of the dispersion curves in Fig. 2a, where the branches involving out-of-plane deformations are predominantly excited. Almost complete superposition can be appreciated between the A 0 and high-frequency modes of the frequency-wavenumber spectrum and that of the Bloch analysis. A small discrepancy arises only in the higher-order modes and in the upper bound of the bandgap, with a downward shift of about 0.05 kHz in the results of the time-transient simulation. From 2 to 4 kHz, additional modes, generated by waves with polarization in the lattice plane (longitudinal or shear), are effectively captured by the finite model. Overall, the results of these two types of analyses reveal good agreement and justify their joint application in the following sections.
The octet unit outlined hitherto, comprising 28 struts of constant, circular cross-section, is denoted as the standard cell. Remarkably, a local resonance-induced bandgap originates in this single-block lattice owing to the refined topology of its design, thus paving the way for a variety of parameters to be leveraged in tuning the dispersion relations. While geometrical and material properties emerge as the primary natural candidate, their influence on the octet band structure has already been assessed 22,58 . The focus here shifts to two novel, so-called, customized octet cells whose parameters allow tailoring the band structure without affecting the original geometry, i.e., cell width and strut thickness. Similarly to the work of Krödel et al., in the customized cell of Fig. 1c the midpoint of the struts is endowed with auxiliary spherical masses, m, of diameter t mass , as this undergoes major displacements when bending deformations are activated. In order to retain a monolithic structure, easier to 3D print with modern techniques, the supplementary mass is made of the same material selected for the lattice and is expressed as a multiple of a reference value, m ref = ρA2t strut . The diameter t mass in Fig. 1c is then derived upon definition of this quantity and its practical feasibility, with respect to the thickness of the struts, is adequately verified in each of the proposed devices. The spheres are added as point masses in the finite element model in COMSOL Multiphysics and the sensitivity analysis is performed on the single customized unit via Bloch theory. Likewise the approach adopted for the standard octet cell, the struts are treated as equivalent fixed-fixed Euler-Bernoulli beams with the distributed mass, ρ A l strut , incremented by the point mass m. Naturally one expects these added masses to provide damping over the whole band structure reducing the frequencies of comparable eigenmodes 12 , as observed in Fig. 5b. The first fundamental frequency, governing the position of the bandgap, has an inverse reliance on the added mass, as already analytically proven by Low et al. 59 . When m increases, the bandgap is instantly shifted to lower frequencies, from the range 3.85-5.08 kHz to 1.73-2.74 kHz in the rightmost inset of Fig. 2c, altering the A 0 and high-frequency modes. Each point in these dispersion curves has an associated phase velocity, v p = ω/k x 5 , which decreases when moving from lighter to heavier masses at constant frequency (red line in Fig. 2c), thereby animating the lattice with slower waves. Although we were inspired by the study of Krödel et al. 12 in adding masses to our standard octet to effectively control the bandgap frequency, it is worth noting that the underlying physics of the two studies differs. Unlike the auxetic microlattice in 12 , the standard octet lattice is already endowed with a bandgap in which Bragg and local resonance coexist, regardless of the added masses.
Similar outcomes arise when exploring the customized cell in Fig. 1d. Following the analogy between octet struts and clamped EB beams, we harness the topology of the connections in the lattice to modify the boundary conditions of the equivalent beams and their natural resonant frequency, in turn affecting the band structure of the architected cell. A node area of length l node , twice as long as t strut , is identified in Fig. 1d, where the diameter of the section, t node , is expressed as a fraction of the original thickness t strut . By decreasing this parameter from t strut to 0.2t strut , clamped-like joints are replaced by quasi-hinged connections that allow almost free rotations between contiguous struts. This yields a substantial reduction of the node bending stiffness, which is proportional to the moment of inertia, I, and consequently to the joint thickness, t 4 strut . The bandgap opening frequency is therefore shifted from the value reported in Eq. (1) to that of a simply supported beam, f pinned = 0.44f clamped , entailing a slope-related change in the low ( A 0 ) and high-frequency branches. Again, the phase velocity of the waves diminishes, leading to slower waves propagating in the hinge-like octet lattice and faster waves being instead transmitted by the stiffer, clamped-based frames.
In designing the cell of Fig. 1d, compliant joints were preferred over decoupling the struts with prefabricated hinge-connectors as they preserve the single-phase material while allowing a gradually varying stiffness of the connections. Although beyond the scope of this study, it is important to point out that this choice may cause a high stress concentration in the joints that must be examined prior to fabrication.
Finally, we must emphasize that the underlying physics of the octet, represented in the bandgap, its generating mechanism and tuning parameters, transcends the choice of finite element model. To validate this remark, we have additionally computed the band structure of the standard octet modeled via 3D solid elements. The outcomes of this analysis, which lies beyond the scope of the current article, prove the physical mechanism described in this section for the Euler-Bernoulli beams to be preserved in the solid model. Consequently, we expect the same to hold for the results reported in the upcoming chapters.
Wave propagation in finite meta-lattices. The computational domain, depicted in Fig. 3a, is a 3 cm thick reticulated plate assembled by spatially tessellating the standard octet cell illustrated in Fig. 1a. A so-called meta-area lies within this structure, where the parameters of the octet, introduced in the previous Section and expressed in terms of supplementary masses or joint stiffness, are progressively tuned to yield exotic attenuation or wave guiding phenomena. Therefore, two types of devices are built, as reported in Fig. 3b www.nature.com/scientificreports/ and metalenses. The metabarriers leverage an arrangement of customized octet units, either tailored at the same or diverse frequencies, to generate an inhibited frequency range, where the displacement field vanishes and ultimately shield a target area. Conversely, the metalens relies upon the spatial variation of lattice properties to modify the refractive index in an inhomogeneous bounded region and manipulate wave trajectories, eventually offering the opportunity to focus and harvest vibration energy.
Wave propagation in lattice metabarriers. Two lattice metabarriers, comprising mass-embedded octets, are implemented within a lattice plate made of 120 by 40 standard cells, arranged as per the configuration in Fig. 3a, and separately investigated in this section. Illustrated in Fig. 3b is the 50 × 20 cells constant metabarrier, where the struts of the octets are equipped with midpoint masses, m, equal to 14m ref ( t mass = 3.5 mm), yielding an attenuation zone from 1.8 to 2.8 kHz (gray region in Fig. 4d). The metastructure is assessed under the excitation of a 60ms Ricker wavelet located at a distance of 45 cells to guarantee a plane wavefront. The spectrum of the source, with dominant component along z, is centered on the bandgap, at 2.1 kHz.
The recorded displacement field u z , filtered inside the bandgap at two sequential instants in Fig. 4a and c (see video in the Supplementary Material), unveils the attenuation potential of the lattice metabarrier, based on Fano-like resonances 61 . At t=7.2 ms, the flexural waves scattered by the resonating mass-struts in the metawedge interfere destructively with the unscattered incoming field, impeding the transmission of energy 62 . The resulting energy gap, originated by the hybridization between the asymmetric propagating Lamb wave and the coalescence of the resonant bending mode of the struts and Bragg scattering, induces splitting of the incident wavefront 63 , as depicted in Fig. 4a, c. This mitigation mechanism, already introduced in the previous Section, is further corroborated by the spectrum of two waveforms, r 1 and r 2 in Fig. 4a, b, located on the masses of a cell belonging to the first row of the metabarrier. The signal of the leftmost mass-strut (blue line in Fig. 4b) exhibits a pronounced Lorentzian symmetric peak centered at the resonant frequency f r (1.8 kHz), whilst the second waveform, r 2 , reveals the typical asymmetric profile, where a sharp enhancement of the spectrum prior to f r is followed by a transmission drop (antiresonance) 61 . This Fano-type interference, limited to a narrow bandwidth if only the first line of cells in www.nature.com/scientificreports/ the metabarrier is considered, intensifies when leveraging an array of octets tuned all at the same resonant frequency. By inspecting the spectral displacement of a point located at the core of the metabarrier, after 25 cells ( r 4 ), this effect becomes evident. The resonances of the mass-equipped struts interact constructively yielding a wider attenuation zone that extends from f r to their antiresonance 64 , as reported by the gray region and red signal in Fig. 4d. Furthermore, by comparing this spectrum with the displacement field in a receiver on the edge of the metastructure ( r 3 ), black line in Fig. 4d, we observe that the majority of the energy with frequency content within the forbidden range is converted into reflections 63 . The wavefront impacting on the metabarrier at 7.2 ms, propagates backward in the plate at t=8.8 ms (Fig. 4c) after interacting with the resonating members. Although effectively hampering wave transmission within the bandgap, the constant metawedge exhibits a fixed frequency range of application stemming from the resonant nature of the attenuation zone that, whilst broadband and easily tuneable, entails a spatially localized phenomenon. A mitigation device, equipped with smoothly increasing masses, is designed in order to overcome this drawback, as these were proven to affect the frequency content of the prohibited range (see Fig. 2c). The graded metabarrier, with masses m 1 to m 15 ranging from 3m ref to 42m ref , respectively with t mass from 2 to 5 mm, is depicted in the right inset of Fig. 3b and tested with an input Ricker source centered at 3 kHz. Following the inverse relationship between added mass and resonant frequency, the heaviest masses ( m 15 ) govern the opening frequency of the total bandgap, while the upper bound depends upon the choice of m 1 (see Fig. 5b). The evolution of the bandgap is also shown in Fig. 5b and reveals how the frequencies within the bandgap decrease spatially as masses increase from grade to grade; waves of a certain frequency propagate through the medium until they reach the spatial location where they exist within a bandgap, at this point they are subsequently trapped and back reflected thus, spatially distillating the spectrum following the Rainbow trapping effect 28,31,65 . As a result, our medium possesses a wider stop band, from 1.14 kHz to approximately 5.1 kHz (see gray region in Fig. 5c). Figure 5a illustrates the displacement field filtered over two distinct frequency ranges (see videos in the Supplementary Material). Between 1.8 and 2.5 kHz, the wave impinging the metabarrier is almost immediately trapped and scattered by the cells with masses equal to 11m ref , owing to the matching with their frequency gap (1.98-3.09 kHz). Similarly, the lower frequency field, corresponding to a pass-band filter between 1.2 and 1.5 kHz, travels undisturbed deeper in the metastructure until it eventually interacts with the cells of mass 31m ref (bottom inset), the exact location these frequencies occur within the varying bandgap.
As validation, the spectral displacement of a point at the right end of the metawedge is compared to that of a receiver at an equivalent position in a bare lattice plate in Fig. 5c. The amplitude reduction, stemming from the presence of the metawedge and approximately equal to 87%, is visibly appreciable, as well as the widening of the bandgap compared to that in Fig. 4d. Notice that, a zero-transmission range exists in both the spectra between 3.8 and 5.1 kHz, which can be attributed to the dispersive behavior of the standard octet cells. While being progressively slowed down within the graded barrier, the impinging energy is trapped by the bending resonant mode Octet-based gradient index lenses. Drawing from the dynamic properties of the octet lattice, two of the most popular GRIN lenses are engineered from this architecture to shape ray trajectories and focus wave energy. The design of these devices relies on the dependency of the phase velocity, extracted from the dispersion relations, upon the frequency and the tuning parameters. By combining this relationship with the refractive index profile, the design variables of the customized octet cells are derived at different locations in the lens. Specifically, the Luneburg and Maxwell lenses are created by gradually varying the supplementary masses embedded in the octets, or the stiffness of their joints, designed to produce the required refractive index profile within a bounded region of radius R, which corresponds to 12 cell layers, each featured by a distinct parameter value (Fig. 3c). The required theoretical refractive index profile is defined as a function of the radial coordinate, r, centered in the middle of the lens, as follows: The effective wavespeeds in the metalenses, v eff , and in the surrounding lattice plate comprising only standard octets, v 0 , are retrieved from the band structure of the cells, computed via Bloch analysis, with the relationship v eff = ω/k. In order to eliminate their frequency dependence, originated by the dispersive nature of the octet, these quantities are evaluated at fixed components, f fix . The frequencies are selected sufficiently far from the stop band, so as to prevent the waves from being localized and amplified by the strut resonances. At 1 kHz and 5.2 kHz, respectively for the Luneburg ( f fix,L ) and Maxwell ( f fix,M ) lens, the phase velocity is calculated for increasing added masses or node thickness and plotted in the top inset of Fig. 6b and c. While both wavespeeds feature a decreasing trend for heavier masses, a significant drop can be observed in the high-frequency branch ( f fix,M ), posing as a suitable candidate for the design of a Maxwell lens, which requires an abrupt variation of the refractive index. The opposite behavior is instead noted when varying the thickness of the nodes yielding a speed-up of the waves that propagate from hinge-based ( t node = 0.3 t strut ) to stiffer lattices ( t node = t strut ). These findings serve as a bridge between the theoretical profiles of GRIN lenses and their actual implementation using reticulated structures, as they relate the octet properties to the refractive index, n, by means of the effective velocity.
After computing the refractive index profiles of the Luneburg and Maxwell lenses with the formula from Eq. (2), as shown in Fig. 6a for increasing radial coordinate, the following expression for the refractive index, www.nature.com/scientificreports/ is used to retrieve the phase velocity distributions illustrated in the bottom inset of Fig. 6a. Finally, an inverse design approach, relying on the relations in top Fig. 6b and c, is adopted to identify the tuning parameters of the customized octets at different r, based on the input effective velocity. The designed point masses and node thicknesses are outlined in Fig. 6b and c for Luneburg and Maxwell lenses. Specifically, embedded masses from 3m ref to 17m ref or node thicknesses ranging from 0.2 to 1 t strut are determined according to the Luneburg lens refractive index (from 1 to 1.41). The fish-eye lens, instead, originates from a stronger variation of the index profile (from 1 to 2), which requires masses from 1 to 34m ref or t node to vary between 0.3 and 1 t strut . These designs, engineered via Bloch theory, are validated by means of time-domain numerical simulations of waves traveling in the finite plate of Fig. 3a, with the Real-ESSI Simulator (see "Methods" section). Both GRIN metalenses are obtained by endowing the octet cells with the masses estimated in Fig. 6b and are excited by sinusoidal sources tuned at the design frequencies f fix,L and f fix,M . As input signal for the Luneburg one, we use a plane wave of 40ms and 1kHz. In the Maxwell case we instead use a point-like source of frequency 5.2 kHz and duration 10ms. The wave field travels through the lenses with minimal losses and reflections, thanks to the gradual variation of properties from the plain octet lattice at r = R to the center of the device that minimizes the impedance mismatch (see videos in the Supplementary Material). In Fig. 7a, a plane wave impinging on the left side of the Luneburg lens at 6.8 ms is progressively bent by the smooth refractive index transition and its trajectory is steered toward a diametrically opposite focal point placed on its surface (t=10.4ms). An array of receivers, positioned along this edge, outlines the significant energy localization generated by the lens. A sharp spike emerges in the out-of-plane displacement u z and its peak value is almost twice as large as its counterpart on the bare plate.
Similarly, two distinct time instants in Fig. 7b validate the Maxwell lens design, where the wave originating from a point source located at the leftmost side of the lens is guided toward an antipodal focus thus, creating the typical fish-eye effect. Notice that, a non-circular wavefront is produced by the fields of the source and focal point www.nature.com/scientificreports/ in the surrounding plate, which instead propagate along preferential directions (approximately at 45 • ) due to the anisotropic character of the standard octet lattice. This behavior is further corroborated by the isofrequency contours of the standard octet lattice at f fix,M , reported in the bottom inset of Fig. 7b. A 70 × 70 plate, complemented with absorbing boundary conditions and a broadband point source placed at its center, is leveraged for the numerical derivation of the isocurves. The resulting wave field is spatio-temporally Fourier transformed and the snapshot corresponding to 5.2 kHz is then extracted. The circular isofrequency contours, typical of an isotropic medium, distort here to a rhombic shape emphasizing the anisotropic nature of the lattice. The wave vectors impinging any of these contours are transmitted along the diagonals of the first Brillouin zone, causing the wave field in the real space to assume the cross-like form that emerges in Fig. 7b. A similar anisotropic trait can be appreciated at lower frequencies ( f fix,L =1 kHz) upon inspection of the corresponding isofrequecy contour (bottom inset of Fig. 7a). In addition, the anisotropy explains the near-zero vertical displacement in the receivers of the bare plate, blue curve in the inset of Fig. 7b, which juxtaposes with the peak of the highly focused signal from the lens (black line). Equivalent results are obtained when the lens design is based on the variation of the node thickness. The only difference lies in the shape of the metadevice, squared instead of circular, to avoid discontinuities in the boundary conditions of the struts.

Discussion
We have demonstrated the potential of octet-based lattice plates for controlling the trajectories of asymmetric Lamb-like waves, be it via attenuation or wave guiding. Starting from the cellular architecture of this structure, we prove the standard octet unit is capable of supporting a Fano-type bandgap (3.85-5.08 kHz) induced by the coalescence of local bending resonances and Bragg frequency, opening the avenue for an easy-to-tune, singlephase material in elastic applications. The design of this reticulated structure is simplified by viewing the struts as Euler-Bernoulli beams clamped at both their ends and leveraging the analytical expression of the fundamental www.nature.com/scientificreports/ flexural mode as a prediction paradigm for the bandgap opening frequency; furthermore, since the octet contains struts of length and direction identical to that of the primitive lattice vectors, it creates a medium in which the clamped-clamped fundamental resonance exists as a standing wave between nodes, hence it is a medium in which local bending resonances and Bragg frequencies coincide, resulting in destructive interferences and a broadband resonant-based bandgap. The two tuning parameters extrapolated from this analogy, midpoint added masses and joint thickness, are proven to control the frequency content of the attenuation zone as well as the slope of the low and high-frequency bands exhibiting out-of-plane polarization. The opening frequency of the bandgap is inversely correlated with the masses and directly correlated with the node thickness; heavier masses drive the forbidden range to lower frequencies and reduce the phase velocity of waves in the lattice, conversely, thicker nodes shift the content of the attenuation zone to higher values and privilege the transmission of faster waves.
Drawing from these intrinsic dynamic behaviors, four devices for controlling ray trajectories are successfully engineered. The devices operate over a range of scales, the metabarriers leverage the bandgap and hence operate within a Bragg-scattering regime; on the other hand, the Luneburg lens is based on spatial variations of the lowfrequency A 0 mode and hence the underlying octet cells behave in a subwavelength manner.
The constant metabarrier confirms the existence of a hybridization bandgap (1.8-2.8 kHz) stemming from the interaction between the propagating flexural field and octet units equipped with auxiliary masses. The wave impinging the resonant members is backscattered and prevented from being transmitted within the metastructure. Whilst the bandwidth of the constant metabarrier is broad, the bandgap and frequency range is spatially localized; we introduce a graded design capable of not only tripling the mitigation width (1.14-5.1 kHz), but also spatially filtering the spectrum via a rainbow trapping mechanism, amplifying a large spectrum of frequencies within its resonating struts. The generating mechanism of the stop band, initially exploited by the lattice metabarriers to match the frequency content of incoming fields and inhibit them, serves to also tune the phase velocity within a bounded region and enable wave steering. Two GRIN lenses, the Luneburg and Maxwell fisheye, are successfully implemented by gradually adjusting the properties of the customized octets according to the theoretical index spatial distribution, thus paving the way for novel applications of this face-centered cubic topology in wave focusing.
Although the dynamic assessment of three-dimensional frame lattices is still at a primordial stage compared to the existing knowledge on their static mechanical properties, the outcomes of our work have shed light on the hidden potential of the octet and tuning parameters. The ease of scaling to different wavelengths that characterizes this structure renders it even more attractive for applications in a variety of domains, including vibration isolation or conversion of focused energy. Indeed, the vibration energy trapped and focused, respectively by the struts and metalenses, could be harvested and converted into other forms of energy, e.g., electrical, by attaching piezoeletric components to the lattice. Future steps of our work will address the development of a more sophisticated, yet computationally expensive, numerical model based on solid elements and its comparison with the simplified Euler-Bernoulli beam model proposed in this study. This would then naturally lead to the experimental validation of the octet-based metadevices upon fabrication using 3D printing. Finally, the complicated matter of stress concentration due to dynamic loading, as well as the extension to different lattice unit cells based both on struts and sheets, such as foams, could be tackled.

Methods
Bloch analysis. The dispersion relation, reported in Fig. 2a (red dotted line), is solved using the finite element method, performed by the software COMSOL Multiphysics 5.6, by sweeping through the required wavevector k determining the permitted frequencies ω = ω(k) as an eigenvalue and the displacement field as an eigenvector. Our choice of fundamental cell (Fig. 1a) is a conventional cubic-like unit cell over a primitive one 19 , which allows us to align the base vectors with the propagation direction of interest. This makes the comparison between the band structure stemming from the FFT analysis performed in the time-transient finite lattice section (Fig. 1b) and the eigenmodal analysis straightforward as we only know the displacement fields along the beams and not everywhere in the cell. With this choice of unit cell the ends of the struts either correspond to a node point or lie on the edges of the cell, we apply either continuity conditions or quasi-periodic Bloch conditions at the strut ends as appropriate; hence the displacement field takes the following well-known form: where k = κ x e x + κ y e y + 0e z denotes the in-plane Bloch wavevector, ω the frequency assuming harmonic timedependence and U(x) is a periodic function satisfying: for any integers n and m. We restrict our attention to waves propagating along the x-axis, hence consider only k along the Ŵ-X portion of the irreducible Brillouin zone (0 -1 2l 1/m) (see Supplementary Material for the full band structure). The struts are modeled as EB beams with 6 degrees of freedom in each node (inset in Fig. 1a), reducing the high computational cost of the simulations that would arise when using solid elements. The EB theory was selected over the Timoshenko theory upon fulfilling the condition EI/0.88AGl 2 strut << 1 , with A and I designating the area and moment of inertia, respectively, and G the shear modulus, as the slenderness of the beams makes shear deformations negligible.
Time-transient simulations. The black portion in the band structure (Fig. 2a)    www.nature.com/scientificreports/ Real-ESSI Simulator 66 . The system is excited by a modulated plane sweep source ranging from 0.1 to 9 kHz with components in all three directions of the space. The asymmetric Lamb waves propagating in the plate are recorded by an array of points lying perpendicularly to the source along the metabarrier center (blue dotted line in Fig. 1b) and the resulting wave field u z is spatio-temporally Fourier transformed to derive the correspondent band structure. The metadevices (Fig. 3) are tested via 3D time-transient analyses conducted with the parallel version of the Real-ESSI Simulator 66 in order to substantially reduce the computational time. The Newmark integration method is employed with γ equal to 0.505. The parameter γ controls the numerical or algorithmic damping, that is inserted into the model to damp out any non-realistic high frequencies stemming from model discretization ( γ equal to 0.5 would correspond to zero numerical damping). It should be noted, here, that this type of damping cannot be used to represent any physical damping mechanism (e.g. viscous damping, hysteretic damping), while improper use could lead to unrealistic attenuation of the real response of the model 67 (i.e. values greater than 0.5 should be used with caution). External software is employed for both pre-and post-processing; specifically, the geometry and mesh (Fig. 3a) are generated with the 3D finite element generator Gmsh and a translator, called gmESSI, is invoked to convert the mesh file from Gmsh to a format readable by Real-ESSI. The mesh of the struts in the lattice consists of three-dimensional EB beam elements. The customized cells in Fig. 1c and d can be easily implemented either by adding point masses to the center of the struts or locally modifying their cross-section in the node area, condensing the degree of complexity to the design stage with Gmsh. Although the Real-ESSI Simulator output results can be visualized with the open-source application ParaView after loading the PVES-SIReader plugin 66 , for this study, an in-house code, written in Python, was employed to post-process and later visualize the wave field in ParaView.
The in-plane size of the full frame plate is specified by the number of cells, N cells,x and N cells,y in Fig. 3a, multiplied by the cell width, l. This domain is tested under vertically polarised excitations with unit amplitude that entail flexural Lamb waves propagating along x. Depending on the assessed metastructure, the input may be a line, for metabarriers and Luneburg lens, or punctual source, for Maxwell lens, and its frequency content span from that of a pulse-like Ricker excitation to a monochromatic source (inset of Fig. 3a). Absorbing layers using quadratically increasing Rayleigh damping (ALID) 60 , defined along the five outer cell arrays of the plate, are implemented to mimic an unbounded domain thus, minimizing undesirable reflections from the boundaries. A total of twenty absorbing layers on every edge is considered, each of thickness l/4, and a mass proportional damping with maximum coefficient c max equal to 2 · 10 4 is assumed. The inner portion of the plate is, instead, left damping-free to avoid compromising the resonance of the struts and given that the selected material, versatile plastic, is characterized by a low quality factor (Q=50).

Data availability
The datasets generated during the current study are available from the corresponding author upon reasonable request.