# Valley-Hall In-Plane Edge States as Building Blocks for Elastodynamic Logic Circuits

###### Abstract

In this work, we investigate theoretically and demonstrate experimentally the existence of valley-Hall edge states in the in-plane dynamics of honeycomb lattices with bi-valued strut thickness. We exploit these states to achieve non-trivial waveguiding of optical modes that is immune to backscattering from sharp corners. We also present how different types of interfaces can be combined into multi-branch junctions to form complex waveguide paths and realize a variety of structural logic designs with unconventional wave transport capabilities. We illustrate this potential with two applications. The first is a direction-selective de-multiplexing waveguide tree featuring a pronounced asymmetric wave transport behavior. The second is an internal waveguide loop along which the energy can be temporarily trapped and periodically released, effectively working as a signal delayer. We provide a rationale for the observed phenomena that blends topological considerations and mechanistic arguments, and we offer a criterion for the proper selection of the junction characteristics that are conducive to non-trivial interface modes. The modal complexity of in-plane elasticity has important consequences on the regime of manifestation of the edge states, as the availability of viable total bandgaps is shifted to higher frequencies compared to the out-of-plane counterpart problem. It also poses additional experimental challenges, associated with proper acquisition and deciphering of the in-plane modes, the solution of which requires a systematic use of in-plane laser vibrometry.

## I I. Introduction

Acousto-elastic metamaterials and phononic crystals are architected structures that are artificially designed to control and manipulate mechanical wave propagation. Some of their most well known properties are mechanical filtering and wave directionality Kushwaha et al. (1993); Liu et al. (2000); Ruzzene et al. (2003); Gonella and Ruzzene (2008), waveguiding (conventional and sub-wavelength) Khelif et al. (2004); Celli and Gonella (2015) and energy trapping Jiao and Gonella (2018); Zhao et al. (2015); Colombi et al. (2014). In recent years, topological insulators (TI) endowed with topologically-protected edge states (TPES) have paved the way to new strategies for wave manipulation in quantum and electronic systems Hasan and Kane (2010); Moore (2010); Qi and Zhang (2010); Chang et al. (2013); Hsieh et al. (2008); Kane and Mele (2005). Inspired by these condensed matter problems, researchers have been able to adapt the working principle of TI to acoustic and elastic systems to achieve analogous wave control capabilities in both the static Kane and Lubensky (2014); Paulose et al. (2015); Rocklin et al. (2017, 2016); Bilal et al. (2017) and dynamic regime Ma et al. (2018); Nash et al. (2015); Souslov et al. (2017); Prodan et al. (2017).

One of the first quantum phenomena to inspire a mechanical analog was the quantum Hall effect (QHE), which allows one-way, non-reciprocal wave propagation by creating chiral edge states that are robust against defects and disorders Nash et al. (2015); Souslov et al. (2017). The QHE can be realized by breaking time-reversal symmetry through the application of an external (e.g. magnetic) field Zhang et al. (2005) that results in the lifting of a degenerate Dirac cone in the band diagram of lattice and the opening of a non-trivial bandgap. In a corresponding mechanical system, it is possible to resort to external active controls, such as spinning rotors Nash et al. (2015) or circulating fluids Souslov et al. (2017), to achieve analogous unidirectionally propagating edge states. However, unlike in electronic systems, the complexity required to apply external mechanical fields partially hinders the practicality of QHE for applications at the device level. In light of these limitations, researchers started exploiting the mechanical analog of the quantum spin Hall effect (QSHE), for which time-reversal symmetry breaking is not required. The QSHE can be realized by manipulating sublattice configurations to create mechanical pseudo-spins and pseudo-spin-dependent effective edges, ultimately achieving topologically protected helical edge modes that are also immune to back-scattering Süsstrunk and Huber (2015); Mousavi et al. (2015); He et al. (2016). While enjoying the inherent advantages of being fully passive, QSHE requires carefully engineered configurations in order to obtain the doubly-degenerate Dirac cones whose lifting is essential to establish helical edge wave conditions. As a result, its practical applicability is also met with difficulties. Mechanical systems involving topologically-protected floppy edge modes Kane and Lubensky (2014); Paulose et al. (2015); Rocklin et al. (2017, 2016); Ma et al. (2018) have also been proposed to obtain asymmetric wave propagation. Although these systems can be easily manufactured and indeed achieve asymmetric wave transport without requiring any active control, the asymmetric behavior only occurs at frequencies lower than the folding region of the acoustic modes Ma et al. (2018).

The most recent efforts in this line of developments have been directed towards creating mechanical analogs of the quantum valley Hall effect (QVHE), whose only requirement is the breaking of space-inversion symmetry (SIS), which can be achieved in a relatively agile fashion in many acoustic and mechanical systems. In acoustics, this effect has been demonstrated for sonic crystals obtained from arrays of rod-like scatterers Lu et al. (2016, 2017). In elastodynamics, the concept was first realized in honeycombs with masses added at lattice nodes Pal and Ruzzene (2017); Vila et al. (2017), and then extended to lattices with curved struts Liu and Semperlotti (2018, 2019) and thin plates endowed with arrays of surface masses Chen et al. (2017); Chaunsali et al. (2018); Zhu et al. (2018). In hexagonal lattices with a symmetry, Dirac points, or valleys are observed at the high-symmetry corners of the Brillouin zone. SIS breaking, which preserves symmetry, results in bandgap opening and in the lifting of the cones. This leaves behind two adjacent valleys, and , that are largely separated in the reciprocal space and feature opposite values of the locally-defined topological index known as valley Chern number. The Chern number dichotomy between the valleys implies low modal compatibility between phonons associated with the and points (which can be interpreted as waves traveling in opposite directions in the lattice), which ultimately manifests as a low degree of back-scattering. A complete and lucid description of the working principle of QVHE is given in Ref. Liu and Semperlotti (2019). The SIS breaking requirement is relatively simple to achieve through a mechanical modulation of the unit cell parameters, making this approach much more amenable to manufacturing and therefore applicable to elastic systems. Most efforts on QVHE to date (with the notable exception of one theoretical work on in-plane elasticity in discrete systems Chen et al. (2018)) have been devoted to the analysis of the out-of-plane flexural (antisymmetric Lamb wave) response of thin elastic plates.

In this work, we shift our attention to the in-plane dynamics of thin lattices in plane-stress conditions. Realizing viable QVHE conditions in in-plane systems features some additional challenges compared to their out-plane counterparts. In out-of-plane problems, a total bandgap can be easily obtained at low frequencies Pal and Ruzzene (2017); Vila et al. (2017); Lu et al. (2016, 2017). This allows a straightforward establishment of edge waves that are not contaminated by any of the bulk modes (as a side note, it is worth noting that edge modes have also been successfully achieved in systems and frequency regimes where SIS breaking only induces partial bandgap opening Liu and Semperlotti (2018); Zhu et al. (2018); this result can be explained by invoking the low modal compatibility and the large difference in density of states between the edge and bulk modes at the frequencies of interests, which makes the establishment of edge modes more favorable). In contrast, in-plane mechanics are characterized by additional modal complexity associated with the co-existence of longitudinal and shear modes even at low frequencies. As a result, identifying total bandgaps that naturally allow for spectrally-isolated edge modes inevitably requires exploring the high-frequency optical regime.

Our structure of choice is a modified hexagonal lattice consisting of beam-like struts with bi-valued thickness. This configuration is characterized by extreme geometric simplicity, which makes it easy to manufacture using conventional cutting or printing techniques. While this structure features multiple bandgaps that are amenable to QVHE, here we focus our attention on the high-frequency regime, where a wide total bandgap is available. Under these conditions, we study numerically and demonstrate experimentally the establishment of non-trivial waveguides and their robustness against back-scattering from sharp corners. Moreover, we leverage this platform to realize complicated paths involving multi-branch junctions and different types of domain walls to achieve a variety of special energy channeling and routing effects. As a first example, we design a direction-selective de-multiplexing waveguide tree that allows one-way wave propagation at high frequencies. We also design a triangular waveguide loop that works as a signal delayer: here a portion of the energy injected in the system is temporarily trapped in the loop, and then released periodically at an output boundary with programmable delays.

## Ii II. Topologically-Protected Edge Modes

To obtain a lattice with a symmetry after breaking the SIS, we start with a regular hexagonal lattice with a symmetry made of 2-cm-long struts with a 7.5:1 in-plane aspect ratio, Fig. 1 (a). The geometric and material properties are chosen to match those of the specimen that we intend to test experimentally, which is manufactured via water jet cutting from a sheet of acrylonitrile butadiene styrene (ABS): length of the struts cm, out-of-plane sheet thickness mm, Young’s modulus, GPa, Poisson’s ratio , density . The phonon dispersion relation, obtained via Bloch unit cell analysis conducted via finite element analysis (FEA), is shown in Fig. 1 (b). The struts are modeled as shear-deformable beams according to Timoshenko beam theory to expedite the calculation with respect to a 2D elasticity model, without compromising the accuracy. In the frequency range of interest (below Hz), dictated by the bandwidth of our actuation equipment, we identify two sets of candidate Dirac-like points at the point. To confirm that these features are indeed Dirac points, we calculate the Berry phase (or Berry flux) around the point:

(1) |

where denotes a close loop around the point and is the Berry connection, which can be obtained from the eigenvectors of a phonon band . For a 2D system that is invariant under time-reversal and two-fold rotation (about ), the Berry flux is generally quantized to or (up to mod ), with being the key signature of a Dirac point.

Numerically, an efficient way to compute the Berry phase is via eigenvector overlap. Let us select wavevector points along a closed loop, labeled , which encapsulate the point of interest and such that the periodic boundary condition is satisfied. The inner product between eigenvectors at two neighboring points can be computed as . When is large and the wavevector points are densely populated along the loop, this inner product is in general a complex number with an absolute value close to unity. Multiplying the inner products computed between all the available pairs yields a complex number whose phase is the Berry phase:

(2) |

As mentioned above, since our system preserves time-reversal and point group symmetry, is quantized to or , and thus the product of eigenvectors appearing on the right-hand-side of Eqn. (2) is a real number. It is then sufficient to look at the sign of this product: if it is positive (negative), the corresponding Berry flux is (), and thus a negative eigenvector overlap product is smoking-gun evidence for the existence of a Dirac point.

To break SIS, we modulate the strut thickness so that the aspect ratio of the upper (lower) half of the unit cell is increased (decreased) by 20%, as shown in Fig. 1 (c). The asymmetry of the lattice opens the Dirac cones and creates band gaps of Hz and Hz for the upper and lower cones, respectively, see Fig. 1 (d). We decide to focus our analysis on the upper gap, which is wider and more tractable. To describe the topology of this gap, we calculate the valley Chern number :

(3) |

where is an arbitrary point in the reciprocal space located in the neighborhood of the valley point, and is the Berry curvature calculated from the eigenvectors of the upper or lower bands near the band gap, which is expressed as:

(4) |

where is the normalized eigenvector, and is the unit cell mass matrix. A color map of the Berry curvature for the upper band is shown in Fig. 2. An analogous Berry curvature map, albeit with opposite signs, can be obtained for the lower band. The numerically calculated Chern number for the valley ( ) point is , approximately tending to of classical QVHE theory.

We now proceed to consider a lattice consisting of two domains (termed lattice type A and B) separated by an interface. The cells above the interface have thin (thick) struts in the top half of the cell and thick (thin) struts in the bottom half, while those below the interface are mirror symmetric, as shown in Fig. 3 (d) [(e)]. To confirm the existence of TPES and to explicitly reveal their paths and morphological attributes in the context of a finite lattice domain, we conduct a analysis. Each supercell consists of a finite strip of 12 cells: the top six cells belong to type A (B), while the bottom to type B (A), as shown in Fig. 3 (a) [(b)]. In this way, the lattice resulting from the repetition of the supercell along the -axis displays a horizontal interface along the zigzag direction. The asymmetry of the unit cell allows two possible types of interface: we name an interface where two thin struts meet a interface, Fig. 3 (a), and the other type interface, Fig. 3 (b). We apply 1D Bloch conditions along the -axis while leaving the top and bottom boundaries free to mimic a finite-width, infinitely long horizontal strip. From the supercell band diagrams obtained for the two different interfaces, superimposed in Fig. 3 (c), we detect the emergence of two new bands within the bandgap. The mode shapes calculated at the valley point () on these new bands, shown in Fig. 3 (d) and (e), respectively, indicate that the emerging bands are indeed interface modes, as the deformation is mostly localized at the interface and decays exponentially into bulk. Interestingly, the supercell mode shape for a thin (thick) interface is antisymmetric (symmetric) about the interface. An antisymmetric (symmetric) point excitation is hence required for a thin (thick) interface to allow wave propagation.

## Iii III. topological protection at multi-branch junctions

An important feature of a TPES is that it allows waves to propagate along sharp corners with negligible back-scattering. While this phenomenon has been effectively explained in the literature for flexural waves, here we attempt to shed further light on the genesis of this phenomenon by providing an intuitive mechanistic explanation. To this end, we need to compute and plot the eigenvector phase rotation for a TPES. Let us note that TPES is, for all intents and purposes, a propagating wave concentrated near the interface boundary. Like any traveling wave, its displacement field is fully characterized in terms of an amplitude and a phase, i.e.,

(5) |

where is the amplitude and is the phase, with and being the real space coordinate and the wavevector, respectively. For a TPES, the amplitude function is concentrated along the interface and decreases exponentially as we move away from the interface. The phase pattern of a TPES contains important information which directly dictates the wave propagation behavior at corners or junctions of interface boundaries.

In Fig. 4, we plot the phase for each node in the FEA model and we monitor the values at the nodes located along the interface defined in Fig. 3 (a). Here, we focus our attention on a TPES with wavevector , i.e., the valley. As previously shown in Fig. 3 (c), this mode has a positive group velocity and thus corresponds to a mode that transports energy from left to right, as denoted by the green arrow in Fig. 4. From the colormap in Fig. 4, as can be seen, along the horizontal interface, the phase angle returns to the same value after every three unit cells following a high-zero-low sequence, which we term “clockwise” (the opposite low-zero-high sequence will be called “counter-clockwise”). The 3-cell periodicity and clockwise pattern are obtained for the selected wavenumber , while the TPES at (not shown) has the same periodicity but counter-clockwise pattern. Now, let us consider other lattice directions along which possible domain interfaces can be realized. Along R1 and R2, which form a angle with the horizontal interface, we observe that the phase also propagates clockwise. In other words, the phase pattern for a right-propagating TPES is fully consistent with TPES modes propagating along R1 and R2.

Phase consistency is a crucial criterion for wave propagation. In general, wave propagation tends to maximize phase consistency, i.e., a wave will usually choose the propagation direction such that the in-coming and out-going waves have the highest phase coherence and therefore the most efficient energy transfer. For example, in a 3D homogeneous medium, the phase pattern of a plane wave is best matched by itself, and thus a wave will choose to maintain its wave-vector such that perfect phase match is preserved. However, if the wave impinges on a 2D interface, because phase coherence cannot be maintained in the whole 3D space anymore, the next best option is selected, i.e., to maintain the phase pattern along the 2D interface, leading to the mechanisms behind the well-known reflection and refraction patterns.

We can apply the same principle to the TPES. First, we consider a multi-segment interface which switches direction from the horizontal to either R1 or R2. At a corner, a wave packet propagating from left to right along the -axis has two options: (a) being reflected towards the left or (b) propagation along R1 or R2. The phase pattern is fully preserved (clockwise in this case) for option (b), while in option (a), the phase will have to change its pattern to counter-clockwise, failing to match that of the incoming wave. Therefore, the wave will prefer traveling around the sharp corner, rather than being scattered back, which explains why a TPES can transmit through such a corner with nearly zero back-scattering.

Let us consider now the propagation path that switches from the horizontal interface to R3 or R4, which form a angle with the horizontal direction, Note that for a 120 corner, the segment of R3 or R4 must involve the other interface type, otherwise the lattice cannot be properly completed. Since, in this case, the landscape of phase rotation is reversed. Naively, we could conclude that R3 and R4 are unfavorable directions. However, it is important to note that the TPES along this new interface features a negative group velocity [see the black curve in Fig. 3 (c)], opposite to that along the horizontal interface. Thus, for a wave packet propagating along R3 or R4, the phase pattern is also clockwise along the group velocity direction, consistent with the scenario of the incoming wave. As a result, waves can still propagate along such a angle, rather than being scattered back, provided that we switch to a different type of interface.

To document the robustness against back-scattering at sharp corners, we model a Z-shape thin interface with two corners, shown in Fig. 5 (a). The interface splits the lattice into two sub-domains with mirror symmetry about the interface. The dynamic response is computed via full-scale dynamic FEA. A burst excitation with frequency Hz with 35 cycles is applied in the form of a horizontal in-plane force at left edge of the interface. As is evident from Fig. 5 (b), the wave indeed propagates along the interface without significant back-scattering until it reaches the right boundary. The back-scattering-immune behavior can also be confirmed by plotting the displacement time history of an arbitrary point located along the interface, as shown in Fig. 5 (g). The presence of a single dominant peak, associated with the incoming wave, indicates that no appreciable reflections are generated at the corners. Since the excitation frequency is within the bulk band gap, no leakage into the bulk is observed.

An alternative strategy to realize a similarly-shaped waveguide is by constructing a trivial interface. This can be done starting from a uniform lattice (say type A) and replacing the cells featuring thickness modulation with regular hexagonal cells along the interface. Let us consider two trivial Z-shaped domain walls with interface widths spanning over one and two cells, as shown in Fig. 5 (c) and (e), respectively. Interestingly, when the trivial interface is one cell in width, as in Fig. 5 (c), the wave remains mostly localized around the excitation point, with little energy propagating along the interface, as shown in Fig. 5 (d). Therefore, a trivial interface with one-cell thickness is not sufficient to effectively close the band gap and establish efficient waveguide transport conditions. By sufficiently widening the interface, as in the case of Fig. 5 (e), involving a two-cell interface, waves can eventually propagate along the trivial waveguide, as shown in Fig. 5 (f). However, strong back-scattering is observed at corners, significantly weakening the signal transmitted and rendering this a very inefficient waveguide. This is best seen by looking at Fig. 5 (h), which again shows the time history of the response sampled at a point along the first segment of the waveguide. A second peak at around s is clear evidence of a strong reflected wave. In conclusion, only the topologically-protected interface guarantees efficient and scattering-free waveguiding.

## Iv IV. Experimental Characterization

We now proceed to perform a series of experiments to confirm the topological attributes of the non-trivial waveguide. To acquire the in-plane wavefields, we scan the lattice with a 3D Scanning Laser Doppler Vibrometer (SLDV, Polytec PSV-400-3D) using one scan point per lattice node. The specimen is a water-jet-cut realization of the aforementioned ABS lattice, comprising unit cells with the same dimensions and parameters used in simulations, see Fig. 6 (a). The green-shaded area denotes the Z-shape interface. For convenience of setup, we rotate the lattice by counter-clockwise and frame it using clamps at the two sides, leaving the top and bottom boundaries free. The excitation is prescribed as an in-plane normal point force applied on the boundary at the opening of the topological interface. The force is exerted using an electrodynamic shaker (Bruel & Kjaer Type 4809, powered by a Bruel & Kjaer Type 2718 amplifier) placed at the bottom of the lattice, as shown in Fig. 6 (b). We prescribe a burst excitation at Hz with 35 cycles in order to activate a very narrow frequency interval falling entirely on the topological edge mode branch within the bandgap, approximately at . From visual inspection of the wavefield, we observe that the wave indeed propagates along the interface without appreciable backscattering, see Fig. 6 (c), consistent with the numerical prediction in Fig. 5 (b). Note that the magnitude is here attenuated along the propagation path due to unavoidable damping in the physical specimen. The robustness against back-scattering is eloquently captured by the spectral plane representation of Fig. 6 (d), obtained by applying a discrete Fourier transform (DFT) to the experimental data sampled along the interface. Indeed, the spectral magnitude map presents a strong asymmetry between the and directions, suggesting negligible reflection from the corner. It can also be appreciated that the spectral response features excellent agreement with the edge mode branch obtained from supercell analysis.

## V V. Multi-interface Junctions for Energy Transport Manipulations

The discussion in Section II offered a rationale for how waves propagate at and corners with or without changes of interface type. These considerations can inspire a plethora of structural logic circuit designs with unconventional energy transport capabilities. Here, we explore the characteristics and opportunities of multi-interface junctions. According to our previous discussion, the lattice types on the two sides of a topologically-protected interface must present a mirror symmetry. Hence, to ensure that all the interfaces converging at a junction are topologically protected, it is necessary that the junction connects an even number of interfaces. Here we illustrate this idea through two examples.

The first design is a direction-selective elastic de-multiplexer intended to allow strongly asymmetric wave transport at high frequencies. The FEA model and experimental specimen are shown in Fig. 7 (a) and (b), respectively. Three “thin” and one “thick” interfaces (according to the terminology introduced before) are connected at point , resulting in three ports on the left side of the lattice and one on the right. The angle between the three thin interfaces is , while the one between the thin and thick interfaces is . A burst excitation with carrier frequency falling inside the bandgap is applied at port on the left, and the simulated and experimentally measured wavefields are shown in Fig. 7 (c) and (d), respectively. It can be seen that the wave prefers to split and steer at the junction and propagate back along and , eventually localizing at the left edge of the lattice without ever reaching the opposite edge. In contrast, the simulated and experimental wavefields obtained exciting at port on the right are shown in Fig. 7 (e) and (f), respectively. As we can see, the wave propagates along the thick interface from right to left and, at junction , splits (with negligible reflections) into two transmitted signals traveling to ports and . The characteristics of the interfaces coalescing at junction ensures that the wave propagates along the two thin interfaces and that form a angle with the thick interface . As a result, transmission form right to left is achieved efficiently. Note that opposite group velocities between the thin and the thick interfaces prevent the wave from propagating straight through the junction, despite the perfect alignment between and .

For our second example, we arrange three thin and one thick interfaces as shown in Fig. 8 (a) to realize a triangular loop embedded in the lattice domain. The idea is to create an energy sink that locally traps the energy injected into the lattice from an edge and delays its release at another edge. When a wave traveling from port impinges onto junction , the energy splits, with a portion propagating along the upper thin interface (which forms a angle with ), and the other following the lower thick interface at a angle. The energy propagating upwards circulates clockwise in the triangular loop, eventually coming back to junction . Here the energy splits again into a packet that enters the loop for a second cycle and a second packet that travels along the thick interface and eventually reaches the lower egde after some delay, see Fig. 8 (c). This process, which is repeated at every cycle, breaks down the burst of energy and releases it at port after finite intervals, with a period that depends on the geometric characteristics of the loop. As a result, in contrast with the input signal, which consists of a single packet, the output signal is broken down into multiple packets that are transmitted over an extended time window, as shown in Fig. 8 (b).

## Vi Conclusion

In this work, we have theoretically analyzed and experimentally demonstrated the in-plane dynamics of topologically-protected interface modes in honeycomb lattices with a bi-valued strut thickness pattern that breaks the space-inversion symmetry. We have documented similarities and differences between topological interfaces formed by connecting thin and thick portions of the unit cell. We have shown the robustness against back-scattering of wave propagation about corners with identical types of interfaces and about with mixed interfaces. Exploiting this knowledge, we have designed multi-interface junctions and used them to realize different structural logic circuits and elastic ports featuring unconventional wave manipulation effects at optical frequencies. Our results extend the applicability of valley-Hall effects to in-plane elastic dynamics and enrich the current understanding of these phenomena in mechanical systems by providing a new interpretation angle that blends topological and mechanistic arguments. The presented modular interface design philosophy can provide inspiration for a broad range of applications for energy management, trapping, routing and harvesting.

## Vii acknowledgment

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

## References

- Kushwaha et al. (1993) M. S. Kushwaha, P. Halevi, L. Dobrzynski, and B. Djafari-Rouhani, Physical Review Letters 71, 2022 (1993).
- Liu et al. (2000) Z. Liu, X. Zhang, Y. Mao, Y. Zhu, Z. Yang, C. T. Chan, and P. Sheng, Science 289, 1734 (2000).
- Ruzzene et al. (2003) M. Ruzzene, F. Scarpa, and F. Soranna, Smart Materials and Structures 12, 363 (2003).
- Gonella and Ruzzene (2008) S. Gonella and M. Ruzzene, Journal of Sound and Vibration 312, 125 (2008).
- Khelif et al. (2004) A. Khelif, A. Choujaa, S. Benchabane, B. Djafari-Rouhani, and V. Laude, Applied Physics Letters 84, 4400 (2004).
- Celli and Gonella (2015) P. Celli and S. Gonella, Applied Physics Letters 107, 081901 (2015).
- Jiao and Gonella (2018) W. Jiao and S. Gonella, Physical Review Applied 10, 024006 (2018).
- Zhao et al. (2015) D.-G. Zhao, Y. Li, and X.-F. Zhu, Scientific Reports 5, 9376 (2015).
- Colombi et al. (2014) A. Colombi, P. Roux, and M. Rupin, The Journal of the Acoustical Society of America 136, EL192 (2014).
- Hasan and Kane (2010) M. Z. Hasan and C. L. Kane, Reviews of Modern Physics 82, 3045 (2010).
- Moore (2010) J. E. Moore, Nature 464, 194 (2010).
- Qi and Zhang (2010) X.-L. Qi and S.-C. Zhang, Physics Today 63, 33 (2010).
- Chang et al. (2013) C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, M. Guo, K. Li, Y. Ou, P. Wei, L.-L. Wang, et al., Science 340, 167 (2013).
- Hsieh et al. (2008) D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nature 452, 970 (2008).
- Kane and Mele (2005) C. L. Kane and E. J. Mele, Physical Review Letters 95, 146802 (2005).
- Kane and Lubensky (2014) C. Kane and T. Lubensky, Nature Physics 10, 39 (2014).
- Paulose et al. (2015) J. Paulose, B. G.-g. Chen, and V. Vitelli, Nature Physics 11, 153 (2015).
- 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).
- Bilal et al. (2017) O. R. Bilal, R. Süsstrunk, C. Daraio, and S. D. Huber, Advanced Materials 29, 1700540 (2017).
- Ma et al. (2018) J. Ma, D. Zhou, K. Sun, X. Mao, and S. Gonella, Physical Review Letters 121, 094301 (2018).
- 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).
- Souslov et al. (2017) A. Souslov, B. C. van Zuiden, D. Bartolo, and V. Vitelli, Nature Physics 13, 1091 (2017).
- Prodan et al. (2017) E. Prodan, K. Dobiszewski, A. Kanwal, J. Palmieri, and C. Prodan, Nature Communications 8, 14587 (2017).
- Zhang et al. (2005) Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 438, 201 (2005).
- Süsstrunk and Huber (2015) R. Süsstrunk and S. D. Huber, Science 349, 47 (2015).
- Mousavi et al. (2015) S. H. Mousavi, A. B. Khanikaev, and Z. Wang, Nature Communications 6, 8682 (2015).
- He et al. (2016) C. He, X. Ni, H. Ge, X.-C. Sun, Y.-B. Chen, M.-H. Lu, X.-P. Liu, and Y.-F. Chen, Nature Physics 12, 1124 (2016).
- Lu et al. (2016) J. Lu, C. Qiu, M. Ke, and Z. Liu, Physical Review Letters 116, 093901 (2016).
- Lu et al. (2017) J. Lu, C. Qiu, L. Ye, X. Fan, M. Ke, F. Zhang, and Z. Liu, Nature Physics 13, 369 (2017).
- Pal and Ruzzene (2017) R. K. Pal and M. Ruzzene, New Journal of Physics 19, 025001 (2017).
- Vila et al. (2017) J. Vila, R. K. Pal, and M. Ruzzene, Physical Review B 96, 134307 (2017).
- Liu and Semperlotti (2018) T.-W. Liu and F. Semperlotti, Physical Review Applied 9, 014001 (2018).
- Liu and Semperlotti (2019) T.-W. Liu and F. Semperlotti, Physical Review Applied 11, 014040 (2019).
- Chen et al. (2017) J.-J. Chen, S.-Y. Huo, Z.-G. Geng, H.-B. Huang, and X.-F. Zhu, AIP Advances 7, 115215 (2017).
- Chaunsali et al. (2018) R. Chaunsali, C.-W. Chen, and J. Yang, New Journal of Physics 20, 113036 (2018).
- Zhu et al. (2018) H. Zhu, T.-W. Liu, and F. Semperlotti, Physical Review B 97, 174301 (2018).
- Chen et al. (2018) H. Chen, H. Nassar, and G. Huang, arXiv preprint arXiv:1802.04404 (2018).