On the topological immunity of corner states in two-dimensional crystalline insulators

A higher-order topological insulator (HOTI) in two dimensions is an insulator without metallic edge states but with zero-dimensional topological boundary modes localized at its corners. In this paper we show that corner modes can nevertheless appear in ordinary, topologically trivial, insulators, but they lack the topological protection characteristic of the midgap modes in a HOTI. By performing"immunity"tests on the corner modes in the breathing Kagome lattice, we demonstrate that these are intensively fragile, and therefore the Kagome lattice does not represent a higher-order topological insulator, contrary to previous theoretical and experimental claims. We also introduce a criterion based on a corner charge-mode correspondence for the presence of topological corner modes in n-fold rotational symmetric chiral insulators that explicitly precludes the existence of a HOTI protected by a threefold rotational symmetry.


I. INTRODUCTION
Topology and geometry find many applications in contemporary physics, ranging from anomalies in gauge theories to string theory. In condensed matter physics, topology is used to classify defects in nematic crystals, characterize magnetic skyrmions, and predict the presence or absence of (anomalous) metallic states at the boundaries of insulators and superconductors 1,2 . For the latter, the "topological" nature of the boundary modes, be they point-like zero modes [3][4][5] , one-dimensional chiral 6,7 and helical states [8][9][10][11][12] , or two-dimensional surface Dirac cones [13][14][15][16][17][18] , resides in their robustness. One can only get rid of these states by a bulk band-gap closing and reopening or by breaking the protecting symmetry, which can be either an internal or a crystalline symmetry. For example, in two-dimensional topological insulators [8][9][10][11][12] one can gap out the helical edge states by introducing a Zeeman term that explicitly breaks the protecting time-reversal symmetry. Similarly, one can move the end states of a Su-Schrieffer-Heeger 19 chain away from zero energy by breaking the chiral (sublattice) symmetry at the edges and/or in the bulk.
New theories exploiting the protecting role of crystalline symmetries have recently led to the discovery of the so-called higher-order topological insulators [20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36] (HOTI): states of matter characterized by the presence of topologically protected modes living at the D − n dimensional boundary of a D dimensional insulator, with n > 1 denoting the order. Thus, a two-dimensional secondorder topological insulators features point-like corner modes, while a three-dimensional second-order topological insulator features helical or chiral edge modes along the one-dimensional edges. The prediction of higherorder topological insulators has triggered an enormous interest in scanning material structures and engineering metamaterials, e.g. electric circuits 28 , exhibiting topological corner or hinge modes. However, in identifying higher-order topological insulating phases different complications arise. First, as in conventional first-order topo-logical insulators, the system can display ordinary in-gap boundary states 37 that are not the prime physical consequence of a non-trivial bulk topology. Second, D − 2 boundary modes of a D-dimensional system can be a manifestation of the crystalline topology of the D − 1 edge rather than of the bulk: the corresponding insulating phases have been recently dubbed as boundaryobstructed topological phases 38 and do not represent genuine (higher-order) topological phases. These complications are particularly severe in second-order higher-order topological insulators in two dimensions, since the corresponding zero-dimensional topological boundary modes fail to possess the anomalous nature chracteristic of, for instance, one-dimensional chiral modes. Singling out a proper second-order topological insulator in two dimensions thus represents a task of exceptional difficulty. In this paper, we prove this assertion by showing that one of the first model system suggested to be a second-order topological insulator in two-dimension 39-42 -the breathing Kagome lattice model -does not host any higherorder topological phase. The corner modes experimentally found in this system [40][41][42] are neither engendered from a second-order bulk topology nor from the edge topology characterizing boundary-obstructed topological phases. They instead are an example of conventional corner modes, with the exact same nature of the edge modes generally appearing in one-dimensional insulating chains 43 . We contrast the fragility of the corner modes in the Kagome lattice with the robustness of the corner modes in chiral-symmetric insulators, which possess the full immunity of the topological corner modes of a HOTI when an additional even-fold rotational symmetry is present. We also formulate a one-to-one correspondence between fractional corner charges 35,44 and corner modes, which predicts the presence or absence of topological zero modes in chiral symmetric insulators with a C n -rotational symmetry, and precludes the appearance of a HOTI phase protected by a threefold rotational symmetry.

II. (NON)-TOPOLOGICAL EDGE STATES IN ONE-DIMENSIONAL INSULATORS
Generally speaking, topological boundary states are a special type of in-gap boundary states, whose presence if ubiquitous in crystalline systems. This can be nicely illustrated using a minimal model for a one-dimensional band insulator: the Rice-Mele model 45 . It is schematically shown in Fig. 1(a). The electrons living on the red -sublattice A (green -sublattice B) sites experience an on-site energy +m (−m), can hop within the unit-cell with hopping amplitude t, and between adjacent unitcells with a hopping amplitude t . Within the bulk energy gap and when |t/t | < 1 we find that the left edge hosts a state at energy +m, whereas the right edge hosts a state at energy −m [c.f. Fig. 1(b)]. On the contrary, when |t/t | > 1 the system fails to exhibit any boundary state [c.f. Fig. 1(c)]. Hence, the Rice-Mele model exhibits boundary states in half of the parameter space, assuming the termination shown in Fig. 1(a) is chosen. In order to establish whether these boundary states are topological in nature, we resort to the main characteristic of topological boundary states: their robustness against smooth perturbations. We notice the following: (i) The energies of the in-gap boundary states change either introducing smooth deformations of the bulk Hamiltonian or applying edge-specific perturbations. An example of this is illustrated in Fig. 1(b).
(ii) The boundary states can completely dissolve into the bulk bands without an accompanying bulkband gap closing, as shown in Fig. 1(c). In the Rice-Mele model this occurs for |t/t | = 1.
(iii) A tailor-made edge potential can even lead to the creation of boundary states out of the bulk bands in the parameter region |t/t | > 1, as explicitly shown in Fig. 1(d).
The boundary states of the Rice-Mele model fail to exhibit any kind of robustness, and can be therefore qualified as ordinary, i.e. non-topological boundary states. This, however, is not yet the end of the story. In fact, in one-dimensional insulating models it is possible to exploit the effect of the internal particle-hole and chiral symmetries. Their existence implies that the model Hamiltonian anticommutes with an (anti)unitary operator that squares to 1 (±1). The Bogoliubov-de Gennes (BdG) Hamiltonian describing a superconductor is, by its very definition, particle-hole symmetric. However, both the particle-hole and the chiral symmetry also play a role outside the superconducting realm. They can in fact arise as approximate symmetries of the effective model Hamiltonian describing an insulator. As an example, in the Rice-Mele chain we may set m = 0, in which case the Hamiltonian anti-commutes with σ z . The resulting model is the well-known Su-Schrieffer-Heeger (SSH) chain.
Using the results above, we find that the SSH chain displays a left and a right boundary state at zero-energy if and only if |t/t | < 1. Importantly, these boundary states now represent truly topological boundary states. Obviously, it would amount to cherry-picking if these states would be dubbed as topological only because they are at zero energy. The rationale for the above is instead based on the following fact: The boundary modes cannot be moved away from zero energy using chiral symmetry preserving bulk or edge perturbations, as long as these do not close the bulk band gap. For example, long-range hoppings between the two sublattices do not move the boundary states since the chiral symmetry is preserved when these processes are included. Hence, the boundary states of the SSH chain are robust zero-energy modes that are protected by the internal chiral symmetry. This robustness can be rigorously proved using the following: let H(λ) parametrize the chiral-symmetric perturbed Hamiltonian (with λ the perturbation strength) of a semi-infinite atomic chain that features a single left end state |Ψ L (0) at zero energy. To show that this edge state will remain pinned at zero energy, we can track the left edge state |Ψ L (λ) , assuming λ is switched on adiabatically. The chiral symmetry guarantees that every state |χ at an energy +E = 0 has a chiral partner σ z |χ at energy −E. Hence, end states can only move away from zero energy in pairs. Thus, we can conclude that the parity of the number of left edge states pinned at zero energy is robust against all continuous perturbations that respect the chiral symmetry. Note that although this conclusion is strictly valid for a semi-infinite chain, it also holds for finite chains, i.e. with a left and a right edge, as long as the decay length of the edge states is small compared to the system size. In particular, this argument cannot be used when the perturbation closes the bulk band gap, since in the latter case the decay length of the boundary states always reaches the system's size. This additionally shows that the presence of the zero-energy edge states is an intrinsic feature of the one-dimensional bulk Hamiltonian.
Since we can discriminate between systems with an odd and even number of left edge states pinned at zero energy, we arrive at the conclusion that one-dimensional chiral-symmetric insulators -but the argument above holds true also for particle-hole symmetric insulators belonging to the class D of the Altland-Zirnbauer table 46,47 -are characterized by the group Z 2 = {0, 1}, with the identity element 0 corresponding to chains featuring an even number of left edge states, and the element 1 corresponding to chains featuring an odd number of left edge states. Physically, the Z 2 group law translates into the fact that if one combines (i.e. hybridizes in a symmetrypreserving way) a chain with an even number of left edge states and a chain with an odd number of left edge states, the combined two-leg atomic chain possesses an odd number of left edge states. On the contrary, hybridizing two chains with both an even or odd number of left edge states results in a two-leg atomic chain with an even number of left edge states. In Fig. 2 we illustrate this by considering a two-leg atomic chain consisting of two SSH chains, each of which featuring a single left edge state, hybridized in two different ways. In the first example, shown in Fig. 2(a) we have hybridized one SSH chain that terminates with an A site, with another SSH chain that terminates with a B site. The result of this hybridization, is that the two left edge states also hybridize, and move away from zero energy, as explicitly shown in Fig. 2(c). This is simply a manifestation of level repulsion. However, the two-leg atomic chain shown in Fig. 2(b) displays a more interesting behavior. There we have hybridized two SSH chains that both terminate with the same sublattice. Despite the hybridization, we find that the two edge states do not move away from zero energy [c.f. Fig. 2(d)]. Even though these two examples provide purely anecdotal evidence, they do suggest that the physics of zero-energy edge states is not fully captured by a Z 2 -invariant. This is indeed the case for chiral-symmetric insulators. In these systems, the isolated zero energy modes must be eigenstates of the unitary chiral symmetry, and can be consequently characterized by their chiral charge. For the specific example of the SSH atomic chain, this also implies that an isolated zero mode is either fully localized on the A-sublattice or on the B-sublattice. Denoting with |Ψ the zero energy end state of a SSH chain and with χ = Ψ| σ z |Ψ the corresponding chiral charge, we have that if χ = 1 the zero-energy mode will be fully localized on the sublattice A, whereas if χ = −1 the zero energy mode will be fully localized on the sublattice B. The perfect localization of the end states clearly implies that having at hand two zero-energy states localized on the same sublattice, and thus with same chiral charge, impedes any level repulsion as this would necessarily break the chiral symmetry. On the contrary, a pair of zero-energy modes localized on different sublattices can be opportunely coupled and moved away from zero energy, in perfect agreement with the features of the two-leg atomic chain shown in Fig. 2(c),(d). This proves that the physics of zero energy states in chiral-symmetric insulators is encoded in the Z-valued chiral charge: where the sum runs over all left edge states |Ψ j . Note, however, that the relation between the chiral charge χ L and the number of left edge states at zero energy does not represent a one-to-one correspondence. With a zero chiral charge, a pair of zero energy boundary states can be still encountered. This thus implies that the absolute value of the left chiral edge charge defines a lower bound for the total number of zero-energy states, and is equal to the number of edge states modulo 2. Hybridizing two insulators with chiral charges n 1 and n 2 , respectively, we will end up with a chain that features at least |n 1 + n 2 | edge states at zero energy. For completeness, we finally note that, being a bulk quantity, the chiral charge χ L can be expressed as a winding number in terms of the one-dimensional Bloch Hamil-tonian using the following formula: For a detailed derivation of the above bulk-boundary correspondence, we refer the reader to Ref. 48.
To summarize this section, we have shown that boundary states are a generic feature of one-dimensional band insulators. However, in the absence of any symmetry these states are not robust, for instance they can dissolve into the bulk bands, or even created as a result of tailor-made edge perturbations. Instead, if one considers systems with a particle-hole or chiral symmetry, we find that the parity of the number of left and right edge states is robust against symmetry allowed perturbations. Moreover, for systems with a chiral symmetry we find that there is a chiral charge associated to each edge, which represents a Z number and is a lower bound for the number of edge states.

III. CORNER MODES IN THE BREATHING KAGOME LATTICE
The preceding discussion on the physical properties of end states in one-dimensional insulators can be applied in an analogous fashion to corner states in twodimensional insulators. This can be nicely illustrated using the breathing Kagome lattice [c.f. Fig. 3(a)], which can be thought of as the two-dimensional cousin of the Rice-Mele atomic chain. For simplicity, we will assume in the remainder that all intra-unit cell hopping amplitudes have an equal magnitude t. We will also make the same assumption for the inter-unit cell hopping amplitudes (magnitude t ). In addition to the inter-and intra-unit cell hopping parameters, we will first allow for different on-site energies in the three sublattices, which we denote with m 1 , m 2 , and m 3 , respectively.
Following Ref. 49 we find that when considering the system in a open-disk geometry respecting the threefold rotational symmetry of the lattice, the lower left, lower right, and upper corner host corner states at energies m 1 , m 2 , and m 3 in half the parameter space, i.e. for |t/t | < 1. This condition being very similar to the one we encountered in the Rice-Mele atomic chain, suggests that also the corner modes of the breathing Kagome lattice correspond to conventional in-gap bound states. And indeed these modes dissolve into the bulk by changing the ratio |t/t |, see, e.g., Ref. 49. Next, let us consider the situation in which the on-site energies of the three sublattices are constrained to be equal, i.e. m 1 ≡ m 2 ≡ m 3 ≡ 0. The fact that the corner modes now reside at zero energy in an extended region of the parameter space could suggest that the corner modes represent in this case genuine topological boundary modes. Contrary to the SSH atomic chain, however, the Kagome lattice model does not possess an internal chiral symmetry that can protect the existence of corner modes pinned at zero energy. This follows from the very simple fact that a chiral-symmetric insulator is incompatible with an odd number of sublattices. Nevertheless, recent studies have suggested that other symmetries could protect the existence of the zero-energy corner modes, qualifying them as topological corner modes and consequently the prime physical consequence of a higher-order non-trivial bulk topology. Specifically, in Ref. 39 it has been suggested that the corner modes are protected by the combination of a three-fold rotational symmetry and a mirror symmetry. Instead, in Refs. 40 and 41 a "generalized" chiral symmetry has been defined in order to prove the topological nature of the corner modes. Such a generalized chiral symmetry is equivalent to requiring that the model Hamiltonian may only be perturbed by hopping processes between different sublattices.
We now show that the opposite is true, and that the corner modes found in the breathing Kagome lattice are nothing but conventional boundary modes even for m 1,2,3 ≡ 0. As mentioned in the preceding section, the defining characteristic of any topological boundary mode is its immunity against perturbations that do not close the insulating band gap and preserve the protecting symmetries. In the context of two-dimensional insulators (without metallic edge states) in general, and of the breathing Kagome lattice in particular, this would mean that the corner modes have to remain pinned at zero energy upon perturbing the two-dimensional bulk, the one-dimensional edges, or the zero-dimensional corners. To see whether the breathing Kagome lattice possesses boundary modes with such robustness, we have considered the effect of applying a local perturbation at the three corners. Specifically, we have introduced longrange hopping processes, with amplitudes s 1 and s 2 , at the three corners, as schematically shown in Fig. 3(b). Note that the corner perturbation fulfills all the symmetry constraints: the three-fold rotational symmetry, the mirror symmetries, and the generalized chiral symmetry are all preserved. Let us now consider the evolution of the corner state energy as we adiabatically "switch on" the perturbation. In Fig. 3(c), we have plotted the evolution of the spectrum assuming the intracell hopping amplitude t ≡ 0 and the corner perturbation hopping amplitudes satisfy s 1 ≡ −s 2 . We immediately find that the corner modes do not remain pinned at zero energy. In fact, we find that upon increasing the perturbation strength the energy of corner mode even crosses the edge and bulk valence bands at energies −t and −2t , respectively. This demonstrates that the corner modes in the breathing Kagome lattice do not possess any topological robustness. Their presence or absence is not in a one-toone correspondence with a topological invariant. Instead they simply constitute ordinary boundary states, as do the ones occurring when the on-site energies m 1,2,3 are different from zero.

IV. CORNER MODES IN CHIRAL-SYMMETRIC INSULATORS
Having established with a concrete microscopic model that the corner modes in the Kagome lattice simply corresponds to fragile ordinary boundary modes, we next introduce a chiral-symmetric insulator featuring robust corner modes. This will also allow us to discuss the different nature and degree of protection provided by a nontrivial "edge" topology as compared to the bulk topology of a higher-order topological insulator.
The microscopic tight-binding model we will consider is schematically depicted in Fig. 4(a). It possesses an internal conventional chiral symmetry and a C 4 fourfold rotational symmetry. When considered in an open-disk geometry that respects both the rotational and the chiral symmetries, the system features four zero-energy states which are completely localized at the corners of the lattice when the intra-cell hopping amplitude t ≡ 0. Precisely as the end modes of the SSH chain, these corner modes can be characterized with a Z number corresponding to their chiral charge χ. More importantly, we find that the the corner modes remain pinned at zero energy under the influence of perturbations that preserve the chiral symmetry, and neither close the edge nor the bulk band gap. This follows from the same reasoning which is behind the stability of zero energy modes in the SSH atomic chain. In particular, it should be stressed that this stability does not rely on the fourfold rotational symmetry.
There is, however, an important distinction between the edge states of the SSH chain and the corner states of a two-dimensional chiral-symmetric insulator. The presence of edge states in the SSH atomic chain is only dependent on the topology of the one-dimensional bulk Hamiltonian. On the contrary, the presence of a corner mode in a two-dimensional chiral-symmetric insulator is dependent on both the topology of the neighboring edges and the topology of the two-dimensional bulk. To illustrate this point, we apply to the model shown in Fig. 4(a) an edge perturbation that can be strong enough to close and reopen the edge gap. Specifically, we introduce a perturbation on the top and bottom edges [c.f. Fig. 4(b)], thereby explicitly breaking the fourfold rotational symmetry. The spectral flow obtained by increasing the strength of the edge perturbation [c.f. Fig. 4(d)] shows that the zero-energy modes remain pinned at zero energy in the weak perturba tion regime. However, in the strong edge perturbation regime, i.e. after the closing and reopening of the edge band gap, the corner modes disappear. Importantly, the edge pertur- Figure 4. (a) C4-symmetric model. Each plaquette is threaded by a π-flux. (b) a strong C4-symmetric edge perturbation is applied, leaving the corner modes intact. However, before the band-closing the corner mode of the lower left corner had a chiral charge χ = +1, after the band-gap reopening its value has changed to χ = −1. (c) If a strong edge perturbation is applied along the vertical edges only, thus breaking the C4-symmetry, we find that the corner modes do not survive the edge band-gap closing and reopening. The spectral flow for going from (a) to (b) is shown in panel (d), whereas the spectral flow in going from (b) to (c) is shown in (e). In both cases the intra-and inter-cell hopping parameters t edge and t edge along the perturbed edges are given by t edge = t sin(θ) and t edge = t cos(θ). bation leaves the bulk of the crystal completely intact independent of its strength. This shows that zero-energy modes carry a halfway protection: they are removable only by edge perturbations causing closing and reopening of the edge band gap. In the language of Ref. 38, the chiral-symmetric insulating model shown in Fig. 4(a), therefore represents a boundary-obstructed topological phase protected by the chiral symmetry. We note that in these phases it is impossible to disentangle the edge topology from the bulk topology. Therefore, the presence of a zero-energy corner state by itself does not provide insights in the bulk topology alone.
Even though the chiral symmetry on its own is insufficient to stabilize the corner modes against strong edge perturbations, we find that the additional presence of the fourfold rotational symmetry does offer this kind of protection. To show this, we consider an additional edge perturbation along the left and right edge, increasing the strength of which eventually leads to the configuration shown in Fig. 4(c) where the C 4 symmetry has been restored. The ensuing spectral flow shown in Fig. 4(e) shows that an additional edge gap closing and reopening point leads to a revival of the zero-energy modes. Therefore, the fourfold symmetric system is characterized by zero energy modes independent of the presence of a strong edge perturbation, qualifying them as the manifestation of the two-dimensional bulk topology. In other words, supplementing the chiral symmetry with the C 4 symmetry turns the model into a second-order topological insulator. This can be also seen using the following argument: in the concomitant presence of the chiral and C 4 rotational symmetry, each corner mode is at the intersection of two adjoined edges related to each other by the fourfold rotational symmetry. The chiral charge χ of this corner mode will change to χ → χ + j, with j an integer upon closing an reopening the band gap along one of the two edges. On the other hand, owing to the fourfold rotational symmetry, the same band gap closing and reopening will occur on the other edge, which will thus contribute with an additional and equal change of the chiral charge. Hence, in total we find that the chiral charge of the corner mode is modified as χ → χ + 2 × j. This therefore implies that, by virtue of the C 4 symmetry, the parity of the chiral charge, and consequently the parity of the number of zero modes per corner ν, is invariant under strong edge perturbations and thus represents a proper Z 2 -invariant of such a two-dimensional insulator. Note that this topology is considerably weaker than the topology of the SSH chain, as the latter is characterized by a Z-number and does not necessitate the presence of a rotational symmetry.

V. CORNER CHARGE-MODE CORRESPONDENCE IN ROTATIONAL SYMMETRIC CHIRAL INSULATORS
We now show that the topological immunity of the zero-energy boundary modes in chiral-symmetric insulators with a fourfold rotational symmetry can be also proved using a corner charge-mode correspondence that can be generalized to chiral insulators with an even-fold rotational symmetry. We first recall that, as shown in Refs. 35 and 44, the crystalline topological indices characterizing rotational symmetric two-dimensional insulators are revealed in the fractional part of the corner charge. Specifically, for corners whose boundaries cross at a maximal Wyckoff position with a site symmetry group that contains the n-fold rotational symmetry, the fractional part of the corner charge is a topological Z n number uniquely determined by the symmetry labels of the occupied Bloch states at the high-symmetry momenta in the Brillouin zone.
Generally speaking, the correspondence between these bulk topological indices and the fractional corner charge is not reflected in the presence or absence of corned modes. However, we now show that a direction relation exist between the fractional part of the corner charge and the parity of zero-energy corner modes in chiralsymmetric insulators featuring a two-, four-, or six-fold rotational symmetry. We will in fact demonstrate that the parity of the number of zero-energy states per corner ν obeys the formula ν = (2 × Q v.b. + N total /n) modulo 2. (1) In the equation above, N total is the total number of sites per unit cell, n, as before, the order of the rotational symmetry, and Q v.b. the corner charge due to the valence bands, which is a quantity quantized to 0 or 1/2 modulo 1. Before deriving the relation above, a few remarks are in order. First, we wish to emphasize that the corner charge Q v.b. should be computed as the total charge within a C n -symmetric corner region that is congruent with unit cell centers whose site symmetry group contains the C n rotational symmetry. This implies that the boundaries of the corner region should be related to each other via the rotational symmetry, and that the boundaries of the region cross at the C n -symmetric unit cell center. Example of these corners for C 2 , C 4 , and C 6 symmetric system are shown in gray in Figs. 5(a)-(c). Second, we remark that the direct relation between corner charge and corner modes is strictly valid for geometric configurations in which the finite-size system can be tiled using an integer number of unit cells. This exclude the presence of fractional unit-cells in the open-disk geometry. Note that the lattice structures shown in Figs. 5(a)-(c) obey this constraint. Finally, we will consider lattices where the positions in the unit cell with a site symmetry group containing the C n symmetry do not host any atomic site. Note that this condition automatically ensures that the total number of atomic sites per unit cell is a multiple of the rotational symmetry order n.
We can now derive the corner charge-mode correspondence of Eq. 1, and first note that at full filling the corner charge is equal to N total /n modulo 2. By inspection of Fig. 5, it can be easily seen that in all three lattice structures N total /n = 1. Next, we note that the corner charge at full-filling can be decomposed into three separate contributions: a valence band Q v.b. , a conduction band Q c.b. , and finally an in-gap state Q in-gap contribution. Here, the valence (conduction) band part accounts for all states whose energies E lie below (above) the band gap, i.e. E ≤ ∆/2 (E ≥ ∆/2), with ∆ the band gap. The in-gap contribution instead is due to states whose energies E lie inside the band gap, i.e., |E| < ∆/2. Furthermore, the presence of chiral symmetry guarantees that the corner charge due to the conduction band is identically equal to the corner charge due to the valence band, i.e. Q v.b. = Q c.b. . Note that the latter is a true equality of number, which includes also the integer part of the charge. As a result, we obtain the following: Finally, we use that the in-gap contribution Q in-gap is equal, modulo 2, to the parity of the number of corner states per corner ν. Hence, upon rearranging terms we arrive at Eq. 1 that proves the correspondence between fractional corner charge and the presence of zero-energy corner modes. As the corner charge is a direct probe of the two-dimensional bulk topology 35,44 , Eq. 1 implies that the parity of the the number of zero-energy states per corner ν itself is a manifestation of this bulk topology. In principle, one may repeat the above analysis for C 3 -symmetric insulators that are also chiral-symmetric. However, an interesting interplay between the chiral symmetry and the three-fold rotational symmetry unfolds. Namely, the chiral symmetry requires that at half-filling the corner charge is a multiple of 1/2, instead the threefold rotation symmetry implies that the corner charge is a multiple of 1/3. Satisfying both conditions leaves only one option, namely a vanishing corner charge. In other words, the chiral symmetry renders the corner charge trivial in a C 3 -symmetric insulator, in the same way that time-reversal symmetry renders the Chern number trivial. Therefore, we find that the mere presence or absence of zero-energy corner states in a C 3 -symmetric insulator does not shine any new light on the two-dimensional bulk topology. In fact, any C 3 -symmetric configuration that can be tiled with an integer number of unit-cells will fail to exhibit an odd number of zero-energy states per corner. To prove this, let us suppose that each of the three corners would host a single zero-energy state. This would then imply that the lattice as a whole exhibits a chiral imbalance. However, this is at odds with the original assumption that the finite geometry can be tiled with an integer number of unit-cells. Ergo, C 3 -symmetric insulators fail to exhibit an odd number of zero-energy modes per corner.

VI. DISCUSSION
To sum up, we have shown that precisely as boundary modes in atomic chains, corner modes are a generic feature of two-dimensional band insulators and are not necessarily a signature of topology. Contrary to previous theoretical and experimental claims, we have indeed proved that the corner modes encountered the breathing Kagome lattice do not exhibit any kind of topological robustness, and have to be instead qualified as ordinary boundary modes. We have also contrasted the fragility of these modes with the robustness of the zero-energy modes appearing in insulators equipped with an internal chiral symmetry. When taken alone, the chiral symmetry provides the halfway topological robustness characteristic of the recently introduced boundary-obstructed topological phases. Furthermore, the presence of a rotational symmetry provides an additional protection mechanism that qualifies the zero energy modes as the prime physical consequence of a higher-order bulk topology.
We have also proved that the immunity of the topological corner modes in the concomitant presence of rotational and chiral symmetry directly follows from a oneto-one correspondence between fractional corner charges, which reveal the crystalline topology of generic insulators, and the parity of the number of zero modes. This one-to-one correspondence only works in crystals possessing an even-fold rotational symmetry, and thus exclude C 3 -symmetric crystals such as the Kagome lattice. We wish to remark, however, that even though the breathing Kagome lattice does not display a genuine "bulk-corner" correspondences, its underlying crystalline topology is still reflected in the fractional charge at corners or other topological defects, such as dislocations 50 .
In closing, we would like to highlight that the topological immunity of zero energy modes in insulators is very different in nature from the one provided by timereversal symmetry in "conventional" first-order topological insulators. This is immediately apparent from the fact that the helical edge states of a two-dimensional topological insulator represent anomalous states: Their anomaly reside in the fact that it is impossible to find a one-dimensional insulator with an odd number of Kramer pairs at the Fermi energy. On the contrary, the zeroenergy fermionic modes encountered in chiral-symmetric insulators do not represent an essential anomaly for the very simple reason that a quantum dot with an odd number of fermionic modes is entirely allowed in nature. This difference is not only interesting per se. In fact, a topological boundary mode that is also anomalous carries an additional degree of protection. Consider for instance a two-dimensional topological insulator whose edges are brought in close proximity to a time-reversal invariant nanowire. In such combined system the edge boundary modes will survive, at odds with what would happen in a Su-Schrieffer-Heeger atomic chain if the termination is changed by an atomic site addition. Moreover, one can also proceed in the opposite direction. For example, one may modify a two-dimensional insulator that is both C 3 -symmetric and chiral-symmetric but fails to exhibit any in-gap mode by attaching single sites to its three corners. Even though this would result in the presence of robust zero-energy corner modes, it should be stressed that these modes are nothing but a microscopic detail of the corner and edges. Put differently, these modes are not in any way informative of the bulk topology.