Open AccessArticle

Effects of the Vertices on the Topological Bound States in a Quasicrystalline Topological Insulator

Simone Traverso

¹,

Niccolò Traverso Ziani

^1,2,* and

Maura Sassetti

^1,2,*

Dipartimento di Fisica, Università Degli Studi di Genova, Via Dodecaneso 33, 16146 Genova, Italy

CNR SPIN, Via Dodecaneso 33, 16146 Genova, Italy

Authors to whom correspondence should be addressed.

Symmetry 2022, 14(8), 1736; https://doi.org/10.3390/sym14081736

Submission received: 27 July 2022 / Revised: 10 August 2022 / Accepted: 16 August 2022 / Published: 19 August 2022

(This article belongs to the Special Issue Symmetry, Topology and Phases of Condensed Matter)

Download

Browse Figures

Versions Notes

Abstract

The experimental realization of twisted bilayer graphene strongly pushed the inspection of bilayer systems. In this context, it was recently shown that a two layer Haldane model with a thirty degree rotation angle between the layers represents a higher order topological insulator, with zero-dimensional states isolated in energy and localized at the physical vertices of the nanostructure. We show, within a numerical tight binding approach, that the energy of the zero dimensional states strongly depends on the geometrical structure of the vertices. In the most extreme cases, once a specific band gap is considered, these bound states can even disappear just by changing the vertex structure.

Keywords:

higher order topological insulators; quasicrystals; bound states

1. Introduction

The quest for quantum states with enhanced robustness to uncontrollable perturbations, such as noise, triggered the study of the topological phases of matter [1]. From the point of view of electronics and spintronics [2,3,4,5,6], the main focus is on non-superconducting systems, and in particular on topological insulators [7,8,9,10]. For superconducting spintronics [11,12,13,14,15] and topologically protected quantum computation [16,17,18,19,20,21,22,23,24,25,26,27,28,29], even through parafermions [30,31,32,33], a prominent role is conversely played by topological superconductors [34].

The physics is similar in the two cases: The systems have a gapped bulk, and bound states located at the physical boundaries with energy in such gaps. The bound states are predicted to enjoy, for topological reasons, enhanced robustness with respect to perturbations and are hence promising for quantum technological applications.

The robustness of such topological states is, however, debated and, in general, appears to be reduced with respect to the theoretical predictions. Prominent examples in this direction deal with the origin of the conductance fluctuations in quantum spin Hall systems, which is not perfectly quantized even when time-reversal symmetry is not explicitly broken [35,36,37,38,39,40,41,42,43], and with the difficulty in establishing the topological (Majorana) nature of the bound states in one-dimensional topological superconductors [44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63].

To make progress, new paradigms for topological protection have hence recently been explored. A very promising one is higher order topology [64,65,66,67,68]. A d-dimensional higher order topological insulator or superconductor is a system in which the topologically protected states live in n dimensions less than the gapped bulk, with

1 < n \leq d

. The index n is referred to as the order of the topological insulator. Very recently, three-dimensional gapped systems with one-dimensional topological states have been realized experimentally [69,70]. On the other hand, proposals for the highly relevant case of two-dimensional bulk and zero-dimensional topological bound states have been put forward, but not experimentally accessed.

Among such proposals, a noteworthy one has been envisioned by S. Spurrier and N. Cooper [71]. The main ingredient is a quasicrystalline two-dimensional lattice, generated by two layers of graphene-like lattices with a 30 degree relative rotation, and possessing 12-fold rotational symmetry. Each layer consists of a Haldane model and the two layers have opposite Chern number. Moreover, the edges of the bilayer structure they consider are of the armchair and bearded type. The result is that the system, thanks to the presence of the rotational symmetry, hosts topological zero-dimensional bound states, which are degenerate, isolated in energy from the others, and localized at the vertices of the structure. The system is hence a higher order topological insulator with

d = 2

and

n = 2

. A subsequent work [72] showed that these properties also characterize samples with armchair-zigzag edges. However, in this case, the energy at which the gap opens drastically changes. While in the armchair-bearded case the gap opens at the charge neutrality point, in the armchair-zigzag case it is substantially moved away from it. This observation poses questions about the sensitivity to details of the model.

In this context, the aim of the present work is to address some questions related to the robustness of the properties of the topological bound states with respect to lattice perturbations at the corners of the structure: Do these corner modes always exist in a given gap? Which is their energy? How strong is the protection granted to their degeneracy?

We will show how the energy, degeneracy and even existence of the predicted corner modes astonishingly depend upon the precise disposition of the lattice sites around the sample vertices. These results seem to narrow the path for a hypothetical experimental realization of the system proposed, which would require site-wise precision in the realization of the sample in order to achieve the HOTI phase close to the Fermi level.

The paper is organized as follows: In Section 2 we introduce the model, both for how it concerns the lattice and for how it concerns the Hamiltonian. In Section 3 we present and discuss our results, answering one by one the questions we addressed above. In Section 4 we draw our conclusions.

2. Model

In order to proceed with the discussion of the results achieved, we must first introduce the model of HOTI on a quasicrystal, which will be the object of our analysis. The present section is organized in two parts: in the first one we introduce the system lattice, while in the latter we present the model Hamiltonian. In what follows, we will consider a square sample instead of a dodecagonal one, as is the case in [71]. Indeed, when dealing with the robustness of the bound states, the choice of a square sample provides major computational advantages, while leaving the physics qualitatively unchanged.

2.1. Lattice

Let us begin by describing the construction of a finite-size square sample of the quasicrystalline lattice under inspection: one starts with the two honeycomb layers superimposed with AA stacking (i.e., each atom of the top layer is right on top of an atom of the bottom layer) and rotates one of the two layers by thirty degrees with respect to the other, taking a rotation axis that goes through the center of two superimposed hexagons (see Figure 1). The resulting lattice is crystallographically equivalent to that of the graphene quasicrystal, which has been shown to possess local bulk

C_{12}

rotational symmetry [73].

We then crop our lattice with a square shape centered in the rotation center, and oriented in such a way that the edges are not jagged. As already noted in [72], the edges of the resulting system—which possesses

C_{4}

global rotational symmetry—can only be given by two possible pairs of the upper and lower layer edges: bearded-armchair or zigzag-armchair. However, even once the combination of the edges has been chosen, there are still two families of square lattices that can be obtained: these differ with each other for the disposition of the lattice sites around the vertices. The zoom of a vertex for each of the four possible configurations just discussed is shown in Figure 2.

In [72], it was shown how the properties of the model for an HOTI on the (spinless) graphene quasicrystal first proposed in [71], depend on the shape of the edges. Here, we want to discuss how, even when the edge combination is fixed, the features of the HOTI phase can still depend on the specific configuration of the sites near the vertices. In order to approach this target, let us first introduce the model Hamiltonian.

2.2. Model Hamiltonian

The HOTI model we analyze [71] is governed by the following tight-binding Hamiltonian

H = t \sum_{〈 i j 〉} c_{i}^{†} τ_{0} c_{j} + λ_{H} \sum_{〈 〈 i j 〉 〉} i ν_{i j} c_{i}^{†} τ_{z} c_{j} + λ_{⊥} \sum_{i j} t_{i j}^{⊥} c_{i}^{†} τ_{x} c_{j},

(1)

where we denote

c_{i} = {(c_{i}^{t}, c_{i}^{b})}^{T}

, with the spinless Fermionic operators

c_{i}^{†, t / b}

creating a fermion on the i-th lattice site of the top/bottom layer of honeycomb lattice. The Pauli matrices

τ_{i}

act in the space of layer pseudospin, coupling the Fermionic operators of the two layers. The Hamiltonian H presents two in-plane hopping terms: the first one describes nearest neighbors (

〈 i j 〉

) hoppings with amplitude t, while the second one involves next-nearest neighbors (

〈 〈 i j 〉 〉

) hoppings of amplitude

λ_{H}

. The latter is associated to a complex phase

i ν_{i j}

(

ν_{i j} = \pm 1

), given by the usual Haldane coupling [74,75]. Also, an inter-layer hopping term is present, parameterized by

t_{i j}^{⊥}

. Its explicit form is borrowed from tight-binding models for real twisted bilayer graphene (TBG) samples [76,77].

t_{i j}^{⊥} = t^{⊥} exp (- \frac{| {\vec{d}}_{i j} | - d_{0}}{δ}),

(2)

with

{\vec{d}}_{i j}

being the vector that goes from the i site on one layer to the j site on the other one. We take the the parameters appearing in Equation (2) as in tight-binding models for real TBG samples as well [76,77]:

δ = 0.184 a

is the decay length of the transfer integral,

d_{0} = 1.362 a

the distance between the two layers and

t^{⊥} = - 0.178 t

a coupling coefficient. All lengths are expressed in units of the in-plane next-nearest neighbor distance a. Finally, with the aim of tuning the interlayer coupling strength, we add the parameter

λ_{⊥}

in Equation (1) (the parameter

t_{⊥}

is taken in such a way that if

λ_{⊥} = 1

, then the value of the coupling for interlayer vertical hoppings matches the value used in realistic tight binding models for TBG [76,77]).

In order to achieve a slightly more intuitive perspective about the model, it is worth noting [71] that without the interlayer hopping term (

λ_{⊥} = 0

), the Hamiltonian actually consists of two Haldane models [74] with opposite Chern number (

\pm 1

in the upper/lower layer respectively) and with a relative

30^{\circ}

twist. This is almost the same as the Kane–Mele model [7,8], just with the spin mapped into layer pseudospin. Indeed, with

λ_{⊥} = 0

the low energy spectrum would be gapless, with the edges of the two planes hosting a pair of counter-propagating modes- just as in the Kane–Mele model. What the interlayer hopping does is to couple together these two Haldane models: more specifically, as was shown in [71], it gaps out the counter-propagating edge modes on the two layers, making it possible for the second order topological insulator (SOTI) phase to occur. Such a phase is characterized by the presence of bound states localized at the corners of the system (corner modes), degenerate among themselves and with energy laying inside the edge-gap.

The phenomenology just discussed is encoded in the effective low energy theory of the model derived in [71]. This theory is only based on the symmetries of the model and it leads to a few possible values for the energy of the corner modes inside the energy gap. For the present case of square sample (

C_{4}

global symmetry), these are given by:

ε_{cm}^{th} = cos (π \frac{p}{4}), p \in Z .

(3)

Having introduced all the necessary ingredients, we can now proceed to presenting and discussing our results.

3. Results

We will now show the results obtained through the numerical diagonalization of Hamiltonian (1) on different square samples of the quasicrystalline lattice. For the construction and diagonalization of the model, the python package pybinding [78] has been used. The fixed parameters of the Hamiltonian were set as described above. We set

λ_{⊥} = 2

, resulting in interlayer hoppings of double strength with respect to real TBG samples. This choice does not change the physics on a qualitative level, but it increases the gap width, reducing the hybridization of the corner modes for the same edge lengths. Also, we set

λ_{H} = 0.3 t

, so that the Haldane bulk gap is maximized [74].

The present section is organized with the aim of separately answering the questions addressed in the introduction about the stability of the corner modes with respect to local geometrical perturbations of the lattice. We start by addressing the issue of the existence of the corner modes in a certain gap.

3.1. Existence

Let us begin by presenting the low energy eigenvalues obtained for four different square samples, one for each of the categories in Figure 2. These are reported in Figure 3, where the plots of the eigenvalues are ordered the same way as the corresponding lattices, exemplified in Figure 2.

Figure 3a,b reports the low-energy eigenvalues for two square samples with bearded-armchair edges and different shapes of the corners (see Figure 2a,b); in both cases we find four degenerate eigenvalues located inside the edge-gap. The existence of the topological bound states in the gap is hence robust with respect to variations in the corner shape. The lattices considered for the diagonalization present an edge with a length of

L \approx 321 a

and

L \approx 347 a

, respectively. The slight difference in the energy of the degenerate quadruplet located inside the edge-gap will be discussed later in the article.

Figure 3c,d reports the low-energy eigenvalues for two square samples with zigzag-armchair edges (in this case, the gap does not open at the charge neutrality point [72]). In this case, the lattices considered for the diagonalization present an edge of length of

L \approx 280 a

and

L \approx 306 a

respectively. Just as before, these two lattices mainly differ for the disposition of the sites near the vertices, as can be understood by comparing Figure 2c,d. The results obtained by considering a lattice such as the one in Figure 2c or the one in Figure 2d are quite different. In the first case, as can be seen in Figure 3c, the quadruplet of degenerate eigenvalues is located right at the middle of the higher energy gap; on the other hand, in the latter case (Figure 3d) the quadruplet is located close to the bottom of the lower energy gap. This result demonstrates that, given a gap, the existence of topological bound states crucially depends on the shape of the vertices.

Before proceeding, it is worth discussing the probability density related to the bound states. This is pictured in Figure 4, for the first in-gap state in the set up of Figure 2a. All other eigenstates behave in a qualitatively analogous way. As expected, the bound states are located at the vertices of the structure. Moreover, their localization length is related to the difference in energy separating them to the lowest closest continuum energy states.

From the analysis just reported, we learn that the existence of the eigenvalues associated to the corner modes in a certain gap is strongly dependent on the disposition of the lattice sites around the vertices in the armchair-zigzag case, while it is not in the bearded-armchair one.

3.2. Energy

By looking at the spectra in Figure 3, it is evident how the energy of the corner modes is strongly influenced by the disposition of the sites at the vertices of the system. In order to be a little more quantitative, we can try to compare the actual results reported above with the predictions of the low energy theory developed in Ref. [71]. If we call m the half-width of the gap and we refer with

E_{mg}

to the midgap energy and with

E_{cm}

to the energy of the corner modes, we can compute the normalized energy of each quadruplet with respect to the corresponding midgap as:

ε_{cm} = \frac{E_{mg} - E_{cm}}{m} .

(4)

The values of such normalized energy for the four spectra in Figure 3 are reported in Table 1.

Recall that the prediction of the low energy theory for the energy of the corner modes was from Equation (3),

ε_{cm}^{th} = cos (\frac{p π}{4})

, with p an integer. Comparing with the values in Table 1, we see that the cases (a) and (c) substantially meet the theoretical value for

p = 3

[71] and

p = 2

respectively. The other two configurations instead, lead to values of the normalized energy which are not compatible with the predictions of the low energy theory. With reference to Figure 2, we can observe that the cases (b) and (d) correspond to samples with lattice sites located just at the vertices. This leads us to the conclusion that the presence of lattice sites just at the corner may act as a sort of on site chemical potential term, causing a shift of the corner mode energies.

Indeed, this statement is corroborated by the results achieved considering the same lattice of Figure 3b, but with the lattice sites located right at the vertices removed, so that the resulting corners are as in Figure 5a. By diagonalizing the Hamiltonian in Equation (1) for this specific lattice configuration, one finds the low energy spectrum reported in Figure 5b. Its comparison with the previous one in Figure 3b shows that the removal of the corner sites modifies the energy of the corner states, which is now visibly higher.

3.3. Degeneracy

In order to assess the degree of protection of the degeneracy of the bound states, we consider a square sample with bearded-armchair edges, with an edge of length

L \approx 347 a

, and we remove the lattice sites located at only one of its corners. In other words, the resulting lattice is just like the ones considered in order to obtain the spectra in Figure 3b and Figure 5b, though with three of the corners as in Figure 2b and with the remaining one as in Figure 5a. Strictly speaking, the system considered would now have lost the original exact

C_{4}

symmetry, which is broken just because of the different disposition of the sites at one of the vertices. On the other hand, the

C_{4}

symmetry is still present concerning the bulk and edges of the system. This should be enough to grant the validity of the low energy theory [71], at least for what it concerns the existence of the bound states. The spectrum resulting from the numerical diagonalization is shown in Figure 6.

One can clearly see that the original quadruplet of Figure 3b splits into a triplet plus a singlet. By carefully examining the precise energies of the eigenvalues, one finds that the triplet is at the energy of the quadruplet in Figure 3b (let us call it

ε

), while the singlet is at the energy of the quadruplet in Figure 5b (let us call it

ε^{'}

). Moreover, one can also plot the probability density of the eigenstates associated with the in-gap eigenvalues (not shown here), finding that corner-modes corresponding to the three degenerate eigenvalues with energy

ε

are localized on the three “untouched” vertices (Figure 2b), while the one with energy

ε^{'}

is localized on the cropped vertex (Figure 5a).

This analysis drives us to two main conclusions: first, the degeneracy of the corner modes is protected just as long as the

C_{4}

symmetry is exact on the whole system (corners included). If we geometrically perturb the system at any vertex causing a local

C_{4}

symmetry breaking, we destroy the degeneracy between the corner modes, though they still survive at different energies inside the gap- if originally present. Second, as we already noted in the previous subsections, the energy of the corner modes inside the edge gap crucially depends on the disposition of the lattice sites at the corners.

4. Conclusions

Our article deals with the quasicrystalline higher order topological insulator predicted in Equation (1). We address three questions related to the properties of the topologically protected zero-dimensional bound states characterizing the model. Such questions are related to the robustness with respect to changes in the structure of the vertices and are the following ones: Do the topological bound states always exist in a given gap? Is their energy fixed? Is their degeneracy protected? We show that the existence of the bound states in a given energy gap is independent of the structure of the vertices if the edges are of the armchair-bearded type, while it is not in the armchair-zigzag case. On the other hand, the energy of the bound states crucially depends on the structure of the vertices in all cases, and the degeneracy is consequently not protected. We hence conclude that if the system is intended to be used as a higher order topological insulator, a precise control of the lattice sites at the vertices is needed.

Author Contributions

Conceptualization, N.T.Z. and M.S.; Programming and calculations S.T. All authors contributed to the interpretation of the results and to the preparation of the draft. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the “Dipartimento di Eccellenza MIUR 2018-2022”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Qi, X.; Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 2011, 83, 1057. [Google Scholar] [CrossRef]
Pesin, D.; MacDonald, A.H. Spintronics and pseudospintronics in graphene and topological insulators. Nat. Mater. 2012, 11, 409. [Google Scholar] [CrossRef] [PubMed]
Michetti, P.; Trauzettel, B. Devices with electrically tunable topological insulating phases. Appl. Phys. Lett. 2013, 102, 063503. [Google Scholar] [CrossRef]
Calvo, M.R.; de Juan, F.; Ilan, R.; Fox, E.J.; Bestwick, A.J.; Mühlbauer, M.; Wang, J.; Ames, C.; Leubner, P.; Brüne, C.; et al. Interplay of Chiral and Helical States in a Quantum Spin Hall Insulator Lateral Junction. Phys. Rev. Lett. 2017, 119, 226401. [Google Scholar] [CrossRef] [PubMed]
Traverso Ziani, N.; Fleckenstein, C.; Dolcetto, G.; Trauzettel, B. Fractional charge oscillations in quantum dots with quantum spin Hall effect. Phys. Rev. B 2017, 95, 205418. [Google Scholar] [CrossRef]
Bendias, K.; Shamim, S.; Herrmann, O.; Budewitz, A.; Shekhar, P.; Leubner, P.; Kleinlein, J.; Bocquillon, E.; Buhmann, H.; Molenkamp, L.W. High Mobility HgTe Microstructures for Quantum Spin Hall Studies. Nano Lett. 2018, 18, 4831. [Google Scholar] [CrossRef]
Kane, C.L.; Mele, E.J. Z₂ topological order and the quantum spin hall effect. Phys. Rev. Lett. 2005, 95, 146802. [Google Scholar] [CrossRef]
Kane, C.L.; Mele, E.J. Quantum spin hall effect in graphene. Phys. Rev. Lett. 2005, 95, 226801. [Google Scholar]
Bernevig, B.; Hughes, T.; Zhang, S. Quantum spin Hall effect and topological phase transition in HgTe quantum wells. Science 2006, 314, 1757. [Google Scholar] [CrossRef]
König, M.; Wiedmann, J.; Brüne, C.; Roth, A.; Buhmann, H.; Molenkamp, L.W.; Qi, X.L.; Zhang, S.-C. Quantum Spin Hall Insulator State in HgTe Quantum Wells. Science 2007, 318, 766. [Google Scholar] [CrossRef]
Linder, J.; Robinson, J.W.A. Superconducting spintronics. Nat. Phys. 2015, 11, 307. [Google Scholar] [CrossRef]
Cayao, J.; Black-Schaffer, A.M. Odd-frequency superconducting pairing and subgap density of states at the edge of a two-dimensional topological insulator without magnetism. Phys. Rev. B 2017, 96, 155426. [Google Scholar] [CrossRef]
Breunig, D.; Burset, P.; Trauzettel, B. Creation of Spin-Triplet Cooper Pairs in the Absence of Magnetic Ordering. Phys. Rev. Lett. 2018, 120, 037701. [Google Scholar] [CrossRef] [PubMed]
Blasi, G.; Taddei, F.; Giovannetti, V.; Braggio, A. Manipulation of Cooper pair entanglement in hybrid topological Josephson junctions. Phys. Rev. B. 2019, 99, 064514. [Google Scholar] [CrossRef]
Keidel, F.; Hwang, S.Y.; Trauzettel, B.; Sothmann, B.; Burset, P. On-demand thermoelectric generation of equal-spin Cooper pairs. Phys. Rev. Res. 2020, 2, 022019. [Google Scholar] [CrossRef]
Kitaev, A. Unpaired Majorana fermions in quantum wires. Phys. Uspekhi 2001, 44, 131. [Google Scholar] [CrossRef]
Ivanov, D.A. Non-Abelian Statistics of Half-Quantum Vortices in p-Wave Superconductors. Phys. Rev. Lett. 2001, 86, 268. [Google Scholar] [CrossRef]
Nayak, C.; Simon, S.H.; Stern, A.; Freedman, M.; Das Sarma, S. Non-Abelian anyons and topological quantum computation. Rev. Mod. Phys. 2008, 80, 1083. [Google Scholar] [CrossRef]
Alicea, J.; Oreg, Y.; Refael, G.; von Oppen, F.; Fisher, M.P.A. Non-Abelian statistics and topological quantum information processing in 1D wire networks. Nat. Phys. 2011, 7, 412. [Google Scholar] [CrossRef]
Fu, L.; Kane, C.L. Josephson current and noise at a superconductor/quantum-spin-Hall-insulator/superconductor junction. Phys. Rev. B 2009, 79, 161408. [Google Scholar] [CrossRef]
Shivamoggi, V.; Refael, G.; Moore, J.E. Majorana fermion chain at the quantum spin Hall edge. Phys. Rev. B 2010, 82, 041405. [Google Scholar] [CrossRef]
Li, J.; Pan, W.; Bernevig, B.A.; Lutchyn, R.M. Detection of Majorana Kramers Pairs Using a Quantum Point Contact. Phys. Rev. Lett. 2016, 117, 046804. [Google Scholar] [PubMed]
Klinovaja, J.; Yacoby, A.; Loss, D. Kramers pairs of Majorana fermions and parafermions in fractional topological insulators. Phys. Rev. B 2014, 90, 155447. [Google Scholar]
Schrade, C.; Zyuzin, A.A.; Klinovaja, J.; Loss, D. Proximity-Induced π Josephson Junctions in Topological Insulators and Kramers Pairs of Majorana Fermions. Phys. Rev. Lett. 2015, 115, 237001. [Google Scholar] [CrossRef] [PubMed]
Keidel, F.; Burset, P.; Trauzettel, B. Tunable hybridization of Majorana bound states at the quantum spin Hall edge. Phys. Rev. B 2018, 97, 075408. [Google Scholar] [CrossRef]
Traverso Ziani, N.; Fleckenstein, C.; Vigliotti, L.; Trauzettel, B.; Sassetti, M. From fractional solitons to Majorana fermions in a paradigmatic model of topological superconductivity. Phys. Rev. B 2020, 101, 195303. [Google Scholar]
Schulz, F.; Plekhanov, K.; Loss, D.; Klinovaja, J. Majorana bound states in topological insulators with hidden Dirac points. Phys. Rev. Res. 2020, 2, 033215. [Google Scholar] [CrossRef]
Laubscher, K.; Chughtai, D.; Chughtai, D.; Klinovaja, J. Kramers pairs of Majorana corner states in a topological insulator bilayer. Phys. Rev. B 2020, 102, 195401. [Google Scholar] [CrossRef]
Fleckenstein, C.; Traverso Ziani, N.; Calzona, A.; Sassetti, M.; Trauzettel, B. Formation and detection of Majorana modes in quantum spin Hall trenches. Phys. Rev. B 2021, 103, 125303. [Google Scholar]
Zhang, F.; Kane, C.L. Time-Reversal-Invariant Z₄ Fractional Josephson Effect. Phys. Rev. Lett. 2015, 113, 036401. [Google Scholar] [CrossRef]
Orth, C.P.; Tiwari, R.P.; Meng, T.; Schmidt, T.L. Non-Abelian parafermions in time-reversal-invariant interacting helical systems. Phys. Rev. B 2015, 91, 081406. [Google Scholar] [CrossRef]
Vinkler-Aviv, Y.; Brouwer, P.W.; von Oppen, F. Z₄ parafermions in an interacting quantum spin Hall Josephson junction coupled to an impurity spin. Phys. Rev. B 2017, 96, 195421. [Google Scholar] [CrossRef]
Fleckenstein, C.; Traverso Ziani, N.; Trauzettel, B. Z₄ parafermions in Weakly Interacting Superconducting Constrictions at the Helical Edge of Quantum Spin Hall Insulators. Phys. Rev. Lett. 2019, 122, 066801. [Google Scholar] [CrossRef] [PubMed]
Sato, M.; Ando, Y. Topological superconductors: A review. Rep. Prog. Phys. 2017, 80, 076501. [Google Scholar] [CrossRef] [PubMed]
Budich, J.C.; Dolcini, F.; Recher, P.; Trauzettel, B. Phonon-Induced Backscattering in Helical Edge States. Phys. Rev. Lett. 2012, 108, 086602. [Google Scholar] [CrossRef] [PubMed]
Schmidt, T.L.; Rachel, S.; von Oppen, F.; Glazman, L.I. Inelastic Electron Backscattering in a Generic Helical Edge Channel. Phys. Rev. Lett. 2012, 108, 156402. [Google Scholar] [CrossRef] [PubMed]
Cheianov, V.; Glazman, L.I. Mesoscopic Fluctuations of Conductance of a Helical Edge Contaminated by Magnetic Impurities. Phys. Rev. Lett. 2013, 110, 206803. [Google Scholar] [CrossRef]
Väyrynen, J.I.; Goldstein, M.; Glazman, L.I. Helical Edge Resistance Introduced by Charge Puddles. Phys. Rev. Lett. 2013, 110, 216402. [Google Scholar] [CrossRef]
Altshuler, B.L.; Aleiner, I.L.; Yudson, V.I. Localization at the Edge of a 2D Topological Insulator by Kondo Impurities with Random Anisotropies. Phys. Rev. Lett. 2013, 111, 086401. [Google Scholar] [CrossRef]
Väyrynen, J.I.; Goldstein, M.; Gefen, Y.; Glazman, L.I. Resistance of helical edges formed in a semiconductor heterostructure. Phys. Rev. B 2014, 90, 115309. [Google Scholar] [CrossRef]
Kainaris, N.; Gornyi, I.V.; Carr, S.; Mirlin, A.D. Conductivity of a generic helical liquid. Phys. Rev. B 2014, 90, 075118. [Google Scholar] [CrossRef]
Väyrynen, J.I.; Pikulin, D.I.; Alicea, J. Noise-induced Backscattering in a Quantum Spin Hall Edge. Phys. Rev. Lett. 2018, 121, 106601. [Google Scholar] [CrossRef] [PubMed]
McGinley, M.; Cooper, N.R. Fragility of time-reversal symmetry protected topological phases. Nat. Phys. 2020, 16, 1181. [Google Scholar] [CrossRef]
Brouwer, P.W.; Duckheim, M.; Romito, A.; von Oppen, F. Probability Distribution of Majorana End-State Energies in Disordered Wires. Phys. Rev. Lett. 2011, 107, 196804. [Google Scholar] [CrossRef] [PubMed]
Brouwer, P.W.; Duckheim, M.; Romito, A.; von Oppen, F. Topological superconducting phases in disordered quantum wires with strong spin-orbit coupling. Phys. Rev. B 2011, 84, 144526. [Google Scholar] [CrossRef]
Pientka, F.; Kells, G.; Romito, A.; Brouwer, P.W.; von Oppen, F. Enhanced Zero-Bias Majorana Peak in the Differential Tunneling Conductance of Disordered Multisubband Quantum-Wire/Superconductor Junctions. Phys. Rev. Lett. 2012, 109, 227006. [Google Scholar] [CrossRef]
Liu, J.; Potter, A.C.; Law, K.T.; Lee, P.A. Zero-Bias Peaks in the Tunneling Conductance of Spin-Orbit-Coupled Superconducting Wires with and without Majorana End-States. Phys. Rev. Lett. 2012, 109, 267002. [Google Scholar] [CrossRef]
Kells, G.; Meidan, D.; Brouwer, P.W. Near-zero-energy end states in topologically trivial spin-orbit coupled superconducting nanowires with a smooth confinement. Phys. Rev. B 2012, 86, 100503(R). [Google Scholar] [CrossRef]
Pikulin, D.I.; Dahlhaus, J.; Wimmer, M.; Schomerus, H.; Beenakker, C. A zero-voltage conductance peak from weak antilocalization in a Majorana nanowire. New J. Phys. 2012, 14, 125011. [Google Scholar] [CrossRef]
Rainis, D.; Trifunovic, L.; Klinovaja, J.; Loss, D. Towards a realistic transport modeling in a superconducting nanowire with Majorana fermions. Phys. Rev. B 2013, 87, 024515. [Google Scholar] [CrossRef]
Roy, D.; Bondyopadhaya, N.; Tewari, S. Topologically trivial zero-bias conductance peak in semiconductor Majorana wires from boundary effects. Phys. Rev. B 2013, 88, 020502(R). [Google Scholar] [CrossRef]
Stanescu, T.D.; Tewari, S. Disentangling Majorana fermions from topologically trivial low-energy states in semiconductor Majorana wires. Phys. Rev. B 2013, 87, 140504(R). [Google Scholar] [CrossRef]
Hui, H.-Y.; Sau, J.D.; Das Sarma, S. Bulk disorder in the superconductor affects proximity-induced topological superconductivity. Phys. Rev. B 2015, 92, 174512. [Google Scholar] [CrossRef]
Klinovaja, J.; Loss, D. Fermionic and Majorana bound states in hybrid nanowires with non-uniform spin-orbit interaction. Eur. Phys. J. B 2015, 88, 62. [Google Scholar] [CrossRef]
Cole, W.S.; Sau, J.D.; Das Sarma, S. Proximity effect and Majorana bound states in clean semiconductor nanowires coupled to disordered superconductors. Phys. Rev. B 2016, 94, 140505(R). [Google Scholar] [CrossRef]
Liu, C.-X.; Sau, J.D.; Stanescu, T.D.; Das Sarma, S. Andreev bound states versus Majorana bound states in quantum dot-nanowire-superconductor hybrid structures: Trivial versus topological zero-bias conductance peaks. Phys. Rev. B 2017, 96, 075161. [Google Scholar] [CrossRef]
Liu, D.E.; Rossi, E.; Lutchyn, R.M. Impurity-induced states in superconducting heterostructures. Phys. Rev. B 2018, 97, 161408(R). [Google Scholar] [CrossRef]
Moore, C.; Zeng, C.; Stanescu, T.D.; Tewari, S. Quantized zero-bias conductance plateau in semiconductor-superconductor heterostructures without topological Majorana zero modes. Phys. Rev. B 2018, 98, 155314. [Google Scholar] [CrossRef]
Fleckenstein, C.; Dominguez, F.; Traverso Ziani, N.; Trauzettel, B. Decaying spectral oscillations in a Majorana wire with finite coherence length. Phys. Rev. B 2018, 97, 155425. [Google Scholar] [CrossRef]
Cao, Z.; Zhang, H.; Lü, H.-F.; He, W.-X.; Lu, H.-Z.; Xie, X.C. Decays of Majorana or Andreev Oscillations Induced by Steplike Spin-Orbit Coupling. Phys. Rev. Lett. 2019, 122, 147701. [Google Scholar] [CrossRef]
Vuik, A.; Nijholt, B.; Akhmerov, A.; Wimmer, M. Reproducing topological properties with quasi-Majorana states. SciPost Phys. 2019, 7, 061. [Google Scholar] [CrossRef]
Kiendl, T.; von Oppen, F.; Brouwer, P.W. Proximity induced gap in nanowires with a thin superconducting shell. Phys. Rev. B 2019, 100, 035426. [Google Scholar] [CrossRef]
Pan, H.; Liu, C.-X.; Wimmer, M.S.; Das Sarma, S. Quantized and unquantized zero-bias tunneling conductance peaks in Majorana nanowires: Conductance below and above 2e²/h. Phys. Rev. B 2021, 103, 214502. [Google Scholar] [CrossRef]
Benalcazar, W.A.; Bernevig, B.A.; Hughes, T.L. Electric multipole moments, topological multipole moment pumping, and chiral hinge states in crystalline insulators. Phys. Rev. B 2017, 96, 245115. [Google Scholar] [CrossRef]
Langbehn, J.; Peng, Y.; Trifunovic, L.; von Oppen, F.; Brouwer, P.W. Reflection-symmetric second-order topological insulators and superconductors. Phys. Rev. Lett. 2017, 119, 246401. [Google Scholar] [CrossRef]
Benalcazar, W.A.; Bernevig, B.A.; Hughes, T.L. Quantized electric multipole insulators. Science 2017, 357, 61. [Google Scholar] [CrossRef]
Schindler, F.; Cook, A.M.; Vergniory, M.G.; Wang, Z.; Parkin, S.S.P.; Bernevig, B.A.; Neupert, T. Higher-order topological insulators. Sci. Adv. 2018, 4, eaat0346. [Google Scholar] [CrossRef]
Fang, C.; Fu, L. New classes of topological crystalline insulators having surface rotation anomaly. Sci. Adv. 2019, 5, eaat2374. [Google Scholar] [CrossRef]
Schindler, F.; Wang, W.; Vergniory, M.; Cook, A.; Murani, A.; Sengupta, S.; Kasumov, A.; Deblock, R.; Jeon, S.; Drozdov, I.; et al. Higher-order topology in bismuth. Nat. Phys 2018, 14, 918. [Google Scholar] [CrossRef]
Noguchi, R.; Kobayashi, M.; Jiang, Z.; Kuroda, K.; Takahashi, T.; Xu, Z.; Lee, D.; Hirayama, M.; Ochi, M.; Shirasawa, T.; et al. Evidence for a higher-order topological insulator in a three-dimensional material built from van der waals stacking of bismuth-halide chains. Nat. Mater. 2021, 20, 1. [Google Scholar] [CrossRef]
Spurrier, S.; Cooper, N.R. Kane-mele with a twist: Quasicrystalline higher-order topological insulators with fractional mass kinks. Phys. Rev. Research 2020, 2, 033071. [Google Scholar] [CrossRef]
Traverso, S.; Sassetti, M.; Traverso Ziani, N. Role of the edges in a quasicrystalline Haldane model. arXiv 2022, arXiv:2207.12996. [Google Scholar]
Ahn, S.J.; Moon, P.; Kim, T.-H.; Kim, H.-W.; Shin, H.-C.; Kim, E.H.; Cha, H.W.; Kahng, S.-J.; Kim, P.; Koshino, M.; et al. Dirac electrons in a dodecagonal graphene quasicrystal. Science 2008, 361, 782. [Google Scholar] [CrossRef] [PubMed]
Haldane, F.D.M. Model for a quantum hall effect without landau levels: Condensed-matter realization of the “parity anomaly”. Phys. Rev. Lett. 1988, 61, 2015. [Google Scholar] [CrossRef]
Bernevig, B.; Hughes, T. Topological Insulators and Topological Superconductors; Princeton University Press: Princeton, NJ, USA, 2013. [Google Scholar]
Moon, P.; Koshino, M. Energy spectrum and quantum hall effect in twisted bilayer graphene. Phys. Rev. B 2012, 85, 195458. [Google Scholar] [CrossRef]
Moon, P.; Koshino, M. Optical absorption in twisted bilayer graphene. Phys. Rev. B 2013, 87, 205404. [Google Scholar] [CrossRef]
Moldovan, D.; Andelković, M.; Peeters, F. pybinding v0.9.5: A Python Package for Tight-Binding Calculations (v0.9.5); Zenodo: Geneva, Switzerland, 2020. [Google Scholar]

Figure 1. Scheme of the lattice construction: one starts from the AA stacking, and then rotates one of the two layers of

30^{\circ}

around an axis passing through the centers of two superimposed hexagons.

Figure 1. Scheme of the lattice construction: one starts from the AA stacking, and then rotates one of the two layers of

30^{\circ}

around an axis passing through the centers of two superimposed hexagons.

Figure 2. A zoom around one vertex for each of the four possible lattice configurations discussed in the main text. Panels (a,b) correspond to lattices presenting bearded-armchair edges, while Panels (c,d) correspond to lattices with zigzag-armchair edges.

Figure 3. Eigenvalues obtained by diagonalizing the Hamiltonian in Equation (1) for square lattices of four different sizes, one for each edge/vertex configuration in Figure 2. We set

λ_{⊥} = 2

and

λ_{H} = 0.3 t

, while all the other fixed parameters are as specified in the main text. The eigenvalues in Panel (a–d) correspond to samples with an edge of approximately

321 a

347 a

280 a

and

306 a

respectively.

Figure 3. Eigenvalues obtained by diagonalizing the Hamiltonian in Equation (1) for square lattices of four different sizes, one for each edge/vertex configuration in Figure 2. We set

λ_{⊥} = 2

and

λ_{H} = 0.3 t

, while all the other fixed parameters are as specified in the main text. The eigenvalues in Panel (a–d) correspond to samples with an edge of approximately

321 a

347 a

280 a

and

306 a

respectively.

Figure 4. Probability density associated to the first in-gap eigenvalue of Figure 2a. On the top-left corner a zoom of the localization pattern of the bound state is reported.

Figure 5. Panel (a): Zoom around the corner of a square lattice such as the one in Figure 2b, but with the corner sites removed. Panel (b): Eigenvalues obtained by diagonalizing the Hamiltonian in Equation (1) for a square lattice with an edge of length

L \approx 347 a

, with the corner sites removed, as depicted in Panel (a). For the diagonalization we set

λ_{⊥} = 2

and

λ_{H} = 0.3 t

as before and all other parameters as specified in the main text.

L \approx 347 a

, with the corner sites removed, as depicted in Panel (a). For the diagonalization we set

λ_{⊥} = 2

and

λ_{H} = 0.3 t

as before and all other parameters as specified in the main text.

Figure 6. Panel (a): Eigenvalues obtained by diagonalizing the Hamiltonian in Equation (1) for a square lattices with an edge of length

L \approx 347 a

, with the sites located at just one of the four corners removed, as in Figure 5a. For the diagonalization we set

λ_{⊥} = 2

and

λ_{H} = 0.3 t

as before and all other parameters as specified in the main text. In Panel (b) a zoom of the spectrum around the gap.

Figure 6. Panel (a): Eigenvalues obtained by diagonalizing the Hamiltonian in Equation (1) for a square lattices with an edge of length

L \approx 347 a

, with the sites located at just one of the four corners removed, as in Figure 5a. For the diagonalization we set

λ_{⊥} = 2

and

λ_{H} = 0.3 t

as before and all other parameters as specified in the main text. In Panel (b) a zoom of the spectrum around the gap.

Table 1. Values of the normalized energy

ε_{cm}

, defined in Equation (4), for the spectra in Figure 3.

Table 1. Values of the normalized energy

ε_{cm}

, defined in Equation (4), for the spectra in Figure 3.

	L	$ε_{cm}$
(a)	$321 a$	−0.708
(b)	$347 a$	−0.772
(c)	$280 a$	+0.010
(d)	$306 a$	−0.430

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Traverso, S.; Traverso Ziani, N.; Sassetti, M. Effects of the Vertices on the Topological Bound States in a Quasicrystalline Topological Insulator. Symmetry 2022, 14, 1736. https://doi.org/10.3390/sym14081736

AMA Style

Traverso S, Traverso Ziani N, Sassetti M. Effects of the Vertices on the Topological Bound States in a Quasicrystalline Topological Insulator. Symmetry. 2022; 14(8):1736. https://doi.org/10.3390/sym14081736

Chicago/Turabian Style

Traverso, Simone, Niccolò Traverso Ziani, and Maura Sassetti. 2022. "Effects of the Vertices on the Topological Bound States in a Quasicrystalline Topological Insulator" Symmetry 14, no. 8: 1736. https://doi.org/10.3390/sym14081736

APA Style

Traverso, S., Traverso Ziani, N., & Sassetti, M. (2022). Effects of the Vertices on the Topological Bound States in a Quasicrystalline Topological Insulator. Symmetry, 14(8), 1736. https://doi.org/10.3390/sym14081736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Effects of the Vertices on the Topological Bound States in a Quasicrystalline Topological Insulator

Abstract

1. Introduction

2. Model

2.1. Lattice

2.2. Model Hamiltonian

3. Results

3.1. Existence

3.2. Energy

3.3. Degeneracy

4. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI