# Edge Modes and Asymmetric Wave Transport in Topological Lattices: Experimental Characterization at Finite Frequencies

###### Abstract

Although topological mechanical metamaterials have been extensively studied from a theoretical perspective, their experimental characterization has been lagging. To address this shortcoming, we present a systematic laser-assisted experimental characterization of topological kagome lattices, aimed at elucidating their in-plane phononic and topological characteristics. We specifically explore the continuum elasticity limit, which is established when the ideal hinges that appear in the theoretical models are replaced by ligaments capable of supporting bending deformation, as observed for instance in realistic physical lattices. We reveal how the zero-energy floppy edge modes predicted for ideal configurations morph into finite-frequency phonon modes that localize at the edges. By probing the lattices with carefully designed excitation signals, we are able to extract and characterize all the features of a complex low-frequency acoustic regime in which bulk modes and topological edge modes overlap and entangle in response. The experiments provide unequivocal evidence of the existence of strong asymmetric wave transport regimes at finite frequencies.

Acousto-elastic metamaterials and phononic crystals are artificially architected materials endowed with the capability to manipulate mechanical waves. This attribute makes them attractive for a variety of emerging technological applications such as acoustic cloaking Cummer and Schurig (2007); Chen and Chan (2007), sound manipulation and control Cummer et al. (2016), smart sensing Chen et al. (2014), imaging Craster and Guenneau (2012), and thermal management Davis and Hussein (2014). An important dynamical property of periodic metamaterials is their inherent directionality, which manifests as a frequency-selective spatial anisotropy of their bulk wave modes. The directionality patterns are dictated by the symmetry of the unit cell. By relaxing the cell symmetry, e.g. through actively reconfigurable microstructural elements, it is possible to alter the global directivity of the medium Celli and Gonella (2014); Celli et al. (2017).

Recent years have seen the advent of topological mechanical metamaterials, which have introduced a new paradigm for wave manipulation enabled by topologically protected wave modes. The concept originally emerged in the realm of quantum physics Haldane (1988); Kane and Mele (2005); Hasan and Kane (2010); Qi and Zhang (2011). For example, time-reversal invariant topological insulators are electronic materials that, in addition to their bulk bandgap behavior, possess topologically protected conducting edge states that are insensitive to defects Hasan and Kane (2010); Qi and Zhang (2011). Interestingly, these edge properties, which are only appreciated in finite domains, are controlled by topological invariants of the unit cell. Inspired by these intriguing properties, significant endeavors have been made to achieve mechanical analogues of these phenomena in the form of static Kane and Lubensky (2014); Paulose et al. (2015a); Rocklin et al. (2017, 2016); Stenull et al. (2016); Bilal et al. (2017) and dynamical Süsstrunk and Huber (2015); Nash et al. (2015); Wang et al. (2015a); Mousavi et al. (2015); Kariyado and Hatsugai (2015); Wang et al. (2015b); Pal et al. (2016); Brendel et al. (2017); Chaunsali et al. (2017, 2018); Prodan et al. (2017) phonon edge modes. A special class of topological phenomena occurs in Maxwell lattices, i.e., frame structures featuring ideal hinges that allow free rotations of the struts Maxwell (1864). In Maxwell lattices, the number of constraints equals that of the degrees of freedom in the unit cell Lubensky et al. (2015). The presence of free rotational mechanisms allows for lattice distortions that do not stretch the struts. This type of mechanism is referred to as a zero-frequency (energy) mode Mao and Lubensky (2018). An example of a Maxwell lattice is the regular kagome lattice (unit cell shown in Fig. 1(a)). While an infinite lattice can be seen as perfectly constrained everywhere, a finite domain with a free boundary necessarily presents zero modes, whose number is proportional to the size of the boundary. Depending on the unit cell, these zero modes can manifest as plane-wave-like features in the bulk or localize at the boundaries as so-called floppy edge modes Sun et al. (2012).

While the behavior of ideal topological lattices has been extensively studied theoretically using ball and spring models (and with a focus on the static response), experimental proofs of concept carried out on physical specimens (realized via fabrication techniques such as cutting, molding, or printing) have been rare Paulose et al. (2015b); Bilal et al. (2017). In this letter, we attempt to bridge this gap in the literature through a suite of experiments conducted on physical topological specimens. The main source of complexity in going from ideal to physical lattices is that the singular hinge connections at the lattice nodes are replaced by finite-thickness ligaments capable of supporting bending deformation. We refer to these conditions as the continuum elasticity limit, to emphasize that the mechanical behavior of the hinges depends on the geometric properties of the ligaments and can be described by continuum elasticity models. Adding bending stiffness lifts the zero-frequency modes to finite frequencies, resulting in a complex acoustic regime where bulk and edge phonons are superimposed and entangled in the response. Our experiments essentially attempt to disentangle the coexisting wave components and interpret them in light of the geometric and topological descriptors of the lattice. In parallel, we visually reconstruct the spatial patterns of the edge modes that localize at the floppy boundaries and we demonstrate the existence of strong asymmetric wave transport regimes at finite frequencies.

Reconstructing the phononic and topological characteristics of kagome lattices requires the capability to obtain in-plane measurements at points on the lattice surfaces. To this end, we use a 3D Scanning Laser Doppler Vibrometer (SLDV, Polytec PSV-400-3D), shown in Fig. 1(e). The lattice specimens are framed on the two sides, leaving the top and bottom edges free. The excitation is prescribed as a point force at the mid point of the bottom edge using an electrodynamic shaker (Bruel & Kjaer Type , powered by a Bruel & Kjaer Type amplifier) placed at the bottom of the structure for in-plane excitation. The specimens are manufactured via water-jet cutting from a sheet of acrylonitrile butadiene styrene (ABS) with the following material parameters: Young’s modulus , Poisson’s ratio , density . The unit cells of kagome lattice specimens corresponding to the configurations of Fig. 1(a) and (b) are shown in Fig. 1(c) and (d), respectively. The regular kagome unit cell is composed of two equilateral triangles connected with a ligament. The topological kagome unit cell contains one equilateral and one isosceles triangles. The side length of all the equilateral triangles for both lattices, and the longer side of the isosceles triangle in the topological lattice are all , while the shorter sides have length . The twist angle of the topological lattice is counter-clockwise, which makes the smaller angle between the triangles . and are the lattice vectors in real space. The ligaments are approximately in width. The specimens consist of arrays of unit cells. In order to excite the topological configuration from both the floppy and the non-floppy edges (which is necessary to probe the topological characteristics), we rotate the specimen by while keeping the excitation point fixed.

Our first preliminary task consists of experimentally reconstructing the bulk wave behavior. Beforehand we conduct a complete numerical characterization of the two unit cells of Fig. 1(c) and (d) using finite element analysis (FEA) and Bloch boundary conditions. The iso-frequency contours of the first phase-constant surface of the dispersion relation are presented in Fig. 2(a) and (b). For completeness, the band diagram for the regular kagome structure is also plotted in Fig. 2(g). Let us recall that band diagrams of ideal kagome lattices feature zero-energy modes involving displacements of lattice sites that occur without any deformation of the lattice Mao and Lubensky (2018). Here, in contrast, as the physical hinges prevent the triangular plates from rotating freely with respect to one another, we do not observe any zero-energy mode (except for the two translational invariants), leaving the phonon spectrum fully gapped except at , where is the wave vector in reciprocal space. Interestingly, the effect of finite-thickness hinges is analogous, in terms of added bending stiffness, to that obtained by adding next-nearest-neighbor interactions in ideal kagome lattices, discussed in Ref Sun et al. (2012) and recalled in Fig. S1. From a mechanics perspective, the rotation-hindering mechanisms are different between the two configurations: next-nearest-neighbor springs introduce additional elastic connections between selected nodes in the structure, while finite-width hinges only provide bending stiffness. However, both produce a shift toward finite frequencies of the isostatic zero-frequency branch in the direction and of the first optical branch at the point. To appreciate the frequency-selective spatial directivity intrinsic to these lattices, we track the evolution of the group velocity vector in the wave-vector plane by taking the gradient of the dispersion surface () with respect to the Cartesian components of the wave vector, i.e., and plotting iso-frequency contours of the gradient surface. Fig. 2(c) and (d) show the group velocity contours for the first phonon mode of the regular and topological kagome lattices, respectively. As expected, the contour of the regular kagome lattice features a strict 6-fold symmetry, while the symmetry is broken in the topological case to reflect the intrinsic asymmetry of the unit cell. In essence, the effects of cell symmetry relaxation manifest in the spatial directivity landscape.

We proceed now to conduct experiments on the finite lattice specimens described above using the setup shown in Fig. 1(e). The laser scans performed for this task are fairly parsimonious, involving one measurement point per triangle. By exciting regular and topological kagome lattices with tone bursts with carrier frequencies falling in the acoustic bulk modes range, we recover the symmetric and asymmetric wave-fields shown in Fig. 2(e) and (f), respectively. Both display wave beaming along the lattice vectors but markedly different degrees of symmetry, in conformance with the corresponding group velocity contours of Fig. 2(c) and (d), respectively. The experimental data also matches the results of full-scale FEA simulations of finite lattices with identical size and characteristics, which are reported in Fig. S2. By exciting the lattice with a broad-band chirp excitation, we can holistically reconstruct the entire phonon band diagram (in the frequency interval of the chirp), as shown in Fig. S3. Dedicated narrow band chirps spanning smaller frequency intervals can also be used to reconstruct partial sectors of the band diagram with higher fidelity. Finally, by performing local fine-scale scans, we can zoom in on selected details of the lattice cells to capture the localized deformation mechanisms of the hinges. For example, for a burst at (which excites the plateau region of the first branch along the direction), in Fig. 2(i) we plot the displacements of the scan points inside a few cells located along the corresponding direction in real space. Indeed, the deformation field presents a combination of transversal and rotational motion, consistent with the mode shape computed via unit cell analysis for the same conditions and shown in Fig. 2(h).

While the bulk wave characteristics are fully captured via Bloch analysis, the topological edge modes, which manifest on the scale of finite lattices, require additional modeling tools. In the context of ideal models, the topological properties are often described by introducing the topological polarization vector , calculated as Kane and Lubensky (2014):

(1) |

where are winding numbers obtained from

(2) |

where is a cycle along the contour of the Brillouin zone connecting and , where is a primitive reciprocal vector (such that ). is the compatibility matrix of the ideal lattice describing the displacement-strain relation in Fourier space. The topological polarization vector reveals the existence of floppy edge modes, and identifies the specific edge(s) where floppy motion is exponentially localized. From Eq. (1), we can see that for the topological cell, and therefore points to the top boundary, thus qualifying it as a floppy edge. Although in this study the lattice connections are realistic ligaments, leading to restrictions on the possible in-plane bending between the plates, we conjecture that the existence of edge modes with floppy characteristics is preserved. The open question is how these modes would manifest in the dynamic response. In essence, we borrow the terminology “floppy modes” from the ideal lattice case and we broaden its significance to describe the more general localization of phonon modes that may occur at low (but finite) frequencies in their realistic counterparts.

Computationally, the fate of the floppy edge modes can be predicted through a supercell analysis. To this end, we consider the finite strip of 6 cells shown in Fig. 3(a), and we apply 1D Bloch conditions along vector , while leaving the top and bottom edges free to mimic an infinitely long topological strip. As can be seen from the band diagram calculated using FEA and shown in Fig. 3 (b), we recognize the emergence of two new branches, both confined in frequency intervals lower than the dispersive region of the lowest bulk modes. The mode shapes shown in Fig. 3(i)-(vi), calculated along the lowest branch, feature indeed deformation localized at one end. This result qualifies these branches as edge modes and characterizes the left edge as the “floppy” boundary, consistent with the polarization vector predictions. These new modes represent the evolution of the floppy modes of ideal topological lattices in the limit of continuum elasticity. The fact that a topological lattice features two edge modes with distinct in-plane mechanisms constitutes an interesting departure from the conventional case of Rayleigh waves arising at the boundary of a 2D half-space, which feature a single mode of deformation Rayleigh (1885).

To prove experimentally the existence of these localized modes, we resort to full scans of the specimen. When we excite at sufficiently low frequencies falling in the range of the edge modes (the two lowest magenta lines in Fig. 3(b)), we capture the emergence of spatial behaviors that are germane to the topological configuration and that differ drastically according to which edge is excited (Fig. 4). Specifically, by exciting from the non-floppy edge, we mostly observe waves that propagate fast in the bulk with very long wavelength. In contrast, by exciting the floppy edge, the energy remains localized at the boundary. These profound differences are reflected in the spectral plane upon discrete Fourier transform (DFT) of the data. Indeed, when we excite from the non-floppy edge, we detect the conventional signature of asymmetric phononic modes (reported for completeness in Fig. S4). In contrast, when we excite from the floppy edge, most of the energy is concentrated in strong spectral features that align perfectly with the branches of the edge modes calculated via supercell analysis, as shown in Fig. 5.

In an effort to quantify the asymmetry of wave transport that has emerged from these tests, we calculate the energy transmissibility, i.e., the ratio between the kinetic energy localized at the edge opposite to the excitation and the energy confined at the edge containing the excitation point. In Fig. 6, we plot this quantity as a function of frequency for two scenarios: one where the excitation is applied at the non-floppy boundary, and the other for excitation applied at the floppy boundary. When we excite from the floppy edge at a frequency of , the transmissibility experiences a dip, which indicates that the energy is indeed localized at the floppy edge with modest propagation in the bulk. In contrast, for the same frequency, the excitation from the non-floppy edge results in a local maximum of the transmissibility, since the wave propagates through the bulk and reaches the opposite (floppy) edge, where it is eventually trapped by floppy mechanisms. However, when we excite above the edge mode frequency cut-off, the two transmissibility curves coalesce and we enter a conventional bulk wave regime.

In conclusion, we have experimentally demonstrated that by working in the limit of continuum elasticity, the phononic characteristics of topological kagome lattices are enriched by new edge modes at finite frequencies which can localize energy at the floppy boundary. This results in strong asymmetric wave transport capabilities in the very low-frequency acoustic regime. Our results help export a rich body of knowledge developed for static lattices to the dynamic regime, and provide innovative protocols for the dynamic experimental characterization of architected and lattice materials at large.

## I acknowledgment

The authors acknowledge the support of the National Science Foundation (NSF grant EFRI-1741618).

## References

- Cummer and Schurig (2007) S. A. Cummer and D. Schurig, New Journal of Physics 9, 45 (2007).
- Chen and Chan (2007) H. Chen and C. Chan, Applied physics letters 91, 183518 (2007).
- Cummer et al. (2016) S. A. Cummer, J. Christensen, and A. Alù, Nature Reviews Materials 1, 16001 (2016).
- Chen et al. (2014) Y. Chen, H. Liu, M. Reilly, H. Bae, and M. Yu, Nature communications 5, 5247 (2014).
- Craster and Guenneau (2012) R. V. Craster and S. Guenneau, Acoustic metamaterials: Negative refraction, imaging, lensing and cloaking, Vol. 166 (Springer Science & Business Media, 2012).
- Davis and Hussein (2014) B. L. Davis and M. I. Hussein, Physical review letters 112, 055505 (2014).
- Celli and Gonella (2014) P. Celli and S. Gonella, Journal of Applied Physics 115, 103502 (2014).
- Celli et al. (2017) P. Celli, S. Gonella, V. Tajeddini, A. Muliana, S. Ahmed, and Z. Ounaies, Smart Materials and Structures 26, 035001 (2017).
- Haldane (1988) F. D. M. Haldane, Physical Review Letters 61, 2015 (1988).
- Kane and Mele (2005) C. L. Kane and E. J. Mele, Physical review letters 95, 226801 (2005).
- Hasan and Kane (2010) M. Z. Hasan and C. L. Kane, Reviews of Modern Physics 82, 3045 (2010).
- Qi and Zhang (2011) X.-L. Qi and S.-C. Zhang, Reviews of Modern Physics 83, 1057 (2011).
- Kane and Lubensky (2014) C. Kane and T. Lubensky, Nature Physics 10, 39 (2014).
- Paulose et al. (2015a) J. Paulose, B. G.-g. Chen, and V. Vitelli, Nature Physics 11, 153 (2015a).
- Rocklin et al. (2017) D. Z. Rocklin, S. Zhou, K. Sun, and X. Mao, Nature communications 8, 14201 (2017).
- Rocklin et al. (2016) D. Z. Rocklin, B. G.-g. Chen, M. Falk, V. Vitelli, and T. Lubensky, Physical review letters 116, 135503 (2016).
- Stenull et al. (2016) O. Stenull, C. Kane, and T. Lubensky, Physical review letters 117, 068001 (2016).
- Bilal et al. (2017) O. R. Bilal, R. Süsstrunk, C. Daraio, and S. D. Huber, Advanced Materials 29 (2017).
- Süsstrunk and Huber (2015) R. Süsstrunk and S. D. Huber, Science 349, 47 (2015).
- Nash et al. (2015) L. M. Nash, D. Kleckner, A. Read, V. Vitelli, A. M. Turner, and W. T. Irvine, Proceedings of the National Academy of Sciences 112, 14495 (2015).
- Wang et al. (2015a) P. Wang, L. Lu, and K. Bertoldi, Physical review letters 115, 104302 (2015a).
- Mousavi et al. (2015) S. H. Mousavi, A. B. Khanikaev, and Z. Wang, Nature communications 6, 8682 (2015).
- Kariyado and Hatsugai (2015) T. Kariyado and Y. Hatsugai, Scientific reports 5, 18107 (2015).
- Wang et al. (2015b) Y.-T. Wang, P.-G. Luan, and S. Zhang, New Journal of Physics 17, 073031 (2015b).
- Pal et al. (2016) R. K. Pal, M. Schaeffer, and M. Ruzzene, Journal of Applied Physics 119, 084305 (2016).
- Brendel et al. (2017) C. Brendel, V. Peano, O. J. Painter, and F. Marquardt, Proceedings of the National Academy of Sciences 114, E3390 (2017).
- Chaunsali et al. (2017) R. Chaunsali, E. Kim, A. Thakkar, P. G. Kevrekidis, and J. Yang, Physical review letters 119, 024301 (2017).
- Chaunsali et al. (2018) R. Chaunsali, C.-W. Chen, and J. Yang, Physical Review B 97, 054307 (2018).
- Prodan et al. (2017) E. Prodan, K. Dobiszewski, A. Kanwal, J. Palmieri, and C. Prodan, Nature communications 8, 14587 (2017).
- Maxwell (1864) J. C. Maxwell, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 27, 294 (1864).
- Lubensky et al. (2015) T. Lubensky, C. Kane, X. Mao, A. Souslov, and K. Sun, Reports on Progress in Physics 78, 073901 (2015).
- Mao and Lubensky (2018) X. Mao and T. C. Lubensky, Annual Review of Condensed Matter Physics (2018).
- Sun et al. (2012) K. Sun, A. Souslov, X. Mao, and T. Lubensky, Proceedings of the National Academy of Sciences 109, 12369 (2012).
- Paulose et al. (2015b) J. Paulose, A. S. Meeussen, and V. Vitelli, Proceedings of the National Academy of Sciences 112, 7639 (2015b).
- Rayleigh (1885) L. Rayleigh, Proceedings of the London Mathematical Society 1, 4 (1885).

Supplemental Materials: Edge Modes and Asymmetric Wave Transport in Topological Lattices: Experimental Characterization at Finite Frequencies

## Appendix A Phonon Band Diagram with Next-nearest-neighbor Interactions

In an ideal regular kagome lattice, when the interaction between lattice nodes is only through the main struts, i.e., when we only consider nearest-neighbor interactions with spring constant , the triangles can freely rotate about their connection points, resulting in an isostatic zero-frequency branch in the direction. This is a zero-energy state, as the rotation of the plates does not result in any elongation or contraction of the struts. In the presence of additional springs with spring constants connecting the next-nearest neighbors (see Fig. S1 (a)), the zero-frequency branch is lifted to a finite-frequency range, as shown in Fig. S1 (b). The effect of is also felt on the first optical branch, where it has been shown that the non-zero optical frequency at the point is proportional to , Sun et al. (2012).

## Appendix B Full-scale finite element simulation

We conducted full-scale finite element simulations of a finite sheet of kagome lattice containing cells. The left and right sides of the lattice are fixed. A burst point excitation with a carrier frequency at is applied at the center of the bottom side of the lattice (see Fig. S2(a)). The red-to-blue color map refers to nodal displacements (large to small). As we can see, the propagation of the wave presents a strong directionality, which conforms with the lattice vector directions. The simulation results are consistent with the experimental observations reported in Fig. 2(e) in the letter. The simulation is repeated for the topological kagome lattice excited at . The results are shown in Fig. S2(b) and are consistent with the experimental ones shown in Fig. 2(f), revealing anisotropic wave beaming and symmetry relaxation.

## Appendix C Band diagram reconstruction via Fourier analysis

The band diagram of a lattice can be reconstructed via Fourier analysis of spatio-temporal data acquired with the laser vibrometer. For this task, the excitation must be engineered carefully to probe the desired frequency intervals in which the branch reconstruction is desired. For the regular kagome lattice case, we collect time histories of velocities at measurement points located along the direction (the problem becomes one-dimensional in space with respect to this fixed direction), thus forming a matrix of velocity data. This matrix is then fed to a 2D discrete Fourier transform (DFT) with respect to time and space, which yields a matrix of spectral velocity amplitudes with respect to frequency and wave number. The DFT is calculated with the following expression:

(S1) |

where and are the amplitude of velocity along/perpendicular to the selected direction in real and Fourier space, respectively; is the position along the selected direction, is the wave number along such direction, is time, is frequency, and and are the total number of evaluation points in space and time, respectively. 2D color maps of display the relation between frequency and wave number in the response. The dominant spectral features in these plots should conform to the branches of the band diagram calculated via unit cell analysis.

We excite the regular kagome lattice with a broadband chirp excitation with frequencies ranging from to . As shown in Fig. S3(a), we can successfully reconstruct the first two phonon bands. To obtain more precise reconstruction and zoom in on the branch details, we use three narrow band chirps, and the reconstructed phonon bands are shown in Fig. S3(b)-(d).

Similarly, we excite the topological kagome lattice from the non-floppy edge with a burst excitation at in order to compare the results with what observed for a similar excitation applied at the floppy edge and reported in Fig. 5(a) in the letter. Here we compare the spectral features obtained via DFT with the branches of the band diagram obtained from unit cell analysis, which capture the bulk modes. We observe that overall, the transversal acoustic mode matches the computational branch more precisely than the longitudinal mode.