# Ultrasound experiments on acoustical activity in chiral mechanical metamaterials

Tobias Frenzel, Julian Köpfler, Erik Jung, Muamer Kadic, Martin Wegener

► To cite this version:

Tobias Frenzel, Julian Köpfler, Erik Jung, Muamer Kadic, Martin Wegener. Ultrasound experiments on acoustical activity in chiral mechanical metamaterials. Nature Communications, 2019, 10, pp.3384. 10.1038/s41467-019-11366-8. hal-02867832

HAL Id: hal-02867832

https://hal.science/hal-02867832v1

Submitted on 27 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-019-11366-8

OPEN

# Ultrasound experiments on acoustical activity in chiral mechanical metamaterials

Tobias Frenzel<fnref n="1" />, Julian Köpfler<fnref n="1" />,², Erik Jung<fnref n="1" />, Muamer Kadic<fnref n="1" />,²,³ & Martin Wegener<fnref n="1" />,²

Optical activity requires chirality and is a paradigm for chirality. Here, we present experiments on its mechanical counterpart, acoustical activity. The notion "activity" refers the rotation of the linear polarization axis of a transversely polarized (optical or mechanical) wave. The rotation angle is proportional to the propagation distance and does not depend on the orientation of the incident linear polarization. This kind of reciprocal polarization rotation is distinct from nonreciprocal Faraday rotation, which requires broken time-inversion symmetry. In our experiments, we spatiotemporally resolve the motion of three-dimensional chiral microstructured polymer metamaterials, with nanometer precision and under time-harmonic excitation at ultrasound frequencies in the range from 20 to 180 kHz. We demonstrate polarization rotations as large as 22° per unit cell. These experiments pave the road for molding the polarization and direction of elastic waves in three dimensions by micropolar mechanical metamaterials.

<footnote n="1">¹Institute of Applied Physics, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany. ²Institute of Nanotechnology, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany. ³Institut FEMTO-ST, CNRS, Université de Bourgogne Franche-Comté, 25000 Besançon, France. Correspondence and requests for materials should be addressed to M.W. (email: Martin.Wegener@kit.edu)</footnote>
An object is geometrically chiral if it lacks centrosymmetry, mirror symmetries, and rotation-reflection symmetries¹⁻³. Chirality, or handedness, has played an important role in physics, optics, chemistry, biology, and medicine for decades already. Examples are optical activity⁴⁻⁷ and circular dichroism⁸,⁹, exclusively right-handed DNA double-helix strands¹⁰, exclusively left-handed neutrinos¹¹, chiral topological fermions¹², stereochemistry¹³, and the effect of drugs depending on the handedness of the underlying molecules¹⁴.

Surprisingly, chirality did not play a major role in the context of mechanics for many years. This has started to change recently. For example, geometrical chirality is important in plant seedpod mechanics¹⁵. Chirality has enabled topological effects at finite frequencies¹⁶,¹⁷, and has opened degrees of freedom beyond Cauchy elasticity in the stationary regime¹⁸,¹⁹. Handedness in shearing auxetics has created rigid and compliant structures²⁰. The reviews²¹⁻²⁵ provide an overview on the status of the field of mechanical metamaterials in general. Chiral phonons in quasi two-dimensional monolayers of noncentrosymmetric tungsten diselenide crystals have only recently been inferred from optical spectroscopy²⁶. However, one of the paradigms of chirality, namely “activity” has so far only been a theoretical possibility for generalized elastic continua²⁷⁻²⁹.

In the presence of optical activity⁴,⁷ or acoustical activity²⁷⁻²⁹, an incident transverse linear polarization of a propagating wave rotates: the eigenpolarizations of the chiral medium are not linear but rather circular, with a lifted degeneracy between left- and right-handed circular modes. Therefore, a linear incident polarization must be decomposed into the two circular-polarization eigenmodes, which propagate with different phase velocities. The resulting phase difference accumulated during propagation leads to a continuous rotation of the linear polarization axis versus the propagation coordinate. Optical and acoustical activity therefore allow, for example, to convert an incident linear polarization into the orthogonal one. Both are reciprocal effects, hence distinct from the nonreciprocal Faraday effect, which requires broken time-reversal symmetry, yet does not need chirality⁷. Regarding mechanics, yet a different kind of mode conversion has recently been observed in two-dimensional centrosymmetric metamaterials³⁰. There, the conversion between in-plane longitudinal

(compression) and in-plane transverse (shear) modes did not require chirality.

In this Letter, we report the experimental observation of the phenomenon of acoustical activity in three-dimensional (3D) chiral polymer metamaterials by direct spatiotemporal imaging at ultrasound frequencies, using cross-correlations of optical microscopy images taken under synchronized stroboscopic illumination. The experimental results agree well with numerical phonon band-structure calculations and simulations for the considered experimental geometry using finite-size samples.

## Results and discussion

**Beyond Cauchy elasticity.** Let us start by recalling that effects of chirality are neglected in mechanics on the level of textbook linear Cauchy elasticity³¹. Plainly speaking, Cauchy continuum mechanics can be seen as “point mechanics”³²⁻³⁴. Therefore, Cauchy elasticity becomes exact if the wavelength of an elastic wave, $\lambda$, is much larger than the crystal lattice constant, $a$, and if the sample size, $L$, is much larger than $a$, too. Therefore, at least one of the two conditions, $\lambda/a \gg 1$ and $L/a \gg 1$, needs to be violated to observe acoustical activity. Violating both conditions simultaneously is expected to yield larger effects. This means that we need to investigate the phonon band structure in the middle between the $\Gamma$-point and the Brillouin-zone edge (by symmetry the bands are degenerate at the Brillouin-zone edge) for beams of finite and infinite lateral extent $L_x = L_y$, with the number of unit cells $N_x = N_y = L_x/a = L_y/a$ in the two directions perpendicular to the wave propagation along $z$.

**Band-structure analysis.** Figure 1a shows a snapshot of a calculated chiral phonon eigenmode (at wave number $k_z = \pi/(4a)$, hence angular frequency $\omega_1 \approx 2\pi \times 107$ kHz) propagating in an infinitely extended 3D chiral micropolar metamaterial with $N_x = N_y = N_z \to \infty$ in the linear elastic regime. For clarity, the displacements are largely exaggerated. The behavior is qualitatively similar for other parameters. Panel b exhibits a simplified representation and panel c the underlying chiral metamaterial crystal unit cell (also see ref. ¹⁹), with geometrical parameters indicated.

<figure><img src="image_2.png" /></figure>

**Fig. 1** Chiral phonon eigenmode and unit cell. **a** Snapshot of a calculated chiral phonon eigenmode with eigenfrequency $\omega_1$ propagating along the z-direction in an infinitely extended 3D chiral micropolar metamaterial crystal. For clarity, all unit cells except for one column along the propagation axis (z-axis) are shown semitransparent and the displacements are largely exaggerated. **b** Simplified representation only showing the behavior of the unit cells' centers of mass, forming a helix that is moving along the z-axis versus time. An animation of **a** and **b** is shown in Supplementary Video 1. **c** Single-unit cell (cf. ref. ¹⁹). Geometrical parameters are cubic lattice constant $a = 250$ µm, $d/a = 0.06$, $r_1/a = 0.32$, $r_2/a = 0.4$, and $\delta = 34.8^\circ$. For the bulk constituent polymer material we have chosen the Young's modulus $E = 4.18$ GPa, mass density $\rho = 1.15$ g cm⁻³, and Poisson's ratio $\nu = 0.4$. The wave number is $k_z = \pi/(4a)$ and the angular frequency is $\omega_1 = 2\pi \times 107$ kHz, equivalent to $a/\lambda_0 = 0.0095$, with the wavelength $\lambda_0$ of pressure waves in the bulk constituent polymer material (cf. Fig. 2).
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Calculated phonon band structures for chiral metamaterial beams. The left vertical scales are in absolute units, the right scales, $a/\lambda_0$, are normalized and scalable to other parameters (cf. Fig. 1). We consider wave propagation with wave number $k_z$ in a beam which is infinitely extended along the z-direction and which contains $N_x \times N_y$ unit cells (with $N_x = N_y$) in the two orthogonal directions. The three panel columns correspond to $N_x = N_y = 1, 3$, and $\infty$, respectively. The upper row is for chiral ($\delta = 34.8^\circ$) metamaterials beams, the lower row for achiral ($\delta = 0$) ones. Transverse or shear-like bands are shown in red (a lifting of degeneracy reflects acoustical activity), pressure-like modes in blue, and twist-like modes in black. The latter are absent for $N_x = N_y = \infty$. For clarity, higher bands not relevant here are plotted in gray. The insets illustrate the corresponding geometries. The spatial mode depicted in Fig. 1 is highlighted by a circle. Parameters are as in Fig. 1</figcaption></figure>

Figure 2 gives an overview on the theoretical expectations for the corresponding band structures $\omega_n = \omega_n(k_z) = \omega_n(-k_z)$, with the angular frequencies $\omega_n$ ($n = 1, 2, ...$) and the wave number $k_z$. All statements concerning the character of the bands made below are based on an investigation of the associated eigenmodes. Here, we assume Bloch-periodic boundary conditions along the z-direction and open boundaries of the beams along the x- and y-direction for finite $N_x = N_y$ (left and center column in Fig. 2). For the limit $N_x = N_y \to \infty$ (right column in Fig. 2), we assume Bloch-periodic boundary conditions along the x- and y-direction as well. For a finite square-shaped cross-section $N_x \times N_y$, four different bands emerge from the Brillouin-zone center. All other bands ($n \ge 5$) are of lesser importance here and are hence plotted in gray for clarity. The lowest two modes ($n = 1, 2$) shown in red correspond to flexural waves with transverse circular polarization. Here, the center of mass of any given unit cell rotates clockwise or counter-clockwise around the wave vector $\mathbf{k}$, respectively, in circles around its rest position (cf. Fig. 1). In other words, these modes are chiral phonons (upper row of Fig. 2). As usual for thin plates and beams³⁵, the flexural modes start from the Γ-point with a parabolic dispersion. The behavior changes toward a linear dispersion if the cross-section of the beam is successively increased.

The absolute and the relative splitting between these two chiral phonon bands tends to zero near the Γ-point and, by symmetry, is also zero at the edge of the first Brillouin zone at $|k_z| = \pi/a$. Thus, maximum splitting occurs roughly around the middle of the one-dimensional Brillouin zone, where $|k_z| = \pi/(2a)$ or $\lambda/a = 4$ with $|k_z| = 2\pi/\lambda$. In analogy to optical activity, for a given frequency $\omega = \omega_1 = \omega_2$, the splitting in wave number, $\Delta k_z$, leads to a polarization rotation angle $\Delta\varphi$ due to acoustical activity over

propagation distance $L_z$ given by $\Delta\varphi = \Delta k_z L_z / 2$. The maximum splitting in Fig. 2 decreases with increasing beam cross-section $N_x = N_y = L_x / a = L_y / a$, but remains finite in the limit $N_x = N_y \to \infty$. This behavior is in line with our above intuitive reasoning starting from Cauchy elasticity as well as with our previous static experiments¹⁹. In the remainder of this Letter, we will focus on these two chiral phonon bands, i.e., on $n = 1, 2$.

Our analysis of phonons based on micropolar cubic-symmetry continua³² yields that these two transverse bands live in a subspace orthogonal to the other two bands ($n = 3, 4$). This statement is true for wave vectors along the three cubic principal directions, e.g., for the z-direction. We, therefore, find several crossings of these bands (rather than avoided crossings). Without chirality, the character of the third (black) and fourth (blue) bands in Fig. 2 would be that of a pure twist and that of a pure compression wave, respectively. In the presence of chirality, their characters are mixed. This band mixing is the direct dynamic counterpart of our recent static experiments on metamaterials with cubic symmetry¹⁹, in which an axial force imposed onto a beam led to a twist of the beam. In the limit $N_x = N_y \to \infty$, the twist mode (black) is absent because it is disallowed for Bloch-periodic boundary conditions along the x-, y-, and z-directions.

**Connection to acoustical activity.** The slopes of all four modes depend on the beam cross-section, hence on $N_x = N_y$. This behavior is the counterpart of the fact that the effective Young's modulus in the stationary case of micropolar materials is no longer a constant "material parameter", but rather depends strongly on the geometry, i.e., on the ratio $N_x = N_y = L_x / a = L_y / a$¹⁸. If we switch off chirality ($\delta = 0$ in Fig. 1a) as a control
calculation and use otherwise identical geometrical and con-
stituent material parameters (cf. Fig. 1), the behavior becomes
ordinary again and the slopes no longer change very much (see
lower row in Fig. 2).

For example, on the basis of the data shown for Nₓ = Nᵧ = 1 in
the upper row of Fig. 2, we estimate from the wave number
splitting at about |kz| = π/(2a) an acoustical activity of 17° per
unit cell. We have confirmed this estimate by independent
numerical calculations for samples with a finite extent along the
z-direction (see Supplementary Fig. 1). Here, we excite a linearly
polarized transverse wave at the bottom of the metamaterial
beam. Nz unit cells upwards, a plate with markers allows for
monitoring the local displacement vectors versus time. Clearly,
the open end of the sample leads to reflections of the elastic
waves. However, on the way back, the polarization rotation is
reversed due to time-reversal symmetry (unlike for the Faraday
effect, see section “Introduction”). Therefore, standing waves
along the beam axis do not change the polarization rotation at the
top of the beam. Standing waves may well influence the amplitude
of the displacements at the top though. As the z-axis has fourfold
rotational symmetry, the polarization rotation does not depend
on the orientation of the incident linear polarization (see
Supplementary Fig. 2).

**Measurement setup.** We have designed our experiments based on these calculations. We fabricate the polymer samples by 3D laser nanoprinting as described previously¹⁹. The dedicated home-built measurement setup has not been described previously: time-harmonic transverse flexural waves are excited by a piezoelectric transducer (Physik Instrumente, PICMA Chip Actuator). The samples are directly glued onto a dedicated holder, which has been 3D printed onto the transducer (see electron micrograph in Fig. 3 and Supplementary Fig. 3). The samples are imaged from the top by a microscope objective lens (Carl Zeiss, LD Achroplan 20×, NA = 0.4). They are stroboscopically illuminated by two infrared (850 nm center wavelength) light-emitting diodes (Vishay, VSLY 3850), the current through which is pulsed with a constant duty cycle of 1.5% (e.g., corresponding to 125 ns pulse duration for a frequency of 120 kHz), and synchronized with respect to the sinusoidal piezoelectric excitation. In this fashion, we take snapshots of the displaced structure versus time. By varying the time delay between excitation and illumination, we acquire slow-motion movies of the markers on top of the plate or at the side of the bottom of the sample. The latter serve to control the incident wave polarization and behavior. Using image cross-correlation analysis³⁶, we detect displacements that are much smaller than the wavelength of the light used for illumination and which are much smaller than one pixel of the complementary metal-oxide-semiconductor black/white camera (FLIR Systems BFLY-PGE-50S5M-C) used to record the images.

**Experimental data.** Example measurements merged into an oblique-view scanning electron micrograph of a metamaterial structure are shown in Fig. 3. The red (blue) scale bar refers to the measured displacement vectors in the middle of the top plate (at the side of the bottom of the sample). These data are shown for five consecutive oscillation periods, demonstrating the reproducibility of the measurements. In this particular example, we derive a polarization rotation of 44° (also see Supplementary Fig. 4). These and further results are summarized in Fig. 4 (full dots) and compared with results from the band-structure calculations (dashed curves) and results from frequency-domain finite-sample calculations (solid curves) already discussed above. The overall agreement between experiment and theory is very good. All

<figure><img src="image_4.png" /><figcaption>**Fig. 3** Oblique-view electron micrograph with overlaid measurement. The sample contains $N_x = N_y = 3$ and $N_z = 12$ unit cells (scale bar: 400 µm). A piezoelectric transducer excites the sample bottom. Top-view optical micrographs are taken under delayed stroboscopic illumination versus time delay. Displacement vectors are extracted by using image cross-correlation analysis. Results for the bottom side of the sample (blue, multiplied by a factor of $5 \times 10^3$) and for the middle of the top of the sample (red, multiplied by a factor of $10^4$) are blended into the electron micrograph. For both cases, five oscillation periods are depicted to emphasize the reproducibility of the experiments (also see Supplementary Fig. 4). From these data, we derive a rotation of linear polarization due to acoustical activity of 44°. Here, the excitation frequency is 160 kHz, all other parameters are as in Figs. 1 and 2. Results from many different similar experiments on different samples are summarized in Fig. 4</figcaption></figure>

qualitative trends are reproduced (as can be seen in Fig. 4a–c), the
rotation angle generally increases with increasing frequency. This
increasing deviation from Cauchy elasticity is expected based on
our reasoning that the fixed unit cell size increases with respect to
the wavelength with increasing frequency, moving away from the
limit of “point mechanics”. This behavior must not be mis-
interpreted in terms of a resonance. It also occurs in generalized
continuum theories²⁷⁻²⁹, which do not account for resonances.
Figure 4a additionally demonstrates that the rotation angle is
maximum for slender beams composed of only a single-unit cell
in the cross-section ($N_x = N_y = 1$). It decreases toward a certain
finite level versus increasing beam cross-section toward a 3D bulk
metamaterial ($N_x = N_y = L/a \gg 1$). The rotary power stays finite
in the bulk (see dashed black curve) because the unit cell size *a* is
finite compared to the wavelength $\lambda$, again deviating from the
limit of “point mechanics”. Figure 4b, c shows that the rotation
angle increases versus beam length or height $N_z$, regardless of the
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Measurements and calculations of acoustical activity. The rotation angle of an incident linear polarization upon propagation of a transverse elastic (flexural) wave is obtained from measurements (full dots) analogous to the example shown in Fig. 3 and plotted versus excitation frequency *f*. Results from band-structure calculations (cf. Fig. 2) are additionally shown as dashed curves, results from frequency-domain finite-sample calculations as solid curves (here $E = E' + E'' = 4.18$ GPa + i0.20 GPa). **a** Rotation angle versus excitation frequency *f* for different beam cross-sections $N_x = N_y$ and fixed beam height $N_z = 12$. **b** Same as **a**, but for fixed beam cross-section $N_x = N_y = 3$ and different heights $N_z$. **c** As **b**, but $N_x = N_y = 1$. The miniatures on the right-hand side, which are colored according to the experimental data points, illustrate the parameter variations. The largest measured rotation is 22° per unit cell for $N_x = N_y = 1$ and $f = \omega/2\pi = 180$ kHz. Achiral control samples show zero polarization rotation within the experimental error (see open blue dots in (b))</figcaption></figure>

beam cross-section (which is different in b and c). This behavior is analogous to the rotation angle in optical activity increasing proportionally with the propagation coordinate.

**Comparison with theory.** The calculations for finite metamaterial samples shown by the solid curves in Fig. 4 exhibit weak wiggles superimposed on a monotonous behavior. We have traced their origin back to mass-spring wobbling resonances, with the mass given mainly by the top plate. The wobbling eigen-frequency decreases with increasing total weight of the top plate, hence it decreases with increasing beam cross-section. Due to the

sample chirality, the wobbling motion induces an additional polarization rotation, which is an artifact in the sense that it does not occur in the band-structure calculations (which refer to $N_z \to \infty$). The behavior near the wobbling resonances depends on whether damping is accounted for. In the calculations shown in Fig. 4, we have added a finite imaginary part to the polymer's Young's modulus according to $E = E' + iE'' = 4.18$ GPa + i0.20 GPa. This complex value has been obtained from independent experiments on tuning forks made by the same 3D laser nanoprinting process under comparable conditions and measured at comparable excitation frequencies (see Supplementary Fig. 5). For comparison, Supplementary Fig. 6 shows computations for $E'' = 0$. Apart from the vicinity of the artifact resonances, the observed polarization rotation does not depend on $E''$. Intuitively, damping reduces the amplitudes of all polarization components alike; therefore, it changes neither the polarization state nor its evolution, making acoustical activity a robust phenomenon.

In conclusion, we have experimentally demonstrated the phenomenon of acoustical activity, the mechanical counterpart of optical activity. Based on generalizations of Cauchy elasticity, acoustical activity in chiral continua has previously been predicted theoretically. Our ultrasound experiments on 3D chiral mechanical metamaterials reveal linear polarization rotations as large as 22° per metamaterial unit cell. These chiral degrees of freedom provide an unprecedented control of the polarization of elastic waves in three dimensions. By scaling the metamaterial lattice constant *a*, a large range of operation frequencies is accessible.

## Methods

**Sample fabrication.** We fabricate the polymer samples by 3D laser nanoprinting (Nanoscribe GmbH, Photonics Professional GT) using a commercially available photoresist (IP-S, Nanoscribe GmbH). During the printing process a 25x objective lens (numerical aperture NA = 0.8, Carl Zeiss) was dipped directly into the photoresist. Underlying 3D models were created in STL-file format using COMSOL Multiphysics (COMSOL Inc.) and further processed using Describe (Nanoscribe GmbH). Therein, the 3D models were split into lines with a horizontal distance between adjacent lines of 0.3 µm and a vertical distance of 0.7 µm between adjacent planes. The scanning of the galvanometric mirrors led to a laser focus speed of 0.1 ms⁻¹, and the laser power was set to 45% of the maximum output power. To remove the unexposed photoresist, the samples were first placed in a bath of mr-Dev 600 (20 min), followed by a bath of acetone (2 min) and finally placed in a bath of isopropanol (20 min). Afterwards, the sample was ultraviolet cured for 20 min.

**Finite-element method (FEM) band-structure and eigenmode calculations.** For the finite-element band-structure calculations $\omega(k_z)$, we solved the eigenvalue problem using COMSOL multiphysics. We used a three-dimensional model typically meshed into a few million tetrahedra and applied second-order ansatz functions (MUMPS solver). In the case of the infinite crystal ($N_x = N_y = N_z = \infty$), the cubic unit cell as shown in Fig. 1c was implemented with Bloch-periodicity on all sides of the cube. In the cases of finite footprints ($N_x = N_y = 1, 2, 3, 5$ and $N_z = \infty$) $N_x \times N_y \times 1$ unit cells were implemented. Bloch boundary conditions were only applied along the z-direction, whereas all other boundaries were left open (i.e., traction free). A linear elastic material was assumed with Young's modulus $E = 4.18$ GPa, a Poisson's ratio $\nu = 0.4$, and a mass density $\rho = 1.15$ g cm⁻³. In order to isolate the structural properties of the unit cell from the ones of the constituent material (Young's modulus, Poisson's ratio, and mass density), the left vertical frequency scale *f* was converted to $a/\lambda_0$ (right vertical scale). Here, $\lambda_0$ is the wavelength of the constituent's longitudinal mode, given by $c/f$, with *c* being the phase velocity of the longitudinal mode $c = \sqrt{\frac{E(1-\nu)}{\rho(1-2\nu)(1+\nu)}}$.

**Finite-structure FEM calculations.** We deduced the polarization rotation angle of the finite structures containing $N_x \times N_y \times N_z$ unit cells by solving the eigenvalue problem for a given frequency $\omega$ using COMSOL multiphysics. As in the experiment, a plate with the same footprint $N_x a \times N_y a$ as the structure and with a thickness of 10 µm was attached to the top. We prescribed a displacement at the bottom of the metamaterial structure. All other boundaries were left open (i.e., traction free). The polarization of the eigenmode was calculated by averaging the displacement field components $u_x$ and $u_y$ over a square with side length $\frac{1}{2} a$ at the center of the top plate. This averaging size is comparable to the distance between the positions of the markers tracked in the experiment. A linear elastic material was
assumed with a complex-valued Young's modulus $E = 4.18$ GPa + i0.20 GPa, a Poisson's ratio $\nu = 0.4$, and a mass density $\rho = 1.15$ g cm$^{-3}$. The only exception are the dots in Fig. S6, where we have used $E = 4.18$ GPa (with zero imaginary part) for comparison.

## Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Received: 24 April 2019 Accepted: 8 July 2019

Published online: 29 July 2019

## References

1. Kelvin, T. W. The second Robert Boyle lecture. *J. Oxf. Univ. Jr. Sci. Club* **18**, 25 (1894).

2. Whyte, L. L. Chirality. *Nature* **180**, 513 (1957).

3. Whyte, L. L. Chirality. *Nature* **182**, 198 (1958).

4. Lindman, K. F. Über eine durch ein isotropes System von spiralförmigen Resonatoren erzeugte Rotationspolarisierung der elektromagnetischen Wellen. *Ann. Phys.* **63**, 621–644 (1920).

5. Lindell, I. V., Sihvola, A. H., Tretyakov, S. A. & Viitanen, A. J. *Electromagnetic Waves in Chiral and Bi-Isotropic*. (Media. Artech House Inc., Norwood, 1994).

6. Serdyukov, A., Semchenko, I., Tertyakov, S. & Sihvola, A. *Electromagnetics of Bi-Anisotropic Materials—theory and Application*. Vol. 11 (Gordon and Breach Science Publisher, Amsterdam, 2001).

7. Hecht, E. *Optics*. 4th ed. (Pearson Education Limited, Harlow, 2013).

8. Sharma, V., Crne, M., Park, J. O. & Srinivasarao, M. Structural origin of circularly polarized iridescence in jeweled beetles. *Science* **325**, 449–451 (2009).

9. Gansel, J. K. et al. Gold helix photonic metamaterial as broadband circular polarizer. *Science* **325**, 1513–1515 (2009).

10. Watson, J. D. & Crick, F. H. C. Molecular structure of nucleic acids. *Nature* **171**, 737–738 (1953).

11. Griffiths, D. *Introduction to Elementary Particles*. (Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2004).

12. Takane, D. et al. Observation of chiral fermions with a large topological charge and associated fermi-arc surface states in CoSi. *Phys. Rev. Lett.* **122**, 076402 (2019).

13. Gal, J. Pasteur and the art of chirality. *Nat. Chem.* **9**, 604–605 (2017).

14. Nguyen, L. A., He, H. & Pham-Huy, C. Chiral drugs: an overview. *Int. J. Biomed. Sci.* **2**, 85–100 (2006).

15. Armon, S., Efrati, E., Kupferman, R. & Sharon, E. Geometry and mechanics in the opening of chiral seed pods. *Science* **333**, 1726–1730 (2011).

16. Süsstrunk, R. & Huber, S. D. Observation of phononic helical edge states in a mechanical topological insulator. *Science* **349**, 47–50 (2015).

17. Ni, X., Weiner, M., Alü, A. & Khanikaev, A. B. Observation of higher-order topological acoustic states protected by generalized chiral symmetry. *Nat. Mater.* **18**, 113–120 (2019).

18. Lakes, R. S. & Benedict, R. L. Noncentrosymmetry in micropolar elasticity. *Int. J. Eng. Sci.* **20**, 1161–1167 (1982).

19. Frenzel, T., Kadic, M. & Wegener, M. Three-dimensional mechanical metamaterials with a twist. *Science* **358**, 1072–1074 (2017).

20. Lipton, J. I. et al. Handedness in shearing auxetics creates rigid and compliant structures. *Science* **360**, 632–635 (2018).

21. Bertoldi, K., Vitelli, V., Christensen, J. & Van Hecke, M. Flexible mechanical metamaterials. *Nat. Rev. Mater.* **2**, 17066 (2017).

22. Bauer, J. et al. Nanolattices: an emerging class of mechanical metamaterials. *Adv. Mater.* **29**, 1701850 (2017).

23. Kadic, M., Milton, G. W., Van Hecke, M. & Wegener, M. 3D metamaterials. *Nat. Rev. Phys.* **1**, 198–210 (2019).

24. Nieves, M. J., Carta, G., Jones, I. S., Movchan, A. B. & Movchan, N. V. Vibrations and elastic waves in chiral multistuctures. *J. Mech. Phys. Solids* **121**, 387–408 (2018).

25. Fernandez-Corbaton, I. et al. New twists of 3D chiral metamaterials. *Adv. Mater.* **31**, 1807742 (2019).

26. Zhu, H. et al. Observation of chiral phonons. *Science* **359**, 579–582 (2018).

27. Portugal, D. L. & Burstein, E. Acoustical activity and other first-order spatial dispersion effects in crystals. *Phys. Rev.* **170**, 673–678 (1968).

28. Bhagwat, K. V. & Subramanian, R. Acoustical activity of crystals: a comparative study of three descriptions. *Acta Cryst.* **A** **44**, 551–554 (1988).

29. Lakhtakia, A., Varadan, V. V. & Varadan, V. K. Elastic wave propagation in noncentrosymmetric, isotropic media: dispersion and field equations. *J. Appl. Phys.* **63**, 5246–5250 (1988).

30. Kweun, J. M., Lee, H. J., Oh, J. H., Seung, H. M. & Kim, Y. Y. Transmodal Fabry-Pérot resonance: theory and realization with elastic metamaterials. *Phys. Rev. Lett.* **118**, 205901 (2017).

31. Sommerfeld, A. *Mechanics of Deformable Bodies: Lectures on Theoretical Physics 2* (Academic Press, New York, 1950).

32. Eringen, A. C. *Microcontinuum Field Theories: I. Foundations and Solids*. (Springer Science & Business Media, New York, 1999).

33. Kadic, M., Frenzel, T. & Wegener, M. Mechanical metamaterials: when scale matters. *Nat. Phys.* **14**, 8–9 (2018).

34. Coulais, C., Kettenis, C. & Van Hecke, M. A characteristic length scale causes anomalous size effects and boundary programmability in mechanical metamaterials. *Nat. Phys.* **14**, 40–44 (2018).

35. Miklowitz, J. *The Theory of Elastic Waves and Waveguides*. 1st ed. (North Holland, Amsterdam, 2015).

36. Bückmann, T. et al. Tailored 3D mechanical metamaterials made by dip-in direct-laser-writing optical lithography. *Adv. Mater.* **24**, 2710–2714 (2012).

## Acknowledgements

We thank Andreas Naber and Alexander Münchinger (both KIT) for discussions. We acknowledge support by the Excellence Cluster “3D Matter Made to Order” (3DMM2O) funded through the German Excellence Strategy via Deutsche Forschungsgemeinschaft (DFG) and by the Carl Zeiss Foundation. We acknowledge additional support by the Helmholtz program “Science and Technology of Nanosystems” (STN) and through the associated “Virtual Materials Design” (VIRTMAT) project as well as by the Karlsruhe Nanostructure Service Laboratory (NSL). T.F. has been supported by the Karlsruhe School of Optics & Photonics (KSOP). T.F. and J. K. have been supported by the Heidelberg Karlsruhe Strategic Partnership (HEiKA) graduate school Functional Materials. J. K. acknowledges support by the Hector Fellow Academy. M. K. acknowledges support by the EIPHI Graduate School (contract “ANR-17-EURE-0002”) and by the French Investissements d'Avenir program, project ISITEBFC (contract “ANR-15-IDEX-03”). We acknowledge support by the KIT-Publication Fund of the Karlsruhe Institute of Technology.

## Author contributions

T.F. and M.W. conceived the experiment. M.W. supervised the work. T.F. conducted the main experiments; J.K. and E.J. performed the tuning fork experiments. T.F. performed numerical simulations with the help of M.K. T.F. analyzed all the experimental and numerical data. T.F. took the scanning electron micrographs. All authors discussed the results. M.W. wrote a first draft of the paper. All authors contributed to the final version of the paper.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-019-11366-8.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Peer review information:** *Nature Communications* thanks Fernando Fraternali and the other anonymous reviewer(s) for their contribution to the peer review of this work.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons

Attribution 4.0 International License, which permits use, sharing,

adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019