Second-order photonic topological insulator with corner states

Second-order photonic topological insulator with corner states

Bi Ye Xie National Laboratory of Solid State Microstructures and Department of Materials Science and Engineering, Nanjing University, Nanjing 210093, China    Hong Fei Wang National Laboratory of Solid State Microstructures and Department of Materials Science and Engineering, Nanjing University, Nanjing 210093, China    Hai-Xiao Wang School of Physical Science and Technology, and Collaborative Innovation Center of Suzhou Nano Science and Technology, Soochow University, 1 Shizi Street, Suzhou 215006, China    Xue Yi Zhu National Laboratory of Solid State Microstructures and Department of Materials Science and Engineering, Nanjing University, Nanjing 210093, China    Jian-Hua Jiang School of Physical Science and Technology, and Collaborative Innovation Center of Suzhou Nano Science and Technology, Soochow University, 1 Shizi Street, Suzhou 215006, China    Ming Hui Lu National Laboratory of Solid State Microstructures and Department of Materials Science and Engineering, Nanjing University, Nanjing 210093, China    Yan Feng Chen National Laboratory of Solid State Microstructures and Department of Materials Science and Engineering, Nanjing University, Nanjing 210093, China

Higher-order topological insulators (HOTIs) which go beyond the description of conventional bulk-boundary correspondence, broaden the understanding of topological insulating phases. Being mainly focused on electronic materials, HOTIs have not been found in photonic systems yet. In this article, we propose a type of two-dimensional second-order photonic crystals with zero-dimensional corner states and one-dimensional boundary states for optical frequencies. All of these states are topologically non-trivial and can be understood based on the theory of topological polarization. Moreover, by tuning the easily-fabricated structure of the photonic crystals, we can realize different topological phases with unique topological boundary states straightforwardly. Our results can be generalized to higher dimensions and provides unprecedented venues for higher-order photonic topological insulators and semimetals.

42.70.Qs, 42.25.Bs, 78.20.Bh

Introduction.—Topological insulators (TIs) and topological semimetals (TSMs) have been theoretically and experimentally studied due to their distinct edge states and transport properties TI1 (); TI2 (). Normally, -dimensional (D) TIs have D gapped bulk states and D gapless boundary states. Recently, the concept of higher-order topological insulator (HOTI) has been put forward to describe those topological insulators (TIs) which have lower-dimensional gapless boundary states HOTI1 (); HOTI2 (); HOTI3 (); HOTI4 (); HOTI5 (); HOTI6 (). Generally speaking, a D TI with D, D, …, D gapped boundary states and D gapless boundary states is called the th-order TI. The HOTIs broaden the family of non-trivial topological insulating phases. Moreover, the HOTIs have unique boundary states which go beyond the conventional bulk-boundary correspondence and are characterized by novel topological invariants HOTI1 (); HOTI2 (); HOTI4 (); HOTI6 ().

However, it is not easy to realize these topological phases in electronic materials. One of the obstacles is that the Fermi levels of electronic materials are not always in the topologically non-trivial band gaps or at the gapless points. The band structures of photonic crystals (PCs) which have no relevance with Fermi level provide us with platforms to study various topological phases such as photonic topological insulators (PTIs) and photonic topological semimetals (PTSMs) without this obstacle PTI1 (); PTI2 (); PTI3 (); PTI4 (); PTI5 (); PTI6 (); PTI7 (); PTI8 (); PTI9 (); PTI10 (); PTI11 (); PTI12 (); PTI13 (); PTI14 (); PTI15 (); PTI16 (); PTI17 (); PTI18 (); PTI19 (); PTI20 (); PTSM1 (); PTSM2 (); PTSM3 (); PTSM4 (); PTSM5 (); PTSM6 (). In terms of the HOTIs, the topological corner and hinge states in PCs can be used to design unprecedented optical devices. So far, the observations of HOTIs are only realized in mechanical metamaterials HOTIEXP1 (), electrical circuits HOTIEXP2 (); HOTIEXP3 (), and coupled optical fiber arrays HOTIEXP4 () which are described predominantly by tight-binding models. The extension of the notion of HOTIs to PCs of which the photonic band structure is dominated by multiple Bragg scattering instead of tight-binding models is still lacking.

In this Letter, we propose a two-dimensional (2D) PTI which is the 2D photonic generalization of the Su-Schrieffer-Heeger (SSH) model SSH (). Previous studies of the 2D SSH model is focused on the zero Berry curvature of the topological phases and gapless 1D edge states 2DSSH0 (); 2DSSH1 (). Here, we demonstrate that there are coexisting edge and corner states when two topologically distinct PCs are placed together to form rectangular boundaries. We reveal that both the bulk polarization, described by the vector , and the edge polarization, described by the nested Wilson-loop eigenvalues and , are quantized by the mirror symmetries and . While the edge states are determined by the bulk topological invariants, the emergence of the corner states are dominated by the edge polarizations. The theoretical predictions and analysis are supported by numerical simulations for all-dielectric PCs at optical frequencies. Our study provides an unprecedented avenue for photonic HOTIs and higher-order topological semimetals.

Second-order topological PCs.—We consider a 2D PC with mirror symmetries as shown in Fig. 1(a). There are four identical dielectric rods in each unit-cell which form isotropic (with symmetry) or anisotropic (with symmetry) PCs, depending on their positions. Origin of the coordinates is set as the center of the unit-cell throughout this Letter where we consider the transverse-magnetic (TM) photonic modes. Photonic band structures of distinct topological indices can be realized by tuning the distance between the four rods. Despite such a simple design, topological edge states and corner states [as illustrated in Fig. 1(b)] can appear in the photonic band gap (PBG) between the first and the second bands (denoted as PBG I).

For isotropic PCs, we found that the intra- and inter-unit-cell distances ( and ) between the neighboring rods control the band topology. and where and are defined in Fig. 1(a). For the PBG I has trivial topology, whereas for the PBG I carries nontrivial topology. A transition appears at (i.e., as in Fig. 1), where the PBG is closed by a quadratic band touching [see Figs. 1(c)-1(e)].

The photonic band structure of our PC, as dominated by multiple Bragg scattering, cannot be quantitatively captured by tight-binding models. Nevertheless, the qualitative feature, particularly the topological properties, can be captured by the tight-binding model depicted in Fig. 1(a) of which the Hamiltonian is


Here , , , , and .

Figure 1: (a) The 2D photonic SSH model and the first Brillouin zone (BZ). We set , and the relative dielectric constant for all the cases. The coupling strength denoted as and . (c)-(e). The band inversion induced by changing . (c) (topologically trivial PBG). (d) (band gap closing). (e) (topologically nontrivial PBG). We choose the values of in order to make maximal overlap between the band gaps in (c) and (e). (b) The eigenmodes of a square supercell where (e)-configuration is inside a box and (c)-configuration is outside the box. There are four degenerate states localized at the four corners in the band gap.

Bulk and edge polarizations.—The tight-binding parameters, , reflect the intra- and inter-unit-cell distances between the neighboring rods, respectively. It is known that the above tight-binding model has nontrivial topology as characterized by the 2D polarization where


where , where , run over all occupied bands, is the periodic Bloch function for the th band. The 2D polarization is connected to the 2D Zak phase 2DSSH0 () via for . 2D Zak phase of the PC equals [i.e., ] for (i.e., ) and reveals that the PC is in topological non-trivial insulating phase, while it equals for (i.e., ) where the PC is in the trivial insulating phase. Beyond the isotropic tight-binding model, the bulk states with nontrivial topology include , , and . We will show that only the last case can lead to both topological edge and corner states, whereas the first two cases can only appear in anisotropic PCs which support nontrivial edge states but not corner states.

Now consider a box boundary between the PC with and the PC with [see Fig. 2(a)]. There are boundaries perpendicular to the and directions which support topological edge states due to the 2D Zak phase. That is, the nontrivial Zak phase for each gives rise to the edge states on the boundaries perpendicular to the direction for the whole edge Brillouin zone . Similarly, there are the edge states on the boundaries perpendicular to the direction due to nontrivial Zak phase .

Furthermore, those 1D edge states themselves are similar to 1D photonic SSH model with nontrivial topology protected by () for edges along the () direction. The topology of the edge states is characterized by the edge polarizations, and , for edges perpendicular to the and directions, respectively. The topological theory of polarization HOTI1 (); HOTI4 () connects the polarization of the edge states to the eigenvalues of the nested Wilson-loops (Mathematical details are given in the Supplemental Materials). In simple terms, the nested Wilson-loops construct the Zak phases of the edge states by projecting the bulk Bloch states into a sector which is topologically equivalent to the edge states using the nested Bloch states, . Here stands for summation over all bands below the PBG, and is the Wilson-loop eigenvector for the Wilson-loop with looping from to at fixed where labels the eigenvalue of the Wilson-loop. With the nested Bloch states, one can calculate the polarization for the edge states perpendicular to the direction,


where . Remarkably, for the PBG considered in this work, there is only one photonic band below the PBG. Hence, there is only one Wilson-loop eigenvalue and the eigenvector . Therefore, and similarly, . Numerical calculations of these topological indices are presented in the Supplemental Materials.

Bulk-edge-corner correspondence.— The topological edge states along the and directions meet at the corner of the box in a square supercell [see Fig. 2(a)]. The topological corner charge is determined by the edge polarizations as [Fig. 2(b)] HOTI4 ()


This quantized corner charge gives rise to a single corner state in each of the four corners [see Fig. 2(a)].

On the other hand, the emergence of the corner modes can also be understood via the tight-binding model. Particularly, when one consider the extreme case where and . In this case, the four corner sites of the tight-binding model is isolated from other parts of the lattice and become zero modes in the band gap. Since the band structure of this extreme case can be adiabatically connected to the case, they are in the same topological class and the corner states will always exist as long as . In terms of the topological invariants, this corresponding to phase.

The above analysis is confirmed by numerical calculation of the photonic eigenmodes for a box geometry illustrated in Fig. 2(a). The spectrum of the eigenmodes is already shown in Fig. 1(b) where four degenerate corner modes are found in the PBG. The electromagnetic fields of the corner modes are strongly localized at each of the four corner [Fig. 2(a) shows a superposition of the four corner states]. These four corner states are protected by the nontrivial topology, as characterized by the bulk and edge polarizations which are quantized by the mirror symmetries and . Therefore, in our photonic systems, there are coexisting edge and corner states protected by the band topology, which is an unique signature of HOTIs.

Figure 2: (a) There are in total four degenerate corner states at the four corners with frequency . The electrical fields are strongly localized at the corners (only one mode is shown here). (b) Schematic of bulk-edge-corner correspondence. The nontrivial topology of the bulk (cyan region) leads to the emergence of the edge states, while the polarization of the edge (indicated by the black arrows) results in the formation of the topological corner state (indicated by the shallow yellow region).

Anisotropic PCs and the topological phase diagram.—Next, we extend our discussions to anisotropic 2D PCs. It turns out that, according to our numerical study, the topological properties of the PCs as well as the bulk-edge-corner correspondence remain valid. The bulk and edge polarizations are determined by the intra- and inter-unit-cell distance between the rods along the and directions, and with . For instance, leads to nontrivial topological indices . There are totally different combinations. In fact, we can treat each combination as the combination of two independent 1D photonic SSH models. If the configurations in two directions are different from each other and are all in the topologically non-trivial phase, there are four degenerate edge states (DESs) and the degeneracy between two sets of 1D DESs are broken with each set only contains two DESs located in the opposite boundaries. We simulate this phenomenon in PC by requiring the same parameters , and as previous isotropic 2D photonic SSH model, however, with a deformed lattice structure as shown in Fig. 3(a). The result is shown in Fig. 3(b) which agrees with the theoretical predictions well. We note that there is also a single corner state for each corner for anisotropic PCs with .

Figure 3: Anisotropic 2D photonic SSH model. (a.) The unit cell of the PC which has the same , and as the isotropic 2D photonic SSH model but with deformed lattice structure characterized by and . Here we choose perfect electric conduction (PEC) as the boundary condition. The eigenstates are solved by COMSOL. (b). and . Non-trivial phase in both -and -directions. Two sets of DESs with different frequencies. For the DESs in the upper figure, the frequency is and for the DESs in the lower figure, the frequency is . (c). and . -direction in topologically non-trivial phase and -direction in topologically trivial phase. There are only 1D DESs in the left and right sides. The frequency is

Another interesting situation is that the anisotropic PC has topologically non-trivial configuration in one direction and topologically trivial configuration in the other direction. Without loss of generality, we assume -direction is the topologically non-trivial direction, i.e., . In this case, there are only two 1D edge states which locate at the left and right boundaries, since only the Zak phase along the direction is nontrivial [see Fig. 3(c)]. Moreover, the polarization of the edge states is trivial, because . Hence, the corner charge is and there is no corner state for the two phases with and .

According to the above discussions, the whole phase diagram for the mirror-symmetric 2D PCs can be classified into four different phases, as depicted in Fig. 4: the topologically trivial phase with for and of which and ; the two phases with nontrivial bulk topology but zero corner charge, i.e., or , and where or ; the phase with both nontrivial bulk topology and corner charge, i.e., and where for and . Therefore, by simply changing the relative distances of nearest rods in anisotropic 2D photonic SSH model, we can achieve various topological phases with different 1D DESs as well as corner states.

Figure 4: Classification of all kinds of generalized 2D photonic SSH model with different configurations in -and -directions. means 2D photonic SSH model is in topologically non-trivial phase along -(-) directions. means 2D photonic SSH model is in topologically trivial phase along -(-) direction. The blue area represents isotropic cases and yellow area represents anisotropic cases. The dashed lines represent there are no 1D edge states and solid lines represent the existence of 1D edge states. The solid circles means there are corner states. The edge (or corner) states of the same color have the same frequency, whereas the states with different colors have different frequencies.

Discussions.—(i). Our study predicts a simple realization of the second-order topological insulator in all-dielectric photonic crystals (PC) with corner states. The exotic corner states can be regarded as the 0D boundary states of the 1D edge of the 2D PC. Moreover, it is topologically protected by mirror symmetries and characterized by the bulk and edge topological polarizations which are coincidentally the same in our PCs. The conventional bulk-edge correspondence is extended to bulk-edge-corner correspondence in this system.

(ii). Besides, we study the anisotropic 2D photonic SSH model and find that the 1D DESs arise due to the 1D structures existing in 2D photonic SSH model. By adjusting the distances between the nearby rods in the -and -directions, the emergence of the edge states and corner states can be controlled straightforwardly. These topologically protected edge and corner states may be valuable for robust waveguides, optical couplers, and lasers.

(iii). If our theory is generalized to 3D photonic SSH model, there can be second-order and third-order topological insulating phases where topological hinge states and corner states can emerge respectively. 3D second-order topological semimetals may be achieved by stacking the 2D photonic SSH model along the direction. The 3D second-order topological semimetals have rich bulk and boundary properties which are yet to be explored HOTSM ().

Note added.— A recent paper showing 0D topological bound state in a 2D PC due to the presence of dislocation (a topological defect) has appeared dislocation (), demonstrating another mechanism for lower-dimensional topological light-trapping.

Acknowledgments.— B.Y.X, H.F.W, X.Y.Z, M.H.L and Y.F.C are supported by the National Key RandD Program of China (Grant No. 2017YFA0303700) and the National Nature Science Foundation of China (Grant No. 51721001) as well as the Academic Program Development of Jiangsu Higher Education (PAPD). H.X.W and J.H.J thank supports from the National Natural Science Foundation of China (No. 11675116).


  • (1) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010).
  • (2) X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
  • (3) F. D. M. Haldane and S. Raghu, Phys. Rev. Lett. 100, 013904 (2008).
  • (4) Z. Wang, Y. Chong, J. D. Joannopoulos, and M. Soljačić, Nature 461, 772-775 (2009).
  • (5) M. Hafezi, E. A. Demler, M. D. Lukin, and J. M. Taylor, Nat. Phys. 7, 907-912 (2011).
  • (6) Y. Poo, R. X. Wu, Z. Lin, Y. Yang, and C. T. Chan, Phys. Rev. Lett. 106, 093903 (2011).
  • (7) A. B. Khanikaev et al. Nat. Mater. 12, 233-239 (2013).
  • (8) M. Hafezi, S. Mittal, J. Fan, A. Migdall, and J. M. Taylor, Nat. Photon. 7, 1001-1005 (2013).
  • (9) M. C. Rechtsman et al. Nature 496, 196-200 (2013).
  • (10) L. Lu, J. D. Joannopoulos, and M. Soljačić, Nat. Photon. 8, 821-829 (2014).
  • (11) W.-J. Chen et al. Nat. Commun. 5, 6782 (2014).
  • (12) S. A. Skirlo, L. Lu, and M. Soljačić, Phys. Rev. Lett. 113, 113904 (2014).
  • (13) W. Gao, M. Lawrence, B. Yang, F. Liu, F. Fang, B. Beri, J. Li, S. Zhang, Phys. Rev. Lett. 114, 037402 (2015).
  • (14) L.-H. Wu and X. Hu, Phys. Rev. Lett. 114, 223901 (2015).
  • (15) T. Ma, A. B. Khanikaev, S. H. Mousavi, and G. Shvets, Phys. Rev. Lett. 114, 127401 (2015).
  • (16) C. He et al. Proc. Natl. Acad. Sci. USA 113, 4924-4928 (2016).
  • (17) D. Leykam, M. C. Rechtsman, and Y. D. Chong, Phys. Rev. Lett. 117, 013902 (2016).
  • (18) F. Gao et al. Nat. Commun. 7, 11619 (2016).
  • (19) L. Xu, H.-X. Wang, Y.-D. Xu, H.-Y. Chen, and J.-H. Jiang, Opt. Express 24, 18059 (2016).
  • (20) X.-C. Sun, C. He, X.-P. Liu, M.-H. Lu, S.-N. Zhu, and Y.-F. Chen, Prog. Quant. Electron. 55, 52 (2017).
  • (21) Xuan Zhu, Hai-Xiao Wang, Changqing Xu, Yun Lai, Jian-Hua Jiang, and Sajeev John, Phys. Rev. B 97, 085148 (2018).
  • (22) Y. Yang, Y. F. Xu, T. Xu, H.X. Wang, J. H. Jiang, X. Hu, Z. H. Hang, Phys. Rev. Lett. 120, 217401 (2018).
  • (23) L. Lu, L. Fu, J. D. Joannopoulos, and M. Soljačić, Nat. Photon. 7, 294-299 (2013).
  • (24) L. Lu, Z. Wang, D. Ye, L. Ran, L. Fu, J. D. Joannopoulos, and M. Soljačić, Science 349, 622 (2015).
  • (25) H.-X. Wang, L. Xu, H.-Y. Chen, and J.-H. Jiang, Phys. Rev. B 93, 235155 (2016).
  • (26) Meng Xiao, Qian Lin, and Shanhui Fan Phys. Rev. Lett. 117, 057401 (2016).
  • (27) H.-X. Wang, Y. Chen, Z. H. Hang, H.-Y. Kee, and J.-H. Jiang, npj Quantum Materials 2, 54 (2017).
  • (28) B. Yang et al. Science 359, 1013-1016 (2018).
  • (29) W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Science 357, 61-66 (2017).
  • (30) Z. Song, Z. Fang, and C. Fang, Phys. Rev. Lett. 119, 246402 (2017).
  • (31) J. Langbehn, Y. Peng, L. Trifunovic, F. von Oppen, and P.W. Brouwer, Phys. Rev. Lett. 119, 246401 (2017).
  • (32) W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Phys. Rev. B 96, 245115 (2017).
  • (33) F. Schindler, A. M. Cook, M. G. Vergniory, Z. Wang, S. S. P. Parkin, B. A. Bernevig, and T. Neupert, arXiv:1708.03636.
  • (34) M. Ezawa. Phys. Rev. Lett. 120 026801 (2018).
  • (35) M. Serra-Garcia, V. Peri, R. Susstrunk, O. R. Bilal, T. Larsen, L. G. Villanueva, S. D. Huber, Nature 555, 342 (2018).
  • (36) C. W. Peterson, W. A. Benalcazar, T. L. Hughes, and G. Bahl Nature 555, 346-350 (2018).
  • (37) S. Imhof, C. Stefan, F. Bayer, J. Brehm, L. Molenkamp, T. Kiessling, F. Schindler, C. H. Lee, M. Greiter, T. Neupert, R. Thomale, arXiv:1708.03647 (2017).
  • (38) J. Noh et al. arXiv:1611.02373 (2016).
  • (39) W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 42, 1698-1701 (1979).
  • (40) F. Liu, and K. Wakabayashi, Phys. Rev. Lett. 118, 076803 (2017).
  • (41) F. Liu, H. Y. Deng, and K. Wakabayashi. Phys. Rev. B 97, 035442 (2018).
  • (42) M. Lin and T. L. Hughes, arXiv:1708.08457.
  • (43) F.-F. Li et al. arXiv:1802.01811 (2018).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description