Chinese Physics B, Volume. 29, Issue 10, (2020)

Symmetry-controlled edge states in graphene-like topological sonic crystal

Zhang-Zhao Yang1, Jin-Heng Chen1, Yao-Yin Peng1, and Xin-Ye Zou1,2、†
Author Affiliations
  • 1Key Laboratory of Modern Acoustics (MOE), Institute of Acoustics, Department of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 20093, China
  • 2State Key Laboratory of Acoustics, Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China
  • show less

    Unique topological states emerged in various topological insulators (TI) have been proved with great application value for robust wave regulation. In this work, we demonstrate the parity inversion related to the definition of the primitive cell in one common lattice, and realize a type of symmetry-controlled edge states confined on the zigzag interfaces of the graphene-like sonic topological crystal. By simply sliding the selected ‘layer’ near the interface, the coupling of the pseudospin states induced by the multiple scattering for the C6v lattice results in the adjustment of the edge states. Based on the physics of the states, we experimentally propose a prototype of acoustic topological filter hosting multiple channels with independent adjustable edge states and realize the selective high transmission. Our work diversifies the prospects for the applications of the gapped edge states in the robust wave regulation, and proposes a frame to design new topological devices.

    Keywords

    1. Introduction

    The concept of topological insulator (TI) has drawn extensive attention since it was proposed. In condensed matter physics, TI has been proved to be an excellent platform for new types of lossless information transmission because of its defects-immune property based on gapless topological-protected edge states.[15] In recent years, optical TIs[616] and acoustic TIs[1726] have also been successfully demonstrated. More recently, the concept of higher-order topological insulator (HOTI) hosting gapped edge states and corner states has shown that the traditional bulk–boundary correspondence is no longer applicable to it.[2734] Different from the topological phases of the traditional TIs, the reported works have shown that the topologies of the HOTIs are characterized by the nonzero bulk polarization, which are strongly related to the point group symmetry (PGS) of the periodic structures.[27,3436] Moreover, some works have demonstrated that HOTIs can also be implemented in two-dimensional elastic waves,[37,38] microwaves,[39] toroidal circuits,[40] optical[4144] and acoustic systems.[4549]

    The previous researches on the HOTIs have always focus on the 0D corner states which are capable of confining the energy on the corners exactly, however, we note that the normally-gapped edge states of the HOTIs also provide unprecedented platforms for the control of wave propagation. In the classical wave systems, the behaviors of the topological states crucially depend on the structure parameters of the material, which implies the possibilities to adjust the bands of the gapped edge states.

    In this work, we propose a graphene-like sonic HOTI and focus on the behaviors of the edge states. We show that the parity inversion of the proposed HOTI depends on the choices of the primitive cells, and the inverted parities of the eigenstates at the high symmetry points near the nontrivial gaps, which indicate nonzero bulk polarization, lead to the generation of the edge states on the interface. Meanwhile, the pseudo-spins carried by the states originating from the graphene-like lattice ensure the robust wave propagation. Importantly, we further demonstrate that the translation of the ‘layer’ leads to the changing of the coupling strength of the multipolar states between the layers, and then results in the expansion or contraction of the gapped edge states. It indicates that the edge states can be adjusted by simply controlling the translation rates. Based on these, a prototype of two-dimensional acoustic topological filter with adjustable edge states is experimentally proposed. By simply sliding the sections of the filter, the bandwidths of the edge states can be controlled, and the propagation of the sound waves at specific frequencies band to the selected ports can be freely manipulated.

    2. The creation of parity inversion

    We start with a regular 2D SC with a C6v symmetry lattice consisting of solid hexagonal columns embedded in the air background as shown in Fig. 1(a). The side length of the column is r0 = 6 mm, and the distance between the nearest-neighbor columns is R = 2.778r0. The lattice constant is a=3R. In the perspective of HOTI, the nontrivial topology is determined by the bulk polarization, which in turn results in the emergence of the edge polarization (Px,Py). A nontrivial edge polarization indicates the existence of the in-gap edge states located on the boundaries of the nontrivial domains, and the robust 0D topological corner state can also emerge on the specific corner as the intersected boundaries are all polarized, which are different from the ordinary defects.[5052] Meanwhile, the topology of the material is strongly related to the point group symmetry, which means that the parity inversion indicating the topological phase transition may be induced simply by selecting different definition of the unit cell in one common lattice.[53] In our work, the solid columns can be treated as hard boundries and the energy is actually propagating in the background air, which is different from the previous works.[24,49] To simplify the model of our work, we only consider two nearest-neighbor hoppings of the aerial meta-atoms respectively along x and y directions, and the structure can then be redefined with two types of rectangular supercells [red dashed lines in Fig. 1(a)]. Figures 1(b) and 1(c) depict the energy band structures of the two supercells, respectively. Due to the same lattice constants ax = a, ay=3a, and the same periods along both directions, the energy band structures of the two supercells shall also be the same. However, the parity of the modal distribution of each band labeled with “+” (even parity) and “–” (odd parity) at the high symmetry points represents distinct topology.

    (a) Schematic of two selected supercells (red dashed lines) based on one common graphene-like lattice, the corners and the edges of the supercells represent the pseudo-atoms and the corresponding hoppings, respectively. The corresponding energy band structures labeled with parities are presented in (b) and (c). (d) and (e) The pressure field distributions of the first two double-degenerate bands from bottom to top at M point of the two supercells, respectively. The opposite parities indicate the parity inversion.

    Figure 1.(a) Schematic of two selected supercells (red dashed lines) based on one common graphene-like lattice, the corners and the edges of the supercells represent the pseudo-atoms and the corresponding hoppings, respectively. The corresponding energy band structures labeled with parities are presented in (b) and (c). (d) and (e) The pressure field distributions of the first two double-degenerate bands from bottom to top at M point of the two supercells, respectively. The opposite parities indicate the parity inversion.

    The topological non-triviality of the parity inversion can be characterized by the integration of the Berry connection over the 1st BZ[27]

    $$ \begin{eqnarray}\begin{array}{lll} & & {P}_{i}=\displaystyle \frac{1}{{(2\pi)}^{2}}\displaystyle \iint {{\rm{d}}}^{2}k{\rm{Tr}}[{A}_{i}],\\ & & {({A}_{i})}_{mn}(k)={\rm{i}}\langle {u}_{m}(k)|{\partial }_{{k}_{i}}|{u}_{n}(k)\rangle,\end{array}\end{eqnarray}$$ (1)

    where i indicates the component of 2D polarization P along the reciprocal vector bi, and |um(k)〉 is the periodic Bloch function for the m-th band. If applying C2 rotation to this system, equation (1) can be rewritten as[36]

    $$ \begin{eqnarray}\begin{array}{lll} & & {P}_{i}=\displaystyle \frac{1}{2}\left(\displaystyle \sum _{m}{q}_{i}^{m}\,\mathrm{mod}\,2\right),\end{array}\end{eqnarray}$$ (2)

    $$ \begin{eqnarray}\begin{array}{lll} & & {(-1)}^{{q}_{i}^{m}}=\displaystyle \frac{{\eta }_{m}(M)}{{\eta }_{m}(\Gamma)},\end{array}\end{eqnarray}$$ (3)

    where ηm (M) represents the parity at M point for the m-th band, and the summation runs over all occupied bands. According to the parities of the modes of the bands at the corresponding high symmetry points, the band gaps of the supercell 1 carry trivial topology, but the gaps of the supercell 2 carry nontrivial topology. To explain it exactly, we present the first four eigenmodes at M point of the two supercells in Figs. 1(d) and 1(e), respectively. The four parities of supercell 1 from bottom to top are characterized as s mode (+), dxy mode (+), and a pair of p modes (–). Contrary, the corresponding parities of supercell 2 are a pair of p modes (–), s mode (+), and dxy mode (+). The topology characterized by the inverted parities demonstrates the emergence of the higher-order topological states on the interfaces between the trivial and nontrivial domains.

    Before introducing the topological states in this crystal, it is important to note that the nontrivial polarization of the domain does not require either the spin-like degrees of freedom or the time-reversal relation.[53,54] However, the pseudo-spins can still be introduced by multiple scattering for the special C6v lattice.[24] We calculate the real-space distributions of the root-mean-square intensity of the energy flow I=1TRe(p)Re(v)dt to describe the pseudo-spins in the two supercells as illustrated in Figs. 1(d) and 1(e) (see details in Appendix A). Crucially, the chiral flows on each corresponding band of the two supercells imply the promotion of the wave propagation between the domains with distinct topology.

    3. Pseudo-spin-protected adjustable edge states

    The bulk–edge–corner correspondence of the HOTIs predicts both the edge and corner states. In this work, we focus on the behaviors of the gapped edge states. Figures 2(a) and 2(b) present the energy band structures of two ribbons, which are periodic along y direction, composed of 30 trivial supercells or 15 trivial and 15 nontrivial supercells along x direction, respectively. Apparently, the predicted edge states vanish on the edge of the complete trivial domain [Fig. 2(a)], but emerge on the nontrivial interface and span the band gaps [Fig. 2(b)]. Meanwhile, due to the strong coupling of the pseudo-spins on the interface, the edge states appear to be gapless (see Fig. A2 in Appendix A).

    Band structures of the ribbons composed of (a) 30 trivial supercells and (b) 15 trivial and 15 nontrivial supercells, being finite along y direction and periodic along x direction.

    Figure 2.Band structures of the ribbons composed of (a) 30 trivial supercells and (b) 15 trivial and 15 nontrivial supercells, being finite along y direction and periodic along x direction.

    We emphasize that the two primitive cells with distinct topology corresponding to the same lattice, which implies that the band inversion can be induced by simply sliding the parts of the structure. As shown in Fig. 3(a), the structure is divided into layers, which can be slid along x direction independently. We define the translation rate βi = ri / ax (i = A, B, C, …), where ri indicates the actual amount of the translation of the i-th layer. We first study one interface between A-layer and B-layer [inset in Fig. 3(a)], and only consider the slide of A-layer. Figure 3(b) illustrates the correlation between the gap of the edge states and βA. As βA increases, the symmetry of the interface reduces and the coupling of the pseudo-spins strengthens, which results in the expansion of the edge states. When βA = 0.5, the anti-symmetry of the interface corresponds to the strongest coupling of the pseudo-spins, which meets the parity inversion exactly. Therefore, the edge states are steep and different from the flat cases without considering pseudo-spins.[36] Figures 3(c)3(e) present the corresponding band structures for βA being 0.125, 0.25, and 0.375, respectively.

    (a) Schematic of the structure divided into independent layers. Inset: Details of the configuration near the interface. (b) βA-dependent band-edge frequencies for the supercells with k = 0 within a translation period. g1 and g2 represent the ranges between the two topological bands. Energy band structures with distinct edge states for (c) βA = 0.125, (d) 0.25, and (e) 0.375, respectively.

    Figure 3.(a) Schematic of the structure divided into independent layers. Inset: Details of the configuration near the interface. (b) βA-dependent band-edge frequencies for the supercells with k = 0 within a translation period. g1 and g2 represent the ranges between the two topological bands. Energy band structures with distinct edge states for (c) βA = 0.125, (d) 0.25, and (e) 0.375, respectively.

    We further consider three independent layers with βA = 0.5, βC = 0.25, and βB = 0. Figure 4(a) shows the energy band structure of the corresponding ribbon. Significantly, there are two pairs of edge states in each gap, which represents the independent topological propagating modes on the two interfaces, respectively. The real-space distributions and the flows of intensity of the pressure field for k = ± 0.05 × 2 π / ax shown in Fig. 4(b) depict the predicted topological modes confined on the specific interface. These results provide the possibility to design a type of multi-channel topological acoustic waveguide to achieve filtering at specific frequencies.

    (a) Energy band structure of the ribbon with two independent interfaces, the two pairs of the edge states are independently emerging on the specific interface. (b) Pressure field distributions and intensity flows of the pseudo-spin-dependent edge states for k = ± 0.05 × 2 π / ax on the interfaces of the ribbon.

    Figure 4.(a) Energy band structure of the ribbon with two independent interfaces, the two pairs of the edge states are independently emerging on the specific interface. (b) Pressure field distributions and intensity flows of the pseudo-spin-dependent edge states for k = ± 0.05 × 2 π / ax on the interfaces of the ribbon.

    4. Experiments

    Here, we construct a prototype of two-dimensional acoustic topological filter based on the crystal proposed above and experimentally measure the corresponding sound wave transmission. As shown in Fig. 5(a), the filter is composed of 3D printed solid hexagonal columns, which can be treated as the acoustic hard boundaries. The whole structure is placed between two hard-boundary waveguide plates to prevent the possible sound energy from leaking into the free space, and is surrounded by sound absorption boundaries to eliminate the interference of the environmental noise. The filter is divided into the four sections labeled with A, B, C, and D, and the four ports are labeled as input, port 2, port 3, and port 4 corresponding to the four inter-paths labeled with P1, P2, P3, and P4, respectively. A broadband sound source is placed at the input, and three transducers are placed at the other three ports to detect the sound energy transmission. We introduce the translation rates βPi (i = 1, 2, 3, 4) to indicate the translation ratio of the corresponding path. If we set βP2 = βP4 = 0 and βP1 = βP3 = 0.25, the symmetry of the structure at P1 and P3 will be broken and two transport paths will be generated, which means that the sound wave can only propagate in one dimension. Figure 5(b) exhibits the high transmission of the sound wave at 6.75–7.39 kHz to port 3 with very low leakage of the energy into the bulk, which corresponds to the energy band structure shown in Fig. 3(d). If βP1 and βP3 are increased to 0.5, the symmetry is the weakest and the sections on both sides are antisymmetric along the path. Hence, the two edge states in the band gaps connect to each other. Figure 5(c) shows the high-transmission of the sound wave at 5.74–7.39 kHz, which corresponds to the energy band structure shown in Fig. 2(b). These results experimentally demonstrate the adjustable edge states hosted by the filter.

    (a) Experimental setup of the two-dimensional acoustic topological filter. (b) Measured bulk (green) and port 3 (red) transmission spectra for βP1 = βP3 = 0.25. (c) Measured bulk (green) and port 3 (red) transmission spectra for βP1 = βP3 = 0.5. The gray regions in both (b) and (c) represent the theoretic range of the edge states.

    Figure 5.(a) Experimental setup of the two-dimensional acoustic topological filter. (b) Measured bulk (green) and port 3 (red) transmission spectra for βP1 = βP3 = 0.25. (c) Measured bulk (green) and port 3 (red) transmission spectra for βP1 = βP3 = 0.5. The gray regions in both (b) and (c) represent the theoretic range of the edge states.

    We further study the two-dimensional transmission of the sound wave in the filter. Figure 6(a) shows the low transmission of the filter’s three ports for βPi = 0 (i = 1,2,3,4). If we set βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25, the four transport paths will all be opened, but only P1 and P3 are antisymmetric. The corresponding transmissions illustrated in Fig. 6(b) indicate that the sound waves from 5.74 kHz to 6.75 kHz can only propagate along P1 and P3, but are hindered within P2 and P4. We show the numerical results of the energy fields for 6400 Hz and 7125 Hz in Figs. 6(c) and 6(d), respectively. Moreover, the total transmission efficiency of the three ports of the filter is almost 100% [dashed line in Fig. 6(b)] if the energy loss of the leakage or viscosity is ignored, which is consistent with the inherent property of TIs. Crucially, the independent translations in two directions indicate that the bandwidth of the edge states on the selected path can be adjusted freely. Therefore, an acoustic topological filter with multiple transmission channels can also be proposed by simply sliding the specific sections. In addition, it is noted that the filter is robust and topological protected, which is distinct from the ordinary defects.[55,56] Simulated results are given in Appendix B. This significant property to freely adjust the bandwidth of the edge states provides excellent ways for the design of adjustable topological acoustic materials.

    (a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    Figure 6.(a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    5. Conclusion

    We have numerically and experimentally demonstrated the adjustable edge states of the two-dimensional graphene-like sonic HOTI. By tuning the translation rate of the layers of the material, we can change the coupling strength of the pseudo-spins on the interface, and then control the gapped edge states. Besides, we have proposed a prototype of multi-channel topological filter with independent adjustable edge states, and experimentally demonstrated its almost lossless transmission of the acoustic waves. These results provide new perspective to control the edge states. Based on this work, we expect to expand the design of the topological insulators, and provide the novel application prospects of acoustic topological materials such as filters, waveguides, and acoustic switches.

    Appendix A: The pseudo-spins in acoustic systems

    For the acoustic waves, the in-plane velocity fields associated with the pressure fields px (r) and py (r) are given by the wave equation[24]v1=?iρω(?px?xx^+?py?yy^)=h1xx^+h1yy^,v2=?iρω(?py?xx^+?py?yy^)=h2xx^+h2yy^.

    Equations (A1) and (A2) can be rewritten as vi=hixx^+hiyy^ (i = 1,2), where x^ and y^ are the unit vectors. Since the wave functions of pressure field px (r) and py (r) are linked to each other by π/2 rotation, one has at origin (h2xh2y)=(h1y?h1x).

    The total velocity fields are obtained from v=?iρω?(px±ipy)=(h1x±ih2x)x^+(h1y±ih2y)y^=(h1x±h1y)(x^?iy^). A ± π / 2 phase shift in x and y components of velocity fields is obtained from the term ± i in front of y^ in Eq. (A4) because of the time harmonic mode exp (iωt). This phase shift represents the circular polarization of the in-plane velocity fields, which corresponds to the angular momentum of the wave function of the pressure fields in the present SC. The circular polarization can be used to simulate acoustic pseudo-spin states. We show the pseudo-spin states in Fig. A1.

    (a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    Figure 6.(a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    (a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    Figure 6.(a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    Appendix B: Edge states vs. ordinary defects

    An important character of the topological states is robust against defects. Here, we simulate the cases with a cavity or a disturbance on the nontrivial interface, compared with the ordinary defects, as illustrated in Figs. B1(a) and B1(b). The corresponding simulated transmission spectra for different cases are present in Fig. B1(c), which depicts the robustness and high transmission of the topological edge states.

    (a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    Figure 6.(a) Measured transmission spectra for port 2 (red), port 3 (green), and port 4 (blue) when βPi = 0. (b) Measured transmission spectra for port 2 (red), port 3 (green), port 4 (blue) and the sum of the three ports for βP1 = βP3 = 0.5 and βP2 = βP4 = 0.25. The green dashed curve representing the total transmission loss exhibits lossless transmission in the two theoretical frequency ranges marked with grey. (c) and (d) The simulated acoustic energy transmission distributions for 6400 Hz and 7125 Hz, respectively.

    [1] C L Kane, E J Mele. Phys. Rev. Lett., 95(2005).

    [2] C L Kane, E J Mele. Phys. Rev. Lett., 95(2005).

    [3] J E Moore. Nature, 464, 194(2010).

    [4] M Z Hasan, C L Kane. Rev. Mod. Phys., 82, 3045(2010).

    [5] X Qi, S C Zhang. Rev. Mod. Phys., 83, 1057(2011).

    [6] F D M Haldane, S Raghu. Phys. Rev. Lett., 100(2008).

    [7] Z Wang, Y D Chong, J D Joannopoulos, M Soljačić. Phys. Rev. Lett., 100(2008).

    [8] M Hafezi, E A Demler, M D Lukin, J M Taylor. Nat. Phys., 7, 907(2011).

    [9] Y Poo, R Wu, Z Lin, Y Yang, C T Chan. Phys. Rev. Lett., 106(2011).

    [10] M C Rechtsman, J M Zeuner, Y Plotnik, Y Lumer, D Podolsky, F Dreisow, S Nolte, M Segev, A Szameit. Nature, 496, 196(2013).

    [11] A B Khanikaev, S H Mousavi, W Tse, M Kargarian, A H MacDonald, G Shvets. Nat. Mater., 12, 233(2013).

    [12] G Q Liang, Y D Chong. Phys. Rev. Lett., 110(2013).

    [13] L Lu, J D Joannopoulos, M Soljačić. Nat. Photon., 8, 821(2014).

    [14] L Wu, X Hu. Phys. Rev. Lett., 114(2015).

    [15] A B Khanikaev, G Shvets. Nat. Photon., 11, 763(2017).

    [16] F Gao, H Xue, Z Yang, K Lai, Y Yu, X Lin, Y Chong, G Shvets, B Zhang. Nat. Phys., 14, 140(2018).

    [17] A B Khanikaev, R Fleury, S H Mousavi, A Alù. Nat. Commun., 6, 8260(2015).

    [18] M Xiao, W Chen, W He, C T Chan. Nat. Phys., 11, 920(2015).

    [19] Z Yang, F Gao, X Shi, X Lin, Z Gao, Y Chong, B Zhang. Phys. Rev. Lett., 114(2015).

    [20] R Fleury, A B Khanikaev, A Alù. Nat. Commun., 7(2016).

    [21] Y Peng, C Qin, D Zhao, Y Shen, X Xu, M Bao, H Jia, X Zhu. Nat. Commun., 7(2016).

    [22] C He, X Ni, H Ge, X Sun, Y Chen, M Lu, X Liu, Y Chen. Nat. Phys., 12, 1124(2016).

    [23] J Lu, C Qiu, L Ye, X Fan, M Ke, F Zhang, Z Liu. Nat. Phys., 13, 369(2017).

    [24] Z Zhang, Q Wei, Y Cheng, T Zhang, D Wu, X Liu. Phys. Rev. Lett., 118(2017).

    [25] J Lu, C Qiu, W Deng, X Huang, F Li, F Zhang, S Chen, Z Liu. Phys. Rev. Lett., 120(2018).

    [26] Y J Ding, Y G Peng, Y F Zhu, X D Fan, J Yang, B Liang, X F Zhu, X G Wan, J C Cheng. Phys. Rev. Lett., 122(2019).

    [27] C Fang, M J Gilbert, B A Bernevig. Phys. Rev. B, 86(2012).

    [28] W A Benalcazar, B A Bernevig, T L Hughes. Phys. Rev. B, 96(2017).

    [29] F Liu, K Wakabayashi. Phys. Rev. Lett., 118(2017).

    [30] J Langbehn, Y Peng, L Trifunovic, F von Oppen, P W Brouwer. Phys. Rev. Lett., 119(2017).

    [31] Z Song, Z Fang, C Fang. Phys. Rev. Lett., 119(2017).

    [32] W A Benalcazar, B A Bernevig, T L Hughes. Science, 357, 61(2017).

    [33] M Ezawa. Phys. Rev. Lett., 120(2018).

    [34] F Liu, H Deng, K Wakabayashi. Phys. Rev. Lett., 122(2019).

    [35] F Liu, H Deng, K Wakabayashi. Phys. Rev. B, 97(2018).

    [36] F Liu, M Yamamoto, K Wakabayashi. J. Phys. Soc. Jpn., 86(2017).

    [37] M Serra-Garcia, V Peri, R Süsstrunk, O R Bilal, T Larsen, L G Villanueva, S D Huber. Nature, 555, 342(2018).

    [38] H Fan, B Xia, L Tong, S Zheng, D Yu. Phys. Rev. Lett., 122(2019).

    [39] C W Peterson, W A Benalcazar, T L Hughes, G Bahl. Nature, 555, 346(2018).

    [40] S Imhof, C Berger, F Bayer, J Brehm, L W Molenkamp, T Kiessling, F Schindler, C H Lee, M Greiter, T Neupert, R Thomale. Nat. Phys., 14, 925(2018).

    [41] J Noh, W A Benalcazar, S Huang, M J Collins, K P Chen, T L Hughes, M C Rechtsman. Nat. Photon., 12, 408(2018).

    [42] B Xie, H Wang, H Wang, X Zhu, J Jiang, M Lu, Y Chen. Phys. Rev. B, 98(2018).

    [43] S Mittal, V V Orre, G Zhu, M A Gorlach, A Poddubny, M Hafezi. Nat. Photon., 13, 692(2019).

    [44] M Li, D Zhirihin, M Gorlach, X Ni, D Filonov, A Slobozhanyuk, A Alù, A B Khanikaev. Nat. Photon., 14, 89(2020).

    [45] H Xue, Y Yang, F Gao, Y Chong, B Zhang. Nat. Mater., 18, 108(2019).

    [46] X Ni, M Weiner, A Alù, A B Khanikaev. Nat. Mater., 18, 113(2019).

    [47] X Zhang, H Wang, Z Lin, Y Tian, B Xie, M Lu, Y Chen, J Jiang. Nat. Phys., 15, 582(2019).

    [48] Z Chen, C Xu, R Al Jahdali, J Mei, Y Wu. Phys. Rev. B, 100(2019).

    [49] Z Zhang, M Rosendo López, Y Cheng, X Liu, J Christensen. Phys. Rev. Lett., 122(2019).

    [50] S Yu, X Ni, Y Xu, C He, N Priyanka, M Lu, Y Chen. Chin. Phys. Lett., 33(2016).

    [51] A Khelif, A Choujaa, S Benchabane, B Djafari-rouhani, V Laude. Appl. Phys. Lett., 84, 4400(2004).

    [52] C Ahmet, A K Olgun, U Bulent. Chin. Phys. B, 22(2013).

    [53] Y Ota, F Liu, R Katsumi, K Watanabe, K Wakabayashi, Y Arakawa, S Iwamoto. Optica, 6, 786(2019).

    [54] O Ryo, H Shin, N Takeshi. Phys. Rev. B, 100(2019).

    [55] X Kong, L Yue, Y Chen, Y Liu. Chin. Phys. B, 21(2012).

    [56] Z Fang, Z Hou. Chin. Phys. Lett., 35(2018).

    Tools

    Get Citation

    Copy Citation Text

    Zhang-Zhao Yang, Jin-Heng Chen, Yao-Yin Peng, Xin-Ye Zou. Symmetry-controlled edge states in graphene-like topological sonic crystal[J]. Chinese Physics B, 2020, 29(10):

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Received: Apr. 11, 2020

    Accepted: --

    Published Online: Apr. 21, 2021

    The Author Email: Xin-Ye Zou (xyzou@nju.edu.cn)

    DOI:10.1088/1674-1056/ab9c0e

    Topics