• Photonics Research
  • Vol. 9, Issue 4, 622 (2021)
Yu Zhang1、2、3, Langlang Xiong1、3, Meng Zhang1、3, and Xunya Jiang1、2、3、*
Author Affiliations
  • 1Engineering Research Center of Advanced Lighting Technology, Fudan University, Ministry of Education, Shanghai 200433, China
  • 2Department of Illuminating Engineering and Light Sources, School of Information Science and Engineering, Fudan University, Shanghai 200433, China
  • 3Institute of Future Lighting, Academy for Engineering and Technology, Fudan University, Shanghai 200433, China
  • show less
    DOI: 10.1364/PRJ.405230 Cite this Article Set citation alerts
    Yu Zhang, Langlang Xiong, Meng Zhang, Xunya Jiang, "Fractal topological band-gap structure induced by singularities in the one-dimensional Thue–Morse system," Photonics Res. 9, 622 (2021) Copy Citation Text show less

    Abstract

    The physical origin of the fractal topological band-gap structure in the one-dimensional Thue–Morse system has been revealed, which is characterized by the evolutions of two types of topological singularities with zero-scattering properties and the paths of phase vortex points, which are the mirrored paths of the first-type singularities. The field distribution of the upper and lower gap-edge states will interchange when the traditional gaps are closed and reopened. The topologically protected edge-states are found at both traditional gaps and fractal gaps. Our work broadens the topological properties of quasicrystals or aperiodic systems and provides potential applications in new optoelectronic devices.

    1. INTRODUCTION

    Topological band-gap (TBG) structures have received extensive attention for electronic and photonic systems [1,2]. Many fantastic topological phenomena are realized in photonic crystals (PhCs) based on the topological photonic band-gap (TPBG) [3,4], such as topological edge states [510], Weyl points and nodal lines [1113], one-way rotating states [14], logic gates [15], and sensors [16]. The TPBGs of one-dimensional (1D) systems are intensively studied, because the systems are simple and can be solved in strict theory. Recently, the connection between the surface impedance of 1D semi-infinite PhC and the geometry phase, also called the Zak phase, was revealed [17]. More importantly, in this original work, a simple method with the self-consistent gauge, based on the reflection coefficient of the semi-infinite 1D PhC models, is established to determine the singularities, which dominate the evolution of TPBG of 1D PhCs. With the introduction of synthetic dimensions or parameter freedoms to the 1D periodic systems, rich topological phenomena, such as the Weyl points and the nodal lines (surfaces) in high-dimensional spaces, are also observed and well studied [1113,18]. In the reflection phase map, there are phase vortex points (PVPs), but only the PVPs whose eigenvector turns to zero are singularities [11,12]. It is revealed that the different evolutions of singularities between bands correspond to different topological phase transitions [12,13].

    Different from the periodic systems, the quasicrystals (QCs) and aperiodic systems, which lack translational symmetry but possess long-range order, display Bragg diffraction spots, and also have complex band-gap structures. Extremely rich transport properties and complex band-gap structures of QCs also imply the rich topology beyond the periodic systems [19]. The first effort to strictly define the topological properties of band-gap structures is done by Kraus et al. [20,21] for 1D QCs, and the topologically protected edge-states are demonstrated experimentally in photonic systems [20,22,23]. The topology of QCs can also be directly retrieved from diffraction experiments [24]. More excitingly, very recently the phenomenon of fractal TBG (FTBG), which is not observed in the periodic systems, has been found in two-dimensional (2D) Floquet QCs and the Bott index originally to describe the topology of random systems is introduced into QCs to characterize the topological phase transition [25]. Characterized by a none-zero spin Bott index, the quantum spin Hall (QSH) effect is also realized in Penrose-type QC lattices, associated with robust edge-states and quantized conductance [26,27]. However, to date, unlike the periodic systems whose topological phase transitions are from the evolution of singularities, any relation between FTBG of QCs and the singularities has never been suggested nor has ever been explored. Without understanding the relation, we cannot answer certain questions, such as “what is the origin of FTBG and how does it evolve with the structural parameters of QCs?” or “what are the similarities and differences between TBG of periodic systems and QCs?”.

    Since there is lack of support from Bloch’s theorem, how waves (electrons, photons, or otherwise) are transmitted through QCs or aperiodic systems is not fully resolved to this day [25]. Generally, researchers have to deal with finite structures for the study of QCs, so that the study of finite-length PhCs is more instructive. Also, inspired by the studies of 1D periodic systems [12,13,18] and with the introduction of synthetic dimensions or parameter freedoms to QCs, more complex band-gap structures in higher dimensional space can possibly open a new window for us to understand the FTBG of QCs. The topological properties of the 1D periodic systems are signed by the singularities of the reflection coefficient, which is from the zero-scattering of each cell, resulting in the non-zero Zak phase or Chern number [12,13,17], even for finite periodic systems [28]. Then it is reasonable for us to make the assumption that the topological number of the 1D QCs or aperiodic systems can also be obtained from the singularities characterized by zero-scattering. Precise theory needs to be established based on the Bott index of 1D QCs and aperiodic systems, which has been proved to be equivalent to the Chern number in 2D translationally invariant systems [29]. In previous works [30,31], it is revealed that the band-gaps of QCs and aperiodic systems can be separated into two classes, the traditional gaps (TGs) and the fractal gaps (FGs), due to the scattering from the periodic interfaces and the multi-scale nested structure of interfaces, respectively. It is natural to think that the topological properties of TGs and FGs may be quite different because of their different scattering origin.

    In this work, we show that the FTBG features of 1D Thue–Morse (TM) systems with one synthetic dimension are generated by the evolutions of two types of singularities and PVPs. The former are characterized by the zero-scattering of the two basic scatterers that compose the TM systems, while the latter are the ordinary transmission resonance from Bragg scattering. Both TGs and FGs of TM systems are closed and reopened when the evolution paths of the first-type singularities pass through them. The second-type singularities only appear when the optical path ratio between two kinds of layers satisfies the ratio of two odd integers, and is indicated by the closing of TGs. A π-phase-jump of the reflection coefficient when a TG or FG is closed and reopened is the sign of the topological phase transition. Before and after topological phase transition, unlike the 1D periodic systems whose upper and lower gap-edge states will interchange both the spatial inversion symmetry (SIS) and the field distribution, only the latter will appear for TM systems because of the absence of SIS. The topologically protected edge-states are also found for both TGs and FGs. Although this work is done for 1D TM systems, we believe the similar generation mechanism of FTBG is also available for other 1D or higher dimensional QCs and aperiodic systems.

    2. STRUCTURE

    The 1D TM system can be generated by the recursive relationship Sn=Sn1+S˜n1, S1=AB, where the S˜n1 is obtained by exchanging A and B in the Sn1. For example, an S3 TM system can be represented as ABBABAAB [30]. The system is composed of the dielectric layers A and B with widths dA and dB, the relative permittivities ϵA=4 and ϵB=6.25, and the relative permeability μA=μB=1, respectively. Without loss of generality, we suppose that the background material outside the finite TM system is the same as the A-kind layer. So all B-kind layers can be thought as scatterers, and they are submerged inside the background material. The electric field in the ith layer can be written as Ei(z)=Ai+eiki(zzi)+Aieiki(zzi). The transfer matrix between layers can be expressed as Mi(Ai+,Ai)T=(Ai+1+,Ai+1)T. Ai+ and Ai are the coefficients of the forward and backward electric field in the ith layer. ki is the wave vector in the ith layer. Mi is the transfer matrix between layers. The coefficients of transmission and reflection and the field distribution can be obtained by the transfer matrix of a finite TM system. The central wavelength of the mth TG is 2m(lA+lB), where lA=ϵAdA and lB=ϵBdB are optical paths of layer A and B, respectively [31]. In this work, lA+lB is set as 600 nm, while the structural parameter is σ=lBlA+lB. Our study is focused on the frequency regions below the third TG.

    It should be noted that the band-gap structure becomes more and more detailed with the increasing order of TM systems. However, when the TM order m is large enough, e.g., m>4, the main band-gap structure of TM systems becomes stable, and the physics studied in this work is not influenced by the order of TM systems. So when we discuss the TGs and FGs, we generally choose the order m (even or odd) of the system larger than four, and it is irrelevant with the topology of band-gap structure of TM systems. The main difference between even-order and odd-order systems is the location of inversion center, which is only important when we use a certain order TM as a super-cell with a periodic boundary condition and study the topology of the band-gap of such periodic systems (see Appendix A).

    3. RESULTS AND DISCUSSION

    (a)–(c) Transmission spectra of the S7 TM system with the structural parameter σ=512,12,712, respectively. The central positions of the first three orders of TGs are labeled by dotted lines. The red dots and blue dots indicate two types of singularities. (d) The refractive index distribution (blue line) and |E(z)| (black line) of the PhC at singularity with frequency f=5×1014 Hz and the same parameters as in (b). (e) and (f) Refractive index distribution (blue lines) and |E(z)| (black lines) of the TM system at two types of singularities with frequencies f1=5×1014 Hz and f2=2.5×1014 Hz, respectively.

    Figure 1.(a)–(c) Transmission spectra of the S7 TM system with the structural parameter σ=512,12,712, respectively. The central positions of the first three orders of TGs are labeled by dotted lines. The red dots and blue dots indicate two types of singularities. (d) The refractive index distribution (blue line) and |E(z)| (black line) of the PhC at singularity with frequency f=5×1014  Hz and the same parameters as in (b). (e) and (f) Refractive index distribution (blue lines) and |E(z)| (black lines) of the TM system at two types of singularities with frequencies f1=5×1014  Hz and f2=2.5×1014  Hz, respectively.

    To understand the field distribution characteristic of singularities, the electric field |E(z)| (black lines) of states at singularities and the refractive index distributions (blue lines) are shown in Figs. 1(d)–1(f) for a 1D PhC and a TM system, respectively. The 1D PhC shown in Fig. 1(d) is composed with the same two kinds of layers as the TM system in Fig. 1(b). It is known that if we take B-kind layers as the scatterers in the background material of A-kind, the zero-scattering occurs at the singularities of 1D PhC [28]. Because of the zero-scattering of B-kind layers, the |E(z)| of a finite 1D PhC in Fig. 1(d) shows two special properties, (i) the flat field profile through the whole system, which means the amplitude of total field is unchanged, (ii) the flat field amplitude in A-kind layers, which means the field in A-kind layers only contains the forward propagating component. Figures 1(e) and 1(f) also show the similar zero-scattering properties at two types of singularities of TM systems. Because of the same frequency condition of the singularities, the |E(z)| shown in Fig. 1(e) has the same properties as those in Fig. 1(d). For the second-type singularities, the |E(z)| is more complex. First, the field profile is still flat through the whole system, which implies the zero-scattering property of the whole TM system. Second, from these A-kind layers with flat field amplitude, we can find that the scattering from S2=ABBA and S2˜=BAAB is zero for the state of the second-type singularities. Because any TM system whose order is greater than two can be decomposed into S2=ABBA and S2˜=BAAB, the scattering of the whole TM system should be zero for all second-type singularities.

    The physical origin and the field distribution of the first-type singularities can be explained by the transfer matrix method. Regarding the A-kind layers as the background material, we can find two kinds of scatterers, (A)B(A) and (A)BB(A) in the TM system. The elements of the transfer matrix for scatterer (A)B(A) are M11=M22*=cos(kBdB)+i2(kAkB+kBkA)sin(kBdB), M12=M21*=i2(kAkBkBkA)sin(kBdB). When the condition of the first-type singularities sin(kBdB)=0 is satisfied, the transfer matrix becomes I or I, where I is the unitary matrix. For scatterer (A)BB(A), we can insert an infinite thin A-kind layer between two B-kind layers, and then the transfer matrix of (A)BB(A) becomes I. Obviously, at such a TM system, the incident field can not sense the scattering of B-kind layers at all.

    The physical origin of the second-type singularities, which appear only when the mth TG is closed, can also been analyzed by similar method. For the central wavelength of the mth TG, the second-order trace map [3234] is x2=[112(ϵAϵB+ϵBϵA)]cos(4mσπ)+1+12(ϵAϵB+ϵBϵA).The sufficient and necessary condition of the mth TG being closed [31] is x2=2, and then we have σ=n2m,nN*. We emphasize that this condition is different from the mth-gap-closed condition of 1D PhC, which is σ=nm and also the condition of nodal line (surface) in the higher-dimensional space with synthetic dimensions [12]. The new condition causes several special properties of the TM system. Firstly, for the mth TG, there are 2m1 closing points with σ changing in the TM system, while there are m1 closing points in 1D PhC. Secondly, for the mth gap of the TM system, besides the m1 closing points which satisfy sin(kBdB)=0, there are m second-type singularities satisfying the condition σ=2n12m,nN*.The physical explanation of the condition for second-type singularities is the zero-scattering of S2=ABBA and S2˜=BAAB, which could be thought of as second-order units of the TM system. From the trace map theory [31], all traces of higher order units also satisfy xn=xn22(xn12)+2=2, which means the TGs are closed for all ordered TM systems when the condition x2=2 is satisfied. It is noted that the first TG can be closed for TM systems while it is always opened for 1D PhCs.

    We note that two types of topological singularities in this work have not been strictly demonstrated so far. The gap-closing-reopening is a phenomenon of topological transition, but could not be used as the judgement of topological transition or the existing of topological singularities. Zero-scattering is a much stronger support of the existing of topological singularities, but still not enough. The theoretical demonstration of two types of singularities could be done by counting the number change of resonant modes in upper and lower bands before and after the gap-closing-reopening, which can be regarded as the strict characteristic of topological phase transition for 1D aperiodic systems [35].

    (a) Transmission map of the S5 TM system in the parameter space {σ, f}. The blue (red) area indicates a gap (band). (b) The reflection phase map of the S5 TM system. The first-type singularities are located at the singular lines (white lines), which are characterized by the π-phase-jump. The black lines are the mirrored lines. Both of them form the basic structure of the FTBG feature. (c), (d) and (e), (f) Local maps of transmission and reflection phase around the first-type and second-type singularities, respectively. A series of PVPs are indicated by the white circles in (d). (g) Transmission distribution error graph of the S5 TM system after 500 times is calculated when the disorder intensity is 0.1.

    Figure 2.(a) Transmission map of the S5 TM system in the parameter space {σ, f}. The blue (red) area indicates a gap (band). (b) The reflection phase map of the S5 TM system. The first-type singularities are located at the singular lines (white lines), which are characterized by the π-phase-jump. The black lines are the mirrored lines. Both of them form the basic structure of the FTBG feature. (c), (d) and (e), (f) Local maps of transmission and reflection phase around the first-type and second-type singularities, respectively. A series of PVPs are indicated by the white circles in (d). (g) Transmission distribution error graph of the S5 TM system after 500 times is calculated when the disorder intensity is 0.1.

    For more details of the FTBG, we carefully studied two special cases of the transmission and reflection phase maps around the two types of singularities, which are shown in Figs. 2(c), 2(d) and Figs. 2(e), 2(f), respectively. In Figs. 2(c) and 2(d), an SL and an ML intersect at the central frequency of the second TG (signed by two white dashed lines). The second TG is closed since the SL cuts through it, and there is a topological transition with a π-phase-jump at two sides of SL. Along the ML represented by the black solid line in Fig. 2(d), we can find a series of PVPs indicated by the white circles which are corresponding to the ordinary resonant transmission for different systems. In the vicinity of the SL and ML intersection, a high-transmission in a wide frequency range is indicated by the white solid line in Fig. 2(c). When σ=12, the relative bandwidth of the transmittance higher than 0.9 can reach 30%. As illustrated as Fig. 2(g), the high-transmission region is also robust against 10% disorder strength introduced into dA, when the optical length lA+lB is kept as constant. Such a robust high transmission case in wide frequency range is quite rare in low-dimension systems, which can be used to design specific optoelectronic devices. In Figs. 2(e) and 2(f), there is a second-type singularity at the central frequency of the first TG when σ=12. Focusing on the frequency range of the first TG (signed by the white dashed lines) and increasing σ continuously near σ=12, a typical topological transition appears such that the first TG becomes thinner, totally closed, and then reopened. Exactly at the point with f=2.5×1014  Hz and σ=12, a π-phase-jump occurs which makes it possible to realize the topological edge-state in the first TG.

    (a) Local transmission map near the second TG of the S7 TM system in the parameter space. The blue (red) area indicates a gap (band). A (B) is the lower (upper) gap-edge state marked by a yellow dot (triangle) when σ=0.4. C (D) is the lower (upper) gap-edge state marked by a green dot (triangle) when σ=0.6. (b) The |E(z)| of state A (B) is indicated by the blue (red) lines. (c) The |E(z)| of state C (D) is indicated by the blue (red) lines. The black lines show the refractive index distribution, s and L means the total length of the system.

    Figure 3.(a) Local transmission map near the second TG of the S7 TM system in the parameter space. The blue (red) area indicates a gap (band). A (B) is the lower (upper) gap-edge state marked by a yellow dot (triangle) when σ=0.4. C (D) is the lower (upper) gap-edge state marked by a green dot (triangle) when σ=0.6. (b) The |E(z)| of state A (B) is indicated by the blue (red) lines. (c) The |E(z)| of state C (D) is indicated by the blue (red) lines. The black lines show the refractive index distribution, s and L means the total length of the system.

    (a) Structure composed of two S7 TM systems which are spliced together. The structural parameters of the left and right TM systems satisfy σl=35 and σr=25, respectively. The black dashed line indicates the interface, and the blue (yellow) areas indicate A-kind (B-kind) layers. (b) The transmission spectrum of the structure shown in (a). (c) The |E(z)| of the topologically protected edge-state with a frequency satisfying f=2.5×1014 Hz. (d) Edge-states in the space of the reflection phase ϕ and frequency f {ϕ, f} when the S5 TM system with σ=0.4 is connected to a reflector with an adjustable phase.

    Figure 4.(a) Structure composed of two S7 TM systems which are spliced together. The structural parameters of the left and right TM systems satisfy σl=35 and σr=25, respectively. The black dashed line indicates the interface, and the blue (yellow) areas indicate A-kind (B-kind) layers. (b) The transmission spectrum of the structure shown in (a). (c) The |E(z)| of the topologically protected edge-state with a frequency satisfying f=2.5×1014  Hz. (d) Edge-states in the space of the reflection phase ϕ and frequency f {ϕ, f} when the S5 TM system with σ=0.4 is connected to a reflector with an adjustable phase.

    In fact, to show the existing of a topologically protected edge-state, we choose any value of σ, as long as ZSR+ZSL=0 is satisfied in the gap. To better show that the edge-state is topological, we propose to connect the S5 TM system with σ=0.4 which is topologically non-trivial and an artificial reflector with adjustable reflection phase ϕ tuned from 0 to π. As illustrated as Fig. 4(d), the edge-states are represented by the blue circles for different ϕ, and the gray areas indicate the band. Obviously, an edge-state crosses the entire gap, because of the topology difference.

    (a) and (b) Maps of transmission and the reflection phase at the initial position of the S8 TM system in the parameter space. (c) Transmission spectrum near the edge-state corresponding to the structure composed of two S8 TM systems with σl=0.5458 and σr=0.6333. (d) The |E(z)| of the topologically protected edge-state with frequency f=4.156907×1014 Hz.

    Figure 5.(a) and (b) Maps of transmission and the reflection phase at the initial position of the S8 TM system in the parameter space. (c) Transmission spectrum near the edge-state corresponding to the structure composed of two S8 TM systems with σl=0.5458 and σr=0.6333. (d) The |E(z)| of the topologically protected edge-state with frequency f=4.156907×1014  Hz.

    Further theoretical and experimental research can be carried out from the following aspects. Firstly, the relationship between the topological number of QCs (the change of Bott index) and the evolution of singularities can be explored in detail. Secondly, the robust edge-states with high-Q can be experimentally demonstrated, and the broadband transparent feature caused by the topological properties of QCs can be widely used on optical or electromagnetic device design. Thirdly, the mechanism based on singularities could be extended to the high-dimensional QCs, such as the study of more complicated singularities and the characteristics of their evolution paths.

    4. CONCLUSION

    In conclusion, we study the FTBG feature of a 1D TM system in the structural parameter-frequency space, which is generated by the evolving paths of two types of topological singularities and the PVPs. The first-type singularities can evolve along π-phase-jump lines which satisfy the condition sin(kBdB)=0, while the paths of PVPs, which satisfy the condition sin(kAdA)=0, are the mirrored lines of the evolving paths of the first-type singularities. The second-type singularities appear in the mth TG when the structural parameter satisfies σ=2n12m,nN*. The self-similarity of the FTBG feature could be understood by the reappearing satisfaction of these conditions. Different from 1D PhCs, there is no symmetry-interchange at the topological transition in a TM system, while the FDI of the upper and lower gap-edge states will take place in a TM system before and after the TG closing, which can be regarded as a special feature of the topological phase transition in the QCs. The topologically protected edge-states with the exponential decay from the interface between two TM systems are realized at both TGs and FGs. These works will open new windows to understand the fractal topological band-gap structures.

    APPENDIX A: THE INVERSION CENTER AND SPACIAL INVERSION SYMMETRY FOR TM SUPER-CELLS WITH PERIODIC BOUNDARY CONDITIONS

    Generally, for a finite TM system, there is no inversion center due to the lack of translational symmetry. But, if we use a certain-order-TM structure as a super-cell and set a periodic condition on the super-cell, we can find two inversion centers for both odd- and even-ordered TM super-cells. We call such systems super-cell-TM systems. For an odd-order S2m+1(mN*) TM super-cell, the two inversion centers are the middle positions of the S2m and S?2m TM super-cells, while these, for even-order or S2m TM super-cells are the start point and the middle position, respectively. For example, the S4 TM system with periodic boundary conditions can be represented as |ABBABAAB|BAABABBA|, and the symbol | means one inversion center. Meanwhile, ?z+zinv=??z+zinv and Ez+zinv=E?z+zinv, where zinv means the coordinate of the inversion center, are automatically satisfied. The super-cell-TM systems have almost the same main gaps as the TM systems, and also show the gap-closing-reopening processes as TM systems with different σ. It has been demonstrated that the reflection phase in the gap range at one inversion center is from 0 to π or ?π to 0 in the binary 1D PhCs [17]. Similarly, under periodic boundary conditions, the reflection phase of a gap at the inversion center in the super-cell-TM system is also from 0 to π or ?π to 0. So the gaps can be divided into two types according to the different range of the reflection phase with different topological properties.

    It has been demonstrated that the changes of the spatial inversion symmetry of the upper and lower gap-edge states, can be regarded as a strong evidence of topologically nontrivial phase transition [12,13,17]. For the super-cell-TM systems, the reflection phase of the gap-edge can only be 0 or π. The total electric field at the inversion center can be represented as (1+eiθr)Ei, where θr and Ei are the reflection phase and incident wave, respectively. Choosing a particular inversion center, the gap-edge state is symmetric (S-state) with θr=0, while that is antisymmetric (A-state) with θr=0.

    Here we choose an S5 super-cell-TM system as an example where σ=0.125, 0.4, 0.6, 0.875, corresponding to four areas in the second TG of the TM system. The |E(z)| of the upper and lower gap-edge states are shown in Fig.?6. For the certain inversion center (red dash line), the symmetry of |E| of the upper and lower gap-edge states is exchanged before and after the gap is closed. For example, the upper (lower) gap-edge state is A-state (S-state) with σ=0.125, while the upper (lower) gap-edge state is S-state (A-state) with σ=0.4. In fact, in all processes of gap closing and reopening of TM systems caused by two types of singularities, whether it is TG or FG, this change of spacial inversion symmetry can be observed in the super-cell-TM systems as long as we select the appropriate inversion center. From the similar studies of quasi-crystals with super-cells [10,43], we can conclude that the gap-closing-reopening caused by two kinds of singularities of TM systems is really the topological phase transition.

    (a) and (b) Maps of transmission and the reflection phase at the initial position of the S5 TM system under the periodic boundary conditions in the parameter space. A–H represent the lower and upper gap-edge states for σ=0.125, 0.4, 0.6, 0.875, respectively. (c)–(j) The |E(z)| of the gap-edge states marked in (a). L is the total length of the S5 TM super-cell.

    Figure 6.(a) and (b) Maps of transmission and the reflection phase at the initial position of the S5 TM system under the periodic boundary conditions in the parameter space. A–H represent the lower and upper gap-edge states for σ=0.125, 0.4, 0.6, 0.875, respectively. (c)–(j) The |E(z)| of the gap-edge states marked in (a). L is the total length of the S5 TM super-cell.

    References

    [1] M. Z. Hasan, C. L. Kane. Colloquium: topological insulators. Rev. Mod. Phys., 82, 3045-3067(2010).

    [2] X.-L. Qi, S.-C. Zhang. Topological insulators and superconductors. Rev. Mod. Phys., 83, 1057-1110(2011).

    [3] L. Lu, J. D. Joannopoulos, M. Soljacic. Topological photonics. Nat. Photonics, 8, 821-829(2014).

    [4] T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zilberberg, I. Carusotto. Topological photonics. Rev. Mod. Phys., 91, 015006(2019).

    [5] Z. Wang, Y. D. Chong, J. D. Joannopoulos, M. Soljačić. Reflection-free one-way edge modes in a gyromagnetic photonic crystal. Phys. Rev. Lett., 100, 013905(2008).

    [6] Y. Poo, R.-X. Wu, Z. Lin, Y. Yang, C. T. Chan. Experimental realization of self-guiding unidirectional electromagnetic edge states. Phys. Rev. Lett., 106, 093903(2011).

    [7] B. Yan, J. Xie, E. Liu, Y. Peng, R. Ge, J. Liu, S. Wen. Topological edge state in the two-dimensional stampfli-triangle photonic crystals. Phys. Rev. Appl., 12, 044004(2019).

    [8] Y. Peng, B. Yan, J. Xie, E. Liu, H. Li, R. Ge, F. Gao, J. Liu. Variation of topological edge states of 2D honeycomb lattice photonic crystals. Phys. Status Solidi (RRL), 14, 2000202(2020).

    [9] Z. Guo, B. Yan, J. Liu. Straight lined and circular interface states in sunflower-type photonic crystals. J. Opt., 22, 035002(2020).

    [10] J. Zhao, S. Huo, H. Huang, J. Chen. Topological interface states of shear horizontal guided wave in one-dimensional phononic quasicrystal slabs. Phys. Status Solidi (RRL), 12, 1800322(2018).

    [11] Q. Wang, M. Xiao, H. Liu, S. Zhu, C. T. Chan. Optical interface states protected by synthetic Weyl points. Phys. Rev. X, 7, 031032(2017).

    [12] Q. Li, Y. Zhang, X. Jiang. Two classes of singularities and novel topology in a specially designed synthetic photonic crystals. Opt. Express, 27, 4956-4975(2019).

    [13] Q. Li, X. Jiang. Singularity induced topological transition of different dimensions in one synthetic photonic system. Opt. Commun., 440, 32-40(2019).

    [14] T. Hou, R. Ge, W. Tan, J. Liu. One-way rotating state of multi-periodicity frequency bands in circular photonic crystal. J. Phys. D, 53, 075104(2019).

    [15] R. Ge, B. Yan, J. Xie, E. Liu, W. Tan, J. Liu. Logic gates based on edge states in gyromagnetic photonic crystal. J. Magn. Magn. Mater., 500, 166367(2020).

    [16] A. Shi, R. Ge, J. Liu. Side-coupled liquid sensor and its array with magneto-optical photonic crystal. J. Opt. Soc. Am. A, 37, 1244-1248(2020).

    [17] M. Xiao, Z. Q. Zhang, C. T. Chan. Surface impedance and bulk band geometric phases in one-dimensional systems. Phys. Rev. X, 4, 021017(2014).

    [18] A. V. Poshakinskiy, A. N. Poddubny, M. Hafezi. Phase spectroscopy of topological invariants in photonic crystals. Phys. Rev. A, 91, 043830(2015).

    [19] A. I. Goldman, R. F. Kelton. Quasicrystals and crystalline approximants. Rev. Mod. Phys., 65, 213-230(1993).

    [20] Y. E. Kraus, Y. Lahini, Z. Ringel, M. Verbin, O. Zilberberg. Topological states and adiabatic pumping in quasicrystals. Phys. Rev. Lett., 109, 106402(2012).

    [21] Y. E. Kraus, O. Zilberberg. Topological equivalence between the Fibonacci quasicrystal and the Harper model. Phys. Rev. Lett., 109, 116404(2012).

    [22] M. Verbin, O. Zilberberg, Y. E. Kraus, Y. Lahini, Y. Silberberg. Observation of topological phase transitions in photonic quasicrystals. Phys. Rev. Lett., 110, 076403(2013).

    [23] M. Verbin, O. Zilberberg, Y. Lahini, Y. E. Kraus, Y. Silberberg. Topological pumping over a photonic Fibonacci quasicrystal. Phys. Rev. B, 91, 064201(2015).

    [24] A. Dareau, E. Levy, M. B. Aguilera, R. Bouganne, E. Akkermans, F. Gerbier, J. Beugnon. Revealing the topology of quasicrystals with a diffraction experiment. Phys. Rev. Lett., 119, 215304(2017).

    [25] M. A. Bandres, M. C. Rechtsman, M. Segev. Topological photonic quasicrystals: fractal topological spectrum and protected transport. Phys. Rev. X, 6, 011016(2016).

    [26] H. Huang, F. Liu. Quantum spin Hall effect and spin Bott index in a quasicrystal lattice. Phys. Rev. Lett., 121, 126401(2018).

    [27] H. Huang, F. Liu. Theory of spin Bott index for quantum spin Hall states in nonperiodic systems. Phys. Rev. B, 98, 125130(2018).

    [28] P. A. Kalozoumis, G. Theocharis, V. Achilleos, S. Félix, O. Richoux, V. Pagneux. Finite-size effects on topological interface states in one-dimensional scattering systems. Phys. Rev. A, 98, 023838(2018).

    [29] Y. Ge, M. Rigol. Topological phase transitions in finite-size periodically driven translationally invariant systems. Phys. Rev. A, 96, 023610(2017).

    [30] X. Jiang, Y. Zhang, S. Feng, K. C. Huang, Y. Yi, J. D. Joannopoulos. Photonic band gaps and localization in the Thue–Morse structures. Appl. Phys. Lett., 86, 201110(2005).

    [31] H. Lei, J. Chen, G. Nouet, S. Feng, Q. Gong, X. Jiang. Photonic band gap structures in the Thue–Morse lattice. Phys. Rev. B, 75, 205109(2007).

    [32] M. Kohmoto, B. Sutherland, K. Iguchi. Localization of optics: quasiperiodic media. Phys. Rev. Lett., 58, 2436-2438(1987).

    [33] N.-H. Liu. Propagation of light waves in Thue–Morse dielectric multilayers. Phys. Rev. B, 55, 3543-3547(1997).

    [34] X. Wang, U. Grimm, M. Schreiber. Trace and antitrace maps for aperiodic sequences: extensions and applications. Phys. Rev. B, 62, 14020-14031(2000).

    [35] T. A. Loring. A guide to the Bott index and localizer index(2019).

    [36] D. Meidan, T. Micklitz, P. W. Brouwer. Topological classification of adiabatic processes. Phys. Rev. B, 84, 195410(2011).

    [37] G. Bräunlich, G. M. Graf, G. Ortelli. Equivalence of topological and scattering approaches to quantum. Commun. Math. Phys., 295, 243-259(2010).

    [38] X. Huang, Y. Yang, Z. H. Hang, Z.-Q. Zhang, C. T. Chan. Geometric phase induced interface states in mutually inverted two-dimensional photonic crystals. Phys. Rev. B, 93, 085415(2016).

    [39] Q. Wang, M. Xiao, H. Liu, S. Zhu, C. T. Chan. Measurement of the ZAK phase of photonic bands through the interface states of a metasurface/photonic crystal. Phys. Rev. B, 93, 041415(2016).

    [40] Y. Wu, C. Li, X. Hu, Y. Ao, Y. Zhao, Q. Gong. Applications of topological photonics in integrated photonic devices. Adv. Opt. Mater., 5, 1700357(2017).

    [41] M. Xiao, G. Ma, Z. Yang, P. Sheng, Z. Zhang, C. T. Chan. Geometric phase and band inversion in periodic acoustic systems. Nat. Phys., 11, 240-244(2015).

    [42] G. Ma, M. Xiao, C. T. Chan. Topological phases in acoustic and mechanical systems. Nat. Rev. Phys., 1, 281-294(2019).

    [43] F. Baboux, E. Levy, A. Lematre, C. Gómez, E. Galopin, L. Le Gratiet, I. Sagnes, A. Amo, J. Bloch, E. Akkermans. Measuring topological invariants from generalized edge states in polaritonic quasicrystals. Phys. Rev. B, 95, 161114(2017).

    [44] A. R. Akhmerov, J. P. Dahlhaus, F. Hassler, M. Wimmer, C. W. J. Beenakker. Quantized conductance at the majorana phase transition in a disordered superconducting wire. Phys. Rev. Lett., 106, 057001(2011).

    [45] I. C. Fulga, F. Hassler, A. R. Akhmerov. Scattering theory of topological insulators and superconductors. Phys. Rev. B, 85, 165409(2012).

    [46] 46On the transmission map, the width, central frequency, and structure of the FGs are different at two sides of the SL, which is the track of the first-type singularity. If we choose two σ very near the SL, which correspond to the two TM systems, the reflection phase of the left (ϕl) and right (ϕr) TM systems can match well, which satisfies ϕl+ϕr=0. But, due to the slow decaying of the field magnitude, it is difficult to observe the topological edge-state between the two TM systems. However, if we choose σ a little farther away from the SL, due to the gap structure difference, it is difficult to satisfy ϕl+ϕr=0 at a certain frequency inside the two gaps.

    Yu Zhang, Langlang Xiong, Meng Zhang, Xunya Jiang, "Fractal topological band-gap structure induced by singularities in the one-dimensional Thue–Morse system," Photonics Res. 9, 622 (2021)
    Download Citation