arxiv: v1 [cond-mat.mes-hall] 26 Aug 2017

Transcription

1 Subwavelength and directional control of flexural waves in plates using topological waveguides arxiv: v1 [cond-mat.mes-hall] 26 Aug 2017 Rajesh Chaunsali, 1 Chun-Wei Chen, 1, and Jinkyu Yang 1 1 Aeronautics and Astronautics, University of Washington, Seattle, WA, USA, (Dated: August 29, 2017) Inspired by the quantum spin Hall effect shown by topological insulators, we propose a plate structure that can be used to demonstrate the pseudo-spin Hall effect for flexural waves. The system consists of a thin plate with periodically arranged resonators mounted on its top surface. We extend a technique based on the plane wave expansion method to identify a double Dirac cone in such a locally resonant plate structure and to realize topological effects at a subwavelength scale. Building on this method, we verify numerically that this topological plate structure not only guides low-frequency flexural waves robustly along a desired path, but also imparts a unique one-way propagation characteristic to them. PACS numbers: n a Cd I. INTRODUCTION A topological insulator has emerged as a new state of matter in condensed matter physics. This is a special type of insulator that conducts electricity only on its surface/boundary. Here topology is relevant because one can predict the boundary properties of these finite materials (i.e., finite-sized lattices) solely by knowing the bulk properties of infinite materials (i.e., infinitely large lattices). Topological framework provides an elegant way to categorize the bulk properties in terms of a topological invariant, and thus, one can expect a topological protection and a degree of robustness for the boundary properties 1,2. It is recent that this whole framework dealing with the flow of electrons has evolved further and influenced other areas such as photonics 3 and acoustics It has also propelled new design paradigm for artificial mechanical structures, so-called topological mechanical metamaterials, to achieve unconventional static (zero frequency) and dynamic (non-zero frequency) responses. In particular, tailoring non-zero frequency responses, i.e., elastic waves in structures, on topological grounds shows tremendous potential to be used for energy harvesting, novel sensing, and impact mitigation purposes 32. One of the most unique topological effects is the quantum spin Hall effect, the underlying phenomenon shown by topological insulators 33,34. These systems are passive in the sense that they do not require any external force/field, but still, possess directional boundary states. This is due to the presence of Kramers partners, i.e., two opposite spins of electron, which travel in the opposite directions on their boundaries, keeping the time reversal symmetry intact. Although mechanical counterparts, being bosonic systems, do not possess these intrinsic spins, one can carefully design the system to have two pseudo-spins by imposing certain symmetries in the lattice, and thus realize pseudo-spin Hall effect 17,20,23,26. While previous studies have successfully reported the feasibility of the pseudo-spin Hall effect in mechanical settings, such realization in continuum plate structures has not been thoroughly explored yet. This is partly due to the complexity in guided stress wave propagation in plates. Also, it is necessary to ingeniously modify the plate design to form a double Dirac cone, which is essential for creating two opposite pseudo-spins of the spin Hall effect. Moreover, the studies on mechanical spin Hall systems so far have relied on the Bragg scattering effect in structures. This means that realizing topological effects at low frequencies would require big structures, and that can restrict their practical usage. Therefore, it would be a significant advancement to the current research trend if one can demonstrate low-frequency topological effects in a continuum mechanical structure such as plate, which is ubiquitous in many engineering disciplines. In order to deal with the aforementioned challenges, we propose a topological plate system that consists of a thin plate with periodically arranged local resonators mounted on its top surface. Inspired by the local resonance mechanism 35, the current system is a subwavelength metamaterial 36. To design such a system, we employ a technique based on the combination of the classical plate theory and the plane wave expansion (PWE) method 36 38, which enables fast and efficient calculation of the wave dispersion relation. Specifically, we invoke the mathematically folding of bands 39 to create a double Dirac cone for low frequencies of the propagating flexural waves. Building on this method, we verify numerically using the finite element method (FEM) that this topological plate structure can guide flexural waves robustly at subwavelength scales in desired directions. Though there have been numerous studies on plate-based waveguides in the past 40 50, we have shown that this locally resonant topological plate system has potential to robustly guide low-frequency flexural wave with a unique one-way propagation characteristic, not found

2 2 FIG. 1: (a) Lattice arrangement of the resonators with a rhombus shaped unit cell in red. Dimensional parameters and lattice vectors are also shown. (b) Detailed view of the unit cell with six resonators mounted on top of a thin plate. in conventional waveguides. The structure of this manuscript is as follows: in section II, we describe the design of the topological plate. In section III, we present the PWE method to calculate dispersion relation. In section IV, we show the mathematical folding of bands and create a double Dirac cone at a low frequency. In section V, we show the formation of a band gap around the double Dirac frequency by perturbing the lattice, which facilitates the system to transition from a topologically trivial state to a non-trivial state. In section VI, we employ the FEM to show the existence of two local modes, each designated by a pseudo-spin, at the interface of topologically trivial and non-trivial lattices. In section VII, we demonstrate the feasibility of guiding low-frequency flexural waves robustly along a path with various bends and having the one-way propagation characteristic. In section VIII, we conclude this manuscript. II. DESCRIPTION OF THE LOCALLY RESONANT TOPOLOGICAL PLATE Our system consists of a thin plate on which multiple local resonators are attached to form a lattice arrangement (Fig. 1). The rhombus shaped unit cell is of length a and consists of six resonators in a hexagonal arrangement (Fig. 1a). Each resonator is at a distance R and rotationally symmetric from the center of the unit cell, showing the C 6 symmetry. a 1 and a 2 are the lattice vectors. We model the resonators as cylindrical heads attached to the plate with a thin neck (Fig. 1b). In order to invoke the topological effects in the system, we will only vary radius R, keeping the C 6 symmetry intact in this unit cell. As a substrate material, we choose an aluminum plate (E = 77.6 GPa, ρ = 2, 730 kg/m 3, ν = 0.352) of thickness h = 1 mm with the unit cell of size a = 45 mm. The resonator neck is made of acrylic plastic (E neck = 3.2 GPa) with h neck = 5 mm and d neck = 2 mm, whereas the resonator head is made of tungsten (ρ head = 19, 260 kg/m 3 ) with h head = 14 mm and d head = 9 mm. The aforementioned material properties are based on nominal values of standard materials. III. CALCULATION OF THE UNIT CELL DISPERSION We first present the PWE method to show dispersion characteristics of the system. To this end, we simplify the resonator design with a lumped mass (m = πρ head d 2 head h head/4) connected to the plate with a linear spring (β = πe neck d 2 neck /4h neck). In this process, we neglect the bending motion of the resonators and only consider their motion in the out-of-plane z-direction. This is because the bending frequencies are away from the Dirac point we will be investigating (see Appendix A). Consequently, the bending motion would not affect the topological behavior emerging around the Dirac point. Because h a, the plate can be assumed to be thin 37, so the transverse motion of the plate is calculated as per the classical plate theory (i.e., Kirchhoff-Love theory) 51. Following the approaches taken by Pal et al. 24, Xiao et al. 38, and Torrent et al. 46, governing equations for the time-harmonic vibration of the unit cell with angular frequency ω can be written as D 4 w(r) ω 2 ρhw(r) = β α [w(r α ) w(r α )]δ(r R α ), (1a) ω 2 m w(r α ) = β[w(r α ) w(r α )], (1b) where D = Eh 3 /12(1 ν 2 ) represents the flexural rigidity of the plate, r = (x, y) denotes the generalized coordinate of the plate, w(r) represents the transverse displacement of the plate, and w(r α ) represents the displacement of the resonating masses attached at points R α. We have α = 1, 2,..., 6 for six different resonating masses per unit cell and δ(r R α ) is a delta function in two dimensions.

3 3 We introduce the following non-dimensional angular frequency Also, the mass of the resonator can be normalized as Ω = ωa 2 ρh D. γ = m ρa c h, where A c = 3 a 2 /2 is the area of a unit cell. We can write the normalized resonance frequency of a resonator as Ω R = a 2 (β/m)ρh/d. Due to periodicity in the system, we employ the PWE method to calculate the dispersion properties of the unit cell. We write the displacement of the plate as a superposition of multiple plane waves such that w(r) = G W (G)e i(k+g) r, (2) where W (G) is a plane wave coefficient, K is the Bloch wave vector, and G denotes the reciprocal lattice vector given by G = pb 1 + qb 2, in which p and q are integers, and b 1 and b 2 are the basis vectors of the reciprocal lattice. We truncate the summation with respect to G by choosing both p and q as M, (M 1),..., 0,..., (M 1), M. Therefore, the reciprocal space is a N N finite grid with N = 2M + 1. The displacement of the plate at the locations where the resonators are attached can be simply deduced as w(r α ) = G W (G)e i(k+g) Rα. (3) Substituting Eq. (2) and Eq. (3) into Eq. (1a) D G K + G 4 W (G )e i(k+g ) r ω 2 ρh G W (G )e i(k+g ) r = β α [ w(r α ) G W (G )e i(k+g ) R α ] δ(r R α ). (5) (4) Multiplying both sides with e i(k+g) r, we obtain [ D ] K + G 4 ω 2 ρh W (G )e i(g G) r = β α G Taking the area integral over the entire unit cell of area A c leads to [ D ] K + G 4 ω 2 ρh W (G ) e i(g G) r dr 2 G = β α A c e i(k+g) r [ w(r α ) G W (G )e i(k+g ) R α ]δ(r R α ). (6) [ w(r α ) G W (G )e i(k+g ) R α ] A c e i(k+g) r δ(r R α )dr 2. (7) We now use the following relations A c e i(g G) r dr 2 = { A c, if G = G 0, otherwise (8a) A c f(r)δ(r R α )dr 2 = f(r α ) (8b) to obtain ] A c [D K + G 4 ω 2 ρh W (G) = β α [ w(r α ) G W (G )e i(k+g ) R α ] e i(k+g) Rα. (9)

4 Using the Bloch s theorem for the resonators, we write w(r α ) = w(0 α )e ik Rα, where w(0 α ) represents the Bloch displacement of the resonator (indexed with α) at the reference unit cell. Thus, we deduce [a 4 K + G 4 Ω 2 ] W (G) = γω 2 R α e ig Rα [ w(0 α ) G W (G )e ig R α ]. (10) Similarly, we simplify the second governing Eq. (1b) (for α = 1, 2,..., 6) as [ ] Ω 2 w(0 α ) = Ω 2 R W (G)e ig Rα w(0 α ). (11) G Given the N N size of the reciprocal space, we arrange Eq. (10) and Eq. (11) in the form of an eigenvalue problem to solve for Ω at a specific Bloch wave vector K and obtain the dispersion relation. Mathematically, this can be expressed as [ ] { } P11 P 12 W (G) P 22 w(0 α ) P 21 { } = Ω 2 W (G). (12) w(0 α ) Here, K + G P 11 = a 4 0 K + G K + G N N 4 G 1 R 1 + γω 2 G 2 [ ] R 2 [ ] R exp i R1 R 2 R 6 exp i G1 G 2 G N N ;.. G N N R 6 G 1 R 1 P 12 = γω 2 G 2 [ ] R exp i R1 R 2 R 6 ; P 21 = Ω 2 R 2 [ ] R exp i G1 G 2 G N N ;.. P 22 = I 6, G N N where exp and I 6 represent the exponential function and the 6 6 identity matrix, respectively. We chose N = 7 for further calculations. R 6 4 IV. BAND FOLDING AND SUBWAVELENGTH DESIGN By using the aforementioned technique, we calculate the dispersion relation for the hexagonal arrangement of resonators, i.e., R = a/3 (Fig. 2a). Torrent et al. 46 showed the existence of a single Dirac cone in such a system. Building on this finding, in the current study, we create a double Dirac cone (two Dirac cone dispersion curves superimposed). For the given configuration, we observe such a double Dirac cone at f d = 0.94 khz at the Γ point (Fig. 2a). This is possible because we have chosen a bigger unit cell consisting of six resonators instead of two (compare the unit cells in Fig. 2b of different colors and corresponding Brillouin zones below). Thus, the dispersion curves fold onto a smaller Brillouin zone 39. Note that the physics is the same in both representations and it is a mathematical folding of bands. However, achieving a double Dirac cone which is a key ingredient of the spin Hall systems guides us to focus on the bigger unit cell and tailor it to realize topological effects (to be further discussed in the Section V). In Fig. 2a, we also plot the analytical parabolic dispersion relation for a bare plate (i.e., the identical plate as the substrate described in Section II, but without local resonators attached). It is given by Ω p = (ak p ) 2, where Ω p and k p are the normalized angular frequency and the wave number of flexural wave, respectively. This gives us an idea of the wavelength λ p (= 2π/k p ) with which the flexural wave would propagate in the bare plate (host medium) at a given frequency. Therefore, the flexural wavelength in the host medium at the Dirac frequency would

5 5 FIG. 2: (a) Dispersion plot of the locally resonant plate, showing a double Dirac cone for the hexagonal lattice with R = a/3 (red curves). Parabolic dispersion of the bare plate is in light gray color. (b) Two unit cell representations to enable mathematical folding of bands (the corresponding Brillouin zones below). The smaller cell (enclosed by blue dashed lines) represents the typical, irreducible unit-cell configuration, while in this study, we consider the bigger unit cell (enclosed by red dashed lines) for topological manipulations. Their sizes are compared with the wavelength in the host medium at the Dirac frequency. be λ d = 2π/ [ ρh(2πf d ) 2 /D ] 1/4. This is indicated by the star marker on the parabolic dispersion curve in Fig. 2a. For the chosen set of design parameters, λ d is approximately 2.3 times longer than the length of the large unit cell (i.e., a), and 4 times longer than the size of the small unit cell (a/ 3). Figure 2b shows the relative sizes of the unit cells compared to this wavelength, marking its subwavelength characteristics. Therefore, as the topological effects will be seen around the Dirac frequency, this opens up new pathways to realizing them for large-wavelengths flexural waves by using a relatively smaller substrate. Note that the Dirac frequency f d is less than the resonant frequency of each resonator f R = (1/2π) β/m = 1.72 khz. Torrent et al. 46 showed analytically that f d f R. By further reducing the resonant frequency, it is possible to shift the Dirac point to even lower frequency regime, thereby making the design deep-subwavelength. However, practical challenges in designing such a system can limit the same, e.g., due to heavy resonating masses leading to a static bending of the thin plate. V. BAND INVERSION AND TOPOLOGY We now vary the radius R and see its effects on the wave dispersion in the system. For R < a/3, as shown in Fig. 3a, there emerges a band gap near the Dirac frequency. Because we still keep the C 6 symmetry intact, two modes on each side (lower or higher side) of the gap are degenerate at the Γ point. Seen in the insets are the corresponding degenerate mode shapes of the plate at the Γ point, which are obtained by the PWE method described in Section III. Here, the lower frequency modes are of p type (p 1 and p 2 as shown in the bottom panel of Fig. 3a), and the higher frequency modes are d type (d 1 and d 2, upper panel in Fig. 3a) as per the analogy to electronic orbital shapes. For R > a/3, however, the band gap still exists, but its topological characteristic is different from the earlier case. As shown in Fig. 3b, the degenerate modes are flipped, i.e., d type modes are at the lower frequency compared to p type modes. This band-inversion as we vary R around R = a/3 indicates a typical topological transition in the system. The validity of this result based on lumped mass model is verified and discussed in Appendix A in comparison with the FEM (using COMSOL), which takes into account all geometrical features in the resonator design. The presence of degenerate modes around the band gap has important implications in realizing pseudo-spin Hall effect. One can take linear combinations of these modes and construct two alternate modes, i.e., pseudo-spin modes, without changing the physics of the system. Let p ± = p 1 ± ip 2 and d ± = d 1 ± id 2 represent such spin modes for these degenerate points. The sign in the middle determines if these are rotating clockwise or counterclockwise. We can interpret the dispersion near the Γ point in terms of the pair of spins by projecting the eigen-states onto the spin basis {p ±, d ± }. Therefore, the effective Hamiltonian of the system around the Γ point reduces to the one for Cd/Te/HgTe/CdTe quantum well 34 and would resemble a mechanical pseudo-spin Hall system. One can show that the bands have non-zero spin Chern number for the case with R > a/3, hence, proving it to be topologically

6 6 FIG. 3: Band inversion when the unit cell is perturbed around R = a/3 configuration. (a) A case with R < a/3 leading to topologically trivial band gap with p type modes having lesser frequency than d type modes (insets). (b) A case with R > a/3 leading to topologically non-trivial band gap with p and d type modes inverted. non-trivial 10,11,39. VI. EMERGENCE OF TOPOLOGICAL BOUNDARY STATE Now that we verified the feasibility of the double Dirac cone formation and the band inversion in the unit-cell level, we move to the investigation of wave guiding characteristics in multi-cell configurations. To account for more complicated geometry and boundary conditions in such a multi-cell setting, we resort to the FEM. According to the bulk-boundary correspondence of topology 1, we expect distinct behaviors on the boundaries of topologically trivial and non-trivial lattices. One way to observe it clearly is to have topologically distinct lattices placed adjacently and investigate their intersections for a non-trivial local response. To this end, we take a supercell, which consists of both topologically trivial (R = 0.8a/3) and non-trivial (R = 1.1a/3) lattices, 10 units of each placed as one strip (Fig. 4a). The periodic boundary condition is introduced in the direction of another lattice vector (at 60 from the horizontal). In this way, such a system provides a quick way to calculate vibration responses at the intersection and monitor their propagation along the periodic direction. In Fig. 4b, we plot eigenfrequencies of the supercell as a function of wave number in the periodic direction. The presence of two modes inside the band gap (shown in blue and red curves) is especially striking, since those have the following non-trivial properties. First, they represent two types of pseudo-spin modes localized at the intersection: one rotates clockwise, while the other rotates counterclockwise. Second, both have opposite group velocities at a given frequency. Figure 4c shows the respective mode shapes corresponding to points S1 and S2 in Fig. 4b, which are excited at 0.91 khz. Opposite spins and group velocities of these modes can be verified by looking at the harmonic evolution of these modes (see Supplementary Movie 1). We also plot in-plane time-averaged mechanical energy flux (I j = σ ij v j, where σ ij and v j are stress tensor and velocity vector, respectively) over a harmonic cycle as black arrows. This further confirms the spin nature of these flexural modes in the locally resonant plate. There is a small gap at the Γ point for these spin modes. This indicates that these pseudo-spin modes are topological but not strictly topologically protected. The protection is only guaranteed by the C 6 symmetry, which we break at the interface of the two lattices. However, we will show in the next section that these modes can still be used to build robust and

7 7 FIG. 4: (a) A supercell made by placing topologically distinct lattices adjacently. (b) Eigenfrequencies of the supercell as a function of the wave number in the periodic direction. The local modes at the interface have been highlighted in red and blue colors. In background are the other modes, i.e., bulk modes and local modes at the extreme left and right end of the supercell. (c) Pseudo-spin modes at the interface corresponding to the points S1 and S2 in (b). The color intensity represents the out-of-plane displacement, and arrows indicate the time-averaged mechanical energy flux confirming their spin nature. directional waveguides. The gap at the Γ point could be reduced by minimizing the effect of the C 6 symmetry breaking, for example, by constructing a graded interface between two topologically distinct lattices (see Appendix B). VII. DIRECTIONAL WAVEGUIDES In the previous section, we have shown that the mechanical spin Hall effect enables us to have two pseudo-spins at one frequency but in the opposite directions. In order to demonstrate how this property can be used to build unconventional waveguides on plates, we combine topologically trivial and nontrivial lattices to form a 2D structure. Figure 5 shows the waveguide (three linear segments with two bends) along the interface of two types of lattices. We give a forced excitation in the z-direction at the center of the plate (indicated by the star symbol) in such a way that we selectively excite spin modes. This could be done, for example, by choosing multiple points in the vicinity but with a phase difference in their forcing. Insets show the three point excitation to a nontrivial unit cell (i.e., with R = 1.1a/3) where the spin is predominately d type. We extract the phase information from the spin modes in the supercell analysis done earlier for 0.91 khz (see the resonators in red and blue colors, representing out-of-phase oscillations in Fig. 4c). We apply this phased excitation as F 1 = F exp(iωt), F 2 = F exp(iωt + 2π/3), and F 3 = F exp(iωt + π). Note that the phase differences in the three excitation points are not equally spaced, but show π/3 and 2π/3 differences between the neighboring ones (i.e., π/3 between F 2 and F 3, and 2π/3 between F 1 and F 2 ). This excitation tactic induces clockwise spin in Fig. 5a and counterclockwise spin in Fig. 5b. We enforce low-reflecting boundary conditions on the plate and perform harmonic analysis using the FEM. We confirm the unique features of this topological waveguide. The clockwise spin mode propagates to the left (Fig. 5a) and the counterclockwise spin propagates only to the right (Fig. 5b). These spin waves propagate robustly along the waveguide interface, and even though there are sharp bends, there is no back scattering as the spins remain intact (see Supplementary Movie 2). These simulation results imply that by using the pseudo-spin Hall effect induced in the locally resonant plate structure, we can propagate flexural waves in selected path and direction direction without resorting to the breakage of time-reversal symmetry. That is, without the necessity of any active components, we simply need to create a topological boundary and perturb the host medium strategically via selective excitation.

8 8 FIG. 5: Pseudo-spin dependent one-way propagation of flexural wave. Star indicates the location of excitation. Multi-point phased excitation in force is used to excite clockwise (a) and counterclockwise (b) spin modes. VIII. CONCLUSIONS We have proposed a locally resonant plate structure to demonstrate the pseudo-spin hall effect for the subwavelength and directional control of flexural waves in plates. We simplify the resonator design with a lumped mass model and employ the plane wave expansion method, which enables us to efficiently investigate the key design parameters responsible for forming a double Dirac cone at low frequencies. Keeping the C 6 symmetry intact, we perturb the unit cell and show the band-inversion process, leading to topologically distinct lattice configurations. When two topologically distinct lattices are placed adjacently, we show the existence of two pseudo-spin modes traveling in the opposite directions along the interface. This unique feature is used to build topological waveguides with multiple bends and robustly guiding the spin-dependent flexural waves in one direction. The finding could be useful in designing compact and robust one-way channels for guiding low-frequency flexural waves in applications such as energy harvesting, sensing, and impact mitigation. Future studies include the optimization of the locally resonant plate configurations by using the proposed numerical techniques, as well as the experimental verification of the waveguiding effects, which will be reported by the authors future publications. ACKNOWLEDGMENTS We gratefully acknowledge fruitful discussions with Krishanu Roychowdhury (Cornell University), Rui Zhu (Beijing Institute of Technology), Raj Kumar Pal (Georgia Institute of Technology), Simon Yves (CNRS), Romain Fleury (EPFL), Cheng He (Nanjing University), Zhiwang Zhang (Nanjing University), and Panayotis Kevrekidis (University of Massachusetts, Amherst). We are grateful for the support from NSF (CAREER and EFRI ). Equally contributed first author 1 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 2 X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 (2011). 3 L. Lu, J. D. Joannopoulos, and M. Solja cić, Nat. Photonics 8, 821 (2014). 4 M. Xiao, G. Ma, Z. Yang, P. Sheng, Z. Q. Zhang, and C. T. Chan, Nat. Phys. 11, 240 (2015). 5 Z. Yang, F. Gao, X. Shi, X. Lin, Z. Gao, Y. Chong, and B. Zhang, Phys. Rev. Lett. 114, (2015). 6 J. Lu, C. Qiu, M. Ke, and Z. Liu, Phys. Rev. Lett. 116, (2016). 7 C. He, X. Ni, H. Ge, X.-C. Sun, Y.-B. Chen, M.-H. Lu, X.-P. Liu, and Y.-F. Chen, Nat. Phys. 12, 1124 (2016). 8 Z. G. Chen and Y. Wu, Phys. Rev. Appl. 5, (2016).

9 9 9 R. Fleury, A. B. Khanikaev, and A. Alù, Nat. Commun. 7, (2016). 10 Z. Zhang, Q. Wei, Y. Cheng, T. Zhang, D. Wu, and X. Liu, Phys. Rev. Lett. 118, (2017). 11 S. Yves, R. Fleury, F. Lemoult, M. Fink, and G. Lerosey, New J. Phys. 19, (2017). 12 C. L. Kane and T. C. Lubensky, Nat. Phys. 10, 39 (2014). 13 J. Paulose, B. G. Chen, and V. Vitelli, Nat. Phys. 11, 153 (2015). 14 D. Z. Rocklin, B. G. G. Chen, M. Falk, V. Vitelli, and T. C. Lubensky, Phys. Rev. Lett. 116, (2016). 15 O. Stenull, C. L. Kane, and T. C. Lubensky, Phys. Rev. Lett. 117, (2016). 16 O. R. Bilal, R. Süsstrunk, C. Daraio, and S. D. Huber, Adv. Mater. 29, (2017). 17 R. Süsstrunk and S. D. Huber, Science 349, 47 (2015). 18 L. M. Nash, D. Kleckner, A. Read, V. Vitelli, A. M. Turner, and W. T. M. Irvine, Proc. Natl. Acad. Sci. 112, (2015). 19 P. Wang, L. Lu, and K. Bertoldi, Phys. Rev. Lett. 115, (2015). 20 S. H. Mousavi, A. B. Khanikaev, and Z. Wang, Nat. Commun. 6, 8682 (2015). 21 T. Kariyado and Y. Hatsugai, Sci. Rep. 5, (2015). 22 Y.-T. Wang, P.-G. Luan, and S. Zhang, New J. Phys. 17, (2015). 23 R. K. Pal, M. Schaeffer, and M. Ruzzene, J. Appl. Phys. 119, (2016). 24 R. K. Pal and M. Ruzzene, New J. Phys. 19, (2016). 25 J. Vila, R. K. Pal, and M. Ruzzene, arxiv: (2017). 26 S.-Y. Yu, C. He, Z. Wang, F.-K. Liu, X.-C. Sun, Z. Li, M.-H. Lu, X.-P. Liu, and Y.-F. Chen, arxiv: (2017). 27 G. Salerno, T. Ozawa, H. M. Price, and I. Carusotto, Phys. Rev. B 93, (2016). 28 C. L. Fefferman, J. P. Lee-Thorp, and M. I. Weinstein, Proc. Natl. Acad. Sci. 111(24), 8759 (2014). 29 R. Chaunsali, F. Li, and J. Yang, Sci. Rep. 6, (2016). 30 R. Chaunsali, E. Kim, A. Thakkar, P. G. Kevrekidis, and J. Yang, Phys. Rev. Lett. 119, (2017). 31 E. Prodan, K. Dobiszewski, A. Kanwal, J. Palmieri, and C. Prodan, Nat. Commun. 8, (2017). 32 S. D. Huber, Nat. Phys. 12, 621 (2016). 33 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, (2005). 34 B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 314, 1757 (2006). 35 Z. Liu, X. Zhang, Y. Mao, Y. Y. Zhu, Z. Yang, C. T. Chan, and P. Sheng, Science 289, 1734 (2000). 36 D.-L. Yu, G. Wang, Y.-Z. Liu, J. Wen, and J. Qiu, Chinese Phys. 15, 0266 (2006). 37 M. M. Sigalas and E. N. Economou, J. Appl. Phys. 75, 2845 (1994). 38 Y. Xiao, J. Wen, and X. Wen, J. Phys. D. Appl. Phys. 45, (2012). 39 L.-H. Wu and X. Hu, Phys. Rev. Lett. 114, (2015). 40 F. L. Hsiao, A. Khelif, H. Moubchir, A. Choujaa, C. C. Chen, and V. Laude, Phys. Rev. E 76, (2007). 41 Y. Pennec, B. Djafari-Rouhani, H. Larabi, J. O. Vasseur, and A. C. Hladky-Hennion, Phys. Rev. B 78, (2008). 42 Y. Pennec, B. Djafari-Rouhani, H. Larabi, A. Akjouj, J. N. Gillet, J. O. Vasseur, and G. Thabet, Phys. Rev. B 80, (2009). 43 T. C. Wu, T. T. Wu, and J. C. Hsu, Phys. Rev. B 79, (2009). 44 M. Oudich, M. B. Assouar, and Z. Hou, Appl. Phys. Lett. 97, (2010). 45 F. Casadei, T. Delpero, A. Bergamini, P. Ermanni, and M. Ruzzene, J. Appl. Phys. 112, (2012). 46 D. Torrent, D. Mayou, and J. Sánchez-Dehesa, Phys. Rev. B 87, (2013). 47 X. P. Wang, P. Jiang, T. N. Chen, and J. Zhu, AIP Adv. 5, (2015). 48 P. Jiang, X. P. Wang, T. N. Chen, and J. Zhu, J. Appl. Phys. 117, (2015). 49 M. Ghasemi Baboly, A. Raza, J. Brady, C. M. Reinke, Z. C. Leseman, and I. El-Kady, Appl. Phys. Lett. 109, (2016). 50 Y. Chen, G. Hu, and G. Huang, J. Mech. Phys. Solids 105, 179 (2017). 51 F. Fahy and P. Gardonio, Sound and Structural Vibration: Radiation, Transmission and Response (Academic Press, Oxford, 2007) Second Edition. APPENDIX A: COMPARISON BETWEEN THE LUMPED MASS MODEL (PWE) AND FULL SCALE MODEL (FEM) We corroborate the results obtained earlier based on the lumped mass model by now using the FEM of the unit cell with all geometrical features described in Fig. 1b. As shown in Fig. 6, we clearly see an excellent agreement of the PWE method (red curves) with the FEM (gray circles). It should be also noted that the p and d type degenerate mode shapes obtained through the lumped mass model comply with those obtained by the full scale model (compare the inset images between Fig. 2 and Fig. 6). In the dispersion relation obtained by the FEM, however, we observe that previously neglected shear-horizontal (SH 0) and shear (S 0) plate modes do appear for a thin plate. In our analysis, it is reasonable to ignore such modes and focus solely on anti-symmetric (A 0) flexural modes for the transverse source excitation at low frequencies. We also see small hybridization band gaps 35 due to the bending of local resonators, which was not accounted for in the lumped mass model. Nevertheless, the band-inversion process around the Dirac frequency is still intact, and it was quickly verified through our lumped mass model.

10 10 FIG. 6: Comparing unit cell dispersion curves obtained from the lumped mass model and the full scale FEM for the cases: (a) R = 0.8a/3, (b) R = a/3, and (c) R = 1.1a/3. Excellent match is observed and band-inversion mechanism is confirmed for flexural modes (A 0). Additional modes emerging in the numerical simulation are due to the bending of resonators, and SH 0 and S 0 guided modes of the plate. APPENDIX B: GRADED INTERFACE BETWEEN TWO TOPOLOGICALLY DISTINCT LATTICES Here we verify the scheme of reducing the gap observed at the Γ point in Fig. 4b. The gap emerges due to the breakage of the C 6 symmetry at the interface, therefore, we minimize the effect of symmetry breakage by modifying the interface. Figure 7a shows a supercell, in which, a topologically non-trivial lattice (R = 1.1a/3) smoothly transitions to a topologically trivial lattice (R = 0.8a/3) via four lattices with a gradient in their radii, i.e., R 1 = 1.05a/3, R 2 = 1.02a/3, R 3 = 0.95a/3, and R 4 = 0.9a/3. The gap reduction is confirmed in Fig. 7b. It is about 3 times lesser than the one observed in Fig. 4b. Spin nature of these interface modes is also intact as shown in Fig. 7c.

11 FIG. 7: The effect of a graded interface on the pseudo-spin modes. (a) A supercell with a smooth interface (radii R 1 > R 2 > R 3 > R 4) between topologically trivial and non-trivial lattices. (b) Eigenfrequencies of the modified supercell. The gap reduction at the Γ point can be seen. (c) Pseudo-spin modes corresponding to S1 and S2 points in the dispersion curve. 11