Local Hilbert space fragmentation and weak thermalization
in Bose-Hubbard diamond necklaces
Eloi Nicolau,1Anselmo M. Marques,2Jordi Mompart,1Ver`onica Ahufinger,1and Ricardo G. Dias2
1Departament de F´ısica, Universitat Aut`onoma de Barcelona, E-08193 Bellaterra, Spain.
2Department of Physics and i3N, University of Aveiro, 3810-193 Aveiro, Portugal.
We study Bose-Hubbard models in a family of diamond necklace lattices with ncentral sites. The
single-particle spectrum of these models presents compact localized states (CLSs) that occupy the up
and down sites of each diamond. By performing an appropriate basis rotation, the fragmentation of
the many-boson Hilbert space becomes apparent in the adjacency graph of the Hamiltonian, showing
disconnected sub-sectors with a wide range of dimensions. The models present a conserved quantity
related to the occupation of the single-particle CLSs that uniquely identifies the different sub-sectors
of the many-boson Hilbert space. Due to the fragmentation of the Hilbert space, the distribution of
entanglement entropies of the system presents a nested-dome structure. We find weak thermalization
through sub-sector-restricted entanglement evolution and a wide range of entanglement entropy
scalings from area-law to logarithmic growth. Additionally, we observe how the distinguishability
between the different domes increases with the number of central sites and we explain the mechanism
behind this fact by analyzing the graph structure of the Hamiltonian.
I. INTRODUCTION
The Eigenstate Thermalization Hypothesis (ETH) pre-
dicts how an excited state of a many-body closed quan-
tum system should thermalize [1–3]. Although most sys-
tems obey this hypothesis, numerous examples of non-
ergodic systems have been found. Perhaps the most
prominent example is integrable systems, where the num-
ber of conserved quantities equals or exceeds the degrees
of freedom of the system, thus exactly determining all
the eigenstates [4]. In many-body localized systems [5],
the interplay between disorder and interactions gives rise
to emergent integrability, which also leads to a strong
violation of the ETH. More recently, it was shown that
the ETH can also be weakly violated by a vanishing sub-
set of non-thermal eigenstates, dubbed Quantum Many
Body Scars (QMBS). They were initially found in one-
dimensional Rydberg arrays [6] with the underlying PXP
model [7,8], and were also discovered in parallel in the
AKLT model [9,10]. Since these initial works, QMBS
have been found in several systems where there is either
a tower of scarred eigenstates [10–23] or an isolated scar
[24–32].
A broader phenomenon that also leads to weak ther-
malization is Hilbert space fragmentation, also known as
Hilbert space shattering or Krylov fracture [33]. The
Hilbert space presents exponentially many dynamically
disconnnected sectors that prevent the system from ther-
malizing completely. Remarkably, this mechanism can
lead both to a weak or a strong violation of the ETH.
This effect can arise in a wide variety of systems, such
as dipole moment or center-of-mass conserving systems
[34–38], the 1D t-Jzmodel [39], the t-Vand t-V1-V2
models [40,41], and models with dipolar interactions
[42]. All the above examples exhibit fragmentation of
the Hilbert space in the product state basis [43], i.e.,
classical fragmentation. Quantum fragmentation, which
occurs in an entangled basis, has been recently shown to
arise in Temperley-Lieb spin chains [43] and in quantum
East models [44]. However, it has yet to be determined if
quantum fragmentation leads to different phenomenology
than its classical analogue.
The fragmentation in the above examples has recently
been referred to as standard Hilbert space fragmenta-
tion, to distinguish it from local Hilbert space fragmen-
tation [45], that arises in models with [21,27,32,46,47]
or without [48] frustration and in flat band models [49].
While standard fragmentation is due to the presence of
non-local conserved quantities, locally fragmented sys-
tems present strictly local conservation laws.
In this work, we report on a family of Bose-Hubbard di-
amond necklaces [50] that exhibit quantum local Hilbert
space fragmentation. Here, the presence of a single-
particle flat band composed of compact localized states
(CLSs) gives rise to the fragmentation of the Hilbert
space when introducing on-site interactions. As a conse-
quence of this fragmentation, one finds a nested distribu-
tion of entanglement entropies, sector-restricted thermal-
ization, and a broad range of sub-sectors of the Hamil-
tonian that range from frozen sub-sectors following area-
law to non-integrable sub-sectors with logarithmic scal-
ing.
The article is structured as follows: in Section II, we
introduce the system and we describe the basis rotation
that reveals the fragmentation of the Hilbert space in Sec.
II A. In Sec. II B, we analyze the conserved quantity that
characterizes the sub-sectors of the Hamiltonian, discuss
the adjacency graphs of the fragmented Hamiltonian, and
demonstrate that the system is strongly fragmented. The
numerical results are discussed in Sec. III, which include
the distribution of entanglement entropies, the entangle-
ment evolution and scaling, the level spacing analysis and
a comparison between the different models of the dia-
mond necklace family. Finally, we summarize our con-
clusions in Sec. IV.
arXiv:2210.02429v2 [cond-mat.stat-mech] 4 Jan 2023