{'fraction_non_alphanumeric': 0.044847667069100206, 'fraction_numerical': 0.012570129962360628, 'mean_word_length': 4.056193895870736, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Second-order topological phases protected by chiral symmetry

We study second-order topological insulators and semimetals characterized by chiral symmetry. We investigate topological phase transitions of a model for construction of the two-dimensional second-order topological insulators protected only by chiral symmetry. By the theory of the phase transitions, we propose a second-order topological semimetal and insulators with flat hinge bands in chiral-symmetric three-dimensional systems. The three-dimensional second-order topological phases can be obtained from the stacked two-dimensional second-order topological insulators with chiral symmetry. Moreover, we show that broken chiral symmetry in the three-dimensional second-order topological phase allows a second-order topological insulator with chiral hinge states. We also demonstrate the second-order topological phases by using a lattice model.


I. INTRODUCTION
Higher-order topological insulators have recently drawn research interest as new topological crystalline phases . Unlike conventional first-order topological insulators, two-dimensional (2D) second-order topological insulators (SOTIs) have topologically protected corner states, and three-dimensional (3D) SOTIs have topological gapless modes on the hinges. In some crystalline SOTIs, the topological corner and hinge states can arise only from the nontrivial bulk topology, when the lattice termination is compatible with the crystal symmetries. The crystalline insulators with the nontrivial second-order topology, which are independent of the lattice termination, are called intrinsic second-order topological phases [6][7][8].
Intriguingly, stacking 2D SOTIs enables the coexistence of topological gapless surface and hinge states while the bulk gap is open. The 3D SOTI can be realized by two anticommuting mirror symmetries [39]. Furthermore, the idea of the crystalline SOTIs has been extended to semimetallic phases, which are called second-order topological semimetals (SOTSMs) [14-16, 20, 24, 39-41]. SOTSMs have not only topological hinge states but also topological gapless nodes in the bulk. The topological hinge states appear between the gapless points projected onto the hinges.
Such nontrivial phases are named extrinsic SOTIs [6][7][8], whose topology of the corner (hinge) states is determined by both the bulk and the edges (surfaces). The SOTI phases are realizable even if we arbitrarily terminate the lattice [6-8, 50, 51, 53]. The topological modes are stable as long as the bulk and the boundary bands are gapped.
However, we need to calculate a complicated boundary Hamiltonian to study extrinsic second-order topology in general. The phase diagram to search for the SOTI is difficult to explore analytically without employing crystal symmetry. Thus, it is useful to give a systematic way to construct extrinsic SOTIs and the explicit topological invariant for the phase diagram. Moreover, topological gapless nodes are characterized by a change of topological invariants for lower-dimensional insulators in momentum space. Therefore, we can produce new 3D second-order topological phases from the phase diagram for the 2D SOTIs.
In this work, we show various second-order topological phases realized by chiral symmetry. We discuss a simple method to construct the chiral-symmetric second-order topological phases. To do so, we develop the previous method [50,51] for the 2D SOTI protected only by chiral symmetry in view of the topological phase transition. From our theory, we also show different 3D second-order topological phases due to chiral symmetry, and one of them is an unconventional topological phase which has a single gapless point on the 2D surface in addition to the topological hinge states. Furthermore, it is found that the 3D SOTI with surface gapless points can become a 3D SOTI with chiral hinge states by broken chiral symmetry.
This paper is organized as follows. In Sec. II, we show the method for the construction of the second-order topological phases due to chiral symmetry. We demonstrate the topological phases by a lattice model constructed from our method in Sec. III. Our conclusion is summarized in Sec. IV.

II. SECOND-ORDER TOPOLOGICAL PHASES CHARACTERIZED BY CHIRAL SYMMETRY
In this section, we study SOTIs and a SOTSM characterized by chiral symmetry in terms of the bulk and edge gap closings for the topological phase transitions. We give a simple Hamiltonian to create the 2D SOTIs only by chiral symmetry and extend it to 3D systems by considering the translation symmetry.

A. Corner states and topological invariant
First of all, we review construction of 2D SOTIs in class AIII according to Ref. 50. We apply the method to 3D systems with chiral symmetry later. We consider the following bulk Hamiltonian: where H i (k i ) (i = x, y) are two Hermitian matrices with chiral symmetry represented as Π i . Here, Π 2 i = 1 i , and 1 i is the identity matrix with the same size as H i (k i ). Because {H i , Π i } = 0 is satisfied, the Hamiltonian in Eq. (1) has chiral symmetry Π = Π x ⊗ Π y . In this paper, we focus on gap closing at zero energy because chiral symmetry is present. We impose open boundary condition (OBC) with a right-angled corner [54]. When the model is gapped at zero energy, it can be characterized by a Z topological invariant given by [50,51] where w i=x,y are conventional winding numbers for H i (k i ). The expression of the winding number is given in the Appendix A. The winding numbers distinguish whether the one-dimensional (1D) bulk Hamiltonians H i (k i ) are in the first-order topological phase by chiral symmetry [55 -58]. As discussed in the next section, ν 2D is unchanged as long as the model is gapped at zero energy in the bulk and on the edges. The system with the nonzero ν 2D is a 2D SOTI with zero-energy corner states. Indeed, the nonzero ν 2D indicates the existence of the corner states when the edges break the translation symmetries in the x and the y directions. We consider a semi-infinite system with one corner for Eq. (1). Because k x and k y are separate in Eq. (1), the real-space Hamiltonian can be described by two real-space Hamiltonians of H x (k x ) and H y (k y ). We denote the Hamiltonians of H i (k i ) as H OBC i . The Hamiltonian under the corner boundary condition (CBC) can be represented as [54] Here, Π OBC y and 1 OBC x are representations of the chiral symmetry and the identity matrix in the terminated system, respectively. To see topological corner states, we assume that the bulk and the edges are gapped at zero energy. If both H x (k x ) and H y (k y ) are topologically nontrivial, i.e. ν 2D = 0, we can see the corner states with the zero energy as follows. Let φ zero i be one of the topological zero-energy eigenvectors of H OBC i at one zerodimensional edge from the nonzero w i . Then, we can find a zero-energy state given by φ zero x ⊗ φ zero y for H CBC in Eq. (3). By assumption, the zero-energy state should be a corner state. Generally, we can obtain |ν 2D | = |w x w y | topological corner states because H OBC i=x,y have |w i | zeroenergy modes. Thus, zero-energy corner states appear when ν 2D is nonzero.

B. Topological phase transitions and gap closing
To grasp the 2D SOTI, we revisit the topological phase transitions in view of the gap closing. We clarify how the topological phases for H(k) change when we continuously deform H i (k i ). In other words, we can see how the band gap necessarily closes in the bulk or the edges when ν 2D changes.
We assume that the system is in a trivial phase with ν 2D = 0 to elucidate the topological phase transitions. H x (k x ) and H y (k y ) need to close the band gap to change ν 2D = w x w y . To begin with, we discuss the bulk gap given by Eq. (1). Because {H y , Π y } = 0, the bulk Hamiltonian satisfies Therefore, if and only if H x (k x ) and H y (k y ) take zerovalued eigenvalues at the same time, the bulk bands close the gap. This gap closing can change both w x and w y , and thus ν 2D . By contrast, ν 2D can change even if H x (k x ) and H y (k y ) do not take zero-valued eigenvalues simultaneously. We consider a trivial phase with (w x , w y ) = (1, 0). From the above discussion, when the system enters a nontrivial phase with (w x , w y ) = (1, 1) from the trivial phase, the bulk Hamiltonian H(k) does not close the band gap. Then, we investigate an edge normal to the x direction. Namely, we see a semi-infinite system with the edge along the y direction. Because we retain periodicity in the y direction, the Hamiltonian under the edge boundary condition (EBC) can be described as Although H EBC x (k y ) satisfies a condition similar to Eq. (4), we note that H y (k y ) determines gap closing on the edge because H x (k x ) has the nonzero w x . H OBC x has the topological zero-mode φ zero x now, and H y (k y ) has Bloch eigenstates ψ ny (k y ) with the eigenvalues E ny (k y ). Therefore, an eigenvector of the H EBC x (k y ) can be obtained from φ zero H y (k y ) can effectively describe the edge states near zero energy thanks to the nontrivial topology of H x (k x ). Hence, the gap closes on the edge when H y (k y ) takes zero-valued eigenvalues. We also analyze a topological phase transition from the phase with (w x , w y ) = (0, 1) to the phase with (1,1). In this case, we show that gap closing happens on the edge normal to the y direction. The Hamiltonian with the edge boundary is Similarly, we can obtain an eigenstate of H EBC because this state is a zero-energy mode of H OBC y . Therefore, the gap closes on the edge through the eigenvalues of H x (k x ).

C. Construction of 3D SOTSMs and SOTIs
Hereafter, we generalize the theory about the 2D SOTI to 3D systems with chiral symmetry by adding translation symmetry in the z direction. Since the wave vector k z is added to Eq. (1) in the 3D systems, the bulk Hamiltonian is given by Thus, we can regard k z as a new parameter leading to the topological phase transition characterized by ν 2D (k z ). The topological phase transition can occur in momentum space, in which gap-closing points can emerge due to the chiral symmetry for the 3D system. Therefore, the phase has the topological gap-closing points and hinge states at zero energy. As a result, we can realize a 3D second-order topological phase by stacking the 2D SOTIs. The positions of the gapless points depend on how H x (k x , k z ) and H y (k y , k z ) give rise to the topological phase transitions. Because ν 2D (k z ) is periodic for k z , point nodes always appear in pairs. If H x and H y change the winding numbers w x and w y simultaneously, a SOTSM phase emerges in the bulk, as shown in Fig. 1(a). The system shows zero-energy hinge states between the nodes projected onto the hinges because of the nonzero ν 2D (k z ). For example, the gapless node appears when (w x , w y ) changes from (0, 0) to (1,1) in the 3D momentum space.
By contrast, the system can have gapless points on the surfaces. While the bulk is insulating, the surface gapless points coexist with topological hinge states. There are two types of topological phase transitions by the change in k z , which results in analogs of surface topological semimetals. We consider the first type where gapless points appear on the same surface. In this type, either w x or w y changes in momentum space while the other one is fixed to a nonzero value, for instance [ Fig. 1(b)]. In the second type, we can obtain gapless points on the different surfaces. The point nodes can be found when w x and w y change in momentum space. We note that this surface gapless structure is unique to the 3D systems because 2D bulk topological semimetals necessarily have topological point nodes in pairs in momentum space [58][59][60][61][62]. As an example, each surface can have a single gapless point if (w x , w y ) changes as (0, 1) → (0, 0) → (1, 0) → (1, 1) → (0, 1), as illustrated in Fig. 1(c). In both types, topological hinge states appear at the zero energy between the gapless points.
Next, we break chiral symmetry in the 3D SOTI with the gapless points on the surfaces. The gap opens at the point nodes because the symmetry protection is absent. Hence, we can discuss a surface Chern insulator due to massive Dirac cones [44,45,48,63]. In the specific case, we can easily diagnose the existence of the gapless hinge states (see Appendix B). Typically, if the pairs become gapped by broken chiral symmetry, they can contribute to the surface Chern numbers. The chiral hinge states appear between the gapped surfaces on the trivial bulk. Consequently, the system can show the 3D SOTI phase with chiral symmetry breaking [ Fig. 1(d)].

III. MODEL
We study a lattice model to demonstrate the secondorder topological phases realized by chiral symmetry, and confirm our theory in the previous section. The 2D lattice model is constructed from the Su-Schrieffer-Heeger (SSH) model [64], which can be transformed to the Benalcazar-Bernevig-Hughes model [1,2].

A. The 2D SOTI model
We construct second-order topological phases and the Hamiltonian using the general method from Eq. (1). We study a 2D tight-binding model described as wherex andŷ are the unit vectors in the x and y directions, respectively. t x , t ′ x , t y and t ′ y are real hopping parameters [see Fig. 2(a)]. We set the lattice constants to unity. The Hamiltonian in momentum space is where σ x,y,z and τ x,y,z are Pauli matrices acting on the sublattices, and σ 0 and τ 0 are the identity matrices. The energy eigenvalues are This model actually consists of the two SSH models represented by and the chiral symmetry is given by Thus, this 2D model has chiral symmetry Π = τ z ⊗ σ z . The topological invariant ν 2D can be calculated from the winding numbers w j for H j (k j ). They are given by Hence, we can obtain the phase diagram for the 2D model from Eq. (2) [ Fig. 2 (b)]. When the corner boundary condition is imposed and ν 2D = 1, the zero mode appears.
In the Appendix C, we show the analytic eigenvector localized at the corner in this model. We note that our model in Eq. (10) is topologically equivalent to the Benalcazar-Bernevig-Hughes model which can have topological corner modes protected by two anticommuting mirror symmetries [1,2]. Alternatively, we can reinterpret the zero-energy corner modes as a manifestation of the extrinsic second-order topology under the corner boundary condition in class AIII [51]. Therefore, the corner modes can remain stable topologically even after the mirror symmetries are broken. We investigate the stability by breaking the crystal symmetries in Appendix D.

B. Stacked SOTI model
We stack the 2D SOTI model to demonstrate 3D second-order topological phases. While we preserve the chiral symmetry, we introduce terms depending on k z . The dependence on k z is determined by how the 2D SO-TIs are stacked in the z direction.
First, we add two terms 2χ 1 cos k z τ x ⊗ σ z and 2χ 2 cos k z τ 0 ⊗ σ x to the Hamiltonian in Eq. (10). These terms effectively alter the hopping terms t x and t y to t x + 2χ 1 cos k z and t y + 2χ 2 cos k z , respectively. Hence, the topological phase transition can happen in the momentum space. When t x + 2χ 1 cos k z = ±t ′ x (|t y /t ′ y | ≤ 1) and/or t y + 2χ 2 cos k z = ±t ′ y (|t x /t ′ x | ≤ 1) can be satisfied, the k z planes have gapless points.
Suppose that t x = t y , t ′ x = t ′ y and χ 1 = χ 2 . In the parameter region, we can obtain a SOTSM with the bulk gapless points. If the parameters change as the red arrow in Fig. 3(a), two fourfold-degenerate points appear at (k x , k y , k z ) = (π, π, ± arccos t ′ x −tx 2χ1 ) in the bulk. The topological hinge states also exist between the projected point nodes, as shown in Fig. 3(b).
Next, we set χ 2 = 0 to see a SOTI with gapless points on the surfaces. Since t y is fixed, only w x can change in momentum space. Therefore, we can realize gapless points on the surface perpendicular to the y direction. For example, we change t x as the green arrow in the diagram. w x and ν 2D change when t x + 2χ 1 cos k z = t ′ x . Figure 3(c) shows the two surface nodes at (k x , k z ) = (π, ± arccos t ′ x −tx 2χ1 ) and the zero-energy hinge states be-  ). The bulk has the point nodes at kz = π/2 and 3π/2 on the line. The top panel in (c) shows bands for the surfaces normal to the y direction. The surface bands have nodal points at (kx, kz) = (π, π/2) and (π, 3π/2). In (d), the blue and the green bands represent electronic structures for the surfaces normal to the y and the x directions, respectively. The gapless point in the blue (green) bands is located at (kx, kz) = (π, 5π/3) ((ky, kz) = (π, 5π/6)).
tween the two points projected to the hinges.
Second, we add the new terms, (2η cos k z + δt x )τ x ⊗ σ z and (2η sin k z + δt y )τ 0 ⊗ σ x , instead of the previous two terms. The hopping t x (t y ) is effectively regarded as t x + δt x + 2η cos k z (t y + δt y + 2η sin k z ). If t ′ x = t ′ y , the modified hoppings form the circle in the diagram, as depicted in Fig. 3(a). Thus, the gap closing can occur on the different surfaces when the circle encloses a phase transition point accompanied by bulk gap closing. We assume that the circle encloses the point t x /t ′ x = t y /t ′ y = 1, and that η/t ′ x > 0. Then, the point nodes appear at (k x , k z ) = (π, − arccos ) and (k y , k z ) = (π, π − arcsin t ′ y −ty−δty 2η ) on the surfaces normal to the y and the x directions, respectively [ Fig. 3 (d)]. The anomalous distribution of the point nodes is allowed because the gap closings both in the bulk and on the surfaces change the hinge topology, similar to a mirror-symmetric topological quadrupolar semimetal [39]. To annihilate the surface point nodes, the gap should close in the bulk.
We emphasize that our Hamiltonian does not need any crystal symmetries to realize the second-order topological phases. The 2D chiral-symmetric model can also be characterized by the two mirror symmetries for quantization of the nested Wilson loop [1,2,39]. However, even if the mirror symmetries are broken, the hinge flat-band states can survive as long as the chiral symmetry exists. We can indeed see the protection of the chiral symmetry by breaking the mirror symmetries. For example, we introduce the term 2t w1 cos k z τ y ⊗ σ z + 2t w2 cos k z τ 0 ⊗ σ y for the case in Fig. 3(c). If t w1,2 ≪ t ′ x,y , the topologi-cal zero-modes are stable because ν 2D does not change [ Fig. 4 (a)].

C. SOTIs with broken chiral symmetry
Finally, we break chiral symmetry to produce a 3D SOTI from the 3D SOTI with the surface gapless points. The system belongs to class A. By using our model, we here show that the SOTI can be realized easily from the surface gapless points in Fig. 3 (c). We add the perturbation 2t m sin k z τ z ⊗ σ z in order to obtain the SOTI without chiral symmetry. Figure 4 (b) shows the gapless hinge states which originate from the zero-energy states.
To understand the topological phase transition, we investigate the surface normal to the y direction from the discussion in Sec.II B. Because H y is topologically nontrivial, we can obtain the surface effective Hamiltonian near the zero energy by expanding H xz = (t x + t ′ x cos k x + 2χ 1 cos k z )τ x + t ′ x sin k x τ y + 2t m sin k z τ z (see also Appendix B). The effective Hamiltonian describes the surface Dirac points with the mass induced by the broken chiral symmetry.
We can calculate the surface Chern number from the surface Hamiltonian. The effective Hamiltonian near the gapless points k ± ≡ (π, ± arccos where q = (q x , q z ) is the wave vector measured from k ± . By the perturbation term, the Dirac cones at k ± obtain the mass gap described as 2t m sin k ± z τ z . Because sin k + z > 0 and χ 1 /t ′ x > 0 in the calculation, the surface has a Chern number of −sgn(t m ). As a result, the system becomes the SOTI phase, which is consistent with the discussion in Appendix B.

IV. CONCLUSION AND DISCUSSION
In the present paper, we have shown that various chiral-symmetric second-order topological insulators and semimetal are realizable. We have suggested the method to construct the second-order topological phases protected by chiral symmetry. The theory reveals that there are several types of chiral-symmetric second-order topological phases in three-dimensional systems, and that the topological phases can be obtained from stacked twodimensional second-order topological insulators with chiral symmetry. The zero-energy hinge states are stable topologically by chiral symmetry. Particularly, we have discovered the second-order topological insulator, which has not only hinge bands but also a single gapless point on the surface thanks to chiral symmetry. In general, this surface gapless point is not allowed in the two-dimensional bulk because both the bulk and surfaces need to contribute topologically. We have also shown that broken chiral symmetry can yield the second-order topological insulator from the surface gapless points in the chiral-symmetric second-order topological insulators. Our method tells us how to give extrinsic second-order topology easily in classes A and AIII.
Importantly, if any Hamiltonian with chiral symmetry can be continuously deformed to Eq. (1) or (8) without the bulk and surface gap closings, we can distinguish whether the system has nontrivial second-order topology because the band topology is unchanged [50,51]. Namely, our theory is useful even though the Hamiltonian does not take the form in Eq. (1) or (8). Thus, our model is available for photonic and phononic crystals and electric circuits [26][27][28][29][30] because their structures are flexibly controllable. Theoretically, corner modes in class BDI are suggested in photonic systems composed of waveguides and of coupled ring resonators with negative hopping amplitudes [47]. We also expect that our secondorder topological phases can be constructed in waveguide arrays, and that the corner modes are detectable through measurement of the outgoing intensity distribution of light injected into a waveguide. Moreover, electric circuits have the possibility of experimental realization of the second-order topological phases protected by chiral symmetry, because capacitors and inductors can provide band structures of the circuit Laplacian [26,65,66]. The nontrivial second-order topological modes in the electric circuits are measurable by the impedance response since the zero modes contribute to large resonances at the boundaries such as corners. Therefore, our model is feasible in systems with the easily controllable structures.
Additionally, topological stability of gapless nodes is characterized locally by the change of a topological invariant in momentum space [58][59][60]. Therefore, the nodal structures in our paper are allowed in other classes with chiral symmetry. Actually, second-order topological phases have been predicted in superconductors [6-8, 10, 67-78]. Some of the topological superconductors are understandable in view of extrinsic second-order topology [77,78]. Recently, robust Majorana hinge states were studied in a second-order topological superconductor as long as the particle-hole symmetry is preserved [78]. Thus, the extrinsic second-order topological nodal phases protected by chiral symmetry are expected to appear in superconductors with time-reversal symmetry. diagonal, the Q matrix can be rewritten as By using the off-diagonal component q(k), the winding number is defined as [56][57][58] The winding number w corresponds to the number of topological edge states in the 1D chiral-symmetric system.
Appendix B: Model for chiral hinge states In this appendix, we give a 3D Hamiltonian which shows chiral hinge states without the protection of any symmetries. The Hamiltonian is beneficial to obtain the 3D extrinsic SOTI in class A. We consider the following bulk Hamiltonian [50]: Here, H xz (k x , k z ) is a Hermitian matrix without symmetry, and 1 x has the same size as H xz (k x , k z ). H y (k y ) is the chiral-symmetric Hamiltonian explained in Sec. II. We assume that the bulk and the surface gaps are open at zero energy. Because Eqs. (B1) and (1) have the same structure, we discuss the topological hinge states in a manner similar to Sec. II by regarding H xz (k x , k z ) as a 2D bulk Hamiltonian. Therefore, we define a Chern number Ch xz for H xz (k x , k z ), which gives the 1D chiral edge states ψ chiral . The Chern number Ch xz is defined for the bulk bands below zero energy since the origin of the energy can be shifted. To see chiral hinge states in the 3D SOTI, we investigate the semi-infinite Hamiltonian with one hinge: We maintain the translation symmetry in the z direction.
In the same way as in Sec. II, we can obtain topological hinge states ψ chiral (k z ) ⊗ φ zero y with the gapless dispersion if Ch xz and w y are nonzero. Thus, we can introduce a Z topological invariant I = Ch xz w y , which characterizes the number of gapless hinge states [50]. Realistically, H y can depend on k z . However, since k z can be interpreted as a parameter in H y (k y ), gapless hinge states can be detected if the model can be adiabatically deformed to Eq. (B1) .
In comparison with our model in Sec. III C, H xz is given by H xz = (t x + t ′ x cos k x + 2χ 1 cos k z )τ x + t ′ x sin k x τ y + 2t m sin k z τ z . The massive Dirac cones on the surface give the nonzero Chern number Ch xz = −sgn(t m ). As a result, we obtain the chiral hinge states from the SOTI with the broken chiral symmetry.
On the other hand, the phase transition from a SOTSM to a SOTI with chiral hinge states was also suggested in 3D systems with parity symmetry [14,40]. The SOTSM hosts nodal lines topologically protected by parity-time symmetry in pairs, giving nearly flat hinge bands between the nodal lines. Time-reversal symmetry breaking in the SOTSM realizes the 3D SOTI, and the nearly flat hinge states become the chiral hinge states in the SOTI. The chiral hinge states are intrinsically Z 2 protected by the parity symmetry [7,9,10,14,40]. By contrast, our topological phase transition stems from broken chiral symmetry, and the induced chiral hinge states are characterized by the extrinsic second-order topology.
that the zero modes are localized at the hinge. When |[t x + f x (k z )]/t ′ x | = 1 and/or |[t y + f y (k z )]/t ′ y | = 1, the hinge flat band merges into the correspondent gapless points.
Appendix D: Topological stability of the corner modes without protection of crystal symmetries Our 2D SOTI Hamiltonian in Eq. (4) does not need protection of crystal symmetries thanks to the chiral symmetry. Here, we show the topological stability of the corner modes by breaking crystal symmetries in the model in Eq. (9). This also describes the robustness of hinge states in the 3D chiral-symmetric second-order nontrivial phases, although the translation symmetry is necessary for the k z dependence. Because the hinge states are characterized locally in momentum space by ν 2D (k z ), we can understand them as corner modes in the 2D chiralsymmetric SOTI parametrized by k z . Therefore, we focus on the corner modes in the 2D SOTIs.
The chiral-symmetric model has twofold-rotational symmetry, and two mirror symmetries with respect to the x and the y axis if t x = t y and t ′ x = t ′ y in addition to time-reversal and particle hole symmetries. To break the symmetries except the chiral symmetry, we add a perturbation term, H p = R αβ [t px (τ y ⊗ σ z ) αβ + t py (τ 0 ⊗ σ y ) αβ ]c † Rα c Rβ + R t p [c † R+xB c RA + c † R+xC c RD + H.c.] in Eq. (9), where the summation of α and β runs over sub-lattices A, B, C, and D. The perturbation preserves the chiral symmetry. Unless the additional term closes the bulk and edge gaps in the nontrivial phase, the corner modes are robust even if the Hamiltonian does not have the form in Eq. (1). Figure 5 shows topological corner modes in the model. As shown in Fig. 5(b), the corner modes survive after the perturbation breaking the crystal symmetries is added. Hence, we can see the topological stability of the corner modes due to the chiral symmetry.