Reconfigurable topological phononic crystal slabs Lu-yang Feng, Hong-bo Huang, Jian-cheng Zhang, Xiao-ping Xie*, and Jiu-jiu Chen* State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, College of Mechanical and Vehicle Engineering, Hunan University, Changsha 410082, China
Abstract We propose two-dimensional reconfigurable phononic crystal slabs that support topologically protected edge states for Lamb waves. It consists of a triangular lattice of triangular-like air holes perforated in a slab and two sandwich cylinders are placed in the matrix of each unit cell. The pseudospin-orbit coupling is achieved by shifting the sandwich cylinders up or down. The robust transport of Lamb waves along reconfigurable interfaces is further demonstrated. The dynamically reconfigurable design provides an excellent platform for tuning the frequency range of topological edge states and steering efficient transport of Lamb waves along arbitrary pathways, which can be employed for elastic-wave communications, signal processing, and sensing. Keywords: topological edge state, reconfigurability, Lamb wave, pseudospin-orbit coupling, robustness
1. Introduction The discoveries of the quantum Hall effect, the quantum spin Hall effect, and topological insulators have revolutionized our scientific cognition of condensed matter physics. Since the gapless edge states in graphene and HgTe quantum wells were discovered [1,2], more studies on analogous states in structures supporting classical waves have been conducted. Topologically protected edge states have been demonstrated in the topological photonic crystals [3-13], which is of great importance in both fundamental science and the practical application. In parallel, these concepts have also inspired a novel field of topological acoustics [14-27]. Achieving topologically protected acoustic-wave propagation depends 1
mainly on three categories of mechanisms. The first is to realize one-way edge states similar to the quantum Hall effect. By introducing a circulating fluid flow to break the time-reversal symmetry [14-16], the topologically protected unidirectional edge states with strong robustness were demonstrated in nonreciprocal systems. The second mechanism exhibits the pairwise topological edge states in analogy with the quantum spin Hall effect. By constructing a pair of pseudospins [19-22], the robust pseudospin-dependent transport has been successfully demonstrated in systems with pseudo-time-reversal symmetry. The third is to achieve valley-polarized topological edge states by emulating the valley Hall effect. The topological valley transport has been successfully observed by breaking the spatial-reversal symmetry [23-26]. Different from the previous studies on acoustics in the fluid domain where acoustic waves are purely longitudinal, the research of topological physics in the solid slabs is extremely challenging due to the mixture of all three polarizations of elastic waves. Recently, topologically protected edge states for Lamb waves have been achieved by using a dual-scale topological phononic crystal (TPC) [19]. However, its topological bandgaps and transport pathways are limited to specific situations that cannot be tuned after machining. Most recently, many studies on the reconfigurability of the TPCs have been performed in the fluid domain due to the merits of tuning the topological bandgaps and transport pathways [22,23], which can enrich the design and use of the TPCs. So far, there is few investigation on the reconfigurability of the TPCs in the solid domain owing to the limitation of the movement of the scatterers. Thus, it is meaningful and challenging to investigate the reconfigurable TPCs in the solid domain. The purpose of this letter is to realize the reconfigurable phononic crystal (PC) slabs supporting topological edge states for Lamb waves. It consists of a triangular lattice of triangular-like air holes perforated in a slab and two sandwich cylinders are inserted in the matrix material of each unit cell. The inversion symmetry breaking caused by moving the sandwich cylinders up or down leads to the pseudospin-orbit coupling. We explicitly demonstrate that the topologically protected transport of 2
Lamb waves is robust against various kinds of defects along reconfigurable interfaces. The tunability of the designed TPC slabs makes it possible to dynamically control the frequency range of topological edge states and reconfigure arbitrary transport pathways for Lamb waves without obvious backscattering.
2. Reconfigurable topological design As shown in Figs. 1(a) and (b), the designed PC slabs are composed of a triangular array of triangular-like air holes perforated in an elastic slab and the host material of each unit cell is inserted by two movable sandwich cylinders, each of which consists of one cylindrical core (blue) in the center and two cylindrical face sheets (green) on both ends. The lattice constant is a = 40 mm and the slab thickness is t = 10.5 mm. The fillet radius of the triangular-like hole is r = 8.3 mm and the distance from center to side of the triangular-like hole is l = 14.85 mm. The diameter of the cylinder, the length of the cylindrical core, and the thickness of the cylindrical face sheet are d = 6.4 mm, h = 10.5 mm, and b = 0.64 mm, respectively; the displacement of the cylinders is denoted as s (suppose that s = 0 mm when the neutral plane of the sandwich cylinders overlaps with that of the slab [middle panel in Fig. 1(b)]). The material of the cylindrical core of the sandwich cylinder is polyetheretherketone (PEEK) which is a common engineering plastics with excellent friction and wear properties as well as mechanical properties, and the material of the remaining structures is steel. The material parameters are chosen as follows: the density ρ = 7800 kg/m3, the Young modulus E = 206 GPa, and the Poisson's ratio σ = 0.3 for steel; ρ = 1285 kg/m3, E = 4.21 GPa, and σ = 0.388 for PEEK. The propagation of Lamb waves
is
governed
by
the
wave
equilibrium
equation
u ( 2 )( u) ( u), where ρ is the material density, u is the displacement vector, and both E /[(1 )(1 2 )] and E /[ 2(1 )] are the Lame’s coefficients. The wave equilibrium equation can be transferred into the eigenvalue problem (K 2 (k )M)U 0 [28], where K and M are the stiffness matrix and mass matrix, respectively, and k represents the wave vector in the 3
irreducible Brillouin zone. The eigenvalue problem is solved by using eigenfrequency analysis of solid mechanics in the commercial finite element software COMSOL Multiphysics. Periodic boundary conditions are imposed on the four external cross section of the unit cell, and the stress-free boundary conditions are set on the other surfaces. When the wave vector k varies along the boundary of the irreducible Brillouin zone (Γ-K(K')-M-Γ) [inset of Fig. 2(a)], the Lamb-wave band structures ω(k) are obtained, as shown in Fig. 1(a), (c), and (d). For a phononic system, it is necessary to increase the degrees of freedom to two-fold states to emulate two pseudospin states [20,22]. Consequently, a four-fold Dirac degeneracy in the phononic band structure is required. For s = 0 mm [middle panel of Fig. 1(b)], both the in-plane C3v and the z-direction symmetries of the PC slabs are maintained, resulting in the four-fold Dirac degeneracy at point K [Fig. 2(a)]. The degenerate Lamb-wave modes are classified by their displacement fields as the symmetric (S) modes (blue lines) and the anti-symmetric (A) modes (red lines). The frequency and group velocity of these two families of modes are matched around the Dirac point, which allows one to use unitary transform of the orthogonal basis (S and A modes) as pseudospin states [19,20]. We then introduce pseudospin-orbit coupling to induce topological phase translation, accomplished by the opening of a complete topological bandgap. When the sandwich cylinders are moved up or down from the equilibrium position [left or right panel of Fig. 1(b)], the z-direction inversion symmetry of the PC slabs is broken, leading to the formation of a complete topological bandgap and two two-fold degeneracies: the lower (upper) A and upper (lower) S modes [Fig. 2(c)]. For all frequencies in the proximity of the original Dirac points where the frequency and group velocity degeneracy are maintained, the A and S modes of the PC slabs become hybridized ( A S )
2 [19], which can be used as two pseudospin states [20]. As
shown in Fig. 2(b), the evolution of displacement fields of the unit cells from s = 0 mm to s ≠ 0 mm also confirms the pairwise hybridizations of the A and S modes. In practice, the z-direction symmetry breaking caused by moving the sandwich cylinders 4
produces an effect analogous to spin-orbit coupling in graphene. To reveal the pseudospin-orbit effect, the effective Hamiltonian for the present PC slabs can be obtained by the perturbation theory as [10,19]:
Hˆ D (ˆz sˆ0ˆ xkx ˆ0 sˆ0ˆ yk y ) mˆz sˆzˆ z ,
(1)
where ˆ i , ˆi , and sˆi (i x, y, z ) are the Pauli matrices acting on orbit, valley, and pseudospin state vectors, respectively; ˆ0 and sˆ0 are the identity matrix, vD is the group velocity near the Dirac point, and m ( m gap 2 , where gap is the topological bandgap) is the effective mass induced by the z-direction symmetry breaking. The last term mˆz sˆzˆ z [Eq. (1)] caused by the z-direction symmetry breaking is identical to the spin-orbit coupling of the Kane-Mele Hamiltonian [1], which demonstrates that the z-direction symmetry breaking leads to the pseudospin-orbit coupling in the present PC slabs. The occurrence of the pseudospin-orbit coupling can induce a topological phase transition. In order to reveal the topological phase transition [29-30], the pseudospin Chern numbers are introduced to characterize various topological phases. For the Hamiltonian given in Eq. (1), the pseudospin Chern numbers of each band can be evaluated as
C sgn (m) [10,19]. For s 0 mm [middle panel of Fig. 1(b)], the pseudospin Chern numbers C are equal to zero due to m 0 , meaning trivial topological phase; while for s 0 mm, the pseudospin Chern numbers C are nonzero due to m 0 , corresponding to the nontrivial topological phase. Obviously, by moving the sandwich cylinders up or down, the PC slabs experience a topological phase transition. Note that although the band structures are identical in the PC slabs with the sandwich cylinders shifted up and down, the sign of the effective mass m is opposite [left or right panel of Fig. 1(b)], indicating the opposite pseudospin Chern numbers and different topological properties. To illustrate this further, we calculate numerically the pseudospin Chern numbers by integrating the Berry curvature (k ) k A(k ) of the same pseudospins throughout the first Brillouin zone, where A(k) is the Berry connection and k ( kx , ky , kz ) [19,20]. The Berry connection of the nth band is 5
defined as An (k ) Im unit cell un (k ) | k | un (k ) dV , where
u n (k )
is the normalized
Bloch wave function that can be obtained through the COMSOL simulation. In this way, we find that the pseudospin Chern numbers for four bands from lower to upper reverse from C = -1, +1, -1, +1 to C = +1, -1, +1, -1, as the sandwich cylinders are moved from up to down. The opposite pseudospin Chern numbers reveal that the two PC slab structures are topologically distinct, which is consistent with the previous discussion. Therefore, the ability to independently shift the sandwich cylinders makes it possible to conveniently control the topological properties of the PC slabs. Next, we investigate the influence of the displacement s of the sandwich cylinders on the topological bandgap. As shown in Fig. 2(d), the topological bandgap enlarges with the increase of |s|. Because the increase of the displacement s leads to the bigger differences of the height of the stubbed surface on both sides of the slabs, the bigger bandgap can form at crossings of the Lamb-wave band structure [31,32]. Therefore, we can tune the topological bandgap by simply moving the sandwich cylinders. The tunability of the topological bandgap, absence in the conventional fixed TPCs, can dramatically enrich the design and use of the TPC slabs.
3. Topological edge states To realize topologically protected edge states for Lamb waves, we construct a ribbon-shaped supercell with one interface which is between two PC slab structures with different pseudospin Chern numbers [Fig. 3(a)]. The sandwich cylinders of the left PC slab structure (type I) are shifted up, while the sandwich cylinders of the right PC slab structure (type II) are shifted down. We calculate the projected band structure by eigenfrequency analysis of solid mechanics in the COMSOL. Periodic (absorbing) boundary conditions are applied to the edges along x (y) direction. When the wave vector k sweeps from - a to a along x direction, the projected band structure is obtained, as shown in Fig. 3(b). It is observed that there are two pairs of topological edge states within the topological bandgap: the forward-propagating edge states (red lines) and the backward-propagating edge states (blue lines). Because the pseudospin 6
Chern numbers differ by 2 across the interface, according to the bulk-boundary correspondence [33-35], there are four topologically protected edge states for the TPC slabs [10,19]. As shown in Fig. 3(c), the displacement fields of the supercell display that two pairs of topological edge states are confined at the interface and decay rapidly into the bulk. In order to further understand the topological edge states, we check the real-space distribution of the Poynting vector P Re[V T ] / 2 ( V is the velocity vector and T is the stress tensor) averaged over a period 2 , which describes the energy flow in the TPC slabs [36]. The energy flow of the forward-propagating edge states [1 and 3 enlarged views in Fig. 3(c)] has the opposite senses of rotation on both sides of the interface, resulting in the directional transport of the Lamb waves [10]. While for backward-propagating edge states, the senses of rotation are reversed [2 and 4 enlarged views in Fig. 3(c)]. It is worth noting that the conjunctive points [green points in Fig. 3(b)] between two topological edge states are transition points where there is no energy flow. Additionally, the frequency range of topological edge states is controllable considering the tunable topological bandgap. To verify topologically protected transmission of edge states, full wave simulation in finite TPC slabs (16×16 unit cells) is carried out via frequency-domain analysis of solid mechanics in the COMSOL. As shown in Fig. 4, this TPC slabs consist of two different PC slab structures (type I and II), separated by a linear interface. Perfectly matched layers are applied around the TPC slabs to suppress the boundary reflections. Lamb waves (51.4 KHz) are excited by a small displacement on the left end (top panel) and right end (bottom panel) of the interface, respectively. The Lamb waves propagate along the linear interface without obvious backscattering, which demonstrates efficient transmission of topological edge states for Lamb waves in the TPC slabs. Considering the unidirectional transmission of each kind of topological edge states [21,22], the propagation of Lamb waves in two directions along the interface [Fig. 4] verifies the presence of two types of edge states in the TPC slabs. Topologically protected edge states are robust owing to the locking of the pseudospins to their respective propagation directions. In order to verify the 7
robustness of topologically protected edge states, we intentionally introduce different defects into TPC slabs and investigate their effects. As shown in bottom panel of Fig. 5, there are three types of common defects (two sharp bends, cavity, and disorder) in a zigzag interface of the TPC slabs. It is observed from the top panel of Fig. 5 that Lamb waves excited by a small displacement at the left end maintain a complete transmission along the interface without obvious backscattering. The result explicitly demonstrates the topological edge states are robust against various kinds of defects, which is a universal characteristic caused by the intrinsic topological feature of the TPC slabs and their topological protection. In order to illustrate the reconfigurability of the TPC slabs, two different transport pathways are successively constructed by shifting the sandwich cylinders up or down in the same TPC slabs. As shown in Fig. 6, during 0 - t1 period, Lamb waves incident from the port 1 transmit to port 4 along the first pathway (green line) with high transmission efficiency; during t1 - t2, the second topological transport pathway (light blue line) is reconfigured rapidly from the first pathway; during t2 - t3 period, Lamb waves effectively propagate from port 3 to port 2 along the second pathway. Therefore, by moving the sandwich cylinders, we can realize topologically protected transmission of Lamb waves along arbitrary pathways from arbitrary input ports to output ports at different time. When each sandwich cylinder is equipped with one electronically controlled linear actuator, the reconfiguration of the topological transport pathways can be completed instantaneously. In other words, the dynamic reconfiguration of the topological transport pathways can be achieved in the TPC slabs, which could be used for steering efficient transport of Lamb waves along arbitrary pathways in real time.
4. Conclusion To sum up, we have presented the PC slabs exhibiting tunable topological edge states and reconfigurable transport pathways for Lamb waves. The pseudospin-orbit coupling is achieved by shifting the sandwich cylinders up or down. With numerical 8
simulations, we unambiguously demonstrate topologically protected propagation of Lamb waves along reconfigurable interfaces and the robustness of topological edge states against defects. In view of the scale invariance of the elastic wave equation, Lamb waves can be manipulated and controlled over a higher frequency range. The dynamically reconfigurable TPC slabs provide an excellent platform for freely steering the topologically protected transport of Lamb waves along arbitrary pathways, which has potential applications in elastic-wave communications, signal processing, and sensing. *Corresponding author. Email address:
[email protected] (X.-P. Xie);
[email protected] (J.-J. Chen).
Acknowledgements The authors gratefully acknowledge financial support from National Science Foundation of China under Grant No.11374093 and Young Scholar fund sponsored by common university and college of the province in Hunan.
9
References [1] C.L. Kane, E.J. Mele, Phys. Rev. Lett. 95 (2005) 226801. [2] B.A. Bernevig, T.L. Hughes, S.C. Zhang, Science 314 (2006) 1757. [3] F.D.M. Haldane, S. Raghu, Phys. Rev. Lett. 100 (2008) 013904. [4] Z. Wang, Y.D. Chong, J.D. Joannopoulos, M. Soljačić, Phys. Rev. Lett. 100 (2008) 013905. [5] Z. Wang, Y. Chong, J.D. Joannopoulos, M. Soljačić, Nature 461 (2009) 772. [6] Y. Poo, R.X. Wu, Z. Lin, Y. Yan, C.T. Chan, Phys. Rev. Lett. 106 (2011) 093903. [7] A.B. Khanikaev, S.H. Mousavi, W.K. Tse, M. Kargarian, A.H. MacDonald, G. Shvets, Nat. Mater. 12 (2013) 233. [8] W.J. Chen, S.J. Jiang, X.D. Chen, B. Zhu, L. Zhou, J.W. Dong, C.T. Chan, Nat. Commun. 5 (2014) 5782. [9] L.H. Wu, X. Hu, Phys. Rev. Lett. 114 (2015) 223901. [10] T. Ma, A.B. Khanikaev, S.H. Mousavi, G. Shvets, Phys. Rev. Lett. 114 (2015) 127401. [11] X.J. Cheng, C. Jouvaud, X. Ni, S.H. Mousavi, A.Z. Genack, A.B. Khanikaev, Nat. Mater. 15 (2016) 542. [12] T. Ma, G. Shvets, New. J. Phys. 18 (2016) 025012. [13] M. Goryachev, M.E. Tobar, Phys. Rev. Appl. 6 (2016) 064006. [14] Z. Yang, F. Gao, X. Shi, X. Lin, Z. Gao, Y. Chong, B. Zhang, Phys. Rev. Lett. 114 (2015) 114301. [15] A.B. Khanikaev, R. Fleury, S.H. Mousavi, A. Alù, Nat. Commun. 6 (2015) 8260. [16] Z.G. Chen, Y. Wu, Phys. Rev. Appl. 5 (2016) 054021. [17] H. Huang, J. Chen, S. Huo, J. Phys. D: Appl. Phys. 50 (2017) 275102. [18] Y.G. Peng, C.Z. Qin, D.G. Zhao, Y.X. Shen, X.Y. Xu, M. Bao, H. Jia, X.F. Zhu, Nat. Commun. 7 (2016) 13368. [19] S.H. Mousavi, A.B. Khanikaev, Z. Wang, Nat. Commun. 6 (2015) 8682. 10
[20] C. He, X. Ni, H. Ge, X.C. Sun, Y.B. Chen, M.H. Lu, X.P. Liu, Y.F. Chen, Nat. Phys. 12 (2016) 1124. [21] J. Mei, Z. Chen, Y. Wu, Sci. Rep. 6 (2016) 32752. [22] Z. Zhang, Q. Wei, Y. Cheng, T. Zhang, D. Wu, X. Liu, Phys. Rev. Lett. 118 (2017) 084303. [23] J. Lu, C. Qiu, L. Ye, X. Fan, M. Ke, F. Zhang, Z. Liu, Nat. Phys. 13 (2017) 369. [24] L. Ye, C. Qiu, J. Lu, X. Wen, Y. Shen, M. Ke, F. Zhang, Z. Liu, Phys. Rev. B 95 (2017) 174106. [25] S.Y. Huo, J.J. Chen, H.B. Huang, G.L. Huang, Sci. Rep. 7 (2017) 10335. [26] J.J. Chen, S.Y. Huo, Z.G. Geng, H.B. Huang, X.F. Zhu, AIP Adv. 7 (2017) 115215. [27] Y.G. Peng, Y.X. Shen, D.G. Zhao, X.F Zhu, Appl. Phys. Lett. 110 (2017) 173505. [28] Y. Wu, K. Yu, L. Yang, R. Zhao, X. Shi, K. Tian, J. Sound Vib. 413 (2018) 101. [29] R. Chen, B. Zhou, Phys. Lett. A 381 (2017) 944. [30] Y. Yang, J. Yang, X. Li, Y. Zhao, Phys. Lett. A 382 (2018) 723. [31] T.T. Wu, Z.G. Huang, T.C. Tsai, T.C. Wu, Appl. Phys. Lett. 93 (2008) 111902. [32] T.C. Wu, T.T. Wu, J.C. Hsu, Phys. Rev. B 79 (2009) 104306. [33] M.Z. Hasan, C.L. Kane, Rev. Mod. Phys. 82 (2010) 3045. [34] R.S. Mong, V. Shivamoggi, Phys. Rev. B 83 (2011) 125109. [35] L. Lu, J D.Joannopoulos, M. Soljačić, Nat. Photonics 26 (2014) 821. [36] B.A. Auld, Acoustic fields and waves in solids, Wiley, New York, 1973.
11
Figure Captions Fig. 1. (a) Perspective view (left panel) of a unit cell and top view (right panel) of the PC slabs formed by a triangular array of triangular-like air holes in a slab. The host material of each unit cell is inserted by two movable sandwich cylinders (b) Sandwich cylinders at different positions in the slab: from left to right, the sandwich cylinders are shifted from the lower to upper position, leading to a topological transition as indicated by the mass term m. The sandwich cylinder consists of one cylindrical core (blue) in the center and two cylindrical face sheets (green) on both ends. The black and green dotted lines represent the position of the neutral planes of the slab and the sandwich cylinders, respectively. Fig. 2. (a) Band structure for s = 0 mm. A four-fold Dirac degeneracy is formed at point K. The red and blue lines indicate the anti-symmetric (A) modes and the symmetric (S) modes, respectively. (b) The evolution of displacement fields of the unit cells around the point K from s = 0 mm to s 0 mm. The four displacement fields from bottom to top on the right or left side correspond to four bands of interest from lower to upper, respectively. (c) Band structure for |s| = 0.64 mm. A complete topological bandgap is induced by shifting sandwich cylinders up or down. The parameters C are the pseudospin Chern numbers of the PC slabs with the sandwich cylinders shifted up. (d) The evolution of the topological bandgap as a function of the displacement |s| of the sandwich cylinders. Fig. 3. (a) Side view (upper panel) and top view (lower panel) of a ribbon-shaped supercell consisting of 16 unit cells on each side of the interface (light blue dotted lines). The red arrow on each side of the interface indicates the moving directions of the corresponding sandwich cylinders. The PC slabs on both sides of the interface are marked as I and II, corresponding to the mass term m > 0 and m < 0. (b) Projected band structure of the supercell. The red and blue dotted lines indicate the forward-propagating and backward-propagating edge states, respectively. The black dotted lines represent the bulk bands. (c) Displacement field and Poynting vector of 12
the supercell at points 1, 2, 3, and 4 in (b). The black and red arrows denote the Poynting vector and the overall energy flow, respectively. Fig. 4. Displacement field (51.4 kHz) of the TPC slabs with a linear interface. The sign I and II represent the PC slabs with the sandwich cylinders shifted up and down, respectively. Lamb waves are excited by a small displacement on the left end (top panel) and right end (bottom panel) of the interface, respectively. The red arrows indicate the input/output ports. Fig. 5. Displacement field (top panel) of the TPC slabs with a zigzag interface including three kinds of defects (bottom panel): sharp bends, cavity and disorder (the disordered sandwich cylinders are marked in red). Lamb waves are excited by a small displacement on the left end of the interface. The red arrows indicate the input/output ports. Fig. 6. Dynamic reconfiguration of topological transport pathways and displacement fields. The green and light blue lines represent two different transport pathways which are successively constructed by shifting the sandwich cylinders up or down in the same TPC slabs. The red and gray arrows indicate the presence and absence of incident or received Lamb waves, respectively.
13
(a)
(b) Fig. 1.
14
(a)
(b)
(c)
(d)
Fig. 2.
15
(a)
(b)
(c) Fig. 3.
16
Fig. 4.
17
Fig. 5.
18
Fig. 6.
19