Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2408.10321 · Physical review. B./Physical review. B

Universal freezing transitions of dipole-conserving chains

J. Classen-Howes, Riccardo Senese, Abhishodh Prakash

Published Aug 20, 2025• 4 citations
arXivHTMLPDF

Abstract

We demonstrate the existence of a universal phase diagram of quantum chains with range- k interactions subject to the conservation of a total charge and its dipole moment. These systems exhibit “freezing” transitions between strongly and weakly Hilbert-space-fragmented phases as the charge filling ν is varied. We show that these continuous phase transitions occur at a critical charge filling of ν c = ( k − 2 ) − 1 of the on-site Hilbert-space dimension d . To this end, we analytically prove that, for any d , any state with ν < ν c hosts a finite density of sites belonging to “blockages,” which we define as subregions of the chain across which transport of charge and dipole moment cannot occur. Some blockages arise from sequences of frozen sites, i.e., sites with an unchanging on-site charge, while others do not involve frozen sites at all. We prove that the presence of blockages implies strong fragmentation of typical symmetry sectors into Krylov subspaces, each of which forms an exponentially vanishing fraction of the total sector. By studying the distribution of blockages we analytically characterize how typical states are subdivided into dynamically disconnected local “active bubbles” and prove that typical eigenstates at these charge fillings exhibit area-law entanglement entropy, while there exist rare eigenstates featuring non-area-law scaling. We also numerically show that for ν > ν c and arbitrary d , typical symmetry sectors are weakly fragmented, with their dominant Krylov sectors constituted of states that are free of blockages. We analytically derive some critical exponents characterizing the transition and numerically determine the density of blockages at ν = ν c , with clear implications for transport at the critical point. Finally, we investigate the properties of special-case models for which no phase transitions occur.

DOI
10.1103/2h1v-yx5l

Related papers

Kinetically constrained freezing transition in a dipole-conserving system

Alan Morningstar, Vedika Khemani, David A. Huse

We study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fractonlike conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase, where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents.

arXiv:2004.00096•Physical review. B./Physical review. B•Jun 19, 2020•123 citations

Freezing transition in the particle-conserving East model

Cheng Wang, Zhi-Cheng Yang

Quantum kinetically constrained models can exhibit a wealth of dynamical phenomena ranging from anomalous transport to Hilbert-space fragmentation (HSF). We study a class of one-dimensional particle-number-conserving systems in which particle hoppings are subjected to an East-like constraint, akin to facilitated spin models in classical glasses. While such a kinetic constraint leads to HSF, we find that the degree of fragmentation exhibits a sharp transition as the average particle density is varied. Below a critical density, the system transitions from being weakly fragmented where most of the initial states thermalize diffusively to strongly fragmented where the dynamics are frozen and the system fails to thermalize. Remarkably, the East model allows for both efficient numerical simulations and analytic solutions of various diagnostics of the phase transition, from which we obtain a set of exact critical exponents. We find that the freezing transition in particle-conserving East models belongs to the same universality class as dipole-conserving fracton systems. Our results provide a tractable minimal model for filling-induced freezing transitions associated with HSF, which can be readily tested on state-of-the-art quantum platforms.

arXiv:2307.01993•Physical review. B./Physical review. B•Oct 31, 2023•10 citations

Structure of many-body quantum dynamics: full counting statistics, fragmentation and integrability

Riccardo Senese

This thesis investigates the dynamics of quantum many-body systems in isolation, in three different setups. We start by considering a paradigmatic quantum chaotic model of magnetism, the Heisenberg antiferromagnet in 2 and 3 spatial dimensions. The focus is on obtaining approximations for the full counting statistics (FCS) of a subsystem order parameter, the staggered magnetization, both in and out of equilibrium (following quantum quenches). We employ mappings of spin variables to bosonic ones, together with self-consistent time-dependent mean-field schemes, to compute the FCS via Gaussian techniques. In particular, we introduce a formula for the FCS of any subsystem observable that is at most quadratic in the bosons. We obtain results both in the presence and absence of long-range magnetic order. We then turn to the first of two ergodicity-breaking scenarios, by considering Hilbert space fragmentation (HSF) in charge- and dipole-conserving quantum chains. We determine the critical charge density at which a "freezing" phase transition from weak to strong fragmentation takes place, for any on-site Hilbert space dimension. To characterize the degree of dynamical connectivity in the system, we introduce a novel framework that enables us to obtain analytic results for the strongly fragmented phase, and numerically explore the weakly fragmented one. We derive some critical exponents, and show that the entanglement entropy of typical eigenstates in the strongly fragmented phase obeys an area law. We also uncover the existence of rare eigenstates whose entropy appears to scale beyond any area law. Finally, we investigate dynamical correlators in the Lieb-Liniger gas, one of the simplest interacting integrable theories solvable by Bethe ansatz. We employ Markov chain Monte Carlo methods to sample matrix elements (form factors) that appear in the Lehmann representation of the finite-temperature single-particle bosonic Green’s function. In the physical regimes we consider, sampling is necessary because a full summation of the relevant form factors is computationally intractable. We benchmark our results against exact solutions available in specific limits, and find them to be very accurate.

•Jan 1, 2025•0 citations

Minimal model for Hilbert space fragmentation with local constraints

Bhaskar Mukherjee, Debasish Banerjee, K. Sengupta +1 more

Interacting quantum systems may exhibit fragmentation of their Hilbert space into dynamically disconnected sectors leading to unusual properties like nonergodicity. Fragmentation usually stems from the simultaneous conservation of charge and dipole moment. The authors demonstrate that an alternative scenario occurs in the absence of dipole conservation via a minimal model of interacting spins where fragmentation, including a novel secondary fragmentation in an entangled basis, emerges. A Floquet version of the model yields detectable stroboscopic freezing of certain states at specific drive frequencies.

arXiv:2106.14897•Physical review. B./Physical review. B•Oct 11, 2021•64 citations

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

Chenrong Liu

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.

arXiv:2605.18119•May 18, 2026•0 citations

Quantum Hilbert Space Fragmentation and Entangled Frozen States

Zihan Zhou, Tian-Hua Yang, Bo-Ting Chen

We find that rank deficiency of the local Hamiltonian in a classically fragmented model is the key mechanism leading to quantum Hilbert space fragmentation. The rank deficiency produces local null directions that can generate entangled frozen states (EFS): entangled states embedded in mobile classical Krylov sectors that do not evolve under Hamiltonian dynamics. When the entangled frozen subspace is non-empty, the mobile classical sector splits into a mobile quantum Krylov subspace and an entangled frozen subspace, and the model exhibits quantum fragmentation. We establish this mechanism in four models of increasing symmetry structure: an asymmetric qubit projector with no symmetry, the $\mathbb{Z}_2$-symmetric GHZ projector, a $\mathbb{Z}_3$-symmetric cyclic qutrit projector, and the Temperley-Lieb model. For the asymmetric and GHZ projector models, we obtain closed-form expressions for irreducible Krylov dimensions, degeneracies, and sector multiplicities. The all-mobile-sector EFS in these two models exhibits a sub-volume-law bipartite entanglement entropy scaling as $S \sim \sqrt{L}$. Further, we introduce the notion of weak and strong quantum fragmentation, the quantum counterpart of the weak-strong distinction in classical fragmentation. After removing the EFS, the mobile quantum Krylov subspace decomposes into irreducible blocks. In the weak case, the number of irreducible blocks remains $O(1)$, each is individually ergodic with Gaussian Orthogonal Ensemble (GOE) level statistics, and the unresolved spectrum follows an $m$GOE distribution. In the strong case, the number of irreducible blocks grows with system size, and the gap-ratio distribution approaches Poisson as $L\to\infty$.

quant-pharXiv:2604.05218•Apr 8, 2026•0 citations
We demonstrate the existence of a universal phase diagram of quantum chains with range- k interactions subject to the conservation of a total charge and its dipole moment. These systems exhibit “freezing” transitions between strongly and weakly Hilbert-space-fragmented phases as the charge filling ν is varied. We show that these continuous phase transitions occur at a critical charge filling of ν c = ( k − 2 ) − 1 of the on-site Hilbert-space dimension d . To this end, we analytically prove that, for any d , any state with ν < ν c hosts a finite density of sites belonging to “blockages,” which we define as subregions of the chain across which transport of charge and dipole moment cannot occur. Some blockages arise from sequences of frozen sites, i.e., sites with an unchanging on-site charge, while others do not involve frozen sites at all. We prove that the presence of blockages implies strong fragmentation of typical symmetry sectors into Krylov subspaces, each of which forms an exponentially vanishing fraction of the total sector. By studying the distribution of blockages we analytically characterize how typical states are subdivided into dynamically disconnected local “active bubbles” and prove that typical eigenstates at these charge fillings exhibit area-law entanglement entropy, while there exist rare eigenstates featuring non-area-law scaling. We also numerically show that for ν > ν c and arbitrary d , typical symmetry sectors are weakly fragmented, with their dominant Krylov sectors constituted of states that are free of blockages. We analytically derive some critical exponents characterizing the transition and numerically determine the density of blockages at ν = ν c , with clear implications for transport at the critical point. Finally, we investigate the properties of special-case models for which no phase transitions occur.
We study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fractonlike conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase, where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents.
Quantum kinetically constrained models can exhibit a wealth of dynamical phenomena ranging from anomalous transport to Hilbert-space fragmentation (HSF). We study a class of one-dimensional particle-number-conserving systems in which particle hoppings are subjected to an East-like constraint, akin to facilitated spin models in classical glasses. While such a kinetic constraint leads to HSF, we find that the degree of fragmentation exhibits a sharp transition as the average particle density is varied. Below a critical density, the system transitions from being weakly fragmented where most of the initial states thermalize diffusively to strongly fragmented where the dynamics are frozen and the system fails to thermalize. Remarkably, the East model allows for both efficient numerical simulations and analytic solutions of various diagnostics of the phase transition, from which we obtain a set of exact critical exponents. We find that the freezing transition in particle-conserving East models belongs to the same universality class as dipole-conserving fracton systems. Our results provide a tractable minimal model for filling-induced freezing transitions associated with HSF, which can be readily tested on state-of-the-art quantum platforms.
This thesis investigates the dynamics of quantum many-body systems in isolation, in three different setups. We start by considering a paradigmatic quantum chaotic model of magnetism, the Heisenberg antiferromagnet in 2 and 3 spatial dimensions. The focus is on obtaining approximations for the full counting statistics (FCS) of a subsystem order parameter, the staggered magnetization, both in and out of equilibrium (following quantum quenches). We employ mappings of spin variables to bosonic ones, together with self-consistent time-dependent mean-field schemes, to compute the FCS via Gaussian techniques. In particular, we introduce a formula for the FCS of any subsystem observable that is at most quadratic in the bosons. We obtain results both in the presence and absence of long-range magnetic order. We then turn to the first of two ergodicity-breaking scenarios, by considering Hilbert space fragmentation (HSF) in charge- and dipole-conserving quantum chains. We determine the critical charge density at which a "freezing" phase transition from weak to strong fragmentation takes place, for any on-site Hilbert space dimension. To characterize the degree of dynamical connectivity in the system, we introduce a novel framework that enables us to obtain analytic results for the strongly fragmented phase, and numerically explore the weakly fragmented one. We derive some critical exponents, and show that the entanglement entropy of typical eigenstates in the strongly fragmented phase obeys an area law. We also uncover the existence of rare eigenstates whose entropy appears to scale beyond any area law. Finally, we investigate dynamical correlators in the Lieb-Liniger gas, one of the simplest interacting integrable theories solvable by Bethe ansatz. We employ Markov chain Monte Carlo methods to sample matrix elements (form factors) that appear in the Lehmann representation of the finite-temperature single-particle bosonic Green’s function. In the physical regimes we consider, sampling is necessary because a full summation of the relevant form factors is computationally intractable. We benchmark our results against exact solutions available in specific limits, and find them to be very accurate.
Interacting quantum systems may exhibit fragmentation of their Hilbert space into dynamically disconnected sectors leading to unusual properties like nonergodicity. Fragmentation usually stems from the simultaneous conservation of charge and dipole moment. The authors demonstrate that an alternative scenario occurs in the absence of dipole conservation via a minimal model of interacting spins where fragmentation, including a novel secondary fragmentation in an entangled basis, emerges. A Floquet version of the model yields detectable stroboscopic freezing of certain states at specific drive frequencies.
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.
We find that rank deficiency of the local Hamiltonian in a classically fragmented model is the key mechanism leading to quantum Hilbert space fragmentation. The rank deficiency produces local null directions that can generate entangled frozen states (EFS): entangled states embedded in mobile classical Krylov sectors that do not evolve under Hamiltonian dynamics. When the entangled frozen subspace is non-empty, the mobile classical sector splits into a mobile quantum Krylov subspace and an entangled frozen subspace, and the model exhibits quantum fragmentation. We establish this mechanism in four models of increasing symmetry structure: an asymmetric qubit projector with no symmetry, the Z2\mathbb{Z}_2Z2​-symmetric GHZ projector, a Z3\mathbb{Z}_3Z3​-symmetric cyclic qutrit projector, and the Temperley-Lieb model. For the asymmetric and GHZ projector models, we obtain closed-form expressions for irreducible Krylov dimensions, degeneracies, and sector multiplicities. The all-mobile-sector EFS in these two models exhibits a sub-volume-law bipartite entanglement entropy scaling as S∼LS \sim \sqrt{L}S∼L​. Further, we introduce the notion of weak and strong quantum fragmentation, the quantum counterpart of the weak-strong distinction in classical fragmentation. After removing the EFS, the mobile quantum Krylov subspace decomposes into irreducible blocks. In the weak case, the number of irreducible blocks remains O(1)O(1)O(1), each is individually ergodic with Gaussian Orthogonal Ensemble (GOE) level statistics, and the unresolved spectrum follows an mmmGOE distribution. In the strong case, the number of irreducible blocks grows with system size, and the gap-ratio distribution approaches Poisson as L→∞L\to\inftyL→∞.