{"id":914171,"date":"2026-11-12T00:26:14","date_gmt":"2026-11-12T00:26:14","guid":{"rendered":"https:\/\/www.newsbeep.com\/us\/914171\/"},"modified":"2026-11-12T00:26:14","modified_gmt":"2026-11-12T00:26:14","slug":"scale-rich-network-based-metamaterials-nature","status":"publish","type":"post","link":"https:\/\/www.newsbeep.com\/us\/914171\/","title":{"rendered":"Scale-rich network-based metamaterials | Nature"},"content":{"rendered":"

We begin by determining the mechanical characteristics of SR microstructures, comparing them against three single-scale microstructures frequently used in engineering applications: the square lattice, the hexagonal lattice and the irregular Voronoi tessellation. By generating finite domains of each microstructure with a sufficiently large number of characteristic features, we confirm that the P(\u03bb), P(\u2113) and P(k) for square and hexagonal lattices result in delta functions, whereas the Voronoi tessellation is characterized by narrow distributions (Fig. 1i\u2013k<\/a>).<\/p>\n

To extend the comparison beyond theoretical differences between architectures, we experimentally and numerically determined the mechanical properties of these metamaterial microstructures, using vat photopolymerization 3D-printing (Methods<\/a>) to fabricate centimetre-scale specimens out of elastic\u2013plastic polymer resin (Fig. 1d<\/a>, left), with a feature resolution of approximately 100\u2009\u03bcm. We focused on SR systems with T\u00a0=\u00a0500\u2014a point where the homogenized mechanical properties show convergence and boundary effects are minimized (Supplementary Information, section 2.3<\/a>). Given this feature resolution and the 500-ligament constraint, the relative densities achievable with vat photopolymerization 3D-printing ranged between \\(20 \\% \\le \\bar{\\rho }\\le 75 \\% \\) (Supplementary Information, section 2.5<\/a>); we can also demonstrate lower fill fractions within the tunable-density phase (\\(\\bar{\\rho } < 20 \\% \\); Fig. 1c<\/a>) through the use of two-photon lithography (Fig. 1d<\/a>, right) with submicrometre minimum feature sizes. Next, we characterized the linear and nonlinear mechanical responses in the quasi-static and dynamic regimes of these SR microstructures.<\/p>\n

Highly tunable linear-elastic responses<\/p>\n

To determine the linear-elastic response of SR metamaterials, we performed computational homogenization using the finite-element method. Imposing affine boundary conditions, we extracted a directional modulus Ed in each direction ed and calculated a non-dimensional modulus \\(\\bar{E}\\equiv {E}_{{\\rm{d}}}\/{E}_{{\\rm{s}}}\\), where Es is Young\u2019s modulus of the constituent material (Methods<\/a>).<\/p>\n

At a relative density of \\(\\bar{\\rho }=40 \\% \\), we performed this analysis on the three single-scale systems (square, hexagonal and Voronoi) and SR geometries, demonstrated by the representative realizations in Fig. 2a<\/a>, via polar representations of \\(\\overline{E}\\) as a function of the probing direction. While the hexagonal and Voronoi structures show a near-constant \\(\\overline{E}\\) in all orientations (that is, an isotropic stiffness response, as expected46<\/a>; Fig. 2a<\/a>(i<\/a>),(ii<\/a>)), the square lattice is highly anisotropic, displaying high stiffness parallel to its struts and minimal stiffness in any other direction (Fig. 2a<\/a>(iii<\/a>)). In contrast to these near extremes, we find the behaviour of the SR structures to be modulated by \u03b1: for \u03b1\u00a0\u2248\u00a00, the stiffness is nearly isotropic (Fig. 2a<\/a>(iv<\/a>)); and as \u03b1\u00a0\u2192\u00a01, a principal orientation of maximum stiffness emerges (Fig. 2a<\/a>(v<\/a>),(vi<\/a>)). Altogether, the full set of SR architectures exhibited highly tunable (an)isotropy exceeding that of the three conventional geometries combined.<\/p>\n

Fig. 2: Linear-elastic behaviour of the SR microstructure.\"Fig.<\/p>\n

a, Directional stiffness plots of the reference hexagonal (i), Voronoi (ii) and square (iii) lattices from finite-element simulation, compared with SR samples with \u03b1\u00a0=\u00a00 (iv), \u03b1\u2009=\u20090.3 (v) and \u03b1\u2009=\u20090.9 (vi) at the same relative density. The right panels show the experimentally measured directional stiffness for SR materials with \u03b1\u00a0=\u00a00 (vii), \u00a0\u03b1\u2009=\u20090.1 (viii) and \u03b1\u2009=\u20090.3 (ix), validating the theoretically predicted stiffness anisotropy. Scale bars, 10\u2009mm. b, Anisotropy\u2013density plot, illustrating the broad range of anisotropy achievable with the SR framework compared with the narrow range of values attained by the reference samples. The black filled circles (vii)\u2013(ix) show the experimentally measured anisotropy. c, As the decay parameter \u03b1 increases, anisotropy arises because the maximum specific directional stiffness increases, whereas the minimum specific directional stiffness is constant. The solid lines represent the mean across all samples, and the shaded regions correspond to a \u00b11\u2009standard deviation from the mean value. d, Phase diagram in (\u03b1,\u00a0\u03bb0) space of the anisotropy, which functions as a design map for a tunable anisotropic response across relative densities \\(0.2\\le \\bar{\\rho }\\le 0.75\\). e, Temporal evolution of ligament strain energy during the SR generation process. f, Snapshots at fixed time t are shown for two example ligaments in an SR network, corresponding to the markers in e. g, Universal scaling of the ligament energy Ui as a function of the degree, ki, thickness, \u03bbi and angle, \u03b8i. For each \u03b1, data are averaged over n = 5 independent realizations, each with T = 5,000 steps.<\/p>\n

To quantify the mechanical implications of this observed range in anisotropy across the full space of SR architectures, we plot the elastic anisotropy ratio \\({\\eta }_{{\\rm{E}}}\\equiv {E}_{\\min }\/{E}_{\\max }\\in (0,1]\\), where \\({E}_{\\text{min}}\\equiv \\text{min}(\\bar{E})\\) and \\({E}_{\\text{max}}\\equiv \\text{max}(\\bar{E})\\), as a function of relative density for the reference geometries and 2,600 distinct SR microstructures with different \u03b1 and \u03bb0 (Fig. 2b<\/a>). We confirm that the hexagonal lattice (0.87\u2009\u2264\u2009\u03b7E\u2009\u2264\u20090.94) and the Voronoi tessellation (0.88\u2009\u2264\u2009\u03b7E\u2009\u2264\u20090.97) remain nearly isotropic across a wide range of relative densities. For the square lattice, the anisotropy decreases from \u03b7E\u00a0=\u00a00.03 at \\(\\bar{\\rho }=20 \\% \\) to \u03b7E\u00a0=\u00a00.55 at \\(\\bar{\\rho }=70 \\% \\). By contrast, we find that the SR materials cover 77% of the available area on the density\u2013anisotropy plot (Fig. 2b<\/a>), reaching values not accessible to the reference systems. Indeed, by changing \u03b1 and \u03bb0, we can attain highly anisotropic systems (\u03b7E\u2009\u2264\u20090.2) across a wide range of densities (\\(17 \\% \\le \\bar{\\rho }\\le 64 \\% \\)) as well as nearly isotropic systems (\u03b7E\u2009\u2265\u20090.85) in the range \\(33 \\% \\le \\bar{\\rho }\\le 72 \\% \\). Simultaneously, we can generate systems of widely varying anisotropies at any fixed relative density: at \\(\\bar{\\rho }=20 \\% \\), we span 0.04\u2009\u2264\u2009\u03b7E\u2009\u2264\u20090.78; at \\(\\bar{\\rho }=40 \\% \\), we span 0.12\u2009\u2264\u2009\u03b7E\u2009\u2264\u20090.88; and at \\(\\bar{\\rho }=70 \\% \\), we span 0.26\u2009\u2264\u2009\u03b7E\u2009\u2264\u20090.93. Directed growth approaches (Supplementary Information, section 2.15<\/a>) further extend the achievable coverage to more than 96% of the density\u2013anisotropy space.<\/p>\n

To elucidate the origin of this tunable anisotropy, we plot the mass-normalized (that is, specific) minimum and maximum stiffness values (\\({E}_{\\min }\/\\bar{\\rho },{E}_{\\max }\/\\bar{\\rho }\\)) as a function of the decay parameter \u03b1 across the entire sample set (Fig. 2c<\/a>). Although the minimum specific stiffness is largely independent of \u03b1, increasing \u03b1 while keeping density constant enables an increase in the maximum achievable specific stiffness, leading to robust tunability in the architecture space. This tunable anisotropy is best captured by the (\u03b1,\u00a0\u03bb0) configuration space (Fig. 2d<\/a>), which shows that anisotropy increases with \u03b1 and decreases with \u03bb0. To validate our finite-element model predictions and confirm the isotropic and anisotropic responses, we performed uniaxial compression experiments across various SR geometries and orientations, confirming that at constant density increasing \u03b1 directly leads to greater anisotropy (Fig. 2a<\/a>(vii<\/a>)\u2013(ix<\/a>), Fig. 2b<\/a>, filled black circles, and Supplementary Information, section 2.7<\/a>).<\/p>\n

A scalar metric that accounts for the alignment of material along a certain loading direction indicates that the tunable directional stiffness of SR materials is the result of the inherent heterogeneity in ligament orientations and sizes. This is best captured by an orientation tensor\u00a0M<\/p>\n

$$M\\equiv \\mathop{\\sum }\\limits_{t=1}^{T}{A}_{t}{{\\bf{d}}}_{t}\\otimes {{\\bf{d}}}_{t},$$<\/p>\n

\n (1)\n <\/p>\n

where At is the area of ligament t and dt is a unit vector along its centre line. M allows us to compute a geometric anisotropy \u03b7G\u00a0\u2261\u00a0\u03bc1\/\u03bc2\u00a0\u2208\u00a0(0,\u00a01], where \u03bc1\u2009\u2264\u2009\u03bc2 are the eigenvalues of M. We find that \u03b7G agrees with the elastic anisotropy \u03b7E (Supplementary Information, section 2.6<\/a>), confirming that the anisotropy of the SR system is a direct consequence of its length-scale heterogeneity.<\/p>\n

Tunable energy-density scaling<\/p>\n

To understand the role of multiscale features on the measured properties, we follow two ligaments, one introduced at step t\u00a0=\u00a020 and one at t\u00a0=\u00a0100, through the full generation process. As Fig. 2e<\/a> shows (see Fig. 2f<\/a> for the network snapshots), these ligament-specific trajectories span more than eight orders of magnitude and are irregular and rather unpredictable. Yet, as we show in\u00a0Supplementary Information, section\u00a03<\/a>, the strain energy distribution across all ligaments scales with the introduction time as i\u2212(3\/4+\u03b1\/2), with the degree as \\({k}_{i}^{3\/2+\\alpha }\\), with the length as \\({{\\ell }}_{i}^{3\/2+\\alpha }\\), and with the thickness as \\({\\lambda }_{i}^{3\/(4\\alpha )+1\/2}\\) (Supplementary Fig. 46<\/a>). Combining these scaling laws, we find that the elastic strain energy follows a single per-ligament law (Fig. 2g<\/a> and Supplementary\u00a0Fig. 47<\/a>)<\/p>\n

$${U}_{i}(t)\\approx {k}_{i}^{3\/2}\\,{\\lambda }_{i}^{1\/2}\\,f({\\theta }_{i}),$$<\/p>\n

\n (2)\n <\/p>\n

where f(\u03b8i) is an orientation factor, with \u03b8i being\u00a0the ligament orientation. This law is, remarkably, independent of \u03b1, confirming that all scales contribute to the total energy. Plotting the per-ligament energy against the combination in equation (2<\/a>) collapses all ligaments onto a single line of slope approximately 1, independent of \u03b1 (Fig. 2g<\/a>), confirming that all ligaments follow the scaling law equation (2<\/a>).<\/p>\n

Integrating the per-ligament law, equation (2<\/a>), over the SR distributions P(k)\u00a0\u2248\u00a0k\u22123 and P(\u03bb)\u00a0\u2248\u00a0\u03bb\u2212(1+1\/\u03b1) allows us to analytically derive the scaling relation between the total elastic energy and relative density, \\({U}_{\\mathrm{tot}}\\approx {\\bar{\\rho }}^{m(\\alpha )}\\) (Supplementary Fig. 49c<\/a>), finding that the scaling exponent m depends on the SR exponent \u03b1 as <\/p>\n

$$m(\\alpha )=\\frac{1-\\alpha \/2+p(\\alpha )}{1\/2-\\alpha },\\qquad \\alpha < 1\/2,$$<\/p>\n

\n (3)\n <\/p>\n

where \\(p(\\alpha )\\) is a fitted exponent describing the residual time dependence of the per-ligament energy. This expression matches the directly fitted slopes within \u00b10.2 (more details in\u00a0Supplementary Information, section 3<\/a>). The SR construction therefore interpolates continuously across the scaling map, with \u03b1 as the only control parameter. These results show that the multiscale mechanics of SR networks can be reduced to a universal ligament-level energy law, which turns the growth exponent \u03b1 into a macroscopic control parameter that continuously tunes the energy-density scaling response, with lower m exponents indicating more efficient mass utilization.<\/p>\n

Delocalized nonlinear deformations<\/p>\n

Next, we characterized the nonlinear response of SR metamaterials by conducting large-deformation uniaxial compression experiments (Methods<\/a>) on the 3 reference geometries and SR systems of \u03b1\u00a0=\u00a00, chosen to minimize the impact of elastic anisotropy and emphasize the effects of the heterogeneity of ligament lengths (Fig. 3a<\/a>; 3 specimens each, all at nominal \\(\\bar{\\rho }=42 \\% \\); experimental results for SR systems with \u03b1\u00a0>\u00a00 are discussed in\u00a0Supplementary Information, section 2.11<\/a>.). Experiments on the square, hexagonal and Voronoi tessellations showed qualitative differences between each architecture (for example, in the densification strain); repeated samples of each architecture showed minimal variation (Supplementary Information, section 2.8<\/a>). By contrast, all realizations of SR architectures (Fig. 3a<\/a>, orange dotted lines) showed a monotonically increasing stress\u2013strain response, coexisting with a wide envelope of responses, elucidating the tunability enabled by the SR framework and breaking the structure\u2013response constraints observed in the single-scale reference geometries.<\/p>\n

Fig. 3: Nonlinear behaviour of the SR microstructure.\"Fig.<\/p>\n

a, Experimental stress\u2013strain curves of SR and reference samples, for \\({\\bar{\\rho }}_{\\mathrm{nom}}=42 \\% \\), under uniaxial compression. For the reference samples, the mean behaviour of three replicate specimens is plotted. For the SR geometries, the solid line represents the mean of all samples with \u03b1\u00a0=\u00a00, whereas the dashed curves are individual samples. The response of SR samples with \u03b1\u00a0>\u00a00 is given in\u00a0Supplementary Fig. 33<\/a>. b, Experimental toughness\u2013density Ashby plot, showing the energy absorption capacity across all relative densities. c, Evolution of localized deformation with increasing macroscopic strain: still images together with computed areal Jacobians, J\u00a0=\u00a0A\/A0, for each cell. The geometries dominated by a characteristic length scale (square lattice, (i); hexagonal lattice, (ii); Voronoi tessellation, (iii)) all develop regions of compaction where localization is observed; only the SR sample (iv) avoids the formation of compaction bands. We tracked images up to a macroscopic strain of 25%, after which compaction precludes the accurate computation of Jacobian values. Scale bars, 10\u2009mm. d, Evolution of the localization parameter \\({\\mathcal{L}}\\) in each geometry together with snapshots of the deformation evolution as recovered from still images. Here, the mean and standard deviation of three unique SR specimens are plotted.<\/p>\n

From each stress\u2013strain response, we computed the toughness, measured as the integral under the stress\u2013strain curve to an engineering strain of 40%\u2014before the onset of complete densification. These toughness metrics, when plotted against relative density (Fig. 3b<\/a>), showcase the ability of SR architectures to match or surpass the mechanical resilience of the three reference architectures.<\/p>\n

To understand the deformation mechanisms that lead to these benefits, we extracted still frames from the experiments and computed the areal Jacobian Ji\u00a0\u2261\u00a0Ai\/Ai,0 of each ith cell (the relative change in void area; Fig. 3c<\/a>) where Ai,0 and Ai denote the undeformed and deformed cell areas, respectively47<\/a>. The square lattice was prone to buckling of vertically aligned ligaments, resulting in the observed local maxima in its stress\u2013strain response (Fig. 3a<\/a>, blue curve). Buckling events began at a single ligament, followed by a horizontal cascade of buckling events that spanned the whole length of the specimen (Fig. 3c<\/a>(i<\/a>)), leading to a sudden drop in its load capacity. This process repeated in a layer-by-layer fashion until the entire specimen buckled. Similarly, the hexagonal lattice (teal; Fig. 3c<\/a>(ii<\/a>)) also developed wide compaction bands that spanned the width of the lattice. Although the Voronoi tessellation did not succumb to well-defined shear or compaction bands owing to its lacking periodicity, consistent with previous experimental findings48<\/a>, it exhibited localized compaction near the bottom, ultimately forming a domain of compacted cells that spanned the width of the specimen (purple; Fig. 3c<\/a>(iii<\/a>)), whereas unit cells near the top remained largely undeformed.<\/p>\n

In contrast, SR architectures lack compaction or shear bands, which can be attributed to the local coexistence of ligaments of different lengths and thicknesses; hence, the failure of the weakest ligament in one domain will not trigger the propagation of contiguous failure pathways (Fig. 3c<\/a>(iv<\/a>)). To understand the origin of this diffuse damage, we computationally evaluated the susceptibility of each unsupported ligament section to buckling under a macroscopic deformation (Supplementary Information, section 2.13<\/a>). We find that the network consists of ligaments characterized by a broad distribution of critical buckling loads, and hence a local buckling event does not trigger a global cascade of instability, explaining the stable macroscopic response. Moreover, selectively removing a range of length scales degrades the load-bearing capacity of the network (Supplementary Information, section 2.14<\/a>), suggesting that length-scale heterogeneity underlies the robustness of the system.<\/p>\n

To further understand this delocalization, we introduce a localization parameter <\/p>\n

$${\\mathcal{L}}\\equiv \\frac{{\\sum }_{i=1}^{N}{({\\mu }_{i}-\\bar{\\mu })}^{2}}{{\\sum }_{i=1}^{N}{({J}_{i}-\\bar{J})}^{2}},$$<\/p>\n

where \\({\\mu }_{i}\\equiv \\frac{1}{M}{\\sum }_{j\\in {{\\mathcal{B}}}_{r}(i)}{J}_{j}\\) is the local average Jacobian of cell i, defined over a neighbourhood \\({{\\mathcal{B}}}_{r}(i)\\) of i containing M cells (Supplementary Information, section 2.10<\/a>), and N is the total number of cells. As a result, \\({\\mathcal{L}}\\in [0,1]\\), where \\({\\mathcal{L}}=0\\) implies that the deformation in every neighbourhood is identical\u2014the ideal case\u2014and \\({\\mathcal{L}}=1\\) indicates that the deformation variation between neighbourhoods is comparable to the deformation variation among all individual cells, corresponding to extreme localization.<\/p>\n

For all samples, \\({\\mathcal{L}}\\) initially increases with increasing strain (Fig. 3d<\/a>); the square and hexagonal lattices exhibit rapidly increasing and high asymptotic \\({\\mathcal{L}}\\) owing to the presence of compaction bands, whereas the disordered Voronoi lattice has lower \\({\\mathcal{L}}\\). Consistently, the SR samples show the lowest value of \\({\\mathcal{L}}\\)\u2014with the localization within \\({\\mathcal{L}}\\le 0.2\\) through 25% macroscopic strain, which is 37% lower than the Voronoi geometry\u2014confirming the material\u2019s ability to uniformly distribute damage and offer a macroscopically stable load-bearing capacity. We can further reduce \\({\\mathcal{L}}\\) by selectively nucleating the ligaments in an SR network (Supplementary Information, section 2.12<\/a>) using a \u2018checkerboarding\u2019 approach, allowing us to constrain \\({\\mathcal{L}}\\le 0.08\\) through 25% strain, representing a 74% reduction over the Voronoi value.<\/p>\n

Efficient dynamic impact absorption<\/p>\n

Looking beyond the quasi-static regime, an ideal impact-absorbing material maximizes the energy dissipated during the impact event while minimizing the peak force49<\/a>. We performed impact experiments on SR and reference architectures at relative density \\(\\bar{\\rho }=40 \\% \\), using a rigid impactor released with controlled initial energy corresponding to an instantaneous strain rate of \\(\\mathop{\\varepsilon }\\limits^{\\cdot }=100\\,{{\\rm{s}}}^{-1}\\) at impact.<\/p>\n

Using a high-speed camera, we analysed frames during the impact event (Fig. 4a<\/a>), finding that similar to the quasi-static experiments, the square lattice fails in a layer-by-layer fashion with fractured ligaments localized within sequentially forming horizontal crushing bands. The square lattice showed a sharp initial peak force that dropped to nearly zero followed by subsequent peaks corresponding to local compactions (Fig. 4b<\/a>, blue curve). The hexagonal lattice also exhibited extreme localization owing to the fracture of diagonally oriented planes in as little as 2\u2009ms post-impact. Although the Voronoi tessellation diffused damage more effectively, it also showed clear localized crush bands, leading to a high maximum force (Fig. 4b<\/a>, purple curve). The SR geometry, by contrast, lacks not only well-defined crushing bands but also diffused damage through the structure, leading to the most stable force\u2013time behaviour (Fig. 4b<\/a>, orange curve), with low initial peak force and a lower compaction force than the Voronoi geometry.<\/p>\n

Fig. 4: Static and dynamic performance of SR and reference geometries.\"Fig.<\/p>\n

a, The SR system efficiently distributes damage under dynamic impact loading, captured by high-speed photography. Scale bars, 10\u2009mm. b, Force\u2013time data reveal a flat, plateau-like force\u2013time response in the SR materials, compared with the reference architectures (square, hexagonal, triangle, sea-sponge53<\/a>, and Kagome lattices; Gilbert55<\/a>, Voronoi, and E-Voronoi48<\/a> tessellations; and isotropic (stretch\u00a0h\u00a0=\u00a01) k-nearest foams54<\/a>). c, This corresponds to simultaneous minimization of peak force and maximization of specific energy dissipation during impact. d, Anisotropy and localization under quasi-static loading. All systems have relative density \\(\\bar{\\rho }=40 \\% \\). Localization corresponds to the experimentally observed value at 25% macroscopic strain and anisotropy is computed from finite-element simulations.<\/p>\n

We extracted the overall peak force, Fmax, and the change in kinetic energy of the impactor before and after the event, representing the energy dissipated by the architecture, allowing us to measure the specific energy dissipation as <\/p>\n

$${W}_{\\mathrm{sp}}\\equiv \\frac{{m}_{\\mathrm{impactor}}({v}_{{\\rm{i}}}^{2}-{v}_{{\\rm{f}}}^{2})}{2{m}_{\\mathrm{specimen}}},$$<\/p>\n

where\u00a0mimpactor\u00a0is the mass of the impactor, mspecimen is the mass of the specimen, and vi and vf represent the impactor velocities before and after the impact event, respectively. In Fig. 4c<\/a>, we plot Wsp against Fmax, finding that the SR system occupies the optimal space with maximum specific energy dissipation Wsp and minimum peak force Fmax among all reference geometries considered in this work, indicating that the advantageous quasi-static properties of SR materials offer clear benefits under dynamic loading as well.<\/p>\n

Steering elastic waves using SR materials<\/p>\n

SR metamaterials also support programmable elastic-wave propagation, owing to their widely tunable elastic anisotropy, fully parameterized with only two parameters (\u03b1,\u00a0\u03bb0). To quantify direction-dependent wave propagation in SR metamaterials, we performed Bloch wave analysis on SR microstructures across the (\u03b1,\u00a0\u03bb0) design space\u2014here reinterpreted as representative volume elements in an infinite tiling (Fig. 5a<\/a> and Supplementary Information, section 2.16<\/a>). Leveraging the tunability of SR metamaterials, we designed acoustic metamaterials out of quasi-isotropic SR architectures. Using the effective refractive index concept from optics neff\u00a0=\u00a0c0\/vs, which compares the wave velocity in the constituent material c0 to the velocity vs of a wave travelling through an SR metamaterial made of the same constituent material, we find that for SR metamaterials, neff spans 1.5\u2009\u2264\u2009neff\u2009\u2264\u20098.3 in the (\u03b1,\u00a0\u03bb0) configuration space, demonstrating a substantial range in achievable refractive indices compared with the hexagonal and square lattices, which span 1.0\u2009\u2264\u2009neff\u2009\u2264\u20091.9 within the same relative-density limits (\\(20 \\% \\le \\bar{\\rho }\\le 75 \\% \\); Supplementary Information, section 2.17<\/a>).<\/p>\n

Fig. 5: Elastic lenses from SR metamaterials.\"Fig.<\/p>\n

a, SR samples with their corresponding directional velocity polar plots, demonstrating a wide range of achievable space of wave velocities and anisotropy. Curves, matching the colour of their respective geometries, denote the quasi-longitudinal wave velocities, and the green curve represents the wave velocity of the monolithic material. Samples in the grey\u2013black colour scale show the breadth of wave anisotropy achieved and the samples in the rainbow scale demonstrate isotropic SR unit cells used to construct the elastic lens in b. We use each sample as a unit cell, adding a frame with small gaps to provide boundary conditions for building the superstructure. The effective refractive index, neff, is defined as the ratio of the green to blue curves along the x direction. b, One-quarter of a Luneburg-inspired lens constructed from SR samples with varying refractive indices, highlighted by colour (blue, high refractive index; orange, low refractive index). c, Schematic illustrations of the target ray-focusing behaviour of an ideal Luneburg lens (top) and the designed refractive index profile of the full pixelated lens (bottom). d, Steady-state wavefield of the SR elastic Luneburg-like lens showing the elastic perturbation plane source (white) and focal point (red). The elastic wave is focused on the side opposite the source, as evidenced by the radial wave pattern emanating from the focal point. e, The normalized displacement intensity along the y-axis cross-section through the lens focal point (white dashed line in d). The Luneburg-like lens (solid blue curve) exhibits a sharp focal peak, whereas the randomized lens (inset, green dashed curve) shows no noticeable focusing. The residual intensity peak in the randomized lens\u2014three orders of magnitude smaller\u2014arises from weak energy trapping at the elastic circular boundary. A suitable excitation frequency was found in the range f\u00a0\u2248\u00a0646.7\u2009Hz, well within the long-wavelength regime of each pixel in the lens, where the refractive indices are still valid.<\/p>\n

The refractive index tunability of SR microstructures enables multiple elastodynamic functionalities such as wave guiding or lensing, with applications in subwavelength imaging and energy harvesting50<\/a>,51<\/a>. We designed a Luneburg-like cylindrical lens of diameter d in a pixelated gradient-index (GRIN) approach52<\/a> (Fig. 5b<\/a> and in Methods<\/a>), and validated it via steady-state dynamic finite-element method, focusing elastic perturbations from an incident plane wave (Fig. 5c<\/a>). Displacement fields reveal that waves traversing the lens focus on the opposite side, producing semi-circular wave fronts at the focal point and indicating that SR lenses qualitatively follow Luneburg patterns50<\/a>. The displacement intensity (squared displacement magnitude) contour along the vertical line through the focal point shows a peak that is an order of magnitude larger than the side peaks, indicating energy concentration into a region approximately 0.1d wide (Fig. 5d<\/a>). For comparison, a randomized reference lens constructed by randomly selecting SR cells at each radius produced a focal intensity roughly three orders of magnitude weaker, with a broader peak and side lobes of comparable magnitude, confirming that focusing arises from the GRIN arrangement of the SR pixels. Spanning such a breadth of neff enables the generation of graded materials with smoothly changing acoustic properties over a much larger range than any individual lattice family, achieved by tuning (\u03b1,\u00a0\u03bb0).<\/p>\n","protected":false},"excerpt":{"rendered":"We begin by determining the mechanical characteristics of SR microstructures, comparing them against three single-scale microstructures frequently used…\n","protected":false},"author":2,"featured_media":914172,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[32],"tags":[191180,1159,31842,139542,1160,79],"class_list":["post-914171","post","type-post","status-publish","format-standard","has-post-thumbnail","category-science","tag-complex-networks","tag-humanities-and-social-sciences","tag-mechanical-engineering","tag-mechanical-properties","tag-multidisciplinary","tag-science"],"_links":{"self":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts\/914171","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/comments?post=914171"}],"version-history":[{"count":0,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts\/914171\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/media\/914172"}],"wp:attachment":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/media?parent=914171"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/categories?post=914171"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/tags?post=914171"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}