Research Article

## Sharkskin-Inspired Magnetoactive Reconfigurable Acoustic Metamaterials

## Abstract

Most of the existing acoustic metamaterials rely on architected structures with fixed configurations, and thus, their properties cannot be modulated once the structures are fabricated. Emerging active acoustic metamaterials highlight a promising opportunity to on-demand switch property states; however, they typically require tethered loads, such as mechanical compression or pneumatic actuation. Using untethered physical stimuli to actively switch property states of acoustic metamaterials remains largely unexplored. Here, inspired by the sharkskin denticles, we present a class of active acoustic metamaterials whose configurations can be on-demand switched via untethered magnetic fields, thus enabling active switching of acoustic transmission, wave guiding, logic operation, and reciprocity. The key mechanism relies on magnetically deformable Mie resonator pillar (MRP) arrays that can be tuned between vertical and bent states corresponding to the acoustic forbidding and conducting, respectively. The MRPs are made of a magnetoactive elastomer and feature wavy air channels to enable an artificial Mie resonance within a designed frequency regime. The Mie resonance induces an acoustic bandgap, which is closed when pillars are selectively bent by a sufficiently large magnetic field. These magnetoactive MRPs are further harnessed to design stimuli-controlled reconfigurable acoustic switches, logic gates, and diodes. Capable of creating the first generation of untethered-stimuli-induced active acoustic metadevices, the present paradigm may find broad engineering applications, ranging from noise control and audio modulation to sonic camouflage.

## 1. Introduction

Acoustic metamaterials with tailored architectures exhibit unconventional capability in controlling acoustic waves [1–5] and have enabled a wide range of previously unachievable applications, such as superlensing [6–12], cloaking [13–15], logic operation [16–20], nonreciprocal propagation [21–25], topological insulation [20, 26–29], and wave guiding [30, 31]. Despite the diverse applications, most of the existing paradigms rely on architected structures with fixed configurations, and thus, their properties cannot be modulated once the structures are fabricated [1–5]. Emerging active acoustic metamaterials highlight a special opportunity to enable on-demand property switching of the fabricated acoustic metamaterials [32, 33]. However, the existing active metamaterials typically require tethered loads, such as mechanical compression or pneumatic actuation [20, 30, 31, 34]. Compared to tethered loads, untethered physical stimuli (such as electromagnetic field, light, and temperature) are more appealing owing to their special advantages including noncontact with the structures, nonlocal modulation, and rapid switching [35, 36]. However, using untethered physical stimuli to on-demand switch property states of acoustic metamaterials remains largely unexplored [35, 36]. Although there are recent reports on using the electromagnetic fields to modulate the effective constitutive parameters of acoustic metamaterials [35, 36], using untethered stimuli to enable reconfigurable acoustic devices such as waveguide, logic gates, or nonreciprocal diodes remains elusive.

Here, we report a class of active acoustic metamaterials whose configurations can be on-demand switched via untethered magnetic fields, thus enabling active switching of acoustic transmission, wave guiding, logic operation, and reciprocity. The key mechanism relies on magnetically deformable Mie resonator pillar (MRP) arrays that can be tuned between vertical and bent states corresponding to the acoustic forbidding and conducting, respectively. The MRPs are made of an iron-filled magnetoactive elastomer and feature wavy air channels to enable an artificial Mie resonance within a designed frequency regime [37–39]. The Mie resonance induces an acoustic bandgap, which is closed when pillars are selectively bent by a sufficiently large magnetic field (e.g., 0.13 T). These magnetoactive MRPs are further harnessed to design stimuli-controlled reconfigurable acoustic switches (i.e., shifting between different propagation pathways), reconfigurable acoustic logic gates (i.e., switching among NOT, AND, and OR gates), and reconfigurable acoustic diodes (i.e., switching between the nonreciprocal diode and reciprocal conductor). Integrating stimuli-responsive smart materials and Mie resonances, the present paradigm highlights a unique and promising avenue for acoustic metamaterials that can reversibly, repeatedly, and on-demand switch acoustic propagation, logic operation, and reciprocity via untethered physical stimuli.

## 2. Sharkskin-Inspired Design Principle

The design principle of magnetoactive reconfigurable acoustic metamaterials is inspired by the sharkskin [40–42]. Fast-swimming sharks feature skin denticles that are shaped like “V” trenches and aligned in the direction of fluid flow (Figure 1(a)). It has been proven that the aligned skin denticles can significantly reduce the flow drag because the V-shaped trenches can guide a turbulent flow to become a laminate flow [43]. It has been further discovered that tilting the skin denticles by a small angle can drastically increase the flow drag [44]. Sharks can smartly switch the skin flow drag by reversibly tilting the skin denticles (Figure 1(a)) [41, 45].

Inspired by this natural paradigm in tuning the flow drag, we here propose a class of magnetoactive MRP arrays that can smartly switch the acoustic transmission by reversibly bending MRPs via a remotely controlled magnetic field (Figure 1(b)–(g)). Each MRP (diameter and height ) has six-section wavy air channels that enable an artificial Mie resonance within a designed frequency regime [37–39]. When the pillar spacing is relatively small (e.g., ), the Mie resonances around adjacent MRPs are coupled to form an acoustic barrier to enable a low acoustic transmission within 9040-9140 Hz (<0.2, Figure 1(b), (c), and (h)). Besides, the pillars are made of a magnetoactive elastomer reinforced by ferromagnetic iron nanoparticles (see Materials and Methods and Fig. S1), and thus, they can be selectively bent via a sufficiently large magnetic field (e.g., 0.13 T) to open a large pillar-to-pillar spacing for the remaining vertical pillars, leading to a high acoustic transmission within 9040-9140 Hz (>0.8, Figure 1(d), (e), and (h)). When the remotely controlled magnetic field is reduced or turned off, the bent MRPs return to the vertical state, and thus, the acoustic transmission turns back to the low level within 9040-9140 Hz (<0.2, Figure 1(f)–(h)). Since the pillar deformation is fully elastic, the pillar bending is rapid, reversible, and repeatable, enabling on-demand switching of the acoustic transmission by tuning the applied magnetic field over multiple cycles (Figure 1(i), Movie S1 with cycle period 5 s, and Movie S2 with cycle period 1 s).

## 3. Mechanism of the MRP Array

Next, we analyze the mechanism of the MRP array. An MRP features six-section wavy air channels with the relative width as (Figure 2(a)). The propagation length of the acoustic wave within the air channel is longer than that in the solid medium. Effectively, the acoustic energy within a certain frequency regime is trapped within the air channels [37]. This point can be verified by a numerical simulation which shows that the acoustic energy is concentrated within the MRP center region around 9100 Hz (Figure 2(b) and S2a). This phenomenon can be explained as an artificial Mie resonance within the MRP structure [37–39, 46, 47]. From the perspective of an effective medium, this Mie resonance is corresponding to negative effective bulk modulus and positive effective density within frequency 9050-9120 Hz (Figure 2(c) and S2a). This frequency regime (9050-9120 Hz) is considered as the Mie resonance frequency. This frequency regime is consistent with the numerically simulated bandgap (Fig. S3). To verify this Mie resonance frequency, we further carry out analytical modeling of the MRP using an equivalent multiphase composite model (Supplementary Information, Fig. S4) [37, 48, 49]. The analytical model shows that the Mie resonance frequency of the MRP with diameter 1.5 cm and 64% solid volume fraction is ~8.9 kHz (Figure S4), which is relatively close to the numerical simulations (Figure 2(c)).

The MRPs can be harnessed to shield the acoustic wave if multiple MRPs are arranged in an array (Figure 2(d) and (e)). When the pillar spacing is relatively small (e.g., ), the Mie resonances in adjacent pillars are coupled to construct an acoustic shielding layer to block the acoustic transmission around the Mie resonance frequency regime (Figure 2(d) and (e)). This acoustic forbidding behavior at 9100 Hz is visualized by a numerical simulation shown in Figure 2(f) (Fig. S2b): the acoustic pressure amplitude drastically decreases to a very small value (<0.05) behind the MRP array. This behavior is validated by the experimentally measured acoustic transmission through the MRP array (Figure 2(g)): the acoustic transmission becomes lower than 0.15 within 9050-9120 Hz. The experimentally measured acoustic transmissions roughly agree with the numerical simulations over 8800-9400 Hz (Figure 2(g)). When the pillar spacing is relatively large (e.g., ), the acoustic wave can escape between MRPs to enable a relatively high acoustic transmission (>0.8, Figure 2(h)–(k) and S2b). To verify this phenomenon, we study a control pillar array with solid pillars without air channels ( cm and ) (Fig. S5). Both experiments and simulations show the acoustic transmissions over 8800-9400 Hz are above 0.8 (Fig. S5), implying that the air channels within the MRPs really play an important role to enable the low acoustic transmission around 9100 Hz. To further verify the phenomenon, we gradually vary the pillar spacing and numerically simulate the acoustic transmission through the MRP array. We find the acoustic transmission increases with increasing pillar spacing within 9050-9150 Hz (Fig. S2c).

In this work, the enlargement of the pillar spacing is enabled by the bending of selected pillars in the pillar array (Figure 1(d) and (e)). Once selected pillars are bent, the remaining vertical pillars construct a pillar array with a larger spacing (i.e., ) that leads to high acoustic transmission.

## 4. Magnetically Induced Bending of MRPs

To achieve precise control of the magnetically induced bending of MRPs, we need to understand the underlining mechanics. We first carry out experiments on the magnetically induced bending of a bottom-fixed MRP under a magnetic field with a tilted angle (Figures 2(l) and S6). With increasing magnetic field , the pillar bends by an angle (Figure 2(l)). When the angle reaches 15-18°, further increasing will suddenly bend and pinch the pillar onto the substrate. This can be understood as a critical condition for the magnetically induced buckling of the MRP. To quantify this phenomenon, we plot the bending angle in a function of the applied magnetic field (Figure 2(m)): after the critical magnetic field , the bending angle abruptly reaches 90°.

We then develop an analytical model to understand the magnetically induced buckling of the MRP. Using a similar analysis in the previously reported study on magnetically induced bulking of tilted beams (Supplementary Information, Fig. S7) [35], we express the critical magnetic field of the buckling as where is a dimensionless parameter dependent on the magnetic field angle and the aspect ratio of the MRP , is the magnetic permittivity of the vacuum, is Young’s modulus of the elastomer, is the second moment of area, is the cross-section area of the solid part, and is the pillar length. is the effective magnetic susceptibility difference between the axial and orthogonal directions and can be estimated as , where is the magnetic susceptibility of the elastomer. To validate Equation (1), we maintain the geometry of the pillar and vary the volume fraction of the ferromagnetic iron particle within the magnetoactive elastomer. We find that the critical magnetic field for the magnetically induced buckling decreases as the volume fraction of the iron particle increases from 3.14% to 18.48% (Figure 2(m)). According to Equation (1), when the geometry parameters and the magnetic field angle maintain constant, the critical magnetic field should scale with . As shown in Figure 2(n), this scaling law agrees with the experimental results for MRPs with various volume fractions of the iron particle (parameters are listed in Table S2).

Note that to enable the magnetic buckling of selected pillars shown in Figure 1(d) and (e), we use pillars with higher iron volume fraction at selected locations. For example, the central two MRPs have an iron volume fraction of 18.48% but the other two have an iron volume fraction of 3.14% in Figure 1(e). When a sufficiently large magnetic field (0.13 T) is applied, only the central two MRPs buckle. It should also be noted that the inhomogeneity of the iron volume fraction among four MRPs does not change the key physics of the Mie resonance shown in Figure 2(a)–(k). The numerical simulations of MRP array with inhomogeneous iron volume fractions display a low acoustic transmission within 9050-9150 Hz (Fig. S8), similar to that of the MRP array with homogeneous iron volume fraction (Figure 2(g)). The numerical results of the MRP array with inhomogeneous iron volume fractions also agree with the corresponding experimental results in Figure 1(h) (Fig. S8).

## 5. Magnetoactive Acoustic Double-Throw Switch

Next, we harness MRPs to design a magnetoactive acoustic double-throw switch (Figure 3). We first use MRPs to design a three-branch channel with two outputs (A and B) and one input (C) (Figure 3(a) and (b)). The channel wall is composed of MRP arrays with a low iron volume fraction (3.14%) and a small pillar spacing () to prevent the acoustic leaking within a designed frequency regime. Two sets of MRPs with a high iron volume fraction (18.48%) and a small pillar spacing () are located on the pathways of C-to-A and C-to-B. At the as-fabricated state, the acoustic transmissions of C-to-A and C-to-B are expected to be blocked by the MRPs around the Mie resonance frequency. This point can be first validated by numerical simulations that show the drastically reduced acoustic pressure at channels A and B (Figures 3(c) and S9a). Experiments show that the acoustic transmissions of C-to-A and C-to-B are both below 0.25 within 8800-8950 Hz (Figure 3(d)). We denote the as-fabricated state as “A off B off” state (Figure 3(a)–(d)). Note that this Mie resonance frequency region (8800-8950 Hz) is slightly different from that shown in Figures 1 and 2 (9050-9120 Hz). This difference can be explained by the geometrical inconsistency among different fabricated MRPs: numerical simulations show that the Mie resonance frequency varies from 8620 Hz to 9300 Hz by varying the pillar diameter by 7% (Fig. S10). This effect of pillar diameter variation can also be verified by the analytical modeling of an MRP which shows that the Mie resonance frequency changes from 9480 Hz to 8300 Hz by varying the pillar diameter from 1.4 to 1.6 cm (Fig. S4, Supplementary Information).

To switch the “A off B off” state to “A on B off” state, we use a magnetic field (0.13 T) to bend the MRPs on the pathway of C-to-A (Figure 3(e) and (f)). Numerical simulations show that the acoustic pressure at channel A increases drastically after bending the MRPs (Figure 3(g) and S9b). Experiments show that the transmission of C-to-A increases to ~0.8 within 8800-8950 Hz, while the transmission of C-to-B remains below 0.25 (Figure 3(h)). Similarly, when the MRPs on the pathway of C-to-B are bent by a magnetic field, the switch transforms to “A off B on” state (Figure 3(i)–(l) and S9c). When the MRPs on the pathway of C-to-A and C-to-B are all bent by respective magnetic fields, the switch transforms to the state of “A on B on” (Figure 3(m)–(p) and S9d). The above four states of the double-throw switch can be on-demand and reversibly modulated via the remotely controlled magnetic fields.

## 6. Magnetoactive Reconfigurable Acoustic Logic Gate

Next, we harness MRPs to design magnetoactive reconfigurable acoustic logic gates (Figure 4). Existing acoustic logic gates primarily rely on designed acoustic metastructures with fixed geometries [16–19], and very few of them can switch logic operators with tethered interventions [20]. Reconfigurable acoustic logic gates that can on-demand switch operators by untethered stimuli have not been explored. Here, we demonstrate a magnetoactive reconfigurable acoustic logic gate with a three-branch channel slightly modified from the acoustic switch shown in Figure 3: we use branches A and B as two inputs and branch C as the output (Figure 4). To design a NOT gate (Figure 4(a)), we place 6 MRPs (a high iron volume fraction 18.48% and a small pillar spacing ) at the central location of the channel. Effectively, 2 rows of MRP arrays are located on the pathway of A-to-C and 2 rows of MRP arrays on the pathway of B-to-C (Figure 4(b)). According to the results shown in Figure 3(d), (h), and (l), the acoustic transmission through one row of MRP array is ~0.35 at 8760 Hz. When input A and input B are both strong signals (normalized acoustic pressure amplitude 1), the output C behind two rows of MRP arrays is expected to feature a normalized pressure of at 8760 Hz (less than 0.5) (Figure 4(c)). We here denote that the normalized pressure equal to or larger than 0.5 as the digital “1” (strong) and otherwise as the digital “0” (weak). Numerical simulations show that the acoustic pressure amplitude at the output channel C is much smaller than the input channels A and B at 8760 Hz (Figure 4(d) and S11a). Experiments show that the normalized pressures in output C are all less than 0.5 whenever the inputs A and/or B are strong with a normalized pressure “1” or weak with a normalized pressure “0.2” at 8760 Hz (Figures 4(e) and S12a).

To switch the NOT gate to an AND gate (Figure 4(f)), one row of MRP array is bent by a magnetic field (i.e., 0.13 T), and thus, only one row of MRP array exists on the pathway of both A-to-C and B-to-C (Figure 4(g)). Since the acoustic transmission (normalized pressure) through one row of MRP array is ~0.35 at 8760 Hz, the normalized pressure in output C is expected to couple the signals from inputs A and B, leading to a normalized pressure of (Figure 4(h)). This point can be validated by numerical simulations which show that the normalized pressure amplitude at output C is around 0.6-0.75 at 8760 Hz (Figure 4(i) and S11b). Experiments show that the normalized pressure at output C is ~0.6 (above 0.5 and marked as “1”) when the inputs A and B are both strong (normalized pressure “1”) at 8760 Hz (Figures 4(j) and S12b). To further switch the AND gate to an OR gate (Figure 4(k)), we use three magnets to bend 6 MRPs in the center of the channel (Figure 4(l) and S13). Numerical simulations show that acoustic wave propagation to the output channel C does not have evident resistance at 8760 Hz (Figure 4(n) and S11c). Experiments show that the normalized pressure at output C is above 0.5 when either the input A or the input B is strong (normalized pressure “1”) at 8760 Hz (Figure 4(o) and S12c). Note that the above three acoustic logic operators can be on-demand and reversibly switched by controlling the untethered magnetic fields. Additional analyses and numerical studies show that the magnetoactive reconfigurable logic gates should work not only at 8760 Hz but also through 8700-8830 Hz and 8930-9070 Hz (Figs. S14-S15, more discussion in Supplementary Information). This point is further confirmed by experiments at 8700 and 9050 Hz which are located at two frequency branches, respectively (Fig. S16).

## 7. Magnetoactive Reconfigurable Acoustic Diode

Finally, we harness the MRPs to design a magnetoactive reconfigurable acoustic diode to smartly switch the reciprocity. Existing acoustic diodes or rectifiers that exhibit nonreciprocal acoustic wave propagation typically rely on acoustic metastructures with time-varying elements or nonlinearity [21–25, 50, 51]. An acoustic diode that can be switched by untethered stimuli has not been reported. Here, we harness untethered magnetic fields to demonstrate on-demand switching between the nonreciprocal acoustic diode and the reciprocal acoustic conductor.

The idea for the acoustic diode is based on a cloak-like waveguide (Figure 5(a)–(d)). When an engineering object is placed on the pathway of an acoustic wave, the acoustic pressure amplitude behind the object is significantly reduced due to the refection or absorption of the object (Figure 5(a), (b), S17a and b). Experimental results show that the acoustic transmission within 8500-8800 Hz is reduced by more than a factor of 2 when the object is placed (Figure 5(b)). When the object is surrounded by 6 MRPs with spacing , the direct interaction between the incoming wave and the object is shielded by the MRPs within a certain Mie resonance frequency (e.g., 8500-8760 Hz), and the wave is guided by a circular pathway (spacing ) between the 6 central MRPs and the outside curved MRP arrays (Figure 5(c)). Numerical simulations show that the acoustic pressure amplitude behind the object increases drastically after the 6-MRP envelope is installed (8650 Hz, Figures 5(c) and S17c). Experiments show that the acoustic transmission of the wave is around 80-90% of that of the reference state without the object over frequency range of 8500-8760 Hz (Figure 5(e)). Note that the waveguide at the cloak-like state is not a perfect acoustic cloak because the incident acoustic mode is disrupted by the structure [52]. Interestingly, when the front two MRPs are bent by a magnetic field (i.e., 0.13 T), the cloak-like state is destroyed because the incoming acoustic wave can have direct interaction with the object, thus leading to a significantly low acoustic transmission within 8500-9000 Hz (Figure 5(d), (e), and S17d). The numerical simulations show that the acoustic amplitude behind the object drastically reduces after bending the front MRPs (Figure 5(d)), with the acoustic pressure field similar to that of the object state shown in Figure 5(b). Figure 5(e) shows the acoustic transmission of four interesting states within 8500-8760 Hz: the cloak-like and the reference states exhibit high acoustic transmission, and the blocked and object states exhibit low acoustic transmission. The experimentally measured acoustic transmissions of the four states within 8500-9000 Hz roughly agree with the respective numerical simulations (Fig. S18). Besides, it should be noted that the cloak-like state (Figure 5(c)) and the blocked state (Figure 5(d)) can be reversibly and cyclically switched by controlling the magnetic field (Figure 5(f), Movie S3).

The acoustic diode (Figure 5(g)) is further designed based on the waveguide shown in Figure 5(a)–(d). The acoustic transmission of the MRP array with an applied magnetic field (i.e., 0.13 T) is relatively low (<0.4) within 8500-9000 Hz in the forward wave direction (Figure 5(h)). However, the acoustic transmission increases drastically to 0.7-0.8 within 8580-9000 Hz in the backward wave direction (Figure 5(i)), because the wave is guided by the MRPs to move around the object. The nonreciprocal wave propagation is visualized by the numerical simulations (Figure 5(h) and (i)) and validated by the experimentally measured acoustic transmissions in two directions (Figure 5(j)). The experimentally measured acoustic transmissions in two directions within 8500-9000 Hz agree with the numerical simulations (Fig. S19). Effectively, the acoustic diode demonstrated here resembles the reported acoustic diodes that are composed of a sonic crystal with its own band structure and a nonlinear medium to destroy the system symmetry [22, 23, 25]. In this system, the MRPs define the band structure and the interactions between the incoming wave and the encapsulated object serve as the nonlinear element. Note that the acoustic mode has been changed during the transmission through the diode structure. It is possible due to the structural complexity, similar to the reported acoustic diode by Liang et al. [23, 25], where the highly nonlinear microbubble medium may also change the acoustic mode.

When the magnetic field is reduced or turned off, the bent MRPs become vertical and the structure turns back to a state of the relatively high acoustic transmission (~0.7) in both directions within 8580-8860 Hz because the waves are guided in both directions (Figure 5(k)–(n)). Therefore, on-demand tuning the applied magnetic field can enable reversible switching between a nonreciprocal acoustic diode (Figure 5(g)–(j)) and a reciprocal acoustic conductor (Figure 5(k)–(n)) within 8580-8860 Hz.

## 8. Discussion

In summary, we report a class of active acoustic metamaterials whose configurations can be on-demand switched via untethered magnetic fields, thus enabling active switching of acoustic transmission, wave guiding, logic operation, and reciprocity. The mechanism primarily relies on synergistic integration of MRP arrays and their large deformation actuated by magnetic fields. The MRP arrays allow a large freedom for constructing various acoustic metadevices (e.g., waveguide, logic gate, and diode) via the judicious design of the MRP layout. The magnetically induced large deformation of MRPs allows the configuration modulation and thus leads to function switching of the acoustic metadevices without fabricating new structures. This study highlights a unique paradigm for applying stimuli-responsive smart materials to acoustic metamaterials and metadevices to enable active control of their acoustic properties [32, 33]. This paradigm may promote the integration between various smart or soft materials and acoustic metamaterials to achieve unprecedented functionalities [53, 54]. The unique paradigm may also promote the study of active acoustic metamaterials in switching of acoustic properties via other untethered stimuli, such as electric field, light, and temperature. Besides, as the first generation of untethered-stimuli-induced active acoustic metadevices (including waveguide, logic gate, and diode), these active acoustic metadevices may be used in broad engineering settings, ranging from on-demand noise control in smart infrastructures, automobile, and aircraft [55] and audio modulation for the next generation sound devices [56] to sonic camouflage [13–15]. Furthermore, the concept of the active acoustic metamaterials enabled by magnetoactive structures may be extended to design various switchable acoustic metadevices. For example, the magnetoactive elastomers can be used to design twisted helicoids to design orbital angular momentum metastructures [57] which can be tuned on and off via magnetically modulating the axial lengths. As another example, magnetoactive acoustic tunneling may be realized by assembling magnetoactive Mie resonator arrays into the tunnel, thus enabling actively switching between the tunneling with supercoupling and regular wave guiding [58]. Likewise, the paradigm may be extended to other active acoustic metadevices, such as superlenses [6–12] and topological insulators [20, 26–29, 59, 60].

## 9. Materials and Methods

### 9.1. Sample Preparation

Mie resonator pillars (MRPs) were fabricated with a molding process usingadditive manufacturing technique (Fig. S1). Molds were made with a fused deposition modeling (FDM) method with either ABS or dissolvable filament. The elastomer compound used for MRP is a mixture of 19.6 ml of liquid silicon rubber with various grams (5-35 g, Mold Max, Smooth-on) of iron nanoparticles (Sigma) (Table S2). After mixing the liquid rubber evenly with the iron nanoparticles for 3-5 minutes, the mixture was slowly filled in the molds. Once the elastomer mixture was cured for 4 to 5 hours, the molds were either dissolved with solvent, d-limonene, or manually dissembled. When the solvent was used, the elastomer structure was dried on the hot plate (60°C) overnight. Then, a thin elastomer strap (1-2 mm) was added around the top of the elastomer pillar to prevent the spreading or opening of MRPs due to their own weight. In each experiment, we attach the MRPs on the testing platform with a superglue (Gorilla).

### 9.2. Magnetic Actuation

We put a magnet under the testing platform, yet 20-30 mm ahead of the MRP. This applies a diagonal magnetic field to the MRP. The magnetic field intensity was changed by moving the magnet upward. When the distance between the magnet and the MRPs becomes small enough, the MRPs pinched on to the substrate. Once the magnet was moved downward or removed, the bent MRPs returned to the upright position.

### 9.3. Measurement of Acoustic Transmission

Acoustic experiments were conducted in testing chambers with acoustic insulators (cotton pads) on the chamber walls. Two loudspeakers (OT19NC00-04, Tymphany) controlled by a function generator (PI-8127, PASCO, USA) were placed in the chambers to provide the acoustic signals. The acoustic signals were collected by microphones (378B02 with 426E01, PCB Piezotronics, USA) and processed by a signal conditioner (482C05, PCB Piezotronics) and displayed by an oscilloscope (TBS1052B, Tektronix, USA). The noise was further removed with Savitzky–Golay filter in MATLAB to calculate the acoustic amplitude. The acoustic transmission was measured as , where and are acoustic amplitudes with and without the metastructures, respectively. The transmission measurement averaged over 2-4 measuring spots along a line normal to the wave propagation direction.

### 9.4. Acoustic Simulation

Numerical simulations were implemented with the acoustic module in COMSOL Multiphysics v5.2—a commercial finite element software. The set-ups for the numerical simulations are illustrated in Figs. S2, S9, S11, and S17. In a typical simulation, three material phases were involved in the numerical models: elastomer with high iron volume fraction (18.48%), elastomer with low iron volume fraction (3.14%), and air. The densities and sound speeds of the elastomers are shown in Table S2. Perfectly matching layers were employed to ensure the open boundary along the acoustic transport direction. The simulation was validated by benchmark calculations, and the mesh accuracy was ascertained through a mesh refinement study. The acoustic transmission was calculated by dividing the averaged acoustic pressure amplitude by the averaged background pressure amplitude over a line normal to the wave propagation direction.

## Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

## Authors’ Contributions

K.H.L. and Q.W. conceived the idea and designed the research. K.H.L. carried out experiments with the technical support of K.Y., A.X., and Z.F.. K.H.L., H.A., and Q.W. developed the analytical models and carried out the numerical simulations. K.H.L. and Q.W. wrote the first manuscript. All authors contributed to interpreting the results and revising the manuscript.

## Acknowledgments

Q.W. acknowledges the funding support from the Air Force Office of Scientific Research Young Investigator Program (FA9550-18-1-0192, program manager: Dr. Ming-Jen Pan) and the National Science Foundation (CMMI-1762567). We thank Dr. Guoliang Huang at the University of Missouri and Dr. Chu Ma from the Universit of Wisconsin-Madison for fruitful discussions.

## Supplementary Materials

Analogy between sharkskin denticles and Mie resonator pillars (MRPs). Modeling of MRPs. Theoretical analysis of magnetically induced buckling. Design principle of the reconfigurable acoustic logic gates. Table S1: experimentally measured materials and geometrical parameters of MRPs. Figure S1: fabrication process of MRPs. Figure S2: simulations of single MRP and MRP array. Figure S3: band structure of an MRP. Figure S4: theoretical analysis of an MRP. Figure S5: control solid pillar array. Figure S6: schematic for the placement of the magnet. Figure S7: schematic for a bottom-filed magnetoactive beam. Figure S8: simulation of the acoustic transmission of MRP arrays at two states. Figure S9: simulation configurations of the magnetoactive double-throw switch. Figure S10: effect of pillar spacing on acoustic transmission. Figure S11: simulation configurations of the magnetoactive logic gates. Figure S12: experimental results for logic gates. Figure S13: schematic for the magnet placement for logic gates. Figures S14-S15: simulations to reveal the mechanism for the acoustic logic gates. Figure S16: experimental results for logic gates at 8700 Hz and 9050 Hz. Figure S17: simulation configurations for the magnetoactive acoustic cloak-like waveguide. Figures S18-S19: simulations for the acoustic diode. Movie S1: magnetically induced cyclic bending of the central two pillars of an MRP array (period 5 s). Movie S2: magnetically induced cyclic bending of the central two pillars of an MRP array (period 1 s). Movie S3: magnetically induced cyclic bending of the front two pillars within an MRP-enabled cloak-like waveguide.* (Supplementary Materials)*

## References

- G. Ma and P. Sheng, “Acoustic metamaterials: from local resonances to broad horizons,”
*Science Advances*, vol. 2, no. 2, article e1501595, 2016. View at Publisher · View at Scopus · View at Google Scholar - S. A. Cummer, J. Christensen, and A. Alù, “Controlling sound with acoustic metamaterials,”
*Nature Reviews Materials*, vol. 1, no. 3, article 16001, 2016. View at Publisher · View at Scopus · View at Google Scholar - M. Yang and P. Sheng, “Sound absorption structures: from porous media to acoustic metamaterials,”
*Annual Review of Materials Research*, vol. 47, no. 1, pp. 83–114, 2017. View at Publisher · View at Scopus · View at Google Scholar - M. I. Hussein, M. J. Leamy, and M. Ruzzene, “Dynamics of phononic materials and structures: historical origins, recent progress, and future outlook,”
*Applied Mechanics Reviews*, vol. 66, no. 4, article 040802, 2014. View at Publisher · View at Scopus · View at Google Scholar - B. Assouar, B. Liang, Y. Wu, Y. Li, J.-C. Cheng, and Y. Jing, “Acoustic metasurfaces,”
*Nature Reviews Materials*, vol. 3, no. 12, pp. 460–472, 2018. View at Publisher · View at Scopus · View at Google Scholar - N. Kaina, F. Lemoult, M. Fink, and G. Lerosey, “Negative refractive index and acoustic superlens from multiple scattering in single negative metamaterials,”
*Nature*, vol. 525, no. 7567, pp. 77–81, 2015. View at Publisher · View at Scopus · View at Google Scholar - S. Zhang, L. Yin, and N. Fang, “Focusing ultrasound with an acoustic metamaterial network,”
*Physical Review Letters*, vol. 102, no. 19, article 194301, 2009. View at Publisher · View at Scopus · View at Google Scholar - A. Spadoni and C. Daraio, “Generation and control of sound bullets with a nonlinear acoustic lens,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 107, no. 16, pp. 7230–7234, 2010. View at Publisher · View at Scopus · View at Google Scholar - J. Zhu, J. Christensen, J. Jung et al., “A holey-structured metamaterial for acoustic deep-subwavelength imaging,”
*Nature Physics*, vol. 7, no. 1, pp. 52–55, 2011. View at Publisher · View at Scopus · View at Google Scholar - C. M. Park, J. J. Park, S. H. Lee, Y. M. Seo, C. K. Kim, and S. H. Lee, “Amplification of acoustic evanescent waves using metamaterial slabs,”
*Physical Review Letters*, vol. 107, no. 19, article 194301, 2011. View at Publisher · View at Scopus · View at Google Scholar - J. Li, L. Fok, X. Yin, G. Bartal, and X. Zhang, “Experimental demonstration of an acoustic magnifying hyperlens,”
*Nature Materials*, vol. 8, no. 12, pp. 931–934, 2009. View at Publisher · View at Scopus · View at Google Scholar - M. Molerón and C. Daraio, “Acoustic metamaterial for subwavelength edge detection,”
*Nature Communications*, vol. 6, no. 1, article 8037, 2015. View at Publisher · View at Scopus · View at Google Scholar - S. Zhang, C. Xia, and N. Fang, “Broadband acoustic cloak for ultrasound waves,”
*Physical Review Letters*, vol. 106, no. 2, article 024301, 2011. View at Publisher · View at Scopus · View at Google Scholar - L. Zigoneanu, B.-I. Popa, and S. A. Cummer, “Three-dimensional broadband omnidirectional acoustic ground cloak,”
*Nature Materials*, vol. 13, no. 4, pp. 352–355, 2014. View at Publisher · View at Scopus · View at Google Scholar - B.-I. Popa, L. Zigoneanu, and S. A. Cummer, “Experimental acoustic ground cloak in air,”
*Physical Review Letters*, vol. 106, no. 25, article 253901, 2011. View at Publisher · View at Scopus · View at Google Scholar - T. Zhang, Y. Cheng, J. Z. Guo, J. Y. Xu, and X. J. Liu, “Acoustic logic gates and Boolean operation based on self-collimating acoustic beams,”
*Applied Physics Letters*, vol. 106, no. 11, article 113503, 2015. View at Publisher · View at Scopus · View at Google Scholar - T. Zhang, Y. Cheng, B.-G. Yuan, J.-Z. Guo, and X.-J. Liu, “Compact transformable acoustic logic gates for broadband complex Boolean operations based on density-near-zero metamaterials,”
*Applied Physics Letters*, vol. 108, no. 18, article 183508, 2016. View at Publisher · View at Scopus · View at Google Scholar - F. Li, P. Anzel, J. Yang, P. G. Kevrekidis, and C. Daraio, “Granular acoustic switches and logic elements,”
*Nature Communications*, vol. 5, no. 1, article 5311, 2014. View at Publisher · View at Scopus · View at Google Scholar - Y. Wang, J. P. Xia, H. X. Sun, S. Q. Yuan, and X. J. Liu, “Binary-phase acoustic passive logic gates,”
*Scientific Reports*, vol. 9, no. 1, article 8355, 2019. View at Publisher · View at Scopus · View at Google Scholar - J. P. Xia, D. Jia, H. X. Sun et al., “Programmable coding acoustic topological insulator,”
*Advanced Materials*, vol. 30, no. 46, article 1805002, 2018. View at Publisher · View at Scopus · View at Google Scholar - B.-I. Popa and S. A. Cummer, “Non-reciprocal and highly nonlinear active acoustic metamaterials,”
*Nature Communications*, vol. 5, no. 1, article 3398, 2014. View at Publisher · View at Scopus · View at Google Scholar - N. Boechler, G. Theocharis, and C. Daraio, “Bifurcation-based acoustic switching and rectification,”
*Nature Materials*, vol. 10, no. 9, pp. 665–668, 2011. View at Publisher · View at Scopus · View at Google Scholar - B. Liang, X. S. Guo, J. Tu, D. Zhang, and J. C. Cheng, “An acoustic rectifier,”
*Nature Materials*, vol. 9, no. 12, pp. 989–992, 2010. View at Publisher · View at Scopus · View at Google Scholar - R. Fleury, D. L. Sounas, C. F. Sieck, M. R. Haberman, and A. Alù, “Sound isolation and giant linear nonreciprocity in a compact acoustic circulator,”
*Science*, vol. 343, no. 6170, pp. 516–519, 2014. View at Publisher · View at Scopus · View at Google Scholar - B. Liang, B. Yuan, and J. C. Cheng, “Acoustic diode: rectification of acoustic energy flux in one-dimensional systems,”
*Physical Review letters*, vol. 103, no. 10, article 104301, 2009. View at Publisher · View at Scopus · View at Google Scholar - R. Fleury, A. B. Khanikaev, and A. Alu, “Floquet topological insulators for sound,”
*Nature Communications*, vol. 7, no. 1, article 11744, 2016. View at Publisher · View at Scopus · View at Google Scholar - Y.-G. Peng, C. Z. Qin, D. G. Zhao et al., “Experimental demonstration of anomalous Floquet topological insulator for sound,”
*Nature Communications*, vol. 7, no. 1, article 13368, 2016. View at Publisher · View at Scopus · View at Google Scholar - C. He, X. Ni, H. Ge et al., “Acoustic topological insulator and robust one-way sound transport,”
*Nature Physics*, vol. 12, no. 12, pp. 1124–1129, 2016. View at Publisher · View at Scopus · View at Google Scholar - Z. Yang, F. Gao, X. Shi et al., “Topological acoustics,”
*Physical Review Letters*, vol. 114, no. 11, article 114301, 2015. View at Publisher · View at Scopus · View at Google Scholar - S. Babaee, N. Viard, P. Wang, N. X. Fang, and K. Bertoldi, “Harnessing deformation to switch on and off the propagation of sound,”
*Advanced Materials*, vol. 28, no. 8, pp. 1631–1635, 2016. View at Publisher · View at Scopus · View at Google Scholar - S. Babaee, J. T. B. Overvelde, E. R. Chen, V. Tournat, and K. Bertoldi, “Reconfigurable origami-inspired acoustic waveguides,”
*Science Advances*, vol. 2, no. 11, article e1601019, 2016. View at Publisher · View at Scopus · View at Google Scholar - M. Kadic, G. W. Milton, M. van Hecke, and M. Wegener, “3D metamaterials,”
*Nature Reviews Physics*, vol. 1, no. 3, pp. 198–210, 2019. View at Publisher · View at Google Scholar - K. Bertoldi, V. Vitelli, J. Christensen, and M. van Hecke, “Flexible mechanical metamaterials,”
*Nature Reviews Materials*, vol. 2, no. 11, article 17066, 2017. View at Publisher · View at Scopus · View at Google Scholar - P. Wang, F. Casadei, S. Shan, J. C. Weaver, and K. Bertoldi, “Harnessing buckling to design tunable locally resonant acoustic metamaterials,”
*Physical Review Letters*, vol. 113, no. 1, article 014301, 2014. View at Publisher · View at Scopus · View at Google Scholar - K. Yu, N. X. Fang, G. Huang, and Q. Wang, “Magnetoactive acoustic metamaterials,”
*Advanced Materials*, vol. 30, no. 21, article 1706348, 2018. View at Publisher · View at Scopus · View at Google Scholar - S. Xiao, G. Ma, Y. Li, Z. Yang, and P. Sheng, “Active control of membrane-type acoustic metamaterial by electric field,”
*Applied Physics Letters*, vol. 106, no. 9, article 091904, 2015. View at Publisher · View at Scopus · View at Google Scholar - Y. Cheng, C. Zhou, B. G. Yuan, D. J. Wu, Q. Wei, and X. J. Liu, “Ultra-sparse metasurface for high reflection of low-frequency sound based on artificial Mie resonances,”
*Nature Materials*, vol. 14, no. 10, pp. 1013–1019, 2015. View at Publisher · View at Scopus · View at Google Scholar - Z. Liang and J. Li, “Extreme acoustic metamaterial by coiling up space,”
*Physical Review Letters*, vol. 108, no. 11, article 114301, 2012. View at Publisher · View at Scopus · View at Google Scholar - Y. Xie, B.-I. Popa, L. Zigoneanu, and S. A. Cummer, “Measurement of a broadband negative index with space-coiling acoustic metamaterials,”
*Physical Review Letters*, vol. 110, no. 17, article 175501, 2013. View at Publisher · View at Scopus · View at Google Scholar - P. Ball, “Engineering shark skin and other solutions,”
*Nature*, vol. 400, no. 6744, pp. 507–509, 1999. View at Publisher · View at Scopus · View at Google Scholar - L. Wen, J. C. Weaver, and G. V. Lauder, “Biomimetic shark skin: design, fabrication and hydrodynamic function,”
*The Journal of Experimental Biology*, vol. 217, Part 10, pp. 1656–1666, 2014. View at Publisher · View at Scopus · View at Google Scholar - B. Dean and B. Bhushan, “Shark-skin surfaces for fluid-drag reduction in turbulent flow: a review,”
*Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences*, vol. 368, no. 1929, pp. 4775–4806, 2010. View at Publisher · View at Scopus · View at Google Scholar - D. W. Bechert, M. Bruse, W. Hage, J. G. T. van der Hoeven, and G. Hoppe, “Experiments on drag-reducing surfaces and their optimization with an adjustable geometry,”
*Journal of Fluid Mechanics*, vol. 338, pp. 59–87, 1997. View at Publisher · View at Scopus · View at Google Scholar - A. G. Domel, M. Saadat, J. C. Weaver, H. Haj-Hariri, K. Bertoldi, and G. V. Lauder, “Shark skin-inspired designs that improve aerodynamic performance,”
*Journal of the Royal Society Interface*, vol. 15, no. 139, article 20170828, 2018. View at Publisher · View at Scopus · View at Google Scholar - J. Oeffner and G. V. Lauder, “The hydrodynamic function of shark skin and two biomimetic applications,”
*Journal of Experimental Biology*, vol. 215, Part 5, pp. 785–795, 2012. View at Publisher · View at Scopus · View at Google Scholar - T. Brunet, A. Merlin, B. Mascaro et al., “Soft 3D acoustic metamaterial with negative index,”
*Nature Materials*, vol. 14, no. 4, pp. 384–388, 2015. View at Publisher · View at Scopus · View at Google Scholar - X. Zhu, B. Liang, W. Kan, Y. Peng, and J. Cheng, “Deep-subwavelength-scale directional sensing based on highly localized dipolar Mie resonances,”
*Physical Review Applied*, vol. 5, no. 5, article 054015, 2016. View at Publisher · View at Scopus · View at Google Scholar - C. Ying and L. Xiao-Jun, “Extraordinary resonant scattering in imperfect acoustic cloak,”
*Chinese Physics Letters*, vol. 26, no. 1, article 014301, 2009. View at Publisher · View at Scopus · View at Google Scholar - X. Hu, K.-M. Ho, C. T. Chan, and J. Zi, “Homogenization of acoustic metamaterials of Helmholtz resonators in fluid,”
*Physical Review B*, vol. 77, no. 17, article 172301, 2008. View at Publisher · View at Scopus · View at Google Scholar - Y.-X. Shen, Y.-G. Peng, D.-G. Zhao, X.-C. Chen, J. Zhu, and X.-F. Zhu, “One-way localized adiabatic passage in an acoustic system,”
*Physical Review Letters*, vol. 122, no. 9, article 094501, 2019. View at Publisher · View at Scopus · View at Google Scholar - X. Zhu, X. Zou, B. Liang, and J. Cheng, “One-way mode transmission in one-dimensional phononic crystal plates,”
*Journal of Applied Physics*, vol. 108, no. 12, article 124909, 2010. View at Publisher · View at Scopus · View at Google Scholar - H. Esfahlani, S. Karkar, H. Lissek, and J. R. Mosig, “Acoustic carpet cloak based on an ultrathin metasurface,”
*Physical Review B*, vol. 94, no. 1, article 014302, 2016. View at Publisher · View at Scopus · View at Google Scholar - Y.-X. Shen, Y. G. Peng, F. Cai et al., “Ultrasonic super-oscillation wave-packets with an acoustic meta-lens,”
*Nature Communications*, vol. 10, no. 1, article 3411, 2019. View at Publisher · View at Scopus · View at Google Scholar - H. Tang, Z. Chen, N. Tang et al., “Hollow-out patterning ultrathin acoustic metasurfaces for multifunctionalities using soft fiber/rigid bead networks,”
*Advanced Functional Materials*, vol. 28, no. 36, article 1801127, 2018. View at Publisher · View at Scopus · View at Google Scholar - R. Martínez-Sala, J. Sancho, J. V. Sánchez, V. Gómez, J. Llinares, and F. Meseguer, “Sound attenuation by sculpture,”
*Nature*, vol. 378, no. 6554, p. 241, 1995. View at Publisher · View at Scopus · View at Google Scholar - Y. Xie, T.-H. Tsai, A. Konneker, B.-I. Popa, D. J. Brady, and S. A. Cummer, “Single-sensor multispeaker listening with acoustic metamaterials,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 112, no. 34, pp. 10595–10598, 2015. View at Publisher · View at Scopus · View at Google Scholar - H. Esfahlani, H. Lissek, and J. R. Mosig, “Generation of acoustic helical wavefronts using metasurfaces,”
*Physical Review B*, vol. 95, no. 2, article 024312, 2017. View at Publisher · View at Scopus · View at Google Scholar - H. Esfahlani, M. S. Byrne, M. McDermott, and A. Alù, “Acoustic supercoupling in a zero-compressibility waveguide,”
*Research*, vol. 2019, article 2457870, 10 pages, 2019. View at Publisher · View at Google Scholar - Y. Ding, Y. Peng, Y. Zhu et al., “Experimental demonstration of acoustic Chern insulators,”
*Physical Review Letters*, vol. 122, no. 1, article 014302, 2019. View at Publisher · View at Scopus · View at Google Scholar - Y.-G. Peng, Y. Li, Y. X. Shen et al., “Chirality-assisted three-dimensional acoustic Floquet lattices,”
*Physical Review Research*, vol. 1, no. 3, article 033149, 2019. View at Publisher · View at Google Scholar