Mathematically rigorous statements on the spectral gap of quantum many-body systems in the thermodynamic limit are notoriously difficult to prove -- yet they are of fundamental importance for classifying quantum phases of matter. Here we prove the existence of a finite excitation gap for a particular Hamiltonian which was proposed in [T. F. Maier et al., PRX Quantum 6, 030340 (2025)] and is motivated by the Rydberg platform. The Hamiltonian exhibits only two-body blockade interactions between two-level systems and has a topologically ordered ground state in the toric code phase. We show that our result also applies to a broader class of blockade Hamiltonians which realize non-Abelian quantum double phases and were proposed in [H. P. Büchler et al., Phys. Rev. B 114, 065113 (2026)]. The proof builds on known gap stability results and exploits the local symmetry of the studied models.
Using density-matrix renormalization-group calculations for infinite cylinders, we elucidate the properties of the spin-liquid phase of the spin-$\frac{1}{2}\phantom{\rule{4pt}{0ex}}{J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the triangular lattice. We find four distinct ground states characteristic of a nonchiral,
Z2 topologically ordered state with vison and spinon excitations. We shed light on the interplay of topological ordering and global symmetries in the model by detecting fractionalization of time-reversal and space-group dihedral symmetries in the anyonic sectors, which leads to the coexistence of symmetry protected and intrinsic topological order. The anyonic sectors, and information on the particle statistics, can be characterized by degeneracy patterns and symmetries of the entanglement spectrum. We demonstrate the ground states on finite-width cylinders are short-range correlated and gapped; however, some features in the entanglement spectrum suggest that the system develops gapless spinonlike edge excitations in the large-width limit.
The effect of a local Kondo coupling and Hubbard interaction on the topological phase of the one-dimensional topological Kondo insulator (TKI) is numerically investigated using the infinite matrix-product state density-matrix renormalization group algorithm. The ground state of the TKI is a symmetry-protected topological (SPT) phase protected by inversion symmetry. It is found that on its own, the Hubbard interaction that tends to force fermions into a one-charge per site order is insufficient to destroy the SPT phase. However, when the local Kondo Hamiltonian term that favors a topologically trivial ground state with a one-charge per site order is introduced, the Hubbard interaction assists in the destruction of the SPT phase. This topological phase transition occurs in the charge sector where the correlation length of the charge excitation diverges while the correlation length of the spin excitation remains finite. The critical exponents, central charge, and the phase diagram separating the SPT phase from the topologically trivial phase are presented.
Using quantum Monte Carlo simulations, we map out the phase diagram of Hamiltonians interpolating between trivial and nontrivial bosonic symmetry-protected topological phases, protected by
Z2 and
Z23 symmetries, in two dimensions. In all cases, we find that the trivial and the topological phases are separated by an intermediate phase in which the protecting symmetry is spontaneously broken. Depending on the model, we identify a variety of magnetic orders on the triangular lattice, including ferromagnetism, $\sqrt{3}\ifmmode\times\else\texttimes\fi{}\sqrt{3}$ order, and stripe orders (both commensurate and incommensurate). Critical properties are determined through a finite-size scaling analysis. Possible scenarios regarding the nature of the phase transitions are discussed.
We revisit the phase diagram of the dimerized XXZ spin-
21 chain with nearest-neighbor couplings which was studied numerically in Phys. Rev. B 106, L201106 (2022). The model has isotropic
XY couplings which have a uniform value and
ZZ couplings which have a dimerized form, with strengths
Ja and
Jb on alternate bonds. We find a rich phase diagram in the region of positive
Ja, Jb. We provide a detailed understanding of the different phases and associated quantum phase transitions using a combination of mean-field theory, low-energy effective Hamiltonians, renormalization group calculations employing the technique of bosonization, and numerical calculations using the density-matrix renormalization group (DMRG) method. The phase diagram consists of two Ising paramagnetic phases called IPM
0 and IPM
π, and a phase with Ising Neel order called IN; all these phases are gapped. The phases IPM
0 and IPM
π are separated by a gapless phase transition line given by
0≤Ja=Jb≤1 which is described by a conformal field theory with central charge
c=1. There are two gapless phase transition lines separating IPM
0 from IN and IPM
π from IN; these are described by conformal field theories with
c=21 corresponding to quantum Ising transitions. The
c=1 line bifurcates into the two
c=21 lines at the point
Ja=Jb=1; the shape of the bifurcation is found analytically using RG calculations. A symmetry analysis shows that IPM
0 is a topologically trivial phase while IPM
π is a time-reversal symmetry-protected topological phase (SPT) with spin-
21 states at the two ends of an open system. The numerical results obtained by the DMRG method are in good agreement with the analytical results. Finally we propose experimental platforms for testing our results.
We perform large-scale density matrix renormalization group simulations of the SU(4) Heisenberg model on the honeycomb lattice to address the long-standing question of its ground state in an unbiased and quantitatively controlled manner. We find reliable numerical evidence that the ground state is a gapped spin-orbital liquid, presumably with a Z_{4} topological order, characterized by a finite topological entanglement entropy close to ln(4), the absence of both SU(4) and lattice symmetry breaking, and a variationally optimized ground-state energy well below the previously proposed π-flux variational state. By exploiting full SU(4) symmetry and keeping up to 12 800 SU(4) multiplets, corresponding to more than 1×10^{6} U(1) states, we achieve unprecedented accuracy for two-dimensional SU(4) quantum magnets. Finite-size scaling of energies and entanglement entropies supports a robust gapped phase in the two-dimensional limit, while a gapless critical state on narrow cylinders is identified as a proximate remnant of a Dirac spin-orbital liquid. Our results find the SU(4) honeycomb Heisenberg model a realization of a gapped topological spin-orbital liquid and provide convincing numerical evidence for topological order in a highly symmetric two-dimensional quantum magnet.