# Local symmetry theory of resonator structures for the real-space control of edge-states in binary aperiodic chains

###### Abstract

We propose a real-space approach explaining and controlling the occurrence of edge-localized gap states between the spectral quasibands of binary tight binding chains with deterministic aperiodic long-range order. The framework is applied to the Fibonacci, Thue-Morse and Rudin-Shapiro chains, representing different structural classes. Our approach is based on an analysis of the eigenstates at weak inter-site coupling, where they are shown to generically localize on locally reflection-symmetric substructures which we call local resonators. A perturbation theoretical treatment demonstrates the local symmetries of the eigenstates. Depending on the degree of spatial complexity of the chain, the proposed local resonator picture can be used to predict the occurrence of gap-edge states even for stronger couplings. Moreover, we connect the localization behavior of a given eigenstate to its energy, thus providing a quantitative connection between the real-space structure of the chain and its eigenvalue spectrum. This allows for a deeper understanding, based on local symmetries, of how the energy spectra of binary chains are formed. The insights gained allow for a systematic analysis of aperiodic binary chains and offers a pathway to control structurally induced edge states.

## I Introduction

Aperiodic systems with deterministic long-range order have long been a subject of intense study, in the endeavor to systematically bridge the gap between crystalline periodicity and complete disorder Maciá (2006). While providing a powerful concept in theoretically modeling the transition to disorder, aperiodic order has become an established property of matter as well. A cornerstone of this was the actual observation of “quasicrystals”—non-periodic but space-filling structures surpassing the crystallographic restriction theorem—by Shechtman Shechtman et al. (1984). In nature quasiperiodicity occurs e.g. in macroscopic constellations such as phyllotaxis Pennybacker and Newell (2013); Maciá (2006). Aperiodically ordered systems even play an important role in material science and technology Maciá (2006, 2012). Owing to their long-range order, they can display interesting physical properties such as a low electrical and thermal conductance Dubois (2012); Maciá (2006), low frictionDubois (2012); Mancinelli et al. (2003) and high hardness Mancinelli et al. (2003). Specific quasicrystalline systems have been shown to enhance solar cells Bauer and Giessen (2013), serve as a catalystYoshimura and Tsai (2002) and could allow for superconductivity Koga (2015); Kamiya et al. (2018).

A general characteristic of aperiodic lattices is the clustering of Hamiltonian eigenvalues into so-called “quasibands” resembling Bloch bands of periodic systems Dal Negro et al. (2016). The corresponding eigenstates generally neither extend homogeneously across the system like Bloch states in regular crystals, nor do they decay exponentially like in disordered systems, and are therefore dubbed “critical” Ostlund and Pandit (1984); Maciá (2014); Kohmoto et al. (1987); Fujiwara et al. (1989); Maciá (2006). In specific cases, quasibands have been shown to originate from the localization of different eigenstates on similar repeated substructures in the system yielding similar eigenenergies de Prunelé and Bouju (2001); de Prunelé (2002); Bandres et al. (2016); Vignolo et al. (2016); Maciá Enrique (2017a). The formation of quasibands typically becomes less distinct with increasing spatial complexity, which in turn can be classified by the structure’s spatial Fourier transform—accordingly altering from point-like to singular continuous to absolutely continuous Maciá (2006); Dal Negro and Boriskina (2012); Negro (2013); Maciá Enrique (2017b). The Fourier spectrum can further be connected to the system’s integrated density of states by the “gap labeling theorem” Johnson and Moser (1982); Delyon and Souillard (1983); Luck (1989); Bellissard et al. (1992); Baake et al. (1993), which assigns characteristic integers to the gaps between quasibands.

As ordered lattice systems are truncated in space into finite setups, they may support the occurrence of eigenstates localized along their edges, energetically lying within spectral gaps. In periodic systems, such edge states (or ‘surface states’ Davison and Steslicka (1996)) may or may not appear depending on how the underlying translation symmetry is broken by the lattice truncation, that is, where in the unit cell the system is cut off Zak (1985). In various types of periodic setups, edge states can also be given a topological origin in terms of nontrivially valued invariants (winding numbers) assigned to the neighboring Bloch bands Bernevig (2013). This has boosted an intensive research activity in the field of topological insulators Kane and Mele (2005); Hasan and Kane (2010); Qi and Zhang (2011) and the quest for interesting novel materials and applications Tian et al. (2017), including e. g. robust lasing via topological edge-states in periodic photonic lattices St-Jean et al. (2017).

Edge states may also be present between quasibands in aperiodic systems, as has been shown for binary 1D systems Zijlstra et al. (1999a); El Hassouani et al. (2006); Lei et al. (2007); Pang et al. (2010); Martínez and Molina (2012) and recently demonstrated for 2D photonic quasiperiodic tilings Bandres et al. (2016). Notably, also here a topological character can be assigned to the edge states in correspondence to the system’s bulk properties. Indeed, a position-space based topological invariant, the so-called Bott index Loring and Hastings (2010), can be applied to aperiodically structured Bandres et al. (2016) or even amorphous systems Agarwala and Shenoy (2017). Moreover, for 1D quasiperiodic systems the winding of edge state eigenvalues across spectral gaps coincide with the gap labels mentioned above Johnson and Moser (1982); Negi and Ramaswamy (2001); Kellendonk and Zois (2005); Kraus and Zilberberg (2012), which have recently also been measured in scattering Baboux et al. (2017) and diffraction Dareau et al. (2017) experiments. Remarkably, edge modes occur also as scattering resonances in open systems with different types of deterministic aperiodic order incorporating long-range couplings between lattice constituents, as demonstrated very recently in terms of the eigenmodes of full vectorial Green matrices Wang et al. (2018).

The ubiquitous presence of edge states in aperiodic systems indicates that it derives primarily from the underlying geometrical structure and not from model-specific assumptions. Departing from periodicity, however, there is no translation symmetry whose breaking (at the boundary) would provide a mechanism for edge state formation. On the other hand, aperiodic systems are imbued with local symmetries, that is, different spatially symmetric substructures are simultaneously present in the composite system which possesses many different domains of local symmetries. Indeed, local “patterns” are known to occur repeatedly in deterministic aperiodic systems, as expressed by Conway’s theorem Gardner (1977). In the specific case of 1D binary lattices, local reflection symmetry is abundantly present and follows, at each scale, a spatial distribution closely linked to the underlying aperiodic potential sequence Morfonios et al. (2014). The encoding of such local symmetries into generic wave excitations have recently been described within a theoretical framework of symmetry-adapted nonlocal currentsMorfonios et al. (2017), which obey generalized continuity equations Morfonios et al. (2017); Röntgen et al. (2017); Spourdalakis et al. (2016); Wulf et al. (2016) and whose stationary form allows for a generalization of the parity and Bloch theorems to locally restricted symmetries Kalozoumis et al. (2014) as well as a classification of perfect transmission Kalozoumis et al. (2013). In the context of finite, aperiodically ordered setups, an appealing question is whether a real-space picture for the formation—and thereby control—of edge states can be brought into connection with local symmetries.

In the present work we propose an intuitive real-space picture of the formation of quasibands and edge states in binary aperiodic tight-binding chains. The approach is based on the analysis of eigenstate profiles in the limit of weak inter-site coupling. In this regime, eigenstates generically fragment, i.e., have non-negligible amplitudes only on a small number of sites, as we show by means of a perturbation theoretical treatment. The amplitudes on these fragments are in almost all cases locally symmetric and can be identified as local resonator modes (LRM), i.e., eigenmodes of local resonators embedded into the full chain. Here, a resonator denotes a substructure that can confine, at certain energies, the wavefunction within its interior. The LRMs can be used to classify states, and those belonging to quasibands are composed of repeated LRMs hosted by resonators within the bulk, while edge states are composed of unique LRMs occurring on the edge. We further investigate the reasons for the formation of quasibands by linking the energy of a state to that of its constituting LRMs, where the energy of an LRM is defined as its energy in the corresponding isolated resonator. From this finding, we see that the multiple occurrence of identical resonator structures automatically leads to the formation of quasibands by their capability of hosting identical (and thus degenerate) LRMs. We further use this energetical insight to move a given edge state into a quasiband by performing tailored changes to the corresponding resonators on the edge. The inference of those properties to moderate inter-site coupling depends on the type of aperiodic order used in the model. We here apply the approach to the prominent representatives of three main classes of structural complexity: Fibonacci, Thue-Morse, and Rudin-Shapiro chains, featuring point-like, singular continuous, and purely singular spatial Fourier spectra, respectively.

The paper is organized as follows. In Sec. II we introduce our setup and show examples of quasibands and edge states in Fibonacci chains. We then develop our approach to edge states based on locally symmetric resonators and apply it to Fibonacci chains in Sec. III, to Thue-Morse chains in Sec. IV.1 and to Rudin-Shapiro chains in Sec. IV.2. In Sec. V we comment on the generality of our framework and on the connection to related work. We conclude our paper and give an outlook in Sec. VI. A perturbative treatment demonstrating the localization onto reflection-symmetric resonators is provided in the appendix, together with further technical details including proofs of major statements, complementary explanations, and further comments.

## Ii Prototype quasiperiodic order: The tight-binding Fibonacci chain

We consider a finite one-dimensional chain of sites with real next-neighbor hoppings described by the Hamiltonian

(1) |

where is the onsite potential of site . In the basis of single site excitations , the above Hamiltonian can be written as a tridiagonal matrix

(2) |

Such a tight-binding chain is used in a plethora of interesting model systems, examples including the Aubry-Andre Aubry and André (1980) model relevant in the study of localizationRoati et al. (2008) and the Su-Schrieffer-Heeger model, a simple prototypical chain supporting a topological phase Asbóth et al. (2016). It also effectively describes one-dimensional arrays of evanescently coupled waveguides Szameit et al. (2012); Efremidis et al. (2010). We here fix the hoppings to a uniform value and restrict the onsite elements to be “binary”, that is, the sites are of two possible types and , and the take on corresponding values and , with the contrast defined as

(3) |

Without loss of generality we will set and throughout, having a single control parameter for a given chain.

In the following, we will investigate the spatial profiles of the eigenvectors of in relation to their eigenvalues , given by

(4) |

for chains with aperiodic order. Note that represents a generic finite tight-binding chain; the choice corresponds, e.g., to the kinetic energy of electrons on a lattice with onsite potential , while (made here) can be used to model the coupling of photonic waveguides Szameit et al. (2012); Efremidis et al. (2010) with propagation constants . Our analysis remains qualitatively unaffected by this choice.

We start by presenting the eigenstates and spectral properties of a finite binary chain following the Fibonacci sequence Negro (2013), a prototypical case of quasiperiodic order. This will serve as an initial point motivating the development of a local resonator approach at high contrast in the next section. Starting with , the sequence is constructed by repeatedly applying the inflation rule , resulting in . This sequence is then mapped onto the onsite elements of the tight-binding chain. The spectrum and eigenvectors of this chain are shown in Fig. 1 for a moderate contrast of and sites. Despite the lack of periodicity, the eigenvalues cluster into so-called quasibands, owing to the long-range order present in the Fibonacci chain Enrique Maciá Barber (2008), and the spectrum attains a self-similar structure of quasibands and gaps in the limit. For presentation reasons, we have here chosen large enough to anticipate this spectral feature, though small enough to visually discern the spatial characteristics of the eigenmodes.

The quasibands are occupied by bulk eigenmodes extending along the interior of the chain. Those are known to follow a power-law decay in space, lying between the exponential decay of modes in a randomly disordered chain and uniformly extending Bloch eigenmodes in periodic chains Kohmoto et al. (1987); Fujiwara et al. (1989). Such modes have recently been shown to consist of locally resonating patterns (i.e., characteristic sequences of amplitudes) which occur on repeating segments of a quasiperiodic structures and are characteristic for a given quasiband Vignolo et al. (2016); de Prunelé (2002); Macé et al. (2016). This is particularly visible for the bulk modes of the uppermost quasiband in Fig. 1. A close inspection reveals that the bulk mode profiles tend to localize into locally reflection-symmetric peaks (see black subregions of high amplitude for a given mode). Those in turn follow the distribution of local symmetry axes (or centers of “palindromes” Droubay (1995)) which are hierarchically present in the Fibonacci chain Morfonios et al. (2014), as seen by comparison with the bar plot on the top. Each bar shows the maximal size of a continuous domain of reflection symmetry centered at position , where can refer here to sites () or to links between sites (). For instance, as the first few characters of are

we have and .

Within the gaps between quasibands there may appear spectrally isolated modes, reminiscent of gap modes localized on defects within a periodic lattice Grundmann (2016); Povinelli et al. (2001). For the example given in Fig. 1, i.e., an unperturbed but finite Fibonacci chain, the gap modes are known Niu and Nori (1990); Zijlstra et al. (1999a); Nori and Rodriguez (1986) to be localized at the edges, decaying exponentially into the bulk.

The control of edge states by local changes in the underlying potential sequence is a central aspect of this work. Our approach is that, due to their localization, the occurrence and spectral position of edge states can be influenced by local modifications on the corresponding edge of the aperiodic lattice. We demonstrate the feasibility of this approach in Fig. 2 by using the following representation Baboux et al. (2017) of the Fibonacci potential sequence:

(5) | |||

(6) |

where is the inverse golden mean and the integer site index runs from to . By continuously varying the so-called “phason” , localized flips are induced at discrete values of , forming a two-dimensional pattern in the plane, see Fig. 2. The finite chain of length constitutes a different segment (or “factor”) of the infinite Fibonacci sequence after each flip Baboux et al. (2017). This allows to investigate different Fibonacci-like configurations while maintaining a constant length . In Fig. 2, we visualize the effect of these flips on the energy spectrum, shown in orange. As one can see, the gap states in the purple rectangle, which are localized on the right edge (not shown here), are influenced only by flips acting on this edge, marked by green circles. From bottom to top, the green flips (i) create the edge state (ii) and (iii) modify its energy and (iv) finally annihilate it. Note that in general for processes of type (ii) and (iii), the energetical change accompanying the change of the edge is stronger for a flip near to the edge than for a flip more distant to the edge.

The occurrence of such edge-localized gap states in a finite 1D quasiperiodic potential was recently very elegantly described within a scattering setting Baboux et al. (2017); Dareau et al. (2017) in a continuous system as a consequence of a resonance condition when varying the phason . At the same time, the connection of the winding of to invariant integers labeling the spectral gaps of the quasiperiodic structure through the so-called “gap labeling” theorem Fu et al. (1997), renders the nature of the 1D edge states topological Kraus and Zilberberg (2012). On the other hand, the flip-induced edge state creation/annihilation demonstrated in Fig. 2 suggests that their origin could also be explained by viewing chain edges as a generalized type of “defects” to the quasiperiodic long-range order. In the following, we will develop this idea in terms of the prototype Fibonacci chain. Our aim is to provide a simple and unifying real-space picture for the appearance of edge states in the energy gaps of non-periodic structures.

## Iii Edge modes from truncated local resonators

The analysis of eigenstates at high contrast [see Eq. 3] is at the heart of our approach, revealing structural information that would remain hidden at lower contrast. Once this information is retrieved, we leverage it to develop a generic framework for the understanding and manipulation of quasibands and edge states in binary tight-binding chains.

In the following, we will focus on a Fibonacci chain, choosing a relatively small size for easier treatment and visualization. The slight modifications needed for the treatment of longer chains are commented on in Appendix E. We split the presentation into three subsections, covering the concept of local resonator modes (Sec. III.1), the structural control of edge states (Sec. III.2), and the behavior at low contrast (Sec. III.3).

### iii.1 Local resonator modes and local symmetry

Our starting point is an analysis of eigenstates at high contrast. Those are shown in Fig. 3 (a) for a th generation Fibonacci chain ( sites) with relatively high contrast . We see that each eigenstate is pinned to a small number of sites where is has non-negligible amplitude, practically vanishing on the remaining sites. This is quite different from the states at low contrast (like in Fig. 1) which are smeared out along the whole chain. An impression of how the transition between those two regimes takes place is given in Fig. 3 (b), showing the amplitudes of a bulk and edge state for varying contrast. When increasing the contrast, the spatial profile of the bulk state becomes gradually fragmented: The amplitudes on most sites become suppressed, and a characteristic remnant of the initial extension appears on the sites with substantial amplitude. The fragmentation with increasing contrast also occurs for the edge state, with the difference that there is only a single fragment remaining; in this case the -site on the left edge.

The fragmentation at high contrast is rigorously predicted within a perturbation-theoretical treatment provided in Appendix C, applying to generic binary tight-binding chains. Its results promote an intuitive picture for the cause of fragmentation, where a chain is viewed as a collection of local resonators. A local resonator is here a discrete substructure which, at high contrast, confines the wavefunction within its interior for a certain eigenenergy. The simplest case consists of a three-site structure , where the vertical lines demarcate the resonator “cavity” (the inner part ) from its “walls” (the outer parts ). The resonator character of this particular substructure is analyzed in more detail in Appendix A. Two such resonators can be combined to form a double resonator , formed by overlapping one wall of each . Note that, for a substructure to function as a local resonator, either (i) the resonator wall and its next-neighboring site in the cavity must be of different type or (ii) the resonator wall must coincide with one of the edges ( or , with ).

We now link the resonator concept to the eigenstate fragmentation seen in Fig. 3. As an example, each fragment of (indicated by orange rectangles in Fig. 3 (b)) is localized on the ’s of the local resonator . We denote this fact as , which represents an eigenmode of the isolated resonator and which we will call a local resonator mode (LRM). The overlines here indicate sites with equally signed and relatively much higher amplitude than non-overlined sites; see Appendix A. At high contrast the state can thus be seen as a collection of identical, non-overlapping LRMs (one on each fragment) with negligible amplitudes on the parts in between. In the same manner, each eigenstate shown in Fig. 3 (a) is composed of identical LRMs. In particular, we notice that all states in a given quasiband are characterized by the same resonator mode, different from that of other quasibands. This is shown on the right side of the figure, where LRMs are depicted schematically. Here, overlines and underlines in an LRM such as denote amplitudes of opposite sign. Contrary to the bulk states of quasibands, edge states feature unique resonator modes which are not repeated elsewhere in the chain, with the underlying resonators located at (one of) the chain edges. We thereby distinguish these two types of LRMs as bulk and edge LRMs (-LRMs and -LRMs, respectively).

The fact that each quasiband is characterized by a single resonator mode can be understood as follows. If a given eigenstate of energy is composed of non-overlapping LRMs such that has very low amplitude on the next-neighboring sites of the corresponding resonators, then each of the energies of those LRMs (that is, their eigenenergies in the isolated underlying local resonator) must fulfill . This statement is proven rigorously in Appendix D. Now, applying the perturbative treatment of Appendix C to the chain of Fig. 3 shows that for any two LRMs to be nearly degenerate energetically they must be identical. Thus, each quasiband—having quasidegenerate modes at high contrast—is characterized by a single LRM.

A similar reasoning explains why bulk states of quasibands are more spatially extended compared to edge states lying in spectral gaps. Indeed, due to the quasiperiodicity of the Fibonacci chain, any local resonator (that is, a binary substructure) in the bulk occurs repeatedly (though not periodically) along the chain—specifically, at spacings smaller than double its size. This is a general result known as Conway’s theorem Gardner (1977). Thus, a -LRM hosted by a given local resonator will also be correspondingly repeated along the chain. If the -LRM has energy , then a state with energy is allowed to simultaneously occupy all copies of this -LRM, and is accordingly spatially extended. Edge states, on the other hand, consisting of -LRMs at high contrast, correspond to local resonators induced by the presence of an edge, which breaks the quasiperiodicity. Due to this truncation at the edge (e. g. of the type with at the right edge, compare Fig. 3 (a)), the -LRM generally does not match the energy of any -LRM, and therefore cannot occupy multiple local resonators in the bulk: The eigenstate is confined to the edge, lying energetically isolated in a gap. This is visualized by the marked edge states in Fig. 3.

A remarkable observation in Fig. 3 (a) is further that each local resonator hosting a -LRM is reflection symmetric, such that all isolated -LRMs have definite parity; see schematic on the right. This means that, at high contrast, the fragments (occupied local resonators) of quasiband eigenstates feature, to a very good approximation, local parity with respect to local reflection symmetries of the chain. The positions and sizes of all such local symmetries are shown in the top panel of Fig. 3 (a). An example is given by the state which is locally symmetric around, e. g., the position , and corresponds to the -LRM , as indicated by the green boxes. This behavior is predicted by the higher-order-degenerate perturbative treatment of Appendix C. In particular, it demonstrates the crucial role of local reflection symmetries in the eigenstate localization profiles of binary aperiodic chains. A promising direction of research would thus be to treat this class of systems within the recently developed theoretical framework of local symmetries Kalozoumis et al. (2014); Morfonios et al. (2017); Röntgen et al. (2017).

### iii.2 Structural control of edge states

Having understood the real-space mechanism for the formation of edge-localized gap states in Fibonacci chains, we can now utilize this insight to systematically manipulate these states. In particular, let us show how structural modifications at the edges of a Fibonacci chain can selectively “annihilate” a given edge mode. Note that whether or not one considers a particular state localized (near or on) the edge to lie in an energetical gap is obviously a question of the scale under consideration. This is due to the fact that any finite chain naturally has a discrete spectrum, for which, strictly speaking, no continuous energy-bands are defined. In the remainder of this work, we will solely consider states as gap-edge ones provided that, at a contrast of , they lie in a clearly visible energetical gap. This simplifies our treatment, and in Appendix E we comment on the extensions of this simplification.

For definiteness, we consider the edge state (orange circle) of the chain in Fig. 4 (a) which simply focuses on states to of Fig. 3 (a). This state corresponds to the the -LRM (see right side of Fig. 3 (a)). The underlying resonator is a left-truncated version of the resonator , which hosts the -LRM characterizing the quasiband below (states to ). Now, as shown in Fig. 4 (b), if we complete the resonator into by attaching a site to the left end of the chain, then the edge can accommodate the -LRM instead of the -LRM . Consequently, the edge mode is replaced by a bulk mode of the quasiband. In other words, the edge state is “absorbed” into a quasiband by converting the -LRM of the former to the -LRM of the latter through a structural modification at the edge. This intuitive procedure can be applied similarly for the other two pronounced edge states (green and magenta circles) in Fig. 4 (a), by completing the corresponding the corresponding edge local resonator into a bulk one. Thus, the selective control of a specific edge state is possible.

Let us note, however, that in most cases such a selective annihilation of one edge state leads to the creation of one (or more) other edge state(s) located elsewhere in the spectrum, as a result of the edge modification. For example, the left edge of the modified chain in Fig. 4 (b) features the resonator , which is a truncated version of hosting the -LRM , thus yielding a new gap-edge state (not shown).

Interestingly, in special cases this issue can be overcome by exploiting the local symmetry of bulk resonators, as we now explain using the example shown in Fig. 4 (c). Here, an site is attached to the right edge, which formerly hosted the -LRMs and (cf. Fig. 4 (a)), corresponding to the edge states (green) and (orange), respectively. In the modified chain, the right edge features a local resonator . The key point now is that this resonator supports two LRMs, and , which are degenerate to the b-LRMs and , respectively, due to the reflection symmetry of the underlying resonator . This symmetry-induced degeneracy is shown rigorously in Appendix B. As a result of their degeneracy, the respective e- and -LRMs can combine linearly to compose quasiband states, as seen in Fig. 4 (c); see states in first and third quasiband from bottom with marked edge resonators. The same procedure can be performed on the left edge by removing an site from it, leaving the edge resonator hosting which is degenerate to (see state in top quasiband with orange marked left edge resonator in Fig. 4 (c)). Note that both (right and left) edge modifications above are consistent with the Fibonacci order: The resulting chain is simply obtained from the former one by a single-site shift to the right along the infinite Fibonacci chain. We thus have a case of finite Fibonacci chain with no edge-localized gap states.

From the above it is clear that edge states in binary quasiperiodic chains can now be understood and manipulated within the framework of local resonators. Structural creation and annihilation represents a first fundamental step in edge state control. Indeed, once an edge state is established, its energetic position within a gap can further be tuned by allowing for non-binary (freely varying) potentials at the edges, while leaving the quasibands intact.

### iii.3 Behavior at low contrast

The local understanding and controllability of edge states at high contrast levels raises the question if these properties are retained also at lower contrast. To address this, in Fig. 5 (a) we show the eigenvalue spectrum of the original Fibonacci chain studied previously [Fig. 3 (a)] for varying contrast . As we see, gap states (localization on edges not shown here) are clearly distinguished for all contrast levels. Figure 5 (b) shows the spectrum of the modified Fibonacci chain of Fig. 4 (c)—where all edge states were annihilated at high contrast—for the same contrast values. Also here the structure of the spectrum is retained with varying . In particular, a real-space analysis (not shown here) confirms that all quasiband states in the modified chain remain extended in the bulk for varying . The effect of lowering the contrast is merely a reduction in the fragmentation of the eigenstate profiles which become more smeared out into regions between the LRMs defined at high contrast.

This finding indicates that the mechanism of edge state formation via truncated local resonators based on an analysis at high contrast, remains valid also for lower contrast, though “hidden” due to the spatial smearing of the states. In other words, the contrast parameter can be used as an intermediate tool to manipulate edge states in binary aperiodic model chains: It is first ramped up to reveal the eigenstate structure in terms of LRMs subject to modifications, and then ramped down again with the bulk/edge-state separation retained.

## Iv Application to non-quasiperiodic chains

Featuring a point-like spatial Fourier spectrum (rendering it, by definition, a quasicrystalMaciá (2006); Dal Negro and Boriskina (2012)), the Fibonacci chain studied above represents the class of lowest structural complexity when departing from periodicity towards disorder, as mentioned in Sec. I. The question naturally arises whether the local resonator framework developed in Sec. III, distinguishing edge from bulk states via LMRs, applies also to other classes of aperiodic chains. In the following, we will demonstrate the generality of our approach by applying it two cases of qualitatively different structural character, the Thue-Morse and Rudin-Shapiro chains. We thereby essentially go though the same analysis steps as in Sec. III—identification of LRMs, edge state control, and low contrast behavior—and assess the particularities of each structural case.

### iv.1 Singular continuous Fourier spectrum: Thue-Morse chain

A well-studied case of aperiodic order which is not quasiperiodic is the Thue-Morse sequence Maciá (2006), produced by the inflation rule yielding .
Its Fourier spectrum is singular continuous, and from this viewpoint it is considered Kolář et al. (1991) closer to the disorder limit (with absolutely continuous spectrum Kroon and Riklund (2004)) than quasiperiodic order (with point-like spectrum).
On the other hand, a subset of eigenstates of the Thue-Morse chain strongly resemble those of periodic chains Riklund et al. (1987).
The eigenstates of a -site Thue-Morse chain ^{1}^{1}1Note that we have not used here a generation of the Thue-Morse sequence (of length Enrique Maciá
Barber (2008) for the -th generation), but the same length as for the Fibonacci chain.This is done for the sake of comparison but also to avoid certain symmetries Noguez and Barrio (2001) of the Thue-Morse sequence generations (spatial mirror symmetry for odd and spectral mirror symmetry around for even ), in favor of the generality of the analysis. are shown in Fig. 6.
Indeed, while some bulk states are more strongly localized into subdomains than in the Fibonacci chain for equal contrast (compare to Fig. 1), others are more extended along the chain.
As we see in Fig. 6, in spite of the quasiband structure being more fragmented, there occur well distinguishable states within gaps which are localized on one of the chain edges.

Local resonator modes. Like in Sec. III, we consider a smaller chain of sites to visually facilitate the detailed spatial analysis. Its eigenstates are shown in Fig. 7 (a) for contrast , together with the distribution of local reflection symmetries in the chain (top). As we see, the bulk state profiles are fragmented in a well-defined manner for different quasibands: Like in the Fibonacci case, each bulk state is composed of copies of a -LRM characterizing the corresponding quasiband, as indicated schematically on the right of the figure. In contrast, the three occurring prominent edge states (marked by colored circles) consist of non-repeated -LRMs at one of the chain ends. Like before, the local resonators underlying the -LRMs can be identified as truncated local resonators underlying the -LRMs. This demonstrates that our LRM-based framework for the formation of edge states applies also for this class of aperiodic order. Notably, also here the -LRMs have definite local parity under reflection, and are present in the eigenstates following the local symmetry axes shown in the bar plot (top of Fig. 7 (a)). This is indeed predicted by the perturbation theory of Appendix C. We thus see that also for the Thue-Morse chain its local symmetries essentially provide the regions of localization of the eigenstate fragments at high contrast.

Edge state control. Following the same principle as in the Fibonacci case, Fig. 7 (b) shows how two edge states (marked by green and orange circles in Fig. 7 (a)) are annihilated by attaching a site to the right end of the original chain. Indeed, those edge states were localized on the truncated resonator which is completed to and can thus host the -LRMs and , so that the edge states are “absorbed” into the corresponding quasibands. However, the right edge of the modified chain now features a new edge state (marked by blue circle) with resonator mode (its previous absence is indicated by a in Fig. 7 (a)). It lies, energetically, in the gap just below the quasiband with -LRM . Note that, as expected from our real-space local resonator picture, the left edge state (orange circle) remains unaffected by the present modification on the right edge of the chain, since it is localized on the opposite edge.

Contrast variation. Finally, we investigate how edge states and quasibands behave for lower contrast in the Thue-Morse chain. Figure 8 (a) shows the spectrum of the chain of Fig. 7 (a) for varying contrast, starting from . As we see, the three edge modes in the spectral gaps are clearly visible also at lower contrast levels. The spectrum of the modified chain (with left-attached site) for varying is shown in Fig. 8 (b). As is highlighted by the black circle and the ellipse, the two former edge states are absorbed into the neighboring quasibands (as shown in Fig. 7 (a)) for all considered contrast levels. Also, the left edge state as well as the modification-induced right edge state (orange and blue in Fig. 7 (b)) remain in their gaps as the contrast is varied. Overall therefore, the impact of the modifications persists at lower contrast levels.

### iv.2 Absolutely continuous Fourier spectrum: Rudin-Shapiro chain

Taking a step towards higher structural complexity, we finally investigate the case of a Rudin-Shapiro chain in terms of our local resonator framework. The Rudin-Shapiro sequence Kroon et al. (2002) is obtained by the inflation rule , , , , yielding for an initial seed . Its Fourier spectrum is absolutely continuous, a property shared with completely disordered chains Kroon and Riklund (2004). Further, there are indications that the tight-binding Rudin-Shapiro chain has both exponentially and weaker-than-exponentially localized eigenstates Maciá (2014); Dulea et al. (1993, 1992a, 1992b), while even extended ones have been shown to exist at low contrast. The different character of the Rudin-Shapiro states compared to the Fibonacci or Thue-Morse chain can be anticipated from the eigenstate map shown in Fig. 9. As we see, at this low contrast () there is now no clear distinction between bulk and edge states. Moreover, no clear energetic clustering into well-defined quasibands is present. Note also that the distribution of local reflection symmetries along the chain (see top of figure) is much less structured than in the Fibonacci or Thue-Morse chains (cf. top of Figs. 6 and 1), with overall smaller symmetry domains present. At the same time, there is clustering of symmetry axes with gaps in between, caused by the occurrence of larger contiguous blocks of single type (up to four or sites in a row) along the sequence. In the following we show that there is still a strong link of the eigenstates and spectral features of the Rudin-Shapiro to the presence of locally symmetric resonators.

Local resonator modes. For the high-contrast analysis, we consider a Rudin-Shapiro chain of sites. The size is now chosen slightly larger in order to better reflect the structural properties of Rudin-Shapiro sequence. Indeed, in accordance with its higher complexity, a given substructure will here repeat at relatively larger distances along the sequence. It may thus occur only once in a too short chain, thereby obscuring its long-range order. Figure 10 (a) shows the eigenstate map of the considered chain at contrast . We see that also here the eigenstates fragment onto locally symmetric substructures, and are again composed of -LRMs corresponding to clustered eigenvalue quasibands, as shown on the right. The difference is now that there are many more different identified -LRMs compared to the Fibonacci and Thue-Morse chains. This is because the increased number of contiguous block sizes allows for a higher diversity of local resonator substructures, with larger resonators additionally hosting a larger number of different LRMs each. In turn, there is a higher possibility that different -LRMs have (nearly) the same energy, since the different resonators may have partially overlapping individual eigenspectra. Therefore, it may now more easily occur that different LRMs participate in the same eigenstate (to which they are quasidegenerate; see Appendix D). An example of this are the states indicated by the green ellipse in Fig. 10 (a): Each of them consists of a on the left and two on the right, consecutively overlapping by one site. The emergence of such modes is explained in detail by means of perturbation theory in Appendix C. Further, edge states appear which localize on corresponding -LRMs. Those are now, however, energetically not as clearly distinct from the clustered eigenvalues of quasibands as in the Fibonacci and Thue-Morse cases. For example, the states marked by blue circles are localized on the left edge, but are composed of the -LRMs (from top to bottom) , , , which are nearly degenerate to the -LRMs of the corresponding quasibands (see right side of figure). Nevertheless, there are also well-distinguished edge states lying in gaps (though close to gap edges) marked by purple circles.

Edge state control. We now manipulate the two states marked by purple color which are localized on the right-edge. These localize on the truncated resonator , and their energy is different from the energy of states localized on the complete resonator which occurs twice in the bulk. In (b) we add a on the right edge of the chain, completing this resonator. Due to this completion, the two former gap-edge states move into the respective energy cluster (or quasiband).

Contrast variation. In Fig. 11 (a) we investigate the eigenvalues of the Rudin-Shapiro chain of Fig. 10 (a) for varying contrast. Compared to the case of the Fibonacci chain presented in Fig. 5 or the Thue-Morse chain presented in Fig. 8, the energetic clusters form only at high contrast values. This already indicates that modifications to the chain done at high contrast can not directly be traced to energetic changes at low contrast as was the case for the Fibonacci and Thue-Morse chain. This can also be seen for the two edge states marked by horizontal lines in Fig. 11 (a). At high contrast, these are caused by a truncated resonator on the right edge. In Fig. 10 (b) we have completed this resonator, causing the two edge states to move (at high contrast) closer to the nearest eigenvalue cluster. As can be seen in Fig. 11 (b), this manipulation is only effective at high contrast. For low contrast, the eigenvalue structure is nearly unchanged compared to the original chain shown in Fig. 11 (a).

In conclusion, we have applied our local-symmetry based resonator strategy to the Thue-Morse and the Rudin-Shapiro chain. The results show that our approach can be used to explain and control gap-edge states of the Thue-Morse chain. At high contrast the gap-edge states of the Rudin-Shapiro chain are likewise explained. However, our approach can not be used to make qualitative statements at low contrast.

## V Applicability and relation to other approaches

Let us briefly comment on the limitations of applicability of the developed framework and its connection to similar approaches in the literature. The presented methodology essentially relies on the fragmentation of eigenstates at high contrast and can thus only be applied onto chains featuring such a behavior of eigenstates. A perturbation theoretical treatment of binary tight-binding chains which serves as basis for our methodology, see Appendix C, indicates that at high contrast the fragmentation of eigenstates is indeed the generic case. However, the necessary conditions for this behavior still need to be determined in order to clarify the range of applicability.

The connection between local resonators and quasibands in quasiperiodic setups has been commented on in Ref. Zijlstra et al., 1999b; Maciá Enrique, 2017a; Bandres et al., 2016. For the Thue-Morse sequence, a similar analysis has been achieved in Ref. Riklund et al., 1987. However, to the best of our knowledge, there is no systematic framework bringing together the three concepts of LRMs, quasibands and edge states into a unified context. An approach related to ours is the renormalization group flow analysis. For the tight-binding chains, this method aims at understanding the energetic behavior of a chain through a series of size reductionsNiu and Nori (1990). At each step, the size of the system is decreased, and the behavior of the decreased one is linked to the bigger one by a renormalization procedure, usually done in terms of perturbation theory. The renormalization group flow is a powerful method, and has been successfully used to explain the fractal nature of the Fibonacci spectrum Kalugin et al. (1986); Piéchon et al. (1995); Liu and Sritrakool (1991); Niu and Nori (1990, 1986). However, it needs to be tailored to the system of interest, and as stated in Ref. Niu and Nori, 1986, finding an appropriate renormalization group flow for a general quasiperiodic chain is not easy. This stands in contrast to the very general method proposed in this work, which was shown to be applicable to a broad range of different setups.

## Vi Conclusions and outlook

We have presented a systematic approach to the analysis of aperiodic binary tight-binding chains regarded as a combination of different resonator-like subsystems rather than a single bulk unit. For low inter-site coupling, each eigenstate is seen to be composed of spatially non-overlapping local resonator modes of these resonator structures. This viewpoint, supported by a rigorous perturbation theoretical treatment, allows for an intuitive explanation of the emergence of both quasibands and gap-edge states in such chains. We demonstrate the power of our approach by applying it to Fibonacci, Thue-Morse and Rudin-Shapiro chains and show how gap-edge states occurring in these chains can be manipulated.

A repeating motif in our analysis of eigenstates at high contrast is the fact that most resonator modes share the local symmetries of the underlying systems. This strong impact of local symmetries is remarkable, especially as it is hidden at lower contrast levels by a substantial background in the eigenstate profiles. In this work we have given an explanation for this finding at high contrast, and we believe that the study of local symmetries in complex setups is a very promising field with rich perspectives and potential applications. The recently established framework of local symmetries Kalozoumis et al. (2015, 2014, 2013); Morfonios et al. (2017, 2017); Röntgen et al. (2017); Zampetakis et al. (2016) provides dedicated tools for this purpose, and extensions of it are of immediate relevance. In this line, our work may enable the local-symmetry assisted design of novel optical devices that support desired quasiband structures and strongly localized edge states at prescribed energies, offering exciting opportunities to control light-matter coupling in complex aperiodic environments.

## Vii Acknowledgments

M.R. gratefully acknowledges financial support by the ‘Stiftung der deutschen Wirtschaft’ in the framework of a scholarship. L.D.N. gratefully acknowledges the sponsorship of the Army Research Laboratory under Cooperative Agreement Number W911NF-12-2-0023. The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of the Army Research Laboratory or the U.S. Government. The U.S. Government is authorized to reproduce and distribute reprints for Government purposes notwithstanding any copyright notation herein.

## Appendix A Discrete resonators

The aim of this appendix is to justify viewing substructures embedded in a larger binary aperiodic lattice as local resonators. To this end, we investigate the behavior of the simplest case of such a structure, , in more detail. Its Hamiltonian is

(7) |

with the (unnormalized) eigenstates

(8) |

where , with . For high contrast , localizes on the central site. The idea now is to view at high contrast as a resonator, where the site effectively plays the role of a cavity, while the outer sites play the role of cavity walls. The resemblance to a resonator becomes clearer for a larger structure with more modes between the resonator walls, like the structure in Fig. 12. As one can see, all but two eigenstates extend nearly exclusively on the internal sites, and the wave-like character of these states is well recognizable. Two states exclusively localize on the outer two sites. The setup thus acts as an extended cavity consisting of sites, with two sites playing the role of the cavity walls. The smaller structure is of the same nature, albeit with a cavity of only a single site . Notationally, we will divide the actual cavity and the cavity walls of a resonator by a vertical line, writing e.g. . Similarly, we also view the “inverse” structure as a resonator with resonator modes of higher energy, assuming . Moreover, closely neighboring resonators of the form

(9) |

can be seen as coupled resonators. To indicate the composite character of such resonators, we omit the inner vertical lines, i.e., .

## Appendix B Symmetry argument for the absence of edge states

Here we explain the absence of edge states in Fig. 4 (c) using the concept of local symmetry. The underlying symmetry concept is very general and not limited to the Fibonacci chain, as we demonstrate in the last paragraph of this appendix. Let us denote an arbitrary sequence of ’s and ’s by , its reverse ordered counterpart by , and by a single site or . Then

(10) |

where denotes the eigenvalue spectrum and the tridiagonal Hamiltonian representing . In words, the eigenvalue spectrum of a resonator is completely contained in that of the reflection-symmetric resonator . For example, if and , then and .

To prove the above statement, we note that the Hamiltonian reads

(11) |

where . The matrices and connect the central site to and , respectively. Now, using the “equitable partition theorem” from Ref. Röntgen et al., 2018, we can transform by a similarity transform into a block-diagonal form

(12) |

The similarity transform conserves , and since is block-diagonal, we have

(13) |

which proves Eq. 10. Moreover, again using the equitable partition theorem, one can show that the eigenvalues of belong to eigenstates of with negative parity with respect to the central site .

Let us now apply the above statement to Fig. 4 (c). Here, for each resonator mode at the edge, there exists one resonator mode within the bulk possessing a similar energy:

(14) | ||||

(15) | ||||

(16) |

where denotes the energy of the resonator mode . In the limit of high contrast, where the resonators present in Eqs. 16, 15 and 14 are disconnected from the remainder of the system, the approximations become equalities, and the edge state eigenenergies are thus “absorbed” into the corresponding quasiband.

In a similar manner, the energetic near-equivalence of resonator modes

at high contrast as occurring in Fig. 10 (a) can be explained.

## Appendix C Perturbation theoretical treatment

In this section, we give an explanation for the fragmentation of eigenstates at high contrast in terms of a perturbation theoretical analysis. This will also show why the dominant entries of the eigenstates are in almost all cases obeying local symmetries. To this end, we write the Hamiltonian Eq. 2 as

(17) |

where solely contains the diagonal part of , i.e., isolated sites, while connects them, i.e., contains the off-diagonal elements of . By means of , an eigenstate of as well as its energy are expanded as

(18) | ||||

(19) |

Inserting Eqs. 19 and 18 into the Schrödinger equation yields the perturbation series which is assumed to converge and thus solved order by order in .

At zeroth order, the perturbation series reduces to the eigenvalue equation for the unperturbed . Since it is binary, the eigenstates of are highly degenerate and form two groups, satisfying

where denote the number of sites with potential , respectively. The so-called “correct” zeroth-order states which fulfill

(20) |

and which occur in Eq. 18 and thus also in the perturbation series are linear superpositions of the . In the following, we will always denote the two sets by and simple call the the zeroth-order states.

At the start of the perturbation theoretical treatment, the are known, but the are usually not, and the can not be directly be determined. However, it can be shownHirschfelder and Certain (1974) that already the knowledge of the is sufficient to obtain a series of particular solutions to the -th order perturbation equation, yielding the energy-corrections as a byproduct. Provided that the degeneracy of a given state is lifted at -th order, then the corresponding correct-zeroth order state can be obtained by diagonalizing a matrix which can be derived from the -th order perturbation equationHirschfelder and Certain (1974). Here, is the number of states which are degenerate with up to order . Then, at order , the state correction can be obtained. Note that for the problem at hand, all degeneracies are guaranteed to be lifted at a finite order, since the eigenvalues of tridiagonal matrices with strictly non-vanishing sub-and superdiagonals (such as the one here) are distinctParlett (1998) (i.e., non-degenerate). Though all degeneracies will eventually be lifted, the order at which this happens is in general different for different states. In many textbooks, all degeneracies are resolved already at first order, and the zeroth-order states are the ones that diagonalize the matrix in the corresponding degenerate subspace. This results in simple expressions for the higher-order corrections for both the states and the energy. For our binary , however, degeneracies are usually resolved only at very high order, and the process becomes complex. For Fibonacci chains, all degeneracies are resolved at -th order for generation , at -th order for , at -th order for but only in -th order for .

In the following, we will first show the feasibility of degenerate perturbation theory by means of the Fibonacci chain, showing that for high contrast already the zeroth-order states are sufficient to explain the fragmentation of states. Next, we will show the process of determining the zeroth-order states in the first three orders, allowing for an intuitive picture of the emergence of fragmentation and locally symmetric amplitudes. We have numerically observed convergence of the perturbation series if the contrast is larger than roughly , depending on the exact chain.

#### Application onto the Fibonacci chain

Figure 13 demonstrates the applicability of degenerate perturbation theory to a -th generation Fibonacci chain [the same as shown in Fig. 3 (a)] at a contrast . In subfigure (a), at each site the difference

is shown. Note that the differences are rather small, and in Fig. 13 (b), a detailed picture is given for the uppermost state . In Fig. 13 (c) and (d), the full state and are shown, respectively. As one can see, already the zeroth-order state matches the fragmentation behavior of the full state quite well, up to the two double resonator modes on the left half of the chain. In Fig. 13 (e), we include the first-order correction . As one can see, the resulting state is very close to the full state shown in Fig. 13 (c). Although we have here only shown the -th state (i.e., uppermost) state in detail, the behavior for all other states is similar. This shows that already the first-order state corrections yield very good results.

If one goes to even higher contrast, already the zeroth-order states are sufficient to get a full picture of the fragmentation of a given state. This is demonstrated in Fig. 14 for a comparatively very high contrast of . Subfigure (a) shows the difference at each site. The subfigures (b) and (c) show the complete state and the zeroth-order state , for which the main features (the resonator modes) are visible very well. Again, this behavior is the same for all other states, indicating that already the zeroth-order states give a good representation of the localization patterns occurring in the full state. Before we explicitly show the computations for the first three orders in degenerate perturbation theory, let us comment on the connection between the symmetry of the underlying potential sequence and that of the non-negligible amplitudes of a given eigenstate by means of Fig. 14 (c). As can be seen, the zeroth-order state is locally parity symmetric individually within the two domains . However, as a whole this state is not locally reflection symmetric w.r.t. an axis denoted by . As we will outline in the following, the reason for this is that the environment of the two domains is different, where environment includes not only next-neighboring sites but also the ones located further away (we will explain the notion of “further away” in more detail below). In Fig. 14 (d), we change the environment of the right domain such that it matches that of the first domain up to the first five neighbors. As a result, the zeroth-order state is now symmetric w.r.t. a reflection through the axis . In the following, we will investigate the connection between local symmetries of the underlying chain and that of the zeroth-order states in more detail.

#### Emergence of localization patterns and their locally symmetric character

We will now show the procedure of finding the zeroth-order states, as can be found e.g. in Refs. Hirschfelder and Certain, 1974; Silverstone, 1971. Since this procedure is quite technical, to help the reader we have visualized the process in a concise form in Fig. 15 for the easily traceable case of .

As stated above, can in general not be determined before its degeneracy is not completely lifted. At higher orders, the states degenerate at zeroth order may split into subsets which are degenerate up to first order, each of which can subsequently split into subsets of states , which are degenerate up to second order, and so on. The determination of the zeroth-order states can be done by means of recursively defined auxiliary statesHirschfelder and Certain (1974); Silverstone (1971)

(21) | ||||

(22) | ||||

(23) | ||||

(24) | ||||

(25) |

appearing on the left-hand side of the above equations, where denotes the -th element of the set and can run from to the number of elements within the set. The index is equal to if and equal to if . Each expansion coefficient is the -th component of the vector defined as

where the matrices are obtained by a recursive processHirschfelder and Certain (1974); Silverstone (1971). Explicitly, for the first three orders they are

where

with and

The above recursive process does the following: At the start, we have . These are then superposed according to Eq. 22, obtaining . Within each degenerate subspace , these are again superposed according to Eq. 23, obtaining . Again, within each degenerate subspace , these are superposed according to Eq. 24, obtaining , and so on. Provided that the degeneracy of a given state is solved at -th order, the degenerate subspace for this state at orders contains only one state, so that naturally andHirschfelder and Certain (1974) .

In the following, we will prove that the simultaneously localize on one or more maximally extended blocks of potentials of the same kind (MEBPS) [Statement 1] and determine on which such blocks a given state can simultaneously localize (Statement 2). Each MEBPS is the cavity of a resonator, thus giving reason for the localization of states on resonators. Statement 1 also shows that the are locally parity symmetric individually on each MEBPS, and Statements 3 and 4 further deal with longer-range symmetries of the zeroth-order states. Out of the many possible choices of (due to its high degeneracy), in the following we chose them such that [] is solely localized on the -th site with potential [] (counted from the left).

###### Statement 1.

Each state simultaneously localizes on one or more maximally extended blocks or of potentials of the same kind (MEBPS) and is locally parity symmetric on each of these blocks.

###### Proof.

The proof is done in three steps. Firstly, we show that is block-diagonal, where each block is related to exactly one MEBPS. Secondly, we show that the eigenvectors of are locally parity symmetric and subsequently the are locally symmetry on each MEBPS. Thirdly, we show that any higher-order states with show this local symmetry as well, and thus the zeroth-order states are locally symmetric as well.

We start by proving the following. For the case that contains only next-neighbor couplings (as is the case here) the become block-diagonal, i.e., can be written as

(26) |

where denote the number of blocks occurring in and each block

(27) |

is a tridiagonal Toeplitz matrix. To prove that is of the above form, we note that by the definition of and , two states are coupled to each other by any of the two matrices provided that (i) the single sites on which they localize are direct neighbors and (ii) they have the same zeroth-order energy, i.e., they must be localized on states with identical on-site potential. If (i) and (ii) are fulfilled for with and , then the corresponding matrix element due to the definition of these states. As a result, for each MEBPS [] containing sites, there is one tridiagonal block of the form Eq. 27 present in [].

We now show that the are locally parity symmetric on each MEBPS. To this end, we use the fact that the eigenvectors of the block-diagonal matrix are

with and where is the vector with identical zero entries and denotes the set of eigenvectors of . All vectors in a given set have non-vanishing components only on one maximally extended block of potentials of the same kind and is parity-symmetric w.r.t. a reflection through the center of this block. The latter is due to the fact that the are real and bisymmetric, and such matrices have definite parityCollar (1962) (in the case of degeneracies, the eigenvectors can be chosen accordingly). A matrix is bisymmetric if it is symmetric both around the main and the anti-diagonal. Since we have ordered the state such that it has non-vanishing amplitude on the -th site with potential , one can easily show that each of the has definite parity on each MEBPS.

For second order degenerate perturbation theory, the states which are degenerate up to first order are superposed to obtain . Now, since all states in a given set have distinct eigenvalues, the states are constructed such that for each set there is at most one state possessing non-vanishing amplitudes on any given MEBPS. Thus, will keep the local parity symmetry, and it is trivial to show that the zeroth-order states are locally parity symmetric on each MEBPS as well. ∎

Due to its maximal extension, each MEBPS is directly neighbored either by potentials of the same kind on one or on both sides, with the former being the case if the MEBPS forms one edge of the chain. Thus, the are seen to localize on resonators. We now show that a given state can only simultaneously localize on resonators fulfilling certain conditions.

###### Statement 2.

A given zeroth-order state can simultaneously localize on a set of MEBPS with individual lengths only if the following conditions are met. (i) All the MEBPS must have potentials of the same kind. (ii) There exist integers such that