pith. sign in

arxiv: 2605.18119 · v1 · pith:RCSPNWQXnew · submitted 2026-05-18 · ❄️ cond-mat.quant-gas

Weak Fragmentation and Thermalization in a Dipole-Conserving Bose-Hubbard Chain

Chenrong Liu This is my paper

Pith reviewed 2026-05-20 00:20 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas
keywords Hilbert-space fragmentationdipole conservationBose-Hubbard chaineigenstate thermalizationquantum chaosnonergodic dynamicsexact diagonalization
0
0 comments X

The pith

Weak fragmentation in a dipole-conserving Bose-Hubbard chain permits thermalization and quantum chaos at moderate interactions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines a one-dimensional Bose-Hubbard model that conserves both particle number and dipole moment. It shows that the Hamiltonian in the Fock basis exhibits weak Hilbert-space fragmentation, splitting the space into many sectors while still allowing mixing within them. An exponentially large family of frozen product states is constructed analytically with explicit bounds. Exact diagonalization on finite chains then tracks half-chain entanglement, density relaxation, and level statistics, revealing a crossover from volume-law entanglement and Wigner-Dyson statistics at weaker interactions to area-law entanglement and Poisson statistics at stronger interactions.

Core claim

In the dipole-conserving Bose-Hubbard chain the Hilbert space undergoes weak fragmentation that does not eliminate quantum chaos or thermalization; instead, increasing the on-site repulsion strength drives a crossover from a regime of eigenstate thermalization to a nonergodic regime while preserving the underlying dipole constraint.

What carries the argument

Weak Hilbert-space fragmentation, which partially decomposes the Fock space into disconnected sectors yet leaves enough connectivity for chaotic mixing at moderate interactions, together with the analytically bounded family of frozen product states.

If this is right

  • Half-chain entanglement entropy follows a volume law in the moderate-interaction regime.
  • Local density profiles relax to a uniform value at moderate interactions but remain frozen at strong interactions.
  • Nearest-neighbor level statistics cross from Wigner-Dyson to Poisson with rising interaction strength.
  • The model supplies a minimal lattice platform in which dipole conservation, weak fragmentation, and ergodicity breaking can be tuned against one another.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same weak-fragmentation mechanism may appear in other lattice models that conserve higher multipole moments.
  • Experimental realization in optical lattices with tilted potentials could test the predicted interaction-driven crossover directly.
  • The exponentially many frozen states suggest a route to prethermal plateaus whose lifetime grows with system size.

Load-bearing premise

Exact diagonalization results on small finite chains capture the qualitative location and nature of the transition from thermalizing to nonergodic behavior in the thermodynamic limit.

What would settle it

A direct measurement, in a larger system or via tensor-network methods, of whether the half-chain entanglement entropy remains volume-law and level statistics remain Wigner-Dyson up to interaction strengths where current small-system data already show Poisson statistics.

Figures

Figures reproduced from arXiv: 2605.18119 by Chenrong Liu.

Figure 1
Figure 1. Figure 1: In contrast to the standard Bose-Hubbard model, single-particle hopping is forbidden in this system due to the conservation of the dipole moment. The kinetic term ˆb † i ˆb 2 i+1 ˆb † i+2 represents a correlated hopping process in which two particles simultaneously tunnel from site i + 1 to the nearest-neighbor sites i and i + 2, a mech￾anism that ensures dipole conservation. For this reason, a single part… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) and (b), for a large value of J/U, the entan￾glement entropies of most eigenstates are located in a relatively narrow band. Near the middle of the energy spectrum, these values are close to the Page value, con￾sistent with ergodic behavior and thermalization. These states come from the dominant Krylov sector. At the same time, a small number of outlier states, including some in the middle of the spectr… view at source ↗
Figure 6
Figure 6. Figure 6: (a) shows, the sharp central peak gradually broad￾ens, and the correlator approaches a Gaussian-like dis￾tribution over the lattice. This behavior indicates the ergodic and ETH features and reveals that the local ob￾servables lose memory of their initial conditions and relax toward values determined by thermal equilibrium. While for small J/U, as shown in [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
read the original abstract

We study Hilbert-space fragmentation and thermalization in a one-dimensional dipole-conserving Bose-Hubbard chain. By analyzing the structure of the Hamiltonian matrix in the Fock basis, we show that the system exhibits weak Hilbert-space fragmentation. We further construct an exponentially large family of frozen product states and derive analytical upper and lower bounds on their number. Using exact diagonalization, we examine the consequences of weak fragmentation for eigenstate half-chain entanglement, density relaxation dynamics, and level statistics. All these quantities reveal a transition from a weak eigenstate thermalization regime to a nonergodic regime with increasing on-site interaction strength. These results show that weak Hilbert-space fragmentation \textit{does not} preclude quantum chaos or thermalization, and provides a minimal platform for studying the interplay of dipole conservation, weak fragmentation, and ergodicity breaking.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript studies Hilbert-space fragmentation and thermalization in a one-dimensional dipole-conserving Bose-Hubbard chain. By analyzing the Hamiltonian matrix structure in the Fock basis, it establishes weak fragmentation. It constructs an exponentially large family of frozen product states and derives analytical upper and lower bounds on their number. Using exact diagonalization on finite chains, it examines half-chain entanglement, density relaxation dynamics, and level statistics, revealing a crossover from a weak eigenstate thermalization regime to a nonergodic regime as the on-site interaction U increases. The work concludes that weak fragmentation does not preclude quantum chaos or thermalization.

Significance. If the central claims hold, the results provide a minimal, tunable platform for investigating the interplay of dipole conservation, weak Hilbert-space fragmentation, and ergodicity breaking. The analytical construction and bounds on frozen states, combined with the ED diagnostics of the transition, add concrete value to the fragmentation literature by showing that subextensive fragment growth can still allow regimes of thermalization and chaos. This strengthens understanding beyond strong-fragmentation models.

major comments (1)
  1. [Numerical diagnostics / exact-diagonalization results] The transition from weak ETH to nonergodicity is diagnosed via exact diagonalization for chains with L up to approximately 16 (as implied by the accessible system sizes in the numerical section). Given that dipole conservation restricts effective Hilbert-space growth to roughly exp(c√L) rather than exp(L), fragment sizes and mixing rates are expected to show pronounced finite-size dependence. No finite-size scaling collapse, extrapolation of the crossover value of U, or L-dependence of the spectral statistics (e.g., average r-ratio) or relaxation timescales is reported. This leaves open whether the reported weak-ETH regime persists or shifts in the thermodynamic limit, which is load-bearing for the claim that weak fragmentation permits thermalization.
minor comments (2)
  1. [Section on frozen product states] The abstract states that analytical upper and lower bounds on the number of frozen states are derived, but the main text does not explicitly compare the tightness of these bounds to the numerically counted frozen states for the system sizes studied.
  2. [Figure captions and legends] In the figures displaying level statistics and entanglement, reference lines for Poisson and Wigner-Dyson distributions (or Page value for entanglement) would improve clarity for readers.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address the major comment regarding the numerical diagnostics and finite-size effects below.

read point-by-point responses

  1. Referee: [Numerical diagnostics / exact-diagonalization results] The transition from weak ETH to nonergodicity is diagnosed via exact diagonalization for chains with L up to approximately 16 (as implied by the accessible system sizes in the numerical section). Given that dipole conservation restricts effective Hilbert-space growth to roughly exp(c√L) rather than exp(L), fragment sizes and mixing rates are expected to show pronounced finite-size dependence. No finite-size scaling collapse, extrapolation of the crossover value of U, or L-dependence of the spectral statistics (e.g., average r-ratio) or relaxation timescales is reported. This leaves open whether the reported weak-ETH regime persists or shifts in the thermodynamic limit, which is load-bearing for the claim that weak fragmentation permits thermalization.

    Authors: We agree that a more detailed analysis of finite-size dependence would strengthen the numerical evidence. In the revised manuscript we will add explicit L-dependence plots (for L = 8, 10, 12, 14, 16) of the average r-ratio and the density relaxation timescale as functions of U. These will show that the crossover from near-Wigner-Dyson statistics and fast relaxation at weak U to Poisson-like statistics and slow relaxation at strong U remains qualitatively stable across the accessible sizes. We will also include a brief discussion of the expected scaling based on our analytical bounds, which establish that the largest fragments grow at most as exp(O(√L)). Because the effective Hilbert-space dimension remains sub-exponential, the observed mixing at weak U is not an artifact of small L. A quantitative scaling collapse or extrapolation of the crossover U_c is beyond the scope of exact diagonalization; however, the combination of the sub-extensive fragment bounds and the consistent numerical trends supports the conclusion that weak fragmentation permits thermalization for sufficiently weak interactions. revision: partial

standing simulated objections not resolved
  • A definitive statement on the location of the crossover in the strict thermodynamic limit would require either substantially larger-scale numerics or an analytical theory of the fragment mixing rates, neither of which is currently available.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper derives its claims through direct inspection of the Hamiltonian matrix elements in the Fock basis to establish weak fragmentation, explicit construction of an exponentially large family of frozen product states together with analytical upper and lower bounds on their number, and standard exact-diagonalization computations of half-chain entanglement entropy, density relaxation, and level statistics as functions of interaction strength U. None of these steps reduces by construction to a fitted parameter renamed as a prediction, a self-definitional loop, or a load-bearing self-citation whose validity is presupposed by the present work. The central statement that weak fragmentation permits a regime of quantum chaos and thermalization is an empirical observation extracted from the finite-size spectra and dynamics rather than an algebraic identity or re-labeling of prior inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard quantum-mechanical description of lattice bosons and the validity of finite-size exact diagonalization for identifying the thermalization transition; no ad-hoc fitted parameters or new postulated entities are introduced.

axioms (1)
  • domain assumption The system is described by a dipole-conserving Bose-Hubbard Hamiltonian on a one-dimensional chain
    This model choice defines the Hilbert space and conservation law used for all fragmentation and dynamics analysis.

pith-pipeline@v0.9.0 · 5664 in / 1285 out tokens · 79256 ms · 2026-05-20T00:20:59.913071+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages · 1 internal anchor

  1. [1]

    and raise several natural questions. Does a minimal dipole-conserving Bose-Hubbard chain exhibit strong or weak fragmentation? How abundant are frozen states in the physically symmetric sectors? Can weak fragmenta- tion coexist with a quantum chaotic energy spectrum and thermalization? Or can the dipole constraint prevent er- godicity? Addressing these qu...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    cnoP6GrsjoC///jeksuSfZ7jJQE=

    (7) Here f0 = (1 + √ 5)/2 is the Golden ratio. Numerically, we can apply a dynamic programming (DP) algorithm to exactly solve this counting problem and get C for a larger system beyond the computation limits of the ED. Details about this algorithm can be seen in AppendixB. The numerical results are presented in Fig. 3. As it shows, the number of frozen s...

    work page 2000

  3. [3]

    In the bulk (i.e., i = −k + 1,

    At the edge sites, n−k and nk may be any integer from 0 to N. In the bulk (i.e., i = −k + 1, . . . , k− 1), the occupation numbers are binary, ni ∈ {0, 1}

    work page

  4. [4]

    Distance-2 constraint: For any sites separated by a distance of 2, the occupations satisfy nini+2 = 0

    work page

  5. [5]

    The configuration satisfies two global conservation laws, ∑ i ni = N and ∑ i ini = 0. An exact expression for the total number of frozen states is difficult to obtain because the counting problem is not purely local: the admissible configurations {ni} are additionally constrained by two global conservation laws, which couple all sites and substantially mod...

    work page

  6. [6]

    Let A(i) be the set of allowed values of ni from the constraints, then if u = 0 , A(i) = {0, 1} for bulk sites and A(i) = {0, 1, · · · , Ntarget − n} for boundary sites

    Allowed ni is determinded by u. Let A(i) be the set of allowed values of ni from the constraints, then if u = 0 , A(i) = {0, 1} for bulk sites and A(i) = {0, 1, · · · , Ntarget − n} for boundary sites. While if u = 1, then only A(i) = {0} is permitted. 12 Algorithm 1: Dynamic programming algorithm Input: L = 2k + 1, Ntarget, Ptarget. Sites i =k, . . . , k...

    work page

  7. [7]

    (B2) We then accumulate Wi(u′, v′; n′, p′) = Wi(u′, v′; n′, p′) + Wi−1(u, v; n, p)

    For each ni ∈ A(i), define binary variblew = si and update n′ = n + ni, p ′ = p + ini, (u′, v′) = ( v, w). (B2) We then accumulate Wi(u′, v′; n′, p′) = Wi(u′, v′; n′, p′) + Wi−1(u, v; n, p). (B3)

    work page

  8. [8]

    Let r = Ntarget − n′ be the number of particles remaining after the update, and note that each future particle can change the dipole moment by at most k

    Feasibility pruning. Let r = Ntarget − n′ be the number of particles remaining after the update, and note that each future particle can change the dipole moment by at most k. Hence, a necessary condition to still be able to reach Ptarget is |Ptarget − p′| ≤ rk. (B4) A partial configuration violating this bound cannot lead to a valid completion and is disc...

    work page

  9. [9]

    cnoP6GrsjoC///jeksuSfZ7jJQE=

    After processing all sites, the desired count is ob- tained by summing over the four terminal window states: count = ∑ (u,v)∈{0,1}2 Wk(u, v; Ntarget, Ptarget). (B5) We initialize at the left boundary (before processing any site) with W−k−1(u = 0, v = 0; n = 0, p = 0) = 1 and all other entries zero. In practice, W is sparse and only a small subset of (n, p...

    work page

  10. [10]

    J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A 43, 2046 (1991)

    work page 2046

  11. [11]

    Srednicki, Chaos and quantum thermalization, Phys

    M. Srednicki, Chaos and quantum thermalization, Phys. Rev. E 50, 888 (1994)

    work page 1994

  12. [12]

    Rigol, V

    M. Rigol, V. Dunjko, and M. Olshanii, Thermalization and its mechanism for generic isolated quantum systems, Nature 452, 854 (2008)

    work page 2008

  13. [13]

    T. c. v. Prosen, Ergodic properties of a generic nonin- tegrable quantum many-body system in the thermody- namic limit, Phys. Rev. E 60, 3949 (1999)

    work page 1999

  14. [14]

    F. H. L. Essler and M. Fagotti, Quench dynamics and relaxation in isolated integrable quantum spin chains, Journal of Statistical Mechanics: Theory and Experi- ment 2016, 064002 (2016)

    work page 2016

  15. [15]

    Nandkishore and D

    R. Nandkishore and D. A. Huse, Many-body localiza- tion and thermalization in quantum statistical mechan- ics, Annual Review of Condensed Matter Physics 6, 15 (2015)

    work page 2015

  16. [16]

    S. D. Geraedts, R. Nandkishore, and N. Regnault, Many- body localization and thermalization: Insights from the entanglement spectrum, Phys. Rev. B 93, 174202 (2016)

    work page 2016

  17. [17]

    Colmenarez, P

    L. Colmenarez, P. A. McClarty, M. Haque, and D. J. Luitz, Statistics of correlation functions in the random Heisenberg chain, SciPost Phys. 7, 064 (2019)

    work page 2019

  18. [18]

    Wang, X.-F

    J.-w. Wang, X.-F. Zhou, G.-C. Guo, and Z.-W. Zhou, Thermalization of quantum many-body scars in kinet- ically constrained systems, Phys. Rev. Res. 7, 043337 (2025)

    work page 2025

  19. [19]

    Shibata, N

    N. Shibata, N. Yoshioka, and H. Katsura, Onsager’s scars in disordered spin chains, Phys. Rev. Lett. 124, 180604 (2020)

    work page 2020

  20. [20]

    A. M. Alhambra, A. Anshu, and H. Wilming, Revivals imply quantum many-body scars, Phys. Rev. B 101, 205107 (2020)

    work page 2020

  21. [21]

    Omiya and M

    K. Omiya and M. Müller, Fractionalization paves the way to local projector embeddings of quantum many-body scars, Phys. Rev. B 108, 054412 (2023)

    work page 2023

  22. [22]

    Hummel, K

    Q. Hummel, K. Richter, and P. Schlagheck, Genuine many-body quantum scars along unstable modes in bose- hubbard systems, Phys. Rev. Lett. 130, 250402 (2023)

    work page 2023

  23. [23]

    Z.-C. Yang, F. Liu, A. V. Gorshkov, and T. Iadecola, Hilbert-space fragmentation from strict confinement, Phys. Rev. Lett. 124, 207602 (2020)

    work page 2020

  24. [24]

    Roy and A

    S. Roy and A. Lazarides, Strong ergodicity breaking due to local constraints in a quantum system, Phys. Rev. Res. 2, 023159 (2020)

    work page 2020

  25. [25]

    Łydżba, P

    P. Łydżba, P. Prelovšek, and M. Mierzejewski, Local inte- grals of motion in dipole-conserving models with hilbert space fragmentation, Phys. Rev. Lett. 132, 220405 (2024)

    work page 2024

  26. [26]

    Głódkowski, F

    A. Głódkowski, F. Peña Benítez, and P. Surówka, Hydro- dynamics of dipole-conserving fluids, Phys. Rev. E 107, 034142 (2023)

    work page 2023

  27. [27]

    Morningstar, V

    A. Morningstar, V. Khemani, and D. A. Huse, Kinetically constrained freezing transition in a dipole-conserving sys- tem, Phys. Rev. B 101, 214205 (2020)

    work page 2020

  28. [28]

    Scherg, T

    S. Scherg, T. Kohlert, P. Sala, F. Pollmann, B. Hebbe Madhusudhana, I. Bloch, and M. Aidelsburger, Observing non-ergodicity due to kinetic constraints in tilted fermi-hubbard chains, Nature Communications 12, 4490 (2021)

    work page 2021

  29. [29]

    Z. Lan, M. van Horssen, S. Powell, and J. P. Garrahan, Quantum slow relaxation and metastability due to dy- namical constraints, Phys. Rev. Lett. 121, 040603 (2018)

    work page 2018

  30. [30]

    Zechmann, J

    P. Zechmann, J. Boesl, J. Feldmeier, and M. Knap, Dy- namical spectral response of fractonic quantum matter, Phys. Rev. B 109, 125137 (2024)

    work page 2024

  31. [31]

    Boesl, P

    J. Boesl, P. Zechmann, J. Feldmeier, and M. Knap, De- confinement dynamics of fractons in tilted bose-hubbard chains, Phys. Rev. Lett. 132, 143401 (2024)

    work page 2024

  32. [32]

    E. Lake, M. Hermele, and T. Senthil, Dipolar bose- hubbard model, Phys. Rev. B 106, 064511 (2022)

    work page 2022

  33. [33]

    Zechmann, E

    P. Zechmann, E. Altman, M. Knap, and J. Feldmeier, Fractonic luttinger liquids and supersolids in a con- strained bose-hubbard model, Phys. Rev. B 107, 195131 (2023)

    work page 2023

  34. [34]

    Lake, H.-Y

    E. Lake, H.-Y. Lee, J. H. Han, and T. Senthil, Dipole condensates in tilted bose-hubbard chains, Phys. Rev. B 107, 195132 (2023)

    work page 2023

  35. [35]

    Y.-T. Oh, J. H. Han, and H.-Y. Lee, Fractonic quantum quench in dipole-constrained bosons, Phys. Rev. Res. 6, 023269 (2024)

    work page 2024

  36. [36]

    Liu, Entanglement and fluctuations of the frac- tonic phases in a one-dimensional dipole-conserving bose- hubbard chain, Phys

    C. Liu, Entanglement and fluctuations of the frac- tonic phases in a one-dimensional dipole-conserving bose- hubbard chain, Phys. Rev. B 112, 245139 (2025)

    work page 2025

  37. [37]

    Adler, D

    D. Adler, D. Wei, M. Will, K. Srakaew, S. Agrawal, P. Weckesser, R. Moessner, F. Pollmann, I. Bloch, and 14 J. Zeiher, Observation of hilbert space fragmentation and fractonic excitations in 2d, Nature 636, 80 (2024)

    work page 2024

  38. [38]

    Kohlert, S

    T. Kohlert, S. Scherg, P. Sala, F. Pollmann, B. Hebbe Madhusudhana, I. Bloch, and M. Aidelsburger, Exploring the regime of fragmentation in strongly tilted fermi-hubbard chains, Phys. Rev. Lett. 130, 010201 (2023)

    work page 2023

  39. [39]

    Moudgalya and O

    S. Moudgalya and O. I. Motrunich, Hilbert space frag- mentation and commutant algebras, Phys. Rev. X 12, 011050 (2022)

    work page 2022

  40. [40]

    Srivastava and S

    A. Srivastava and S. Dutta, Hierarchy of degenerate sta- tionary states in a boundary-driven dipole-conserving spin chain, SciPost Phys. 18, 111 (2025)

    work page 2025

  41. [41]

    P. Frey, D. Mikhail, S. Rachel, and L. Hackl, Probing hilbert space fragmentation and the block inverse partic- ipation ratio, Phys. Rev. B 109, 064302 (2024)

    work page 2024

  42. [42]

    Jeyaretnam, T

    J. Jeyaretnam, T. Bhore, J. J. Osborne, J. C. Halimeh, and Z. Papić, Hilbert space fragmentation at the origin of disorder-free localization in the lattice schwinger model, Communications Physics 8, 172 (2025)

    work page 2025

  43. [43]

    Francica and L

    G. Francica and L. Dell’Anna, Hilbert space fragmenta- tion in a long-range system, Phys. Rev. B 108, 045127 (2023)

    work page 2023

  44. [44]

    P. Sala, T. Rakovszky, R. Verresen, M. Knap, and F. Pollmann, Ergodicity breaking arising from hilbert space fragmentation in dipole-conserving hamiltonians, Phys. Rev. X 10, 011047 (2020)

    work page 2020

  45. [45]

    Lehmann, P

    J. Lehmann, P. S. de Torres-Solanot, F. Pollmann, and T. Rakovszky, Fragmentation-induced localization and boundary charges in dimensions two and above, SciPost Phys. 14, 140 (2023)

    work page 2023

  46. [46]

    Khudorozhkov, A

    A. Khudorozhkov, A. Tiwari, C. Chamon, and T. Neu- pert, Hilbert space fragmentation in a 2D quantum spin system with subsystem symmetries, SciPost Phys. 13, 098 (2022)

    work page 2022

  47. [47]

    Adler, D

    D. Adler, D. Wei, M. Will, K. Srakaew, S. Agrawal, P. Weckesser, R. Moessner, F. Pollmann, I. Bloch, and J. Zeiher, Observation of hilbert space fragmentation and fractonic excitations in 2d, Nature 636, 80 (2024)

    work page 2024

  48. [48]

    Russomanno, M

    A. Russomanno, M. Fava, and R. Fazio, Nonergodic be- havior of the clean bose-hubbard chain, Phys. Rev. B 102, 144302 (2020)

    work page 2020

  49. [49]

    Saad, Analysis of some krylov subspace ap- proximations to the matrix exponential operator, SIAM Journal on Numerical Analysis 29, 209 (1992) , https://doi.org/10.1137/0729014

    Y. Saad, Analysis of some krylov subspace ap- proximations to the matrix exponential operator, SIAM Journal on Numerical Analysis 29, 209 (1992) , https://doi.org/10.1137/0729014

    work page doi:10.1137/0729014 1992

  50. [50]

    Y. Y. Atas, E. Bogomolny, O. Giraud, and G. Roux, Distribution of the ratio of consecutive level spacings in random matrix ensembles, Phys. Rev. Lett. 110, 084101 (2013)

    work page 2013

  51. [51]

    Oganesyan and D

    V. Oganesyan and D. A. Huse, Localization of interacting fermions at high temperature, Phys. Rev. B 75, 155111 (2007)

    work page 2007

  52. [52]

    Z. Wang, Y. Wang, and B. Wu, Quantum chaos and phys- ical distance between quantum states, Phys. Rev. E 103, 042209 (2021)

    work page 2021

  53. [53]

    D. J. Luitz, N. Laflorencie, and F. Alet, Many-body local- ization edge in the random-field heisenberg chain, Phys. Rev. B 91, 081103 (2015)

    work page 2015

  54. [54]

    E. R. Castro, J. Chávez-Carlos, I. Roditi, L. F. Santos, and J. G. Hirsch, Quantum-classical correspondence of a system of interacting bosons in a triple-well potential, Quantum 5, 563 (2021)

    work page 2021

  55. [55]

    Nakerst and M

    G. Nakerst and M. Haque, Chaos in the three-site bose- hubbard model: Classical versus quantum, Phys. Rev. E 107, 024210 (2023)

    work page 2023

  56. [56]

    Kullback and R

    S. Kullback and R. A. Leibler, On Information and Suf- ficiency, The Annals of Mathematical Statistics 22, 79 (1951)

    work page 1951

  57. [57]

    weak fragmentation and ther- malization in a dipole-conserving bose-hubbard chain

    C. Liu, Data for the paper “weak fragmentation and ther- malization in a dipole-conserving bose-hubbard chain” (2026), doi: 10.5281/zenodo.20137063

    work page doi:10.5281/zenodo.20137063 2026