Topological crystalline protection in a photonic system

Topological crystalline protection in a photonic system

Jian-Xiao Zhang    Mikael C. Rechtsman    Chao-Xing Liu Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802-6300, USA
July 12, 2019

Topological crystalline insulators are a class of materials with a bulk energy gap and edge or surface modes, which are protected by crystalline symmetry, at their boundaries. They have been realized in electronic systems: in particular, in SnTe. In this work, we propose a mechanism to realize photonic boundary states topologically protected by crystalline symmetry. We map this one-dimensional system to a two-dimensional lattice model with opposite magnetic fields, as well as opposite Chern numbers in its even and odd mirror parity subspaces, thus corresponding to a topological mirror insulator. Furthermore, we test how sensitive and robust edge modes depend on their mirror parity by performing time dependent evolution simulation of edge modes in a photonic setting with realistic experimental parameters.

Introduction - Symmetry and topology are two fundamental mathematical tools in the classification of states of matter in condensed matter physics. Recently, intensive research interests have been focused on the role of symmetry in the classification of topological states, ever since the discovery of time reversal invariant topological insulatorsQi and Zhang (2010); Moore (2010); Qi and Zhang (2011); Hasan and Kane (2010); Kane and Mele (2005); Bernevig and Zhang (2006); Moore and Balents (2007); Fu and Kane (2007); Roy (2009), in which time reversal plays an essential role. A large variety of symmetry protected topological states have been identified theoretically for different symmetry classes and dimensions Ryu et al. (2010); Schnyder et al. (2008) In contrast, the corresponding material realization in experimentally feasible systems of these new topological states has only been limited in several symmetry classes, mainly for time reversal invariant topological insulatorsBernevig et al. (2006); König et al. (2007); Hsieh et al. (2009a); Zhang et al. (2009); Chen et al. (2009); Hsieh et al. (2009b); Chen et al. (2009); Hsieh et al. (2009c); Xia et al. (2009); Liu et al. (2008) , the quantum anomalous Hall insulatorsChang et al. (2013); Yu et al. (2010) , topological superconductorsMourik et al. (2012); Fu and Kane (2008); Das et al. (2012); Lutchyn et al. (2010); Sau et al. (2010); Alicea (2010) and topological mirror insulatorsFu (2011); Teo et al. (2008); Dziawa et al. (2012); Hsieh et al. (2012); Xu et al. (2012); Tanaka et al. (2012). Therefore, searching for new topological systems for symmetry protected topological states is vital for the development of this field.

Topological crystalline insulator (TCI) phasesFu (2011), topological phases that are protected by crystalline symmetry, can exist in a large number of crystal structures with different space group symmetries and a classification of TCIs in different point symmetry groups and space symmetry groups have recently been the focus of researchTanaka et al. (2012); Dziawa et al. (2012); Xu et al. (2012); Hsieh et al. (2012); Jadaun et al. (2012); Fang et al. (2013); Dong and Liu (2015); Liu et al. (2014) . Compared to quantum spin-Hall-type topological insulators, which can only occur for spinful fermions, topological crystalline insulators can exist in both fermionic and bosonic systems. This fact opens up the possibility to realize TCIs in various bosonic systems, including photonic, phononic, magnonic and cold atom systems. As we describe below, photonic systems in particular have generated a great deal of interest as a probe for topological physics.

The realization by Haldane and RaghuHaldane and Raghu (2008); Raghu and Haldane (2008) that gyromagnetic photonic crystals could exhibit non-trivial topological invariants opened the door to the new field of Òtopological photonicsÓLu et al. (2014) Ð in which the propagation of photons in a dielectric structure is protected in a similar sense to electrons in a crystal lattice. The first experimental realization of this phenomenon was made in the group of SoljačićWang et al. (2008, 2009) for the microwave regime. However, scaling the wavelength down to the optical regime (in order to realize topological states in optical devices) was not possible using this mechanism due to weak magnetic response in that frequency regime. Other mechanisms were proposedUmucalılar and Carusotto (2011); Hafezi et al. (2011); Fang et al. (2012); Khanikaev et al. (2013), and finally experimental demonstrations were madeRechtsman et al. (2013a) in a system based on evanescently-coupled helical waveguides; as well as in two-dimensional coupled ring resonatorsHafezi et al. (2013). While photonic topological protection is conceptually similar to that of electronic topological protection (after all, this phenomenon comes down to non-interacting wave dynamics), photonics offers unique advantages and potential applications. To name a few examples, photonic systems can be designed directly by fabrication (allowing any desired lattice structure to be realized); experiments can be carried out at room temperature, meaning any emergent devices can be brought to application more realistically than those that require very low temperatures; and the robustness associated with topological protection could be of use in an array of devices that rely on the flow of light (e.g., sensors, optical interconnects, electrooptic modulators, isolators, among others).

Very recently, a prediction was made that topological photonic crystals with surface states protected by the glide symmetry could be realized for microwave photonsLu et al. (2015) in a macroscopic ferrimagnetic structure. In this work, we propose the optical realization of crystalline symmetry protected boundary modes in photonic crystals phase using an entirely distinct mechanism from the microwave. The systems is quasi-one-dimensional, where an auxiliary parameter, , of which the Hamiltonian is a function, is used in place of a second spatial dimension. This is reminiscent of Ref. Kraus et al., 2012, in which a family of one-dimensional systems (defined by a pump parameter) were used to realize a photonic topological edge state and pump that are mapped to the integer quantum Hall effect. We note that another work Wu and Hu (2015) has predicted an analogue of a quantum spin Hall system that requires crystalline symmetry; however the edge always breaks that symmetry and so the topological edge states are not protected (and are gapped). The realization of crystalline symmetry protected topological states in the optical regime allows for a novel paradigm in exploring topological crystalline phases, opening up the possibility to explore the relationship between topological photonic crystals and nonlinear/interacting bosonic systems and quantum optical effects such as multi-photon quantum walks.

Tight-binding model - We start from a simple tight-binding model to illustrate our main idea and then simulate the system in a more realistic situation. We consider a quasi-one-dimensional (1D) chain along the direction with each unit cell consisting of four sites, denoted by , as marked by “” in Fig. 1A. The positions of four sites are chosen to preserve -directional mirror symmetry with respect to the line denoted by “” in Fig. 1A. Bound states can be induced by lowering on-site energy, forming a potential well on each site. We denote these states as , and the corresponding creation and annihilation operators as and , respectively. Thus, the tight-binding Hamiltonian for this system is given by

Figure 1: (A) Schematics for the quasi-1D chain model. There are four confining potentials, denoted by , in one unit cell, labeled by “”. Four sites in one unit cell are symmetric with respect to the mirror line “”. The intra-unit-cell hopping and inter-unit-cell hopping are labeled by “” and “”. (B) Schematics for the 3D visualization of the waveguide chain. The light propagates along direction. (C) Schematics of the quasi-2D lattice model which can be mapped from the 1D model. The mirror plane (orange) is parallel with the plane. The complex hopping terms along and directions between unit cells are symbolized by purple arrows. is conjugate with the parameter . (See Hamiltonian Eq. 6).

where and denoting the unit cell index. describes the on-site energy and the hopping between sites within the -th unit cell. is for the hopping between two nearest neighbor unit cells and . In the basis, the Hamiltonians and are given by


where is the on-site energy of each site in the -th unit cell, chosen to be identical. denotes the hopping strength along the direction and are hopping parameters between two unit cells, with . Only vertical and horizonal nearest-neighbor(NN) intra-site hoppings, denoted by “” in Fig. 1A, as well as horizonal and diagonal NN inter-site hoppings, denoted by “” in Fig. 1A, are considered.

The above Hamiltonian preserves the -direction mirror operator , given by


which interchanges the sites and ( and ). We set , , . Let us neglect , , and , and set for the moment. In this limit, the Hamiltonian in (Eq. 2) can be diagonalized and the eigen-energies of four eigenstates are shown in Fig. 2. We find that two degenerate states and with eigenenergy are well separated from other two states and with energies , where we use () to represent odd (even) mirror parity of .

Figure 2: The eigen energies and eigen wavefunctions for the Hamiltonian in one unit-cell with four sites. White and black colors indicate the sign for four different eigen wavefunctions (real).
Figure 3: (A) Energy spectrum for a finite chain () of our tight-binding model. Blue and red colors represent odd and even mirror parities of eigen wavefunctions, respectively. There are four bands in energy spectrum. The top and bottom bands come from the states and , while two bands in the energy range originate from the states . A mini-gap, as indicated by “”, exist within the bands. (B) Zoomed-in graph of energy spectrum around the energy range . Boundary modes , , , and , as indicated by “”, “”, “” and “”, are found within the mini-gap “”. “” (“”) denotes the crossing between “” and “” (“” and “”). The light red and light blue shading indicate the mini-gap region “” of two corresponding parities. (C) A gap opening for the crossing “” can be induced by a mirror-symmetry-breaking term , leading to an anti-crossing “”.

Thus, in the limit , we can focus on two degenerate eigenstates and , and project the total Hamiltonian onto the subspace spanned by these two states. The eigenstates of two degenerate states are given by , , and let us denote and to be the annihilation and creation operators for the states on the -th unit cell, respectively. The effective Hamiltonian after projection is given by


where the summation is over the parity subspaces. The parameters and are related to the original parameters , , and by and . If choose , , we find correspondingly. As a result, the Hamiltonian (Eq. 4) is nothing but two copies of Aubry-Andre-Harper (AAH) Hamiltonian Aubry and André (1980); Ganeshan et al. (2013); Kraus et al. (2012) with effective flux strength and in the odd and even subspaces, respectively. The parameter is also known as the pump parameter. Here shifts the index of the first unit cell, introducing a degree of freedom which will be discussed later. As a consequence of opposite mirror parities under , these two copies of Hamiltonian are decoupled from each other once mirror symmetry is preserved. It has been shown that by correctly choosing parameters, boundary modes can exist at the end of a finite chain described by the AAH model, within fractal sets of band gapsHofstadter (1976). Thus, we perform a direct calculation of energy spectrum of a finite chain for the model (Eq. 1). Here for demonstration, we choose the parameters to be , , , and . An extra offset between and is set by taking and . This offset shifts the relative energy level for odd and even bands, in order to emphasize the crossings of edge states inside the mini-gap described later. The four energy levels in one unit cell , , and expand into four energy bands as a function of the parameter , as shown in Fig. 3A. The and bands overlap with each other in the energy regime to the lowest order approximation, which is about for the parameters listed above. Each band is split into several sub-bands with mini-gap. We focus on the mini-gap in the energy range [] denoted as a in Fig. 3A. There are other mini-gaps, for example in the range of [] and []. Here we show the zoomed-in spectrum in Fig. 3B of the mini-gap region “”. We find four energy levels in the mini-gap, denoted by , , , and , or “”, “”, “” and “” in Fig. 3B, respectively, where L and R denote the position where the wave function of the state is localized in the chain, and indicates the parity of the boundary modes, shown by blue (odd) and red (even) in Fig. 3B. A more realistic calculation of the wave function in this system will be given in the next section, based on simulating the full continuum problem. The crossing between the boundary modes and , marked by “” in Fig. 3B, is topologically protected due to opposite mirror parities between them, while the crossing between and , marked by “” in Fig. 3B, is gapless because these two states are located at opposite boundaries of the chain. As long as the chain size is large compared with the penetration length of boundary mode, which is around twice the lattice constant for the parameters listed above, the overlapping between wavefunctions at opposite boundaries is negligible.

To test the symmetry protection by adding a mirror-symmetry breaking term onto , where


and , we find a gap opening between two edge modes (“”) and (“”), as marked by “” in Fig. 3C. Thus, we conclude that these boundary modes are stable only when mirror symmetry is present in the system (Eq. 1).

We would like to emphasize that although our model is written in one dimension, it can be mapped to a two-dimensional (2D) lattice model with complex hopping terms, corresponding to a 2D topological mirror insulator. This is a similar mapping as that made in Ref. Kraus et al., 2012. We may introduce a fictitious dimension and extend our 1D tight-binding model to the plane, forming a 2D square lattice as shown in Fig. 1C. Similar to the unit cell defined in the 1D model, each unit cell at the integer coordinates consists of 4 sites, each having one state described by the creation operators where . The 2D Hamitlonian can be written as


where is the hopping matrix between sites labeled by , in the adjacent unit cells along the () direction, and is the hopping matrix among sites in one unit cell. The detailed forms of and are given in the Appendix.

This Hamiltonian is also invariant under the mirror operation (Eq. 3) about the -axis. We define the basis creation operators in the same way as in Eq. 4. where are defined for the subspaces with odd and even mirror parities as before. In this new basis, the Hamiltonian takes the form


where the parameter represents the phase shift during the hopping along the direction. We notice that the form of corresponds to the Landau gauge for a magnetic field with the flux in one unit cell. Thus, a Chern number can be defined in each mirror subspace for a non-zero , as described in the appendix. The opposite signs for the phase shift in the mirror even and odd subspaces indicate that the total Chern number is canceled for the whole system. However, a mirror Chern number, as defined in Appendix, can characterize non-trivial topological property of our 2D systemHsieh et al. (2012); Kraus et al. (2012). Therefore, this mapping suggests that our tight-binding model provides a realization of 2D topological mirror insulators.

Numerical Simulations for photonic systems - We have now established topological mirror insulator phases in our simple 1D AAH type of tight-binding model and its relation to 2D topological mirror insulators is also illustrated above. The realization of this model in a photonic system requires more sophisticated numerical simulation, which includes continuum degrees of freedom and are based in Maxwells’ equations. Next, we describe the detailed experimental setup and perform a numerical simulation of photonic lattice systems.

To realize the tight-binding model in a realistic photonic system, we consider an array of evanescently coupled elliptical waveguides in fused silica glass. The difference in refractive index inside and outside the waveguide is utilized to confine the light in the plane and serve as a potential well, as shown in Fig. 1B. We follow the standard paraxial approximation for this type of waveguides, in which the Maxwells’ equation of light can be simplified to a Schrödinger-type equation, namely:


where is the envelope function of electric field, is the wavenumber of ambient light in the medium, and are respectively the background and the waveguide deviation from background refraction index. In comparison with the Schrödinger equation, one can see that the direction can be regarded as time, therefore the diffraction of light through the waveguides is equivalent to the time evolution of a particle in a potential determined by . This equation has been utilized to describe other topological phases in photonic latticesRechtsman et al. (2013a), as well as a wide variety of effects in linear and nonlinear opticsYariv (1997), including (but not limited to) the predictionChristodoulides and Joseph (1988) and observation Eisenberg et al. (1998); Fleischer et al. (2003) of lattice solitons, stable photorefractive solitons; Anderson localization in opticsSchwartz et al. (2007); conical diffractionPeleg et al. (2007); and optical pseudomagnetismRechtsman et al. (2013b). The waveguides used in this work can potentially be fabricated using the femtosecond direct laser writing technique; and are inspired by those described in Ref. Szameit and Nolte, 2010.

Figure 4: (A) and (C) show energy spectra for our waveguide chain I and II, respectively. Boundary modes, depicted by “” and “” in (A) and (C), exist in the mini-gap and the corresponding wave functions, as shown in (B) and (D), are localized at the left boundary and possess even and odd parities respectively. The light red and light blue bar to the right indicate the range of interested mini-gap of the respective parities mapped to Fig. 3B of the tight-binding model. A tight-binding model simulation using the same parameters are shown in the Appendix for comparison.

As shown in Fig. 1B, four elliptical waveguides reside on the corners of a rectangle to simulate one unit cell in our tight-binding model. We place these unit cells one by one to form a chain, as described in the tight-binding model. The confining potential is provided by the deviation of the refraction index inside the waveguides from the background, , and the hopping term is determined by the overlap between two nearby sites, which can be controlled by adjusting the distance between them. Details of the relation between distance of different sites and energy level positions are discussed in the Appendix.

Figure 5: Numerical simulation of time-dependent evolution for the injection of the light beam into the waveguide chains configuration I and II. The corresponding boundary mode parities and input beam parities are listed in the second and third column, respectively. Evolution time is expressed in unit of -direction coordinate. Density profiles for wavefunctions at 1, 25, 100 and 200 mm are presented. The dark parts in first three columns of the images are cropped.

For this waveguide chain with a finite number of unit cells (), we first discuss our numerical simulation of the energy spectrum, which is shown in Fig. 4A and C for two waveguide chains with different parameter sets (See appendix for details of the choice of parameters). For simplicity, we label the waveguide chain for the energy spectrum in Fig. 4A and C as numbers I and II, respectively. Here we only focus on energy bands that correspond to bands in our tight-binding model. The energy spectrum of our continuum also reveals a sub-band structure with mini-band gaps for both parity bands and a detailed comparison between continuum model and tight-binding model is shown in the Appendix. In short, we find that the continuum and tight-binding models strongly agree. For both the bands with even (red color) and odd (blue color) parities, they are split into three mini-bands, separated by mini-gaps. Within the mini-gap, boudary modes can exist for both parity bands. Examples of boudary modes are marked by a for even parity bands in Fig. 4A (waveguide configuration I) and b for odd parity bands in Fig. 4C (waveguide configuration II), both at the left edge of the chain. For these two states, the corresponding wave functions are shown in Fig. 4B and D, respectively, from which one can clearly see that both modes are highly localized with a penetration length around twice of the lattice constant. We emphasize that we need to carefully choose the value of in Fig. 4A (Fig. 4C), and at that value, only one boudary mode “” (“”) at that edge is supported and this is important for the experiments of parity selected time evolution of edge modes described below.

We next inject a wave packet into one end of the waveguide chain and study time evolution of wave packets with different parities in this chain by solving the time-dependent Schrödinger equation, Eq. 8, numerically. Technique details of our simulation are included in the appendix. The initial wave packet is chosen to be in Gaussian form, given by


where is a normalization factor, is the center of Gaussian form and chosen to be the center of the first two unit cells. The function describes the parity of the injected light beam and chosen to be for even parity mode and for odd parity mode. The size of the envelope function is tuned to cover the majority of the mode wavefunction. We inject an even beam into the waveguide chain I with (the value of for the mode “” in Fig. 4A), and the corresponding time evolution is shown in the first row of Fig. 5, from which one can see that the wave packet keeps localized after a long time evolution. In contrast, when we inject into the same waveguide configuration, the light spreads into the bulk of the waveguide chain, as shown in the second row of Fig. 5. For chain II with (the value of for the mode “” in Fig. 4C), the light will get delocalized with initial wave packet , but localized with initial wave packet , as shown in the last two rows of Fig. 5. Therefore, we conclude that the localization of edge modes at the boundary sensitively depends on their mirror parities, and our numerical simulation also provides an approach to probe topological crystalline protection of boundary modes in realistic experiments.

I Discussion and conclusion

In conclusion, based on a tight-binding model and continuum numerical simulation, we propose an experimental setup of a photonic lattice to realize topological crystalline protection of boundary modes in photonic systems. Numerical simulation also suggests a possible experimental configuration to detect topological edge modes, of which the edge localization behavior depends on the parities of both the injected wave packets and localized edge modes. In addition, we would like to emphasize that two degenerate states realized in our system can also be regarded as “pseudo-spin” and thus provide a natural platform to construct SU(2) Landau levels Li et al. (2013) and spin-orbit coupling in future studies. The fact that the model was simulated with both tight-binding and more experimentally realistic continuum simulations suggests that can be implemented in a scheme similar to previously realized time-reversal-broken topological phases in photonic lattices. This will open the door to the exploration of topological crystalline protection in interacting systems based on the nonlinear optical response of the ambient medium (giving rise to the nonlinear Schrödinger / Gross-Pitaevskii equation), and entangled quantum walks in such phases (based on injecting entangled photons into the structure). A question of central importance in TCI physics is the question of disorder: since edge state protection is achieved using a symmetry that is easily broken (i.e., crystalline/mirror symmetry) by disorder, what is the nature of the robustness of these states? Will all protection simply break down, or is it preserved in an ensemble-sense? These are fundamental questions where photonics provides a unique and versatile path forwards.


We would like to acknowledge helpful discussions with Ling Lu, Marin Soljačić, Rui-Xing Zhang, Jiabin Yu and Yang Ge. C.-X.L. acknowledges the support from Office of Naval Research (Grant No. N00014-15-1-2675). M.C.R. acknowledges the support of the National Science Foundation under grant number ECCS-1509546.

Appendix A Mapping to the 2D topological mirror insulator model

In this section of the appendix, we will describe the detailed form of 2D topological mirror insulator model, which is mapped from the 1D topological mirror insulator model in the main text. The 2D tight-binding Hamiltonian is given by



satisfying . Here plays a similar role as that in the 1D model. ensures the and states are away from the bands. Here are indices for four sites in one unit cell. is the identity matrix. We will show that is linked to a virtual external gauge field, and the parameters and here match the definition in Eq. 2.

We can project the Hamiltonian onto the subspace of basis, which is written as


where is for odd and even subspaces, respectively. The parameter could be interpreted as the phase shift due to an external gauge field, where the field directions for the odd and even subspaces are opposite. We may choose the phase factor as for two parity subspaces, which correponds to the Landau gauge of a magnetic field in the plane in both the even and odd subspaces.

Substituting the Fourier transform along the direction and , one can get


We find that takes the same form of Eq. 4. From the Fourier transform, the phase parameter in the original model plays a role of the momentum in the extra dimension .

Appendix B Numerical Simulation

In this section, we will describe the details of our numerical simulation of time dependent Schrödinger equation Eq. 8, as well as the parameter setup of our waveguide chains I and II. In the numerical simulation, the lattice constant is set to 28 m. The separations between center of sites in the same unit cell is roughly 12 m in and 13 m in direction. The other system parameters are the wavelength (nm); the change in refractive index, ; and (for silica). The system is discretized into a grid of pixels for each unit cell on average.

The shape of the potential well for each site is given by


where m, m. The parameter determines the depth of the potential well. Our numerical test shows that the on-site energy is linearly proportional to this depth within the energy range that we are interested in.

The overlap between the wave functions for neighboring sites, which contributes to the hopping parameter, drops exponentially as a function of distance, as shown in Fig. 6. From this numerical test, we find the hopping between neighboring sites dominates over the hopping along the off-diagonal direction (shown in the inset of Fig. 6), and thus we can neglect the off-diagonal hopping term between the sites 1 and 3 (2 and 4) and obtain the Hamiltonian Eq. 2. By controlling carefully the distance between different sites, four elliptical waveguides in one unit cell can indeed reproduce four eigen states with the desired energy levels as in our tight-binding model. By interpolation on Fig. 6, we can establish the relation between the parameters described in the tight-binding model and the realistic parameters, such as the position of the potential well and refraction index in each site, in the continuum simulation model. Thus we are able to map the configuration in numerical simulation to the parameters in our tight-binding model. The parameters of the two models fall in the same range if we take the unit of hopping parameters in the tight-binding model to be mm.

Figure 6: Hopping strength as a function of direction separation between two sites. The diagonal hopping strength is also calculated while keeping the hopping strength along direction to be by modulating direction separation. The magnitude of is negligible. The inset (A) shows the two hopping terms inside a unit cell.

If we focus on the structure of in Eq. 2 with the condition , the commutation relationship of suggests that there will be multiple boundary states within the band expanded by and . For an arbitrary , there could exist more than one boundary mode with different energies, as shown in Fig. 7. A similar argument can be applied to numerical simulations and results in the excitation of multiple modes when the input beam is injected on the edge. To find a specific where there is only one boundary mode occupying a boundary, an integer setting the index of the leftmost unit cell is introduced into the tight-binding model and it does not change the topological properties of the systemKraus et al. (2012). By carefully tuning , we are able to find the configurations with only one boundary mode of one mirror parity.

Parameter Value
/mm -0.8
/mm 0.1
/mm -0.275
/mm -0.3
/mm -0.13
/mm -0.062
22 for chain I, 0 for chain II
Table 1: Parameters in tight-binding model to generate Fig. 7. The values are fitted from the mapping between the tight-binding model and the continuum simulation.

With the mapping of parameters between tight-binding model and continuum simulations, we can also identify the required wave-guide chain configurations for our continuum simulations. The dispersion spectra for the chains I and II have been given for continuum model in Fig. 4A and C. As a comparison, we also show the spectra of the chains I and II in tight-binding model in Fig. 7. The positions of boundary modes for both parities match with the continuum simulations, although band widths and energies reveal minor difference. This is due to neglecting the long-range hopping and the correction of from the change of separations within nearby sites. Since we are only interested in boundary modes, this discrepancy is not important for our purpose. The detailed parameters used for the tight-binding model is listed in Table 1.

The normalized input beam Eq. 9 is treated as the initial wave function at time=0. For each timestep in the simulation, the real-space wave function is multiplied by the time evolution operator , where the potential part of Hamiltonian is discretized from the sum of the cubic-Gaussian potential well functions Eq. 13 and the Fourier transform of the kinetic part. The reflection of light on the boundary of the grid is manually reduced to emulate experimental conditions. The timestep is set to 5um/).

Figure 7: Tight-binding simulation of chain I (A) and chain II (B). Compared with Fig. 4A and C, the positions of boundary modes appear to be the same. The light red and light blue bar to the right indicate the range of interested mini-gap of the respective parities, corresponding to the gap in Fig. 3B.

Appendix C Topological invariant

In this section, we will describe how to define mirror Chern number in our system, which follows the definition of Chern number for a quasicrystal system in Ref. Kraus et al., 2012. For the quasicrystal system with one real dimension and one quasi-periodic term, one can define a periodic approximated Hamiltonian in a quasi-1D lattice with number of sites as


where the lattice constant is set to 1. is obtained by making the substitution in in Eq. 2 by , where is the floor function. is essentially an rational approximation of the parameter . Introducing implies a twisted boundary condition along direction.Niu and Thouless (1987) In the limit and , coincides with in Eq. 1. It is shown that for a finite large the topological properties of is the same as those of .

Next, we follow the Ref. Kraus et al., 2012 to define topological invariants of this system through the projector operator for parity , given by


where is the eigenstate of with energy , and the summation is performed over all the bands with parity below the energy . is the center energy of a gap. The Chern number for the bands with parity is defined by Avron et al. (1983); Kraus et al. (2012)




The above expression is equivalent to the famious Thouless-Kohmoto-Nightingale-Nijs formula Thouless et al. (1982) . With this definition of Chern number, the mirror Chern number can be defined as


which characterizes topological mirror insulator phase for a system with mirror symmetry.


  • Qi and Zhang (2010) X.-L. Qi and S.-C. Zhang, Physics Today 63, 33 (2010).
  • Moore (2010) J. E. Moore, Nature 464, 194 (2010).
  • Qi and Zhang (2011) X.-L. Qi and S.-C. Zhang, Reviews of Modern Physics 83, 1057 (2011).
  • Hasan and Kane (2010) M. Z. Hasan and C. L. Kane, Reviews of Modern Physics 82, 3045 (2010).
  • Kane and Mele (2005) C. L. Kane and E. J. Mele, Physical review letters 95, 146802 (2005).
  • Bernevig and Zhang (2006) B. A. Bernevig and S.-C. Zhang, Physical review letters 96, 106802 (2006).
  • Moore and Balents (2007) J. E. Moore and L. Balents, Physical Review B 75, 121306 (2007).
  • Fu and Kane (2007) L. Fu and C. L. Kane, Physical Review B 76, 045302 (2007).
  • Roy (2009) R. Roy, Physical Review B 79, 195322 (2009).
  • Ryu et al. (2010) S. Ryu, A. P. Schnyder, A. Furusaki,  and A. W. Ludwig, New Journal of Physics 12, 065010 (2010).
  • Schnyder et al. (2008) A. P. Schnyder, S. Ryu, A. Furusaki,  and A. W. Ludwig, Physical Review B 78, 195125 (2008).
  • Bernevig et al. (2006) B. A. Bernevig, T. L. Hughes,  and S.-C. Zhang, Science 314, 1757 (2006).
  • König et al. (2007) M. König, S. Wiedmann, C. Brüne, A. Roth, H. Buhmann, L. W. Molenkamp, X.-L. Qi,  and S.-C. Zhang, Science 318, 766 (2007).
  • Hsieh et al. (2009a) D. Hsieh, Y. Xia, D. Qian, L. Wray, F. Meier, J. Dil, J. Osterwalder, L. Patthey, A. Fedorov, H. Lin, et al., Physical review letters 103, 146401 (2009a).
  • Zhang et al. (2009) H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang,  and S. Zhang, CrossRef, CAS, Web of Science® Times Cited 801 (2009).
  • Chen et al. (2009) Y. Chen, J. Analytis, J.-H. Chu, Z. Liu, S.-K. Mo, X.-L. Qi, H. Zhang, D. Lu, X. Dai, Z. Fang, et al., Science 325, 178 (2009).
  • Hsieh et al. (2009b) D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. Dil, J. Osterwalder, F. Meier, G. Bihlmayer, C. Kane, et al., Science 323, 919 (2009b).
  • Hsieh et al. (2009c) D. Hsieh, Y. Xia, D. Qian, L. Wray, J. Dil, F. Meier, J. Osterwalder, L. Patthey, J. Checkelsky, N. Ong, et al., Nature 460, 1101 (2009c).
  • Xia et al. (2009) Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. Grauer, Y. Hor, R. Cava, et al., Nature Physics 5, 398 (2009).
  • Liu et al. (2008) C. Liu, T. L. Hughes, X.-L. Qi, K. Wang,  and S.-C. Zhang, Physical review letters 100, 236601 (2008).
  • Chang et al. (2013) C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, M. Guo, K. Li, Y. Ou, P. Wei, L.-L. Wang, et al., Science 340, 167 (2013).
  • Yu et al. (2010) R. Yu, W. Zhang, H.-J. Zhang, S.-C. Zhang, X. Dai,  and Z. Fang, Science 329, 61 (2010).
  • Mourik et al. (2012) V. Mourik, K. Zuo, S. Frolov, S. Plissard, E. Bakkers,  and L. Kouwenhoven, Science 336, 1003 (2012).
  • Fu and Kane (2008) L. Fu and C. L. Kane, Physical review letters 100, 096407 (2008).
  • Das et al. (2012) A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum,  and H. Shtrikman, Nature Physics 8, 887 (2012).
  • Lutchyn et al. (2010) R. M. Lutchyn, J. D. Sau,  and S. D. Sarma, Physical review letters 105, 077001 (2010).
  • Sau et al. (2010) J. D. Sau, R. M. Lutchyn, S. Tewari,  and S. D. Sarma, Physical review letters 104, 040502 (2010).
  • Alicea (2010) J. Alicea, Physical Review B 81, 125318 (2010).
  • Fu (2011) L. Fu, Physical review letters 106, 106802 (2011).
  • Teo et al. (2008) J. C. Teo, L. Fu,  and C. Kane, Physical Review B 78, 045426 (2008).
  • Dziawa et al. (2012) P. Dziawa, B. Kowalski, K. Dybko, R. Buczko, A. Szczerbakow, M. Szot, E. Łusakowska, T. Balasubramanian, B. M. Wojek, M. Berntsen, et al., Nature materials 11, 1023 (2012).
  • Hsieh et al. (2012) T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil,  and L. Fu, Nature communications 3, 982 (2012).
  • Xu et al. (2012) S.-Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian, I. Belopolski, J. Denlinger, Y. Wang, H. Lin, L. Wray, et al., Nature communications 3, 1192 (2012).
  • Tanaka et al. (2012) Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, T. Takahashi, K. Segawa,  and Y. Ando, Nature Physics 8, 800 (2012).
  • Jadaun et al. (2012) P. Jadaun, D. Xiao, Q. Niu,  and S. K. Banerjee, arXiv preprint arXiv:1208.1472  (2012).
  • Fang et al. (2013) C. Fang, M. J. Gilbert, S.-Y. Xu, B. A. Bernevig,  and M. Z. Hasan, Physical Review B 88, 125141 (2013).
  • Dong and Liu (2015) X.-Y. Dong and C.-X. Liu, arXiv preprint arXiv:1507.00241  (2015).
  • Liu et al. (2014) C.-X. Liu, R.-X. Zhang,  and B. K. VanLeeuwen, Physical Review B 90, 085304 (2014).
  • Haldane and Raghu (2008) F. Haldane and S. Raghu, Physical review letters 100, 013904 (2008).
  • Raghu and Haldane (2008) S. Raghu and F. Haldane, Physical Review A 78, 033834 (2008).
  • Lu et al. (2014) L. Lu, J. D. Joannopoulos,  and M. Soljačić, Nature Photonics  (2014).
  • Wang et al. (2008) Z. Wang, Y. Chong, J. D. Joannopoulos,  and M. Soljačić, Physical review letters 100, 013905 (2008).
  • Wang et al. (2009) Z. Wang, Y. Chong, J. Joannopoulos,  and M. Soljačić, Nature 461, 772 (2009).
  • Umucalılar and Carusotto (2011) R. Umucalılar and I. Carusotto, Physical Review A 84, 043804 (2011).
  • Hafezi et al. (2011) M. Hafezi, E. A. Demler, M. D. Lukin,  and J. M. Taylor, Nature Physics 7, 907 (2011).
  • Fang et al. (2012) K. Fang, Z. Yu,  and S. Fan, Nature photonics 6, 782 (2012).
  • Khanikaev et al. (2013) A. B. Khanikaev, S. H. Mousavi, W.-K. Tse, M. Kargarian, A. H. MacDonald,  and G. Shvets, Nature materials 12, 233 (2013).
  • Rechtsman et al. (2013a) M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. Podolsky, F. Dreisow, S. Nolte, M. Segev,  and A. Szameit, Nature 496, 196 (2013a).
  • Hafezi et al. (2013) M. Hafezi, S. Mittal, J. Fan, A. Migdall,  and J. Taylor, Nature Photonics 7, 1001 (2013).
  • Lu et al. (2015) L. Lu, C. Fang, L. Fu, S. G. Johnson, J. D. Joannopoulos,  and M. Soljačić, arXiv preprint arXiv:1507.00337  (2015).
  • Kraus et al. (2012) Y. E. Kraus, Y. Lahini, Z. Ringel, M. Verbin,  and O. Zilberberg, Physical review letters 109, 106402 (2012).
  • Wu and Hu (2015) L.-H. Wu and X. Hu, Physical review letters 114, 223901 (2015).
  • Aubry and André (1980) S. Aubry and G. André, Ann. Israel Phys. Soc 3, 18 (1980).
  • Ganeshan et al. (2013) S. Ganeshan, K. Sun,  and S. D. Sarma, Physical review letters 110, 180403 (2013).
  • Hofstadter (1976) D. R. Hofstadter, Physical review B 14, 2239 (1976).
  • Yariv (1997) A. Yariv, Optical electronics in modern communications, Vol. 1 (Oxford university press, 1997).
  • Christodoulides and Joseph (1988) D. Christodoulides and R. Joseph, Optics letters 13, 794 (1988).
  • Eisenberg et al. (1998) H. Eisenberg, Y. Silberberg, R. Morandotti, A. Boyd,  and J. Aitchison, Physical Review Letters 81, 3383 (1998).
  • Fleischer et al. (2003) J. W. Fleischer, M. Segev, N. K. Efremidis,  and D. N. Christodoulides, Nature 422, 147 (2003).
  • Schwartz et al. (2007) T. Schwartz, G. Bartal, S. Fishman,  and M. Segev, Nature 446, 52 (2007).
  • Peleg et al. (2007) O. Peleg, G. Bartal, B. Freedman, O. Manela, M. Segev,  and D. N. Christodoulides, Physical review letters 98, 103901 (2007).
  • Rechtsman et al. (2013b) M. C. Rechtsman, J. M. Zeuner, A. Tünnermann, S. Nolte, M. Segev,  and A. Szameit, Nature Photonics 7, 153 (2013b).
  • Szameit and Nolte (2010) A. Szameit and S. Nolte, Journal of Physics B: Atomic, Molecular and Optical Physics 43, 163001 (2010).
  • Li et al. (2013) Y. Li, S.-C. Zhang,  and C. Wu, Physical review letters 111, 186803 (2013).
  • Niu and Thouless (1987) Q. Niu and D. J. Thouless, Physical Review B 35, 2188 (1987).
  • Avron et al. (1983) J. Avron, R. Seiler,  and B. Simon, Physical review letters 51, 51 (1983).
  • Thouless et al. (1982) D. Thouless, M. Kohmoto, M. Nightingale,  and M. Den Nijs, Physical Review Letters 49, 405 (1982).
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