Quasiparticle localization and ergodicity breaking in flat-band lattice models

Abstract

In lattice models, destructive interference can quench single-particle motion and produce (nearly) dispersionless Bloch bands, known as flat bands. Suppressing kinetic energy enhances the role of interactions and can lead to pairing, constrained dynamics, and localization-like behavior even in translationally invariant settings. This thesis studies attractive Hubbard-type interactions in such flat-band systems using an isolated-band (projected) description, which yields an effective Hamiltonian in a Wannier-orbital basis. The π-flux dice lattice is a convenient focus because its two lowest bands are perfectly flat and degenerate, and their subspace admits compactly localized Wannier functions. This makes the projection straightforward to implement and yields a projected Hamiltonian with a simple, short-range structure. A one-dimensional reference model is used alongside the dice lattice to separate and benchmark the effects of different processes. In Publication I, we construct a class of exact many-body eigenstates of the projected Hamiltonian of generic lattice models with flat bands in which localized unpaired quasiparticles coexist with a background of paired particles in the zero-momentum state. This general result motivates the study of a specific model, the π-flux dice-lattice. Using exact diagonalization and realtime evolution from few-particle initial states, we show that two-body bound states behave ergodically. At the same time, the density and spin distribution retain strong memory of the initial quasiparticle positions at long times. These results demonstrate localization-like dynamics in a translationally invariant flat-band model, where memory arises from projection-induced constraints rather than quenched disorder. In Publication II, we introduce a one-dimensional on-site/bond singlet model whose terms mirror those of the projected dice-lattice Hamiltonian, with tunable couplings controlling their relative strengths. Exact diagonalization shows that, for some couplings, the spectrum and dynamics of the on-site/bond-singlet model closely reproduce those of the projected dice lattice. We then use spectral statistics to probe how the different dynamical regimes manifest in energy-level correlations. In Publication III, we apply spectral statistics to the projected dice lattice using higher-order spacing and gap-ratio diagnostics. We find a robust parity dependence: even particle-number sectors show multi-block GOE behavior within fixed symmetry sectors, while odd particle-number sectors exhibit exact doublets and inter-doublet correlations that are GUE at even orders. The appearance of distinct random-matrix ensembles (GOE and GUE) within the same physical model is an unusual and interesting feature, which calls for a deeper understanding of the hidden symmetry structure of the Hamiltonians obtained by projecting onto the flat-band subspace in the dice lattice and other lattice models with flat bands realized with ultracold gases in optical lattices and in the newly discovered class of moiré materials.