Next Article in Journal
Magneto-Viscoelastic Materials: Memory Functionals and Rate Equations
Previous Article in Journal
Excellent Thermoelectric Performance of 2D CuMN2 (M = Sb, Bi; N = S, Se) at Room Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inverse Design of Valley-Like Edge States of Sound Degenerated Away from the High-Symmetry Points in a Square Lattice

1
College of Aerospace Engineering, Chongqing University, Chongqing 400044, China
2
Nanjing Research Institute, Chongqing University, Nanjing 211800, China
3
School of Mechanical Engineering, and National Engineering Research Center for Technological Innovation Method and Tool, Hebei University of Technology, Tianjin 300401, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(19), 6697; https://doi.org/10.3390/ma15196697
Submission received: 30 August 2022 / Revised: 16 September 2022 / Accepted: 22 September 2022 / Published: 27 September 2022

Abstract

:
Robust edge states of periodic crystals with Dirac points fixed at the corners or centers of the Brillouin zones have drawn extensive attention. Recently, researchers have observed a special edge state associated with Dirac cones degenerated at the high symmetric boundaries of the first irreducible Brillouin zone. These nodal points, characterized by vortex structures in the momentum space, are attributed to the unavailable band crossing protected by mirror symmetry. By breaking the time reversal symmetry with intuitive rotations, valley-like states can be observed in a pair of inequivalent insulators. In this paper, an improved direct inverse design method is first applied to realize the valley-like states. Compared with the conventional strategy, the preparation of transition structures with degeneracy points is skipped. By introducing the quantitative gauge of mode inversion error, insulator pairs are directly obtained without manually tuning the structure with Dirac cone features.

1. Introduction

Topological materials have drawn wide attention in the fields of condensed matter physics and engineering [1,2,3,4]. Valley materials are first proposed in graphene structures [5,6,7,8]. Similar two-dimensional materials are observed in other hexagonal lattices [9,10,11,12,13,14,15,16,17], such as bilayer graphene [15,16,17] and transition metal dichalcogenides [13,14], which are typical valleytronic materials.
Similar to electron spin, topological valley transport, which is usually induced by breaking the inversion symmetry of the system, has the potential to transfer information and energy. Inspired by the intriguing ability of supplying a topologically protected transport with high transmission rate and good robustness against backscattering, tremendous progress has been made in the study of non-trivial states in different physical systems, e.g., electromagnetic [18,19,20,21,22,23,24], acoustics [25,26,27,28,29,30,31,32,33,34] and mechanics [35,36,37,38,39,40,41].
Current studies mainly focus on the Dirac cones [42] or double Dirac cones [43,44,45] fixed at the high-symmetric points of the Brillouin zones in hexagonal lattices. However, phononic insulators beyond typical graphene-like lattices are relatively rare. Existing studies only focus on the realization of valley states or valley-like states through conventional procedures by: (1) finding a structure with Dirac cones; (2) splitting the Dirac cone by introducing a disturbance, which breaks the inversion symmetry or mirror symmetry of the base structure; (3) verifying the physical properties of the outcomes after disturbance. From a practical point of view, the desired wild bandgap (i.e., operation bandwidth) may not be obtained through the intuitive procedure mentioned above. Due to the existence of the base structure, the design space for breaking the symmetries is very limited. In general, designers can only choose from some common operations, such as rotation [28], offset [46], and zooming [47]. The trial-and-error procedure lacks mathematical optimality and is often time consuming. To overcome the aforementioned drawbacks, inverse design methods have been proposed to design topological materials with Dirac cones [48,49]. Furthermore, topological insulator pairs have been successfully obtained at certain frequencies by optimizing the transmission path [50,51] or the field distribution of eigenmodes [52,53].
In this paper, a direct inverse design method is adopted to realize valley-like edge states. The geometries of potential insulators are described by the B-spline method. By setting control points in an equiangular distribution, an arbitrary geometry can be determined by a finite set of parameters. In addition, by introducing symmetry restrictions, mirror symmetry breaking and control parameter reduction are achieved simultaneously, which are beneficial to gradient-free algorithms. By utilizing band inversion, mode inversion error is set as a quantitative gauge that is used as a constraint in the design process. Additionally, maximizing the width of common bandgap is set as the objective function. Then, the inverse design is converted into an optimization problem, which can be solved by a genetic algorithm. Compared with the design in the literature [34], the relative bandgap of the inverse design is significantly improved by around 70%.

2. Description of the Geometry of Scatters

According to existing studies [34,54], valley-like states can be realized with only one scatter in a square lattice. Thus, it is appropriate to describe the geometry of scatter with a B-spline curve. As shown in Figure 1, an arbitrary geometry can be described with well-constructed control points. Through introducing extra constraints on the distance between the control point and the origin point, different symmetries can be effectively realized. Moreover, a reference angle is introduced to improve the representation of scatters. After many trials, we found that 14 control points were sufficient for the geometry description, which can be described by a B-spline formulated as:
C ( u ) = i = 0 n + 1 N i , 2 ( u ) P i , ( u l o w e r u u u p p e r )
where n is the number of independent control points (i.e., n = 12) and Pi represents the coordinates of control points, which can be expressed as:
P i = ( r i cos ( ( i 1 ) π 6 + θ ) , r i sin ( ( i 1 ) π 6 + θ ) ) T , i = 1 , , n
with r 1 , , r n denoting the distances from the origin point to the control points and θ denoting the reference angle between the x-axis and the line along the origin point and P 1 . In order to have a closed B-spline curve, two independent control points are introduced, i.e., P 0 = P 12 and P 13 = P 1 . By merging the first two control points with the last two control points, the continuity of the closed curve is guaranteed. In addition, the second order B-spline basis function N i , 2 can be obtained with recurrence formulas as follows:
N i , 0 ( u ) = 1 ,     i f   u i u < u i + 1 0 ,   o t h e r w i s e
N i , p ( u ) = u u i u i + p u i N i , p 1 ( u ) + u i + p + 1 u u i + p + 1 u i + 1 N i + 1 , p 1 ( u )
Since the first p and last p control points are wrapped together (p = 2), the value of knot u i and the numerical range of parameter u in Equation (1) should be in the following form:
u i = i n + p + 1 ,                       i = 0 , 1 , , n + p + 1 u l o w e r = p n + p + 1 u u p p e r = n + 1 n + p + 1
Then, the closed B-spline curve can be determined by a vector V = ( r 1 , r 2 , , r 12 , θ ) T accurately. By flexibly setting these parameters, the group symmetries can be inherently realized. For example, if we set r i = r i + 3 = r i + 6 = r i + 9 ( i = 1 , 2 , 3 ) and θ = 0 , the geometry of scatters will have C4v point-group symmetry. Furthermore, a previous study has shown that the emergency of Dirac points used in this paper is only protected by mirror symmetry [55]. Thus, mirror symmetry breaking can be easily realized by introducing C4 symmetry in the geometry of scatters when we set r i = r i + 3 = r i + 6 = r i + 9 ( i = 1 , 2 , 3 ) . Inspired by the reference design [34], chiral constraints can also be introduced in our design procedure to describe a pair of insulators and reduce the design variables. Then, the number of independent variables is reduced to 4, i.e., V = ( r 1 , r 2 , r 3 , θ ) T . The calculation cost of gradient-free algorithms such as the genetic algorithm is acceptable when applied to these highly reduced design variables.

3. Inverse Design

3.1. Quantitative Measures of Band Inversion

Studies have shown [28,48] that the emergency of localized topological states is accompanied by band inversion, i.e., parts of the relevant two bands will exchange when nontrivial states occur. In addition, the relevant mode fields in a phononic crystal with different chiralities will be reversed compared with the other crystals. The quality of band inversion can be quantified by the difference of the modes of two chiral phononic crystals around Dirac cones. To achieve such requirements, the mode inversion error of two chiral crystals about their i-th and (i + 1)-th bands is introduced as
E i , i + 1 = Ω ( ψ ¯ A , i ψ ¯ B , i + 1 + ψ ¯ B , i ψ ¯ A , i + 1 ) d Ω Ω ( max ( ψ ¯ A , i , ψ ¯ B , i + 1 ) + max ( ψ ¯ B , i , ψ ¯ A , i + 1 ) ) d Ω
where Ω denotes the intersection area of the solid part of these two crystals, subscript A (B) represents the chiral crystal A (B), and ψ ¯ is the normalized amplitude field defined by the orthonormal property of the basis function [56], i.e., ( 2 π ) Ω c e l l ψ i , k 0 * ( r ) 1 B r ( r ) ψ j , k 0 ( r ) d r = δ i j . Considering the difference between two chiral crystals and the computational accuracy, the mode inversion error should be very small but never drop to zero. It is therefore recommended to use this error as a constraint or multi-objective function during optimization.

3.2. Common Bandgap Width

Considering the edge states appeared in the nontrivial common bandgap of a pair of insulators, we adopt the common bandgap width as an objective function, which is formulated as:
G i , i + 1 = min ( ( f A , i + 1 f A , i ) , ( f B , i + 1 f B , i ) , ( f A , i + 1 f B , i ) , ( f B , i + 1 f A , i ) )
where subscript f A , i represents the frequency of the i-th band of the crystal A. It should be noted that, though the value of i used in this paper is four, the value can be different depending on specific situations. For example, parameter i can be six in a 2D photonic crystal with a triangular lattice [55] or one in a 2D sonic crystals with a triangular lattice [28].

3.3. Problem Formulations

Inverse design of valley-like states will be achieved by applying the objective optimization in the following form:
F o r   k   o n   X M : F i n d :   V = ( r 1 , r 2 , r 3 , θ ) T M a x i m i z e :   G i , i + 1 = min ( ( f A , i + 1 f A , i ) , ( f B , i + 1 f B , i ) , ( f A , i + 1 f B , i ) , ( f B , i + 1 f A , i ) ) S u b j e c t   t o :   [ 1 ρ r ( r ) p ] = ω 2 c 2 p B r ( r )                                             E i , i + 1 = Ω ( | | ψ | ¯ A , i | ψ | ¯ B , i + 1 | + | | ψ | ¯ B , i | ψ | ¯ A , i + 1 | ) d Ω Ω ( max ( | ψ | ¯ A , i , | ψ | ¯ B , i + 1 ) + max ( | ψ | ¯ B , i , | ψ | ¯ A , i + 1 ) ) d Ω E p r e s e t                                             ( 0.1 × a ) r i ( 0.4 × a ) , i = 1 , 2 , 3                                             0 θ π 6
where a represents the lattice constant ( a = 24 mm ); the density and the sound speed of wave medium (air) are 1.25 kg/m3 and 343 m/s; ρ r ( r ) and B r ( r ) are the relative mass density and bulk modulus, respectively; E p r e s e t is a positive number, which is close to zero; r i represents the distance from the origin point to the control point P i . Additionally, the value of r i is restricted to [ 0.1 × a , 0.4 × a ] in order to avoid self-interpretation.

4. Valley-Like Vortex States

The geometry of phononic crystals has C4 point-group symmetry, which leads to the mismatch in mirror symmetries between the square scatter and lattice. Since our method is a direct method, the construction procedure of the specific Dirac cones is skipped and the transition structure remains unknown. The outcomes of the imaginary structure after a so-called open–close–reopen operation [28,48] are illustrated in Figure 2, which clearly shows that the eigenmodes of the fourth and fifth bands are inversed near the Dirac cones, i.e., a valley Hall phase transition is induced. This phenomenon can be regarded as a band inversion process, which is similar to the valley Hall phase transition in electronic systems.
Moreover, the acoustic valley-like states possess chiral characteristics. Figure 2 shows the notable features of vortices of power flows in these two different crystals. The forth eigenmode of crystal A exhibits the same vortex chirality as the fifth eigenmode of crystal B. The two eigenmodes of crystal A (or B) exhibit the opposite vortex chirality, while the forth eigenmode of crystal B shows the same vortex chirality with the fifth eigenmode of crystal A. In other words, the inversion of the acoustic valley-like states results in an analogy of the valley Hall phase inversion.

5. Valley-Like Edge States

For a supercell consisting of a 20 × 1 unit cell with rigid boundaries at the top and bottom and Bloch periodic conditions at the left and right boundaries, the array of the unit cells can be regarded as a quasi-one-dimensional structure. Through adjusting the configuration of crystals and comparing relevant energy bands, the edge bands will be observed efficiently. As shown in Figure 3, the supercell with 20-unit cells of the same structure (crystal A) can be regarded as a trivial material, and the two overlapped edge bands emerge inside the bulk bandgap. By observing the mode fields, we can find that one is the top edge mode while the other is the bottom edge mode. No interface mode occurs in this trivial configuration.
We design a new supercell containing 20 × 1-unit cells with different structures, i.e., the upper 10 × 1-unit cells consist of crystal A while the lower 10 × 1-unit cells consist of crystal B. Different crystals at the opposite sides of the interface show different valley Hall phases. Four edge states can be observed in the bulk bandgap, as shown in Figure 1. These edge bands represent four different edge modes. Two of these modes are edge modes located at the top and bottom boundary which are similar to those of the supercell consisting of the same unit cells, while the other two modes are the interface modes. The acoustic valley-like edge states accompanied with clockwise and counterclockwise vortices can robustly lead to the propagation of energy flow in opposite directions.

6. Observation of Energy Propagation against Defects

Figure 4a shows a perfect structure consisting of a pair of valley-like insulators. Two different defects, i.e., cavities and disorders, are considered in Figure 4b,c, which are very common during the processing procedure. The cavity defects are induced by removing several unit cells at the interface between two different structures. Compared with that in Figure 4e, the disturbance can seriously affect the propagation ability of the band-gap-guiding transmission consisting of trivial insulators, but has little influence on the nontrivial valley-like design. The second type of defects is the disorders waveguide, which is induced by the random rotated scatters along the path of wave propagation. The amplitude of rotation is 20% of the upper limit of the reference angle, i.e., the design parameter θ in Equation (2). For any trivial guided edge mode shown in Figure 4, these defects can cause strong backscattering and even block wave propagation. Since cavities and disorders are not spin-mixing defects, the valley-like edge states along the interface will not be broken by these “nonmagnetic” impurities [57,58]. The simulated acoustic pressure fields shown in Figure 4b,c indicate that the wave incident from the left side can bypass all these defects at the valley-like protected interface and then efficiently propagate to the right side without significant attenuation. The valley-like protection has suppressed the backscattering of acoustic waves and ensured the robust transmission against various types of defects, including cavities and disorders.
As shown in Figure 5, we also investigate the transmission properties of the valley-like edge states. The black line in Figure 5b represents the transmission spectrum of the perfect configuration with no defect. The two typical defects mentioned in Figure 4b,c are represented by the red and blue lines, respectively. As a comparison, the transmission spectrum of a trivial structure formed by a single crystal (crystal A) is given in Figure 5b (the green line). The excitation sources are the same for all of these cases. By comparison, the transition of the trivial structure in the bulk bandgap is very low. In contrast, the valley-like insulators with a perfect configuration exhibit very high transmission. Moreover, high-transmission spectra can be observed in two imperfect interfaces with defects, suggesting that our valley-like structure has robust edge transmission properties against “non-magnetic” defects.

7. Conclusions

In conclusion, this paper shows an inverse design method for valley-like edge states. By quantifying the mode inversion error and maximizing the common bandgap, the inverse design can be transformed into an optimization problem. Instead of constructing Dirac cones at the momentum away from the high-symmetric points of the square lattice, the products are directly obtained after the break of mirror symmetries. Such a design has robust transmission properties against defects and can guide sound waves through the twisted path despite suffering two sharp corners (bent by 135° as shown in Figure 6). Compared to the literature [34], this method does not require trial-and-error design and achieves a larger operation bandgap (70% wider than the original design). Since the inverse strategy does not require certain physical fields, it can theoretically be applied to other classical systems such as optical systems and mechanical systems.

Author Contributions

J.Y. conceived and wrote the paper; Y.L., D.S., N.H. and H.N. revised the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Chinese National Natural Science Fund (Grant No.: U1864208), National Science and Technology Major Project (2017-VII-0011-0106), Science and Technology Planning Project of Tianjin (20ZYJDJC00030), the Key Program of Research and Development of Hebei Province (202030507040009), the Fund for Innovative Research Groups of Natural Science Foundation of Hebei Province (A2020202002) and the Key Project of Natural Science Foundation of Tianjin (S20ZDF077).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schaibley, J.R.; Yu, H.; Clark, G.; Rivera, P.; Ross, J.S.; Seyler, K.L.; Yao, W.; Xu, X. Valleytronics in 2D Materials. Nat. Rev. Mater. 2016, 1, 16055. [Google Scholar] [CrossRef]
  2. Ozawa, T.; Price, H.M.; Amo, A.; Goldman, N.; Hafezi, M.; Lu, L.; Rechtsman, M.C.; Schuster, D.; Simon, J.; Zilberberg, O. Topological Photonics. Rev. Mod. Phys. 2019, 91, 15006. [Google Scholar] [CrossRef]
  3. Zhang, X.; Xiao, M.; Cheng, Y.; Lu, M.; Christensen, J. Topological Sound. Commun. Phys. 2018, 1, 97. [Google Scholar] [CrossRef]
  4. Ma, G.; Xiao, M.; Chan, C.T. Topological Phases in Acoustic and Mechanical Systems. Nat. Rev. Phys. 2019, 1, 281–294. [Google Scholar] [CrossRef]
  5. Rycerz, A.; Tworzydło, J.; Beenakker, C. Valley Filter and Valley Valve in Graphene. Nat. Phys. 2007, 3, 172–175. [Google Scholar] [CrossRef]
  6. Xiao, D.; Yao, W.; Niu, Q. Valley-Contrasting Physics in Graphene: Magnetic Moment and Topological Transport Phys. Rev. Lett. 2007, 99, 236809. [Google Scholar] [CrossRef]
  7. Semenoff, G.W.; Semenoff, V.; Zhou, F. Domain Walls in Gapped Graphene. Phys. Rev. Lett. 2008, 101, 87204. [Google Scholar] [CrossRef]
  8. Gorbachev, R.V.; Song, J.; Yu, G.L.; Kretinin, A.V.; Withers, F.; Cao, Y.; Mishchenko, A.; Grigorieva, I.V.; Novoselov, K.S.; Levitov, L.S. Detecting Topological Currents in Graphene Superlattices. Science 2014, 346, 448–451. [Google Scholar] [CrossRef]
  9. Liu, Z.; Feng, W.; Xin, H.; Gao, Y.; Liu, P.; Yao, Y.; Weng, H.; Zhao, J. Two-Dimensional Spin–Valley-Coupled Dirac Semimetals in Functionalized Sbas Monolayers. Mater. Horiz. 2019, 6, 781–787. [Google Scholar] [CrossRef]
  10. Martin, I.; Blanter, Y.M.; Morpurgo, A.F. Topological Confinement in Bilayer Graphene. Phys. Rev. Lett. 2008, 100, 36804. [Google Scholar] [CrossRef] [Green Version]
  11. Qiao, Z.; Jung, J.; Niu, Q.; MacDonald, A.H. Electronic Highways in Bilayer Graphene. Nano Lett. 2011, 11, 3453–3459. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, F.; MacDonald, A.H.; Mele, E.J. Valley Chern Numbers and Boundary Modes in Gapped Bilayer Graphene. Proc. Natl. Acad. Sci. USA 2013, 110, 10546–10551. [Google Scholar] [CrossRef] [PubMed]
  13. Xu, X.; Yao, W.; Xiao, D.; Heinz, T.F. Spin and Pseudospins in Layered Transition Metal Dichalcogenides. Nat. Phys. 2014, 10, 343–350. [Google Scholar] [CrossRef]
  14. Mak, K.F.; McGill, K.L.; Park, J.; McEuen, P.L. The Valley Hall Effect in Mos2 Transistors. Science 2014, 344, 1489–1492. [Google Scholar] [CrossRef] [PubMed]
  15. Ju, L.; Shi, Z.; Nair, N.; Lv, Y.; Jin, C.; Velasco, J.; Ojeda-Aristizabal, C.; Bechtel, H.A.; Martin, M.C.; Zettl, A. Topological Valley Transport at Bilayer Graphene Domain Walls. Nature 2015, 520, 650–655. [Google Scholar] [CrossRef]
  16. Li, J.; Wang, K.; McFaul, K.J.; Zern, Z.; Ren, Y.; Watanabe, K.; Taniguchi, T.; Qiao, Z.; Zhu, J. Gate-Controlled Topological Conducting Channels in Bilayer Graphene. Nat. Nanotechnol. 2016, 11, 1060–1065. [Google Scholar] [CrossRef]
  17. Li, J.; Zhang, R.; Yin, Z.; Zhang, J.; Watanabe, K.; Taniguchi, T.; Liu, C.; Zhu, J. A Valley Valve and Electron Beam Splitter. Science 2018, 362, 1149–1152. [Google Scholar] [CrossRef]
  18. Deng, F.; Sun, Y.; Wang, X.; Xue, R.; Li, Y.; Jiang, H.; Shi, Y.; Chang, K.; Chen, H. Observation of Valley-Dependent Beams in Photonic Graphene. Opt. Express 2014, 22, 23605–23613. [Google Scholar] [CrossRef]
  19. Dong, J.; Chen, X.; Zhu, H.; Wang, Y.; Zhang, X. Valley Photonic Crystals for Control of Spin and Topology. Nat. Mater. 2017, 16, 298–302. [Google Scholar] [CrossRef]
  20. Bleu, O.; Solnyshkov, D.D.; Malpuech, G. Quantum Valley Hall Effect and Perfect Valley Filter Based On Photonic Analogs of Transitional Metal Dichalcogenides. Phys. Rev. B 2017, 95, 235431. [Google Scholar] [CrossRef] [Green Version]
  21. Ma, T.; Shvets, G. All-Si Valley-Hall Photonic Topological Insulator. New J. Phys. 2016, 18, 25012. [Google Scholar] [CrossRef]
  22. Noh, J.; Huang, S.; Chen, K.P.; Rechtsman, M.C. Observation of Photonic Topological Valley Hall Edge States. Phys. Rev. Lett. 2018, 120, 63902. [Google Scholar] [CrossRef]
  23. Gao, F.; Xue, H.; Yang, Z.; Lai, K.; Yu, Y.; Lin, X.; Chong, Y.; Shvets, G.; Zhang, B. Topologically Protected Refraction of Robust Kink States in Valley Photonic Crystals. Nat. Phys. 2018, 14, 140–144. [Google Scholar] [CrossRef]
  24. Chen, X.; Zhao, F.; Chen, M.; Dong, J. Valley-Contrasting Physics in All-Dielectric Photonic Crystals: Orbital Angular Momentum and Topological Propagation. Phys. Rev. B 2017, 96, 20202. [Google Scholar] [CrossRef]
  25. Xie, B.; Liu, H.; Cheng, H.; Liu, Z.; Chen, S.; Tian, J. Acoustic Topological Transport and Refraction in a Kekulé Lattice. Phys. Rev. Appl. 2019, 11, 44086. [Google Scholar] [CrossRef]
  26. Lu, J.; Qiu, C.; Ke, M.; Liu, Z. Valley Vortex States in Sonic Crystals. Phys. Rev. Lett. 2016, 116, 93901. [Google Scholar] [CrossRef] [PubMed]
  27. Ye, L.; Qiu, C.; Lu, J.; Wen, X.; Shen, Y.; Ke, M.; Zhang, F.; Liu, Z. Observation of Acoustic Valley Vortex States and Valley-Chirality Locked Beam Splitting. Phys. Rev. B 2017, 95, 174106. [Google Scholar] [CrossRef]
  28. Lu, J.; Qiu, C.; Ye, L.; Fan, X.; Ke, M.; Zhang, F.; Liu, Z. Observation of Topological Valley Transport of Sound in Sonic Crystals. Nat. Phys. 2017, 13, 369–374. [Google Scholar] [CrossRef]
  29. Lu, J.; Qiu, C.; Deng, W.; Huang, X.; Li, F.; Zhang, F.; Chen, S.; Liu, Z. Valley Topological Phases in Bilayer Sonic Crystals. Phys. Rev. Lett. 2018, 120, 116802. [Google Scholar] [CrossRef]
  30. Zhang, Z.; Tian, Y.; Wang, Y.; Gao, S.; Cheng, Y.; Liu, X.; Christensen, J. Directional Acoustic Antennas Based On Valley-Hall Topological Insulators. Adv. Mater. 2018, 30, 1803229. [Google Scholar] [CrossRef] [Green Version]
  31. Zhang, Z.; Tian, Y.; Cheng, Y.; Wei, Q.; Liu, X.; Christensen, J. Topological Acoustic Delay Line. Phys. Rev. Appl. 2018, 9, 34032. [Google Scholar] [CrossRef]
  32. Wang, M.; Ye, L.; Christensen, J.; Liu, Z. Valley Physics in Non-Hermitian Artificial Acoustic Boron Nitride. Phys. Rev. Lett. 2018, 120, 246601. [Google Scholar] [CrossRef] [PubMed]
  33. Ni, X.; Gorlach, M.A.; Alu, A.; Khanikaev, A.B. Topological Edge States in Acoustic Kagome Lattices. New J. Phys. 2017, 19, 55002. [Google Scholar] [CrossRef]
  34. Xia, B.; Zheng, S.; Liu, T.; Jiao, J.; Chen, N.; Dai, H.; Yu, D.; Liu, J. Observation of Valleylike Edge States of Sound at a Momentum Away from the High-Symmetry Points. Phys. Rev. B 2018, 97, 155124. [Google Scholar] [CrossRef]
  35. Vila, J.; Pal, R.K.; Ruzzene, M. Observation of Topological Valley Modes in an Elastic Hexagonal Lattice. Phys. Rev. B 2017, 96, 134307. [Google Scholar] [CrossRef]
  36. Chen, J.; Huo, S.; Geng, Z.; Huang, H.; Zhu, X. Topological Valley Transport of Plate-Mode Waves in a Homogenous Thin Plate with Periodic Stubbed Surface. Aip Adv. 2017, 7, 115215. [Google Scholar] [CrossRef]
  37. Huo, S.; Chen, J.; Huang, H.; Huang, G. Simultaneous Multi-Band Valley-Protected Topological Edge States of Shear Vertical Wave in Two-Dimensional Phononic Crystals with Veins. Sci. Rep. 2017, 7, 10335. [Google Scholar] [CrossRef]
  38. Liu, T.; Semperlotti, F. Tunable Acoustic Valley–Hall Edge States in Reconfigurable Phononic Elastic Waveguides. Phys. Rev. Appl. 2018, 9, 14001. [Google Scholar] [CrossRef]
  39. Pal, R.K.; Ruzzene, M. Edge Waves in Plates with Resonators: An Elastic Analogue of the Quantum Valley Hall Effect. New J. Phys. 2017, 19, 25001. [Google Scholar] [CrossRef]
  40. Yan, M.; Lu, J.; Li, F.; Deng, W.; Huang, X.; Ma, J.; Liu, Z. On-Chip Valley Topological Materials for Elastic Wave Manipulation. Nat. Mater. 2018, 17, 993–998. [Google Scholar] [CrossRef]
  41. Lera, N.; Torrent, D.; San-Jose, P.; Christensen, J.; Alvarez, J.V. Valley Hall Phases in Kagome Lattices. Phys. Rev. B 2019, 99, 134102. [Google Scholar] [CrossRef]
  42. Lu, J.; Qiu, C.; Xu, S.; Ye, Y.; Ke, M.; Liu, Z. Dirac Cones in Two-Dimensional Artificial Crystals for Classical Waves. Phys. Rev. B 2014, 89, 134302. [Google Scholar] [CrossRef]
  43. Chen, Z.; Ni, X.; Wu, Y.; He, C.; Sun, X.; Zheng, L.; Lu, M.; Chen, Y. Accidental Degeneracy of Double Dirac Cones in a Phononic Crystal. Sci. Rep. 2014, 4, 4613. [Google Scholar] [CrossRef] [PubMed]
  44. Li, Y.; Wu, Y.; Mei, J. Double Dirac Cones in Phononic Crystals. Appl. Phys. Lett. 2014, 105, 14107. [Google Scholar] [CrossRef]
  45. Dai, H.; Liu, T.; Jiao, J.; Xia, B.; Yu, D. Double Dirac Cone in Two-Dimensional Phononic Crystals Beyond Circular Cells. J. Appl. Phys. 2017, 121, 135105. [Google Scholar] [CrossRef]
  46. Kong, X.; Zhou, Y.; Xiao, G.; Sievenpiper, D.F. Spin-Momentum Locked Modes on Anti-Phase Boundaries in Photonic Crystals. Opt. Express 2020, 28, 2070–2078. [Google Scholar] [CrossRef]
  47. Li, J.; Wang, J.; Wu, S.; Mei, J. Pseudospins and Topological Edge States in Elastic Shear Waves. Aip Adv. 2017, 7, 125030. [Google Scholar] [CrossRef]
  48. Nanthakumar, S.S.; Zhuang, X.; Park, H.S.; Nguyen, C.; Chen, Y.; Rabczuk, T. Inverse Design of Quantum Spin Hall-Based Phononic Topological Insulators. J. Mech. Phys. Solids 2019, 125, 550–571. [Google Scholar] [CrossRef]
  49. Chen, Y.; Meng, F.; Jia, B.; Li, G.; Huang, X. Inverse Design of Photonic Topological Insulators with Extra-Wide Bandgaps. Phys. Status Solidi Rapid Res. Lett. 2019, 13, 1900175. [Google Scholar] [CrossRef]
  50. Christiansen, R.E.; Wang, F.; Sigmund, O. Topological Insulators by Topology Optimization. Phys. Rev. Lett. 2019, 122, 234502. [Google Scholar] [CrossRef] [Green Version]
  51. Christiansen, R.E.; Wang, F.; Sigmund, O.; Stobbe, S. Designing Photonic Topological Insulators with Quantum-Spin-Hall Edge States Using Topology Optimization. Nanophotonics 2019, 8, 1363–1369. [Google Scholar] [CrossRef]
  52. Du, Z.; Chen, H.; Huang, G. Optimal Quantum Valley Hall Insulators by Rationally Engineering Berry Curvature and Band Structure. J. Mech. Phys. Solids 2020, 135, 103784. [Google Scholar] [CrossRef]
  53. Dong, H.; Zhao, S.; Zhu, R.; Wang, Y.; Cheng, L.; Zhang, C. Customizing Acoustic Dirac Cones and Topological Insulators in Square Lattices by Topology Optimization. J. Sound Vib. 2021, 493, 115687. [Google Scholar] [CrossRef]
  54. Xia, B.; Wang, G.; Zheng, S. Robust Edge States of Planar Phononic Crystals Beyond High-Symmetry Points of Brillouin Zones. J. Mech. Phys. Solids 2019, 124, 471–488. [Google Scholar] [CrossRef]
  55. He, W.; Chan, C.T. The Emergence of Dirac Points in Photonic Crystals with Mirror Symmetry. Sci. Rep. 2015, 5, 8186. [Google Scholar] [CrossRef] [PubMed]
  56. Mei, J.; Wu, Y.; Chan, C.T.; Zhang, Z. First-Principles Study of Dirac and Dirac-Like Cones in Phononic and Photonic Crystals. Phys. Rev. B 2012, 86, 35141. [Google Scholar] [CrossRef]
  57. Wang, Z.; Chong, Y.; Joannopoulos, J.D.; Soljačić, M. Observation of Unidirectional Backscattering-Immune Topological Electromagnetic States. Nature 2009, 461, 772–775. [Google Scholar] [CrossRef]
  58. He, C.; Sun, X.; Liu, X.; Lu, M.; Chen, Y.; Feng, L.; Chen, Y. Photonic Topological Insulator with Broken Time-Reversal Symmetry. Proc. Natl. Acad. Sci. USA 2016, 113, 4924–4928. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) An illustration of the potential scatter described by a B-spline curve with 14 control points. (b) A C4-symmetric design with reduced independent parameters, i.e., r i = r i + 3 = r i + 6 = r i + 9 ( i = 1 , 2 , 3 ) . (c) A C3-symmetric design with constraints r i = r i + 4 = r i + 8 ( i = 1 , 2 , 3 , 4 ) .
Figure 1. (a) An illustration of the potential scatter described by a B-spline curve with 14 control points. (b) A C4-symmetric design with reduced independent parameters, i.e., r i = r i + 3 = r i + 6 = r i + 9 ( i = 1 , 2 , 3 ) . (c) A C3-symmetric design with constraints r i = r i + 4 = r i + 8 ( i = 1 , 2 , 3 , 4 ) .
Materials 15 06697 g001
Figure 2. An illustration of band inversion of a pair of valley-like insulators after the break of mirror symmetries near the original Dirac point. Accompanied with the reopen of the degenerate point, the fourth eigenmode of crystal A is exchanged with the fifth eigenmode of crystal B. The typical vortex feature can be observed in the time-averaged Poynting vectors marked with additional green arrows.
Figure 2. An illustration of band inversion of a pair of valley-like insulators after the break of mirror symmetries near the original Dirac point. Accompanied with the reopen of the degenerate point, the fourth eigenmode of crystal A is exchanged with the fifth eigenmode of crystal B. The typical vortex feature can be observed in the time-averaged Poynting vectors marked with additional green arrows.
Materials 15 06697 g002
Figure 3. (a) The energy band of a supercell consisting of 20 × 1-unit cells with the same structure (crystal A). (b) The eigenmodes (acoustic pressure field) localized at the top boundary with clockwise vortex and at the bottom with counterclockwise vortex. (c) The energy band of a supercell consisting of 20 × 1-unit cells with different structures (crystal A and crystal B). (d) The acoustic pressure field shows the edge states localized at the top boundary and the bottom boundary, as well as the interface states between the two different crystal structures. In (a,c), the red dotted line and the black solid line represent the edge states localized at the top boundary and the bottom boundary, respectively; the green hollow circle and the blue solid line indicate the two interface states induced by the valley-like insulators; the gray dotted lines represent the bulk bands. In (b,d), the arrows in the acoustic eigenmode fields represent the Poynting vectors; the additional green arrows indicate the location of vortex cores.
Figure 3. (a) The energy band of a supercell consisting of 20 × 1-unit cells with the same structure (crystal A). (b) The eigenmodes (acoustic pressure field) localized at the top boundary with clockwise vortex and at the bottom with counterclockwise vortex. (c) The energy band of a supercell consisting of 20 × 1-unit cells with different structures (crystal A and crystal B). (d) The acoustic pressure field shows the edge states localized at the top boundary and the bottom boundary, as well as the interface states between the two different crystal structures. In (a,c), the red dotted line and the black solid line represent the edge states localized at the top boundary and the bottom boundary, respectively; the green hollow circle and the blue solid line indicate the two interface states induced by the valley-like insulators; the gray dotted lines represent the bulk bands. In (b,d), the arrows in the acoustic eigenmode fields represent the Poynting vectors; the additional green arrows indicate the location of vortex cores.
Materials 15 06697 g003
Figure 4. (a) Simulated acoustic pressure field in the valley-like structures without defects. (b,c) Simulated acoustic pressure field with cavities and disorders, where the incident frequency of the acoustic wave is 16 kHz (within the bulk bandgap). (d) Simulated acoustic pressure field of the trivial band-gap-guiding structures without defects. (e,f) Simulated acoustic pressure field with cavities and disorders. The 4 green circles in subfigure (b,c,e,f) mark the location of defects.
Figure 4. (a) Simulated acoustic pressure field in the valley-like structures without defects. (b,c) Simulated acoustic pressure field with cavities and disorders, where the incident frequency of the acoustic wave is 16 kHz (within the bulk bandgap). (d) Simulated acoustic pressure field of the trivial band-gap-guiding structures without defects. (e,f) Simulated acoustic pressure field with cavities and disorders. The 4 green circles in subfigure (b,c,e,f) mark the location of defects.
Materials 15 06697 g004
Figure 5. (a) Schematic diagram for the investigation of transmission properties. (b) Transmission spectra of different structures. The green line represents the trivial design consisting of only one type of insulator (crystal A). The black, red, and blue curves represent the structures consisting of a pair of different insulators with no defects, with cavities and with disorders, respectively.
Figure 5. (a) Schematic diagram for the investigation of transmission properties. (b) Transmission spectra of different structures. The green line represents the trivial design consisting of only one type of insulator (crystal A). The black, red, and blue curves represent the structures consisting of a pair of different insulators with no defects, with cavities and with disorders, respectively.
Materials 15 06697 g005
Figure 6. Simulated acoustic pressure field along a curved path with sharp corners (bent by 135°).
Figure 6. Simulated acoustic pressure field along a curved path with sharp corners (bent by 135°).
Materials 15 06697 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, J.; Liu, Y.; Sun, D.; Hu, N.; Ning, H. Inverse Design of Valley-Like Edge States of Sound Degenerated Away from the High-Symmetry Points in a Square Lattice. Materials 2022, 15, 6697. https://doi.org/10.3390/ma15196697

AMA Style

Yang J, Liu Y, Sun D, Hu N, Ning H. Inverse Design of Valley-Like Edge States of Sound Degenerated Away from the High-Symmetry Points in a Square Lattice. Materials. 2022; 15(19):6697. https://doi.org/10.3390/ma15196697

Chicago/Turabian Style

Yang, Jishi, Yaolu Liu, Dongyang Sun, Ning Hu, and Huiming Ning. 2022. "Inverse Design of Valley-Like Edge States of Sound Degenerated Away from the High-Symmetry Points in a Square Lattice" Materials 15, no. 19: 6697. https://doi.org/10.3390/ma15196697

APA Style

Yang, J., Liu, Y., Sun, D., Hu, N., & Ning, H. (2022). Inverse Design of Valley-Like Edge States of Sound Degenerated Away from the High-Symmetry Points in a Square Lattice. Materials, 15(19), 6697. https://doi.org/10.3390/ma15196697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop