Abstract
Adding interactions to many-body Hamiltonians of geometrically frustrated lattices often leads to diminished subspaces of localized states. In this paper, we show how to construct interacting many-body Hamiltonians, starting from the non-interacting tight-binding Hamiltonians, that preserve or even expand these subspaces. The methods presented involve modifications in the one-body network representation of the many-body Hamiltonians which generate new interacting terms in these Hamiltonians. The subspace of many-particle localized states can be preserved in the interacting Hamiltonian, by projecting the interacting terms onto the subspace of many-body extended states or by constructing the interacting Hamiltonian applying origami rules to the network. Expanded subspaces of localized states are found if interacting terms that mix subspaces with different number of particles are introduced. Furthermore, we present numerical methods for the determination of many-body localized states that allows one to address larger clusters and larger number of particles than those accessible by full diagonalization of the interacting Hamiltonian. These methods rely on the generalization of the concept of compact localized state in the network. Finally, we suggest a method to determine localized states that use a considerable fraction of the network.
Similar content being viewed by others
Introduction
Flat-band systems have been heavily studied in the past two decades1,2,3,4,5,6,7,8,9. The attention gathered by these systems is due to the interesting phenomenon they manifest, such as fractional quantum Hall effect10 or high-temperature superconductivity11, among others. Flat bands are associated to subspaces of one-body localized states, and such states are a consequence of destructive interference between different classical trajectories of particles and are associated with particular lattice geometry and particular choice of the hopping parameters of the respective tight-binding (TB) Hamiltonian (hopping terms between nearest neighbors in the simplest cases, but longer range hoppings can also allow the existence of localized states)12,13,14.
From the flat-band Bloch eigenstates, applying the inverse Fourier transform, and properly choosing the relative phase between the flat-band Bloch states, it is possible to write maximally localized Wannier states15,16. These Wannier states may be in certain lattices compact, but in general they will have small tails in the neighboring unit cells. A flat-band basis of compact localized (CL) states usually still exists in the latter case, if one drops the condition of an orthogonal basis. This is case of the Lieb lattice where the most compact one-body localized states are described by standing waves that travel along one plaquette, and the particle wavefunction amplitude is zero in the rest of the lattice12,13 [see Fig. 1(a)]. The previous example of the Lieb lattice allows one to conclude that interacting many-body localized states may still exist in the presence of interacting terms that require overlap of the flat-band states if the number of interacting particles is less than the number of plaquettes and the overlap can be avoided. These are states that are insensitive to the interacting terms. In simple cases (for example in the case of the Lieb lattice or diamond chain), one-body CL states are described by standing waves with normalized amplitudes (+1, 0, −1, 0, +1, ...) along a short path. However, in general, the flat-band subspace is perturbed by the introduction of interacting terms and it is not obvious if compact (many-body) localized states survive and how to determine them. In the case of interacting many-body localized states, the localized particles may be quite far apart in the lattice and the notions of standing wave and CL state need to be generalized.
Recently, we have suggested a different approach to tackle this problem17. It is possible to describe many-body localized states through the network generated from the adjacency matrix associated with the respective Hamiltonian. In this representation, a many-body interacting Hamiltonian can be interpreted as a one-body Hamiltonian with possibly long range hoppings and local potentials in the many-body Wannier base, where the single particle moves in the network defined by the finite Hamiltonian matrix elements between the many-body Wannier states (nodes of the network) [see17 for details]. In the network representation, localized many-body states may still be described by independent standing waves [see Fig. 1(b)], the majority of them with support in a small compact region of the network. These compact many-body states correspond to interlocked one-body standing waves in the network, which we designate as bubbles [see Fig. 1(c)]17. The dimension of these bubbles is related to the average number of links per node (degree of the node) in the network, which is determined by the number of particles.
In this paper, we show how to construct interacting many-body Hamiltonians that allow the existence of CL states and expand some of the ideas that were introduced in reference17. We present methods that involve modifications in the real lattice, in the network or directly in the Hamiltonian, such that the subspace of many-body localized states of the non-interacting Hamiltonian is diminished, expanded or preserved. In addition, since it becomes computational demanding to diagonalize many-body Hamiltonians and impractical to draw the corresponding networks (the matrices become too large even for a few number of particles), we suggest some methods that allow one to work with parts of a lattice or network, at a time, reducing the size of the matrices to diagonalize and the size of the networks to draw. We also present an algorithm that allows one to determine many-body localized states that are neither compact in the real lattice or in the network. We apply these methods to determine full sets of localized states of Lieb clusters and diamond chains.
Preserving CL states: projecting interactions out of the subspace of localized states
It is possible to add any kind of interacting terms to a Hamiltonian such that the subspace of states admits exactly the same many-body localized states (states such that all particles are localized) of the non-interacting Hamiltonian, as long as the interactions are projected out of the subspace of localized states (in other words, the interaction is projected onto the subspace of extended states). Since the subspace of localized states is orthogonal to the subspace of extended states, one has \({\hat{P}}_{ext}| loc\rangle =0\) and \({\hat{P}}_{loc}| ext\rangle =0\), where \({\hat{P}}_{ext}\) and \({\hat{P}}_{loc}\) are projector operators that act only in, respectively, extended and localized states. Note that states \(| loc\rangle \) have only localized particles and \(| ext\rangle \) has at least one particle in an extended state (but one or more particles may be localized). This is different from another possible approach where: (i) the single particle annihilation and creation site operators are written as a linear combination of Wannier states of the set of bands of the non-interacting model; (ii) a projection of the interacting Hamiltonian onto extended states is obtained discarding in the single particle annihilation and creation site operators, the part in the linear combination corresponding to the Wannier states of the flat band (see14,18 for an application of this method). The latter is a more severe projection method while the one we present here can be described as the projection closest to the identity operator that applied to the interacting network recovers the many-body localized states of the non-interacting network.
Adding terms to the Hamiltonian implies modifications in the network, so how can the network change and still be associated to a Hamiltonian with the same localized states? In order to answer this question, for simplicity, we discuss this method first in the case of the diamond chain, shown in Fig. 2(a). The diamond chain has three types of sites, A, B and C. The tight-binding Hamiltonian for the diamond chain with spinless fermions, assuming zero on-site potentials and periodic boundary conditions, is given by
where the index j labels the unit cells, Nc is the number of unit cells, and t is the hopping constant. The first term corresponds to intra-unit cell hoppings terms and the second term corresponds to inter-unit cell hopping terms. From the diagonalization of HTB, one obtains three electronic energy bands. Two of those bands are dispersive (non-flat) with energy \({\epsilon }_{k,\pm }=\pm 2t\sqrt{2}\cos (\sqrt{2}{k}_{x}/2).\) The third band is a zero energy (dispersionless) flat-band and it is Nc-fold degenerated. The one-particle CL states associated with the flat band arise from the anti-bonding combination of B and C orbitals, \(| loc;j\rangle =({B}_{j}^{\dagger }-{C}_{j}^{\dagger })/\sqrt{2}| vac\rangle .\) These states result from the linear combination of the doubly degenerate opposite momentum Block states k = π/2 and k = −π/2 of a diamond plaquette [see Fig. 2(b)]. The flat band has only support in the B and C sublattices, and admits Nc orthogonal localized states (the orthogonality of the localized states subspace basis makes the determination of the projection operators easier).
In Fig. 2(c), it is shown the network associated to the TB Hamiltonian of a diamond chain with four plaquettes, with one spin up and one spin down. Each node of the network corresponds to a position pair (i↑, j↓) where i↑ is the position of the spin-up and j↓ is the position of the spin-down according to the site indexing shown in Fig. 2(a). The Hamiltonian admits sixteen two-particle localized states where each particle is in a localized state (dropping the spin index) \(| loc;j\rangle \). These two-particle localized states in the network are linearly independent standing waves that form squared bubble structures (with twelve sites), with finite particle wave-function amplitude in the (red) nodes that correspond to having both spin up and spin down localized at sites B or C. Adding interacting terms to the previous Hamiltonian may lead to a loss of localized states. For instance, in the presence of Hubbard interactions, the states in which the spin-up and the spin-down were localized in the same plaquette become extended (note that the Hubbard interaction becomes a local potential in the one-body network of many-particle Wannier states).
Nonetheless, if the Hubbard terms (U∑ini↑ni↓, with \({n}_{i,\sigma }={X}_{i,\sigma }^{\dagger }{X}_{i,\sigma }\), X = A, B or C, and where the index i labels the sites of the lattice) are projected onto the subspace of extended states, one finds that the interacting Hamiltonian admits exactly the same number of localized states as the non-interacting version of this Hamiltonian (we choose U/t = 10, but one could choose any finite value for U). In this case, the network is slightly different from the network shown in Fig. 2(c), having some modifications in the off-diagonal bubble like structures. The new network includes links that connect the red nodes in Fig. 2(c) that correspond to having both spin occupying sites B or C of the same plaquette. These modifications are shown in Fig. 2(c) at the right, for the top right corner of the network (in the real lattice this correspond to having the spins occupying the top plaquette of the diamond chain shown in Fig. 2(a)). The new edges (drawn in green) correspond to new terms: long range hoppings, t′, [edges linking (2, 3) to (2, 2) to (3, 2) to (3, 3) to (2, 3)], a pair tunneling term, Δ, [edge from (2, 2) to (3, 3)] and a spin exchange term, J, [edge from (3, 2) to (2, 3)]. Besides these modifications, the new Hamiltonian includes long range spin interactions, V, between spins located in sites B and C of the same plaquette. Note also the lowering of the Hubbard interaction values for sites B and C, U′<U (sites 2 and 3 in this example). The fractions α/U (where \(\alpha =t{\rm{{\prime} }}=J=V\)), Δ/U and U′/U are equal to, respectively, 1/8, 3/8 and 5/8, and we found that these ratios remain the same for any value of U.
This means that it is possible to construct a Hamiltonian that combines different kinds of interactions in such a way that the effect of the Hubbard interacting terms is not felt by the subspace of CL states. Note that superconducting interactions and exchange interactions are often assumed to be competing interactions19,20, but in this case the different interactions are cooperating in order to generate the many-particle CL state subspace. The interacting Hamiltonian can be written, symbolically, as a combination of such interactions in the following way
The term (\(\hat{U}{\rm{{\prime} }}+\hat{V}+\hat{t}{\rm{{\prime} }}+\hat{\Delta }+\hat{{\boldsymbol{J}}}\)) does not affect the subspace of localized states and thus, the last two terms of Eq. 2 can be seen as perturbation terms in what concerns the subspace of localized states.
Let us now illustrate this method in more detail in the case of the Lieb lattice. The Lieb lattice is a square bipartite lattice characterized by a unit cell with three types of sites: A (at the vertices), B (at the vertical edges) and C (at the horizontal edges) [see Fig. 3(a)]. The tight-binding Hamiltonian for spinless fermions in the Lieb lattice, assuming zero on-site potentials, is given by
where periodic boundary conditions are implicit and where Lx (Ly) correspond to the number of unit cells in the x (y) direction, and t is again the hopping parameter. The first and second terms correspond, respectively, to intra-unit cell and inter-unit cell hopping terms. From the diagonalization of HTB, one obtains two dispersive (non-flat) bands and one (dispersionless) flat band. The energy of the dispersive bands is given by \({\epsilon }_{k,\pm }=\pm 2t\sqrt{{\cos }^{2}({k}_{x}/2)+{\cos }^{2}({k}_{y}/2)},\) where kα = 2πnα/Lα with nα = 0, 1, . . . , Lα − 1 and α ∈ x, y. The flat band is a zero energy band LxLy-fold degenerated, and the one-particle localized states associated with this band can be written as \(| loc;x,y\rangle =({B}_{x,y}^{\dagger }-{C}_{x,y}^{\dagger }+{B}_{x+1,y}^{\dagger }-{C}_{x,y+1}^{\dagger })/2| vac\rangle .\) The flat band has only support in the B and C sublattices, or, in other words, flat band localized states have zero amplitudes at A sites. The most compact one-particle localized states of Lieb lattices are localized in a single plaquette [see Fig. 1], and are described by standing waves that results from the linear combination of doubly degenerate opposite momentum Bloch states k = π/2 and k = −π/2 of a single plaquette [see Fig. 3(b)]. Two side-by-side CL states overlap at one site (site at the middle of the rung). This implies that the set of one-particle maximally CL states generates a non-orthogonal basis of the flat band subspace, but they are still linearly independent, and thus the number of one-hole localized states is equal to the number of plaquettes, Nplaq. A global many-body basis of the Lieb lattice can be written as \(\{| loc,j\rangle \}\cup \{| ext,i\rangle \}\) where \(\{| ext,i\rangle \}\) is the orthogonal basis of many-particle extended states subspace (again, recall that some of the particles in these states may be localized) and \(\{| loc,j\rangle \}\) is the non-orthogonal basis of compact many-particle localized states subspace.
Due to the fact that the CL states subspace basis, \(\{| loc,j\rangle \}\), is non-orthogonal, it follows that \({\hat{P}}_{loc}\ne {\sum }_{j}| loc,j\rangle \langle loc,j| \), and instead the projector operator is given by (see Supplementary Material and ref. 21 for more details)
where m runs over \(\{| loc,m\rangle \}\), and i runs over the global many-body basis of the Lieb lattice. In the case of a non-orthogonal basis set {i} the overlap matrix obeys the condition \({S}_{ij}=\langle i| \,j\rangle \ne {\delta }_{ij}\) for some i and j. If the basis \(\{| i\rangle \}\) is normalized, i.e., \(\langle i| i\rangle =1\) ∀ i, then S is symmetric (in our case S will also be real since we work with real Hamiltonians). The matrix S can be broken down into blocks in the following way
where the block \({\hat{{\boldsymbol{I}}}}_{ext}\) corresponds to the overlap between extended states, and, assuming an orthogonal basis for the subspace of extended states, it is just an identity matrix. The matrix elements of the anti-diagonal blocks are equal to zero due to the subspace of localized states being orthogonal to the subspace of extended states. The block \({\hat{S}}_{loc}\) refers to the overlap between localized states.
The four localized states of the TB Hamiltonian for the Lieb lattice with two plaquettes with one spin up and one spin down (with tij = t ∀ ij) are shown in Fig. 3(c). The overlap matrix in the basis \(\{| loc,j\rangle \}\cup \{| ext,i\rangle \}\) is in this case
This overlap matrix is easy to understand looking at the states in Fig. 3(c). The values displayed in the 4x4 block (\({ {\hat{S}} }_{loc}\)) were obtained doing the inner product between the four localized (non-orthogonal) eigenstates, of the non-interacting Hamiltonian. That is, this block is \({S}_{ij}^{loc}=\langle lo{c}_{i}| lo{c}_{j}\rangle \), with i, j = 1, 2, 3, 4. Note that the particular order chosen for the states in Fig. 3(c) is irrelevant. Since these states have 16 equal finite absolute amplitudes in the network (red nodes), for each shared red node between two states of Fig. 3(c), one has a 1/16 contribution to the overlap (all contributions have the same sign). Thus, the overlap is 1, 1/4 and 1/16 when the localized states share respectively 16, 4 and 1 nodes of finite amplitude.
In Fig. 3(d), it is shown the network associated to the interacting Hamiltonian,
with Hubbard terms projected onto the subspace of extended states, where one introduced the spin-1/2 variable σ (and the respective sum over the spin degrees of freedom in the Hamiltonian), X = A, B, or C, and \({n}_{i\sigma }={X}_{i\sigma }^{\dagger }{X}_{i\sigma }\) and the sum over nearest neighbors ⟨ij⟩ for the Lieb lattice is given in Eq. 3. The matrix representation of the projection operator \({\hat{P}}_{ext}\) in the many-body Wannier basis can be determined numerically faster from \({\hat{P}}_{ext}=\hat{I}-{\hat{P}}_{loc}\) (where the matrix representation of \(\hat{I}\) is the identity matrix and \({\hat{P}}_{loc}\) is given above) since in general the dimension of the subspace of localized states is much smaller than the dimension of the subspace of extended states. This network is very different from the network associated with the non-interacting Hamiltonian, having new edges connecting all (red) nodes corresponding to electronic configurations with the spin particles occupying sites B and C, including nodes that correspond to having spins trapped in different plaquettes. This means that jumps of particles (or pair of particles) between plaquettes are allowed, and this is a consequence of the non-orthogonal nature of the CL states subspace basis. Despite these modifications, the new network is associated with an interacting Hamiltonian that holds exactly the same localized states of the Hamiltonian that lead to the network in Fig. 3(c), and the reasoning for this is the same as in the case of the diamond chain. The terms that result from the projection of the Hubbard terms combine and negate the effect of the Hubbard interaction. In this case, two new kind of terms emerge from the projection of the interacting Hubbard terms: long range hopping terms that allow particles to jump between sites of different plaquettes, and an interaction that promotes tunneling of displaced pairs of spins to other B or C sites of the lattice.
Preserving CL states: Origami rules in the network
There is a complete set of modifications that can be applied to a geometrically frustrated lattice that preserve the single-particle CL states. This set of rules12 has been proposed for non-interacting geometrically frustrated lattices and clusters, but they remain valid for many-body interacting systems, and better, these rules can be applied to the network generated from the adjacency matrix associated with a many-body non-interacting Hamiltonian and the respective modifications of the network structure lead to modifications in the Hamiltonian, in the form of interactions. As an example, one of the Origami rules states that one may add or drop links between nodes of the network with zero many-body wavefunction amplitude and the resulting state is still a CL eigenstate of the modified many-body Hamiltonian. Thus, it is possible to add interacting terms to a non-interacting Hamiltonian that admits compact many-body localized eigenstates, in such a way that the localized states of the non-interacting Hamiltonian are preserved.
In Fig. 4(b), it is shown the single localized state, in the network representation, admitted by the TB Hamiltonian for a single diamond plaquette (Fig. 4(a)) with one spin up and one spin down particles, assuming that any additional plaquettes will connect to sites 2 and 4. As an example of application of the Origami Rules, one may introduce a term \({\sum }_{ijmn}{c}_{i,\uparrow }^{\dagger }{c}_{j,\downarrow }^{\dagger }{c}_{m,\uparrow }{c}_{n,\downarrow }\) in the TB Hamiltonian, where the indices run over all sites of the network region shown in Fig. 4(b). This term correspond to the first Origami rule that says that an additional set of hopping bonds may be introduced if they do not lower the local symmetry of the network region where the localized state has finite amplitudes and in this case the CL state is preserved. In the network, the term \({c}_{i,\uparrow }^{\dagger }{c}_{j,\downarrow }^{\dagger }{c}_{m,\uparrow }{c}_{n,\downarrow }\) becomes “hopping” terms between two nodes of the network. This term links all sites equally, keeping the symmetry of the network [see Fig. 4(c)] and it is equivalent to a \({c}_{k=0,\uparrow }^{\dagger }{c}_{k=0,\downarrow }^{\dagger }{c}_{k=0,\uparrow }{c}_{k=0,\downarrow }\) term (which is a density operator product in the k-space of one plaquette). The CL state of the plaquette remains the same since it is a linear combination of k = π/2 and k = −π/2 Bloch states in the diamond plaquette. We opted to chose k = 0 but other k values are possible as long as they are different from k = π/2 and k = −π/2 (this is also valid for diamond chains with more plaquettes, and since all localized states can be written as a linear combination of CL states, one may use any linear combination of k values of different plaquettes too).
Diminishing the number of CL states: adding interactions
It it possible to add interacting terms to an Hamiltonian such that the new Hamiltonian still preserves some of the localized states of the original one, lifting partially the degeneracy of the subspace of localized eigenstates. For instance, as discussed previously in this paper, adding a Hubbard term to a TB Hamiltonian of a Lieb lattice (or diamond chain) with one spin up and one spin down particles leads to a diminished subspace of localized states. Note that, if the number of repulsive spins particles is larger than the number of plaquettes then no localized states survives.
However note that adding interactions may also have the opposite effect, i.e, it is possible to add interactions to a Hamiltonian that, originally, does not admit localized states and find localized states, or expand the subspace of localized states of a Hamiltonian that admits localized states in the first place. For the latter case, as example, we show in the next section that the subspace of localized states of the TB Hamiltonian of diamond chains with spin up particles can be expanded by adding terms to the Hamiltonian that mix states that involve components with different spin up particle number.
Generating new CL states: mixed filling states
One can generalize the argument given in section 3 to other types of interactions that do not involve the linear combination of k = π/2 and k = −π/2 plaquette Bloch states corresponding to the CL state of the diamond chain. For example, one can construct Hamiltonians that admit localized states that involve components with different particle number (with links connecting the nodes of the Hamiltonian networks of different fillings). We illustrate this in the case of diamond plaquettes with different number of spin-up particles. In Fig. 5(b), the network associated to the tight-binding Hamiltonian of a single diamond plaquette (Fig. 5(a)) for all possible fillings of spin-up particles is shown. The finite amplitudes of the single CL state in the network admitted by this Hamiltonian are shown in red. This state can be described as a one-body standing wave in the network, with finite particle amplitude in the nodes that correspond to having a single particle occupying the sites 1 and 3 (this agrees with the fact that, due to Pauli exclusion, localized states with more than one spin up particles in a single plaquette can not exist). Again we assume that additional plaquettes will be connected to sites 2 and 4.
In order to construct many-body CL states that involve different particle number we added the interacting pair tunneling terms c1c2 + c1c4 + c2c3 + c4c3 + H. c. to the non-interacting Hamiltonian (recall that the CL state in the diamond plaquette is the anti-bonding combination of the 1 and 3 Wannier states). In Fig. 5(c), the network associated with the new Hamiltonian is shown, as well as the amplitudes of the two CL states that are now found. In the left part of the network Wannier states of one and three particles are linked, while in the right part of the network Wannier states of four, two and zero particles are connected. In the localized state shown at the top, the many-body particle wavefunction amplitude is finite for the single-particle Wannier states 1 and 3 (just as in the case of the localized state of the original, non-interacting Hamiltonian). In the bottom localized state of Fig. 5(c), the many-body particle wavefunction amplitude is finite for the node (1, 3) which corresponds to having two spin-up particles occupying the 1 and 3 sites of the real lattice and for the node corresponding to the vacuum state. Longer range Hubbard-like type interactions may also lead to new localized states in a spinless diamond chain (such localized states have been previously described in22 and are due to the interplay of geometric frustration and interacting terms).
This method can be generalized for lattices with any number of plaquettes (and in particular, for any geometrically frustrated lattice). In Fig. 6(b), the network associated with the non-interacting Hamiltonian in the case of a diamond chain with two plaquettes [see Fig. 6(a)] is shown, as well as the amplitudes of the three CL states. There are two one-body CL states, each corresponding to a standing wave traveling along one of the plaquettes, and a two-particle CL state with each spin-up particle occupying a different plaquette. In Fig. 6(c), the network for the same lattice is shown but now for the interacting version of the TB Hamiltonian where interacting terms c1c2 + c1c4 + c2c3 + c4c3 + H. c. were introduced (these terms act only in the left plaquette and we dropped the spin index). The interacting Hamiltonian admits now four localized states and they involve components with different number of particles. The first localized state occurs in the network with Wannier states of one and three particles (recall that the interaction acts only in the left plaquette). A single spin-up is localized in the left plaquette, having finite amplitude in the nodes 1 and 3. The second localized state occurs in the network that links Wannier states of two and four particles and the two-particle localized state has finite amplitudes in nodes that correspond to having the spins localized at different plaquettes. The third localized state involves a network of Wannier states of zero, two and four particles, with finite amplitude in nodes that correspond to having both spins in the left plaquette. The last localized state appears in the network of Wannier states of one and three particles, with finite amplitude in the nodes 5 and 7 which correspond to having a single spin located in the right plaquette and in the nodes (1, 3, 5) and (1, 3, 7) which correspond to having two spins in the left plaquette and one in the right plaquette.
We emphasize that the above mixed-filling compact many-body states are found only if the interacting term acts on a single plaquette. If the interacting term acts in all lattice sites, this leads to a reduction of the number of localized states because the interaction leads to weight transfer between plaquettes, generating an extended version of some CL states.
Finding many-particle CL states
In the previous sections, we have described how to construct interacting Hamiltonians that admit many-particle CL states. The form of these states was known since it was the same as that of the non-interacting Hamiltonian or because the Hamiltonian was defined in such a way that a particular CL eigenstate was present. Below (with the exception of subsection 6.5), we discuss how to find many-particle CL states of an interacting Hamiltonian without any prior knowledge about its possible form or existence.
Selecting parts of the network
For large number of particles in a lattice, the diagonalization of the interacting Hamiltonian in order to obtain the full set of localized states becomes numerically demanding. However, since most many-body localized states have support in a small region of the network, it is possible to determine the majority of the localized states by diagonalizing the Hamiltonian submatrix associated with that small region of the network17.
In Fig. 7, at the bottom, it is shown one cycle of our method to draw parts of the network and identify CL states. The diagrams 1–5 in Fig. 7 sketch the growth of a small part of the network associated with the Lieb-1 × 2 cluster with one hole and one flipped spin in a ferromagnetic background. Each node of the network corresponds to a position pair (iº, j↓) where iº is the position of the hole and j↓ is the position of the flipped spin according to the site indexing shown in Fig. 7. The method works in the following way. First, one arbitrarily chooses one node of the network (a state of the many-body Wannier basis) from which one starts to construct the network, by adding nodes with links to that initial node according to the elements of the adjacency matrix (see diagram 1). Next, one adds nodes with links to the current cluster of nodes and repeat this process (diagrams 2,3,4) until a bubble-like structure forms in the network (diagram 5), which may indicate the existence of a localized state (in order to verify if the bubble is associated to a localized state, one should determine the eigenstates of the Hamiltonian for this small region, and check if any of the states corresponds to a standing wave with zero amplitudes at the nodes that connect the bubble to the rest of the network). At the end of each step, one checks if any bubble-like structure emerge in the network and if it does, one diagonalizes the partial Hamiltonian. After such structure is found, one selects (randomly) another node, not present in the cluster grown so far, and starts to grow another part of the network, repeating the process until all bubble-like structures are grown.
Besides allowing to obtain localized states without having to diagonalize the full Hamiltonian, working with smaller adjacency matrices also always gives CL states while the diagonalization of the full network can generate “extended” versions of the localized states. Note that some many-body localized states correspond to standing waves that use large portions of the network, and thus can not be obtained by using this method. Further below in this section we discuss how to determine those states.
Drop sites in the lattice
Another way to work with small matrices is to drop sites of the real lattice. Again, one-body (and some many-body) CL states are described by standing waves that travel only along a small part of the whole lattice, and so it is possible to drop one or a set of sites and still determine the localized states that have support in the remaining lattice. For instance, in the case of one-hole localized states of a Lieb ladder (or diamond chain), one can remove any number of outer-most plaquettes and still find CL states of the Hamiltonian associated with the original lattice. If one determines the subset of localized states for each set of dropped points, it is possible to obtain the full set of one-body localized states. This method can be also applied to the many-body network, but it is trickier to choose which nodes to drop.
In Fig. 8(a), it is shown a cropped version of the Lieb-1 × 3 ladder, with the site in the middle of the rung between the first two plaquettes being dropped. The network associated to the Hamiltonian for that lattice with two holes and one flipped spin (in a ferromagnetic background) is shown in Fig. 8(b), where the set of amplitudes given corresponds to the single localized state (out of eighteen of the original Hamiltonian). This state uses substantial fraction of the network. The form of this state becomes clearer when one drops the nodes of the network with zero particle wavefunction amplitude that are not connected to red nodes and finds that the cropped network obtained this way [see Fig. 8(c)] is not much different from the whole network.
Lanczos algorithm
Some many-body compact localized states require a considerable fraction of the network and can not be obtained by using the methods described in the Section 6, and full diagonalization of the corresponding Hamiltonians is too time consuming. In order to determine such states we suggest the application of a Lanczos algorithm23 after using the previous methods. This algorithm is very efficient when applied to sparse Hamiltonian matrices, and in the cases presented in this paper, they are definitively sparse since we consider systems with short-range tight-binding terms.
Next, we explain how to apply this method to determine hole localized states of the TB Hamiltonian for the Lieb lattice (or diamond chain), in the infinite coupling limit (U∕t = ∞), that use a considerable part of the corresponding networks. In this case, the localized states are not the states with the lowest energy, but since they are zero energy state, they are the lowest energy states of the squared Hamiltonian. Note however that in order to lift possible degeneracies associated with extended states one has first to apply local potentials (for example, to sites A in the case of the Lieb lattice and diamond chain). Thus, the energy of the localized states is no longer equal to zero, and in order to shift that energy to zero one has to add a potential (diagonal) term to the Hamiltonian. After the previous steps, a Lanczos method should be applied to the Hamiltonian matrix squared in order to obtain its ground state. Since we are interested to determine the localized states that use a substantial fraction of the network (CL states that have support in small network regions can be obtained by the methods previously described), one needs to add an energy cost to any finite projection of our state onto the already determined CL states (to increase their energy). However, if there are more than one localized state that uses a considerable fraction of the network, the Lanczos algorithm will only select one of them. In this case, one may use the Jacobi-Davidson algorithm24 which is a variational method good at dealing with degeneracies, and is still very efficient when applied to sparse Hamiltonian matrices.
Symmetries and local spin rotation
If a given Hamiltonian with many-particle localized states has symmetries, then, as usual, the full set of CL states may be obtained from a subset of CL states, applying these symmetries. The many-hole Hubbard Hamiltonian of the Lieb lattice in the infinite coupling limit (U∕t = ∞) admits localized states such that the holes are localized in small regions of the lattice, making possible to treat the rest of the lattice as an independent system, and since this rest is half-filled, it has a local SU(2) spin symmetry25,26,27. On the other hand, the compact regions where holes localize have the usual global SU(2) spin symmetry of the Hubbard model. This implies that given a CL state, many others can be found applying any combination of the local spin-raising and spin-lowering operators at any site of the rest of the lattice and of spin-raising and spin-lowering operator defined within the region where the CL particle has support (a Lieb plaquette in the case of the Lieb lattice). So, as mentioned in the Introduction, it is possible to determine part of localized states subspace relying on symmetry arguments.
Note that one can add interactions to the Hubbard Hamiltonian that change its symmetries and one may still have localized states. For instance, one may add a nearest-neighbor interacting term to the infinite coupling Hubbard model on the Lieb lattice and compact one-hole localized states will still exist. This results from the fact that a hole CL state in the Lieb lattice does not involve components with different number of nearest-neighbor pairs of particles. Also one may break the SU(2) symmetry adding a Zeeman term (which only cares about the number of particles with up and down spin) and again compact one-hole localized states will still exist.
Equal amplitudes
In the case of small systems, one may still use the full diagonalization of the Hamiltonian in order to determine subspaces of localized states, since it may have a computing time of the order of the previous methods. However, the usual process of diagonalizing the Hamiltonian to obtain eigenstates may not give the most CL states, giving instead eigenstates that are a linear combination of the most CL states. In the Introduction, we mention that when Fourier transforming the flat band Bloch states in order to obtain maximally localized Wannier states, one may generate the most compact form of the localized states. The relative phases of the Bloch states have to be adjusted with that goal in mind and furthermore this method does not work if the most CL states are not orthogonal. Here, we present an easier procedure that allows one to obtain CL states without having to first determine the maximally localized Wannier states of real flat-band Hamiltonians. The method works in the following way. First, one finds the eigenvalues and eigenstates as usual (direct diagonalization of the Hamiltonian), and identifies the localized states, ψloc, and the respective energies, εloc. Then, one uses the εloc to determine again localized states, using the eigenvalue equation
but assuming that the finite amplitudes of the localized states, \({\psi }_{loc}^{{\prime} }\), are all equal in absolute value. This will avoid the linear combination of degenerate CL states that have overlap. The result of Eq. 8 is a set of CL states linearly independent from each other (or a linear combination of CL states with a clear separation of the finite amplitude regions). This method should be applied after using the previous methods.
Conclusion
The physics of flat bands is currently an hot topic in condensed matter studies. The experimental realization of flat bands has been achieved in, artificial, optical28,29,30,31, photonic32,33,34,35,36 and metallic37,38 geometrically frustrated lattices. Detecting flat bands becomes more challenging in the case of real materials39,40,41, but recently flat bands have been observed by angle-resolved photoemission spectroscopy experiments in bilayer graphene42, and in the case of twisted bilayer graphene, under electrostatic doping, zero-resistance states (associated to a flat band), were found43.
Flat bands are usually associated to subspaces of one-body localized states. In this paper we addressed many-body CL states. We presented methods of construction of interacting Hamiltonians such that the subspaces of CL states of the respective non-interacting TB Hamiltonians are preserved, diminished or expanded. We showed that localized states survive the addition of interactions if those interactions are projected out of the subspace of localized states. The terms that result from this projection combine and negate the effect of the original interacting terms in the subspace of localized states. Other modifications to the network may lead to diminished or expanded subspaces of localized states. We showed that the latter can be achieved by constructing Hamiltonians that mix states with different fillings (by adding links that connect networks associated to Hamiltonians with different particle number). Numerical methods for the determination of many-body localized states were discussed. Since most of these states are compact in the many-body network, the full diagonalization of the Hamiltonian is not necessary and larger cluster and larger number of particles can be addressed.
References
Huber, S. D. & Altman, E. Bose condensation in flat bands. Phys. Rev. B 82, 184502 (2010).
Imada, M. & Kohno, M. Superconductivity from flat dispersion designed in doped Mott insulators. Phys. Rev. Lett. 84, 143 (2000).
Neupert, T., Santos, L., Chamon, C. & Mudry, C. Fractional quantum Hall states at zero magnetic field. Phys. Rev. Lett. 106, 236804 (2011).
Peotta, S. & Törmä, P. Superfluidity in topologically nontrivial flat bands. Nat. Commun. 6, 8944 (2015).
Rhim, J.-W. & Yang, B.-J. Classification of flat bands according to the band-crossing singularity of Bloch wave functions. Phys. Rev. B 99, 045107 (2019).
Sun, K., Gu, Z., Katsura, H. & Sarma, S. D. Nearly flatbands with nontrivial topology. Phys. Rev. Lett. 106, 236803 (2011).
Tang, E., Mei, J.-W. & Wen, X.-G. High-temperature fractional quantum Hall states. Phys. Rev. Lett. 106, 236802 (2011).
Urban, D. F., Bercioux, D., Wimmer, M. & Häusler, W. Barrier transmission of Dirac-like pseudospin-one particles. Phys. Rev. B 84, 115136 (2011).
Wang, Y.-F., Gu, Z.-C., Gong, C.-D. & Sheng, D. Fractional quantum Hall effect of hard-core bosons in topological flat bands. Phys. Rev. Lett. 107, 146803 (2011).
Tsui, D. C., Stormer, H. L. & Gossard, A. C. Two-dimensional magnetotransport in the extreme quantum limit. Phys. Rev. Lett. 48, 1559 (1982).
Heikkilä, T. T. & Volovik, G. E. Flat bands as a route to high-temperature superconductivity in graphite.In Basic Physics of Functionalized Graphite, 123–143 (Springer, 2016).
Dias, R. G. & Gouveia, J. D. Origami rules for the construction of localized eigenstates of the Hubbard model in decorated lattices. Sci. Rep. 5, 16852 (2015).
Derzhko, O., Richter, J. & Maksymenko, M. Strongly correlated flat-band systems: The route from Heisenberg spins to Hubbard electrons. Int. J. Mod. Phys. B 29, 1530007 (2015).
Tovmasyan, M., Peotta, S., Törmä, P. & Huber, S. D. Effective theory and emergent SU(2) symmetry in the flat bands of attractive Hubbard models. Phys. Rev. B 94, 245149 (2016).
Marzari, N., Mostofi, A. A., Yates, J. R., Souza, I. & Vanderbilt, D. Maximally localized Wannier functions: Theory and applications. Rev. Mod. Phys. 84, 1419 (2012).
Marzari, N. & Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands. Phys. Rev. B 56, 12847 (1997).
Santos, F. D. R. & Dias, R. G. Hole-localized states in interacting geometrically frustrated systems. Phys. Rev. B 99, 125152 (2019).
Tovmasyan, M., van Nieuwenburg, E. P. L. & Huber, S. D. Geometry-induced pair condensation. Phys. Rev. B 88, 220510 (2013).
Berk, N. & Schrieffer, J. Effect of ferromagnetic spin correlations on superconductivity. Phys. Rev. Lett. 17, 433 (1966).
Rusinov, A. Ferromagnetism in superconducting alloys (1964).
Soriano, M. & Palacios, J. J. Theory of projections with nonorthogonal basis sets: partitioning techniques and effective Hamiltonians. Phys. Rev. B 90, 075128 (2014).
Lopes, A. A. & Dias, R. G. Interacting spinless fermions in a diamond chain. Phys. Rev. B 84, 085124 (2011).
Lanczos, C. An Iteration Method For The Solution Of The Eigenvalue Problem Of Linear Differential And Integral Operators (United States Governm. Press Office Los Angeles, CA, 1950).
Sleijpen, G. L. & Van der Vorst, H. A. A Jacobi-Davidson iteration method for linear eigenvalue problems. SIAM Rev. 42, 267–293 (2000).
Essler, F. H., Korepin, V. E. & Schoutens, K. Complete solution of the one-dimensional hubbard model. Phys. Rev. Lett. 67, 3848 (1991).
Östlund, S. & Mele, E. Local canonical transformations of fermions. Phys. Rev. B 44, 12413 (1991).
Carmelo, J. M. P., Östlund, S. & Sampaio, M. J. Global SO(3) × SO(3) × U(1) symmetry of the Hubbard model on bipartite lattices. Ann. Phys. (N. Y.) 325, 1550–1565 (2010).
Bercioux, D., Urban, D., Grabert, H. & Häusler, W. Massless Dirac-Weyl fermions in a T 3 optical lattice. Phys. Rev. A 80, 063603 (2009).
Goldman, N., Urban, D. & Bercioux, D. Topological phases for fermionic cold atoms on the Lieb lattice. Phys. Rev. A 83, 063601 (2011).
Wu, C., Bergman, D., Balents, L. & Sarma, S. D. Flat bands and Wigner crystallization in the honeycomb optical lattice. Phys. Rev. Lett. 99, 070401 (2007).
Wu, C. & Sarma, S. D. p x,y -orbital counterpart of graphene: Cold atoms in the honeycomb optical lattice. Phys. Rev. B 77, 235107 (2008).
Endo, S., Oka, T. & Aoki, H. Tight-binding photonic bands in metallophotonic waveguide networks and flat bands in kagome lattices. Phys. Rev. B 81, 113104 (2010).
Mukherjee, S. et al. Observation of a localized flat-band state in a photonic Lieb lattice. Phys. Rev. Lett. 114, 245504 (2015).
Vicencio, R. A. et al. Observation of localized states in Lieb photonic lattices. Phys. Rev. Lett. 114, 245503 (2015).
Xia, S. et al. Demonstration of flat-band image transmission in optically induced Lieb photonic lattices. Opt. Lett. 41, 1435–1438 (2016).
Zong, Y. et al. Observation of localized flat-band states in Kagome photonic lattices. Opt. Express 24, 8877–8885 (2016).
Kajiwara, S., Urade, Y., Nakata, Y., Nakanishi, T. & Kitano, M. Observation of a nonradiative flat band for spoof surface plasmons in a metallic Lieb lattice. Phys. Rev. B 93, 075126 (2016).
Nakata, Y., Okada, T., Nakanishi, T. & Kitano, M. Observation of flat band for terahertz spoof plasmons in a metallic kagome lattice. Phys. Rev. B 85, 205128 (2012).
Chen, H. et al. Prediction of two-dimensional nodal-line semimetals in a carbon nitride covalent network. J. Mater. Chem. A 6, 11252–11259 (2018).
Jiang, W., Liu, Z., Mei, J.-W., Cui, B. & Liu, F. Dichotomy between frustrated local spins and conjugated electrons in a two-dimensional metal-organic framework. Nanoscale 11, 955–961 (2019).
Su, N., Jiang, W., Wang, Z. & Liu, F. Prediction of large gap flat Chern band in a two-dimensional metal-organic framework. Appl. Phys. Lett. 112, 033301 (2018).
Marchenko, D. et al. Extremely flat band in bilayer graphene. Sci. Adv. 4, eaau0059 (2018).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43 (2018).
Acknowledgements
This work is funded by the FEDER funds through the COMPETE 2020 Programme and National Funds throught FCT - Portuguese Foundation for Science and Technology under the projects UIDB/50025/2020 & UIDP/50025/2020 and under the project PTDC/FIS-MAC/29291/2017. F. D. R. Santos acknowledges the financial support from the FCT through the grant PD/BD/135009/2017. R. G. Dias appreciates the support by the Beijing CSRC.
Author information
Authors and Affiliations
Contributions
F.D.R. Santos and R.G. Dias contributed to every stage of the work, including preparing the figures and writing and reviewing the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
Supplementary material
In this appendix, we discuss how to define projector operators for a generic non-orthogonal basis, following reference21. Let us consider an Hilbert space of dimension N and a generic normalized but non-orthogonal basis set \(\{| i\rangle \}\). The overlap beween any two states of this basis defines the overlap matrix \(\langle i| \,j\rangle ={S}_{ij}\). If this basis is the eigenbasis of a real Hamiltonian, then it can be chosen real and that means that S is a real symmetric matrix.
A generic one-body operator \(\hat{{\rm{A}}}\) can be written in terms of this basis as
where \({\tilde{\rm{A}}}\) is the “nuclear” matrix of the operator21. In the particular case of the identity operator we may write
where X is an unknown matrix. Since
and
this implies the matrix equation S ⋅ X ⋅ S = S and consequently X = S−1, that is,
(for simplicity we use the notation \({S}_{ij}^{-1}\) for the elements of the inverse matrix \({({{\rm{S}}}^{-1})}_{ij}\)).
The matrix elements of a generic one-body operator \(\hat{{\rm{A}}}\) are easily obtained
and therefore the “nuclear” matrix of the operator is \({\tilde{\rm{A}}}={S}^{-1}A{S}^{-1}\).
For a general basis set with a overlap matrix S there exist a basis set dual to this (direct) with the property \(\langle i| \,{j}^{\ast }\rangle ={\delta }_{ij}\). Only for orthogonal basis sets both direct and dual sets coincide. The dual representation of a generic operator \(\hat{{\rm{A}}}\) does not give the matrix elements, instead it gives the matrix elements of the nuclear matrix.
In order to obtain the matrices elements Aij in the direct representation, \(\hat{{\rm{A}}}\) can be written in term of the dual basis sets
The definition for \(\hat{{\rm{A}}}\) written in the direct basis in Eq. 9 can be recovered from Eq. 16, by using the expressions
and
Furthermore, since \(\langle i| \hat{{\rm{I}}}| j\rangle ={S}_{ij}\) the identity operator can be written as
Let us now describe a case directly related to the topic addressed in this paper. Assume that in the non-orthogonal basis \(\{| i\rangle \}\), only a few states \(\{| \neg \perp ,m\rangle \}\) have non-zero overlap and we wish to define a projection operator into the subspace defined by these states. This means that \(\{| i\rangle \}=\{| \perp ,n\rangle \}\cup \{| \neg \perp ,m\rangle \}\) and a general state can be written as
We will start by defining the operator \({\widehat{{\rm{P}}}}_{i}\) such that
It is easy to show that
Note that this operator is not a projection operator according to the usual definition (projection operators are unitary and Hermitian). Nonetheless, the usual Hermitian projection operator \(| i\rangle \langle i| \) can be obtained by using both \({\hat{{\rm{P}}}}_{i}\) and its adjoint \({\hat{{\rm{P}}}}_{i}^{\dagger }\) in the following way \({\hat{{\rm{P}}}}_{i}{\hat{{\rm{P}}}}_{i}^{\dagger }=| i\rangle \langle i| .\)
Projecting onto the subspace of basis \(\{| \neg \perp ,m\rangle \}\) is possible through the use of the generalized operators
and
that allow a projection of bras and kets onto the non-orthogonal subspace
The restriction of a operator to the non-orthogonal subspace is therefore obtained using both projectors in the following way
Note that
therefore
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Santos, F.D.R., Dias, R.G. Methods for the construction of interacting many-body Hamiltonians with compact localized states in geometrically frustrated clusters. Sci Rep 10, 4532 (2020). https://doi.org/10.1038/s41598-020-60975-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-60975-7
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.