Topological phase transition in a symmetric blockade structure

AI summary

Overview

Research area: Quantum many-body physics — topologically ordered phases of matter in Rydberg-blockade-type lattice models, studied with large-scale tensor-network numerics.

Technical level: Advanced. The paper assumes familiarity with topological order, the toric code, lattice gauge symmetries, and density matrix renormalization group (DMRG) methods, though its central claims are stated in fairly physical terms.

Scope (one sentence): The paper maps out the full phase diagram of a $\mathbb{Z}_2$-symmetric blockade Hamiltonian on the honeycomb lattice using infinite-cylinder DMRG, establishing the stability of its toric-code topological order up to Rabi drives comparable to the system's characteristic energy scale.

What This Paper Is About

Previous rigorous work (Maier et al., PRX Quantum 6, 030340 (2025)) proved that a class of "fully symmetric" blockade Hamiltonians — transverse-field Ising models with specially engineered antiferromagnetic blockade graphs — host a topologically ordered ground state, but only in the limit of weak Rabi driving. This paper asks what happens beyond that limit: how far the topological phase survives as the Rabi drive $\Omega$ is increased, and what kind of transition (if any) ends it. The authors build a simplified but closely related version of that blockade structure, which needs fewer atoms per plaquette, and simulate it numerically to answer this question.

Key Contributions

1. A reduced, numerically tractable blockade model. The authors introduce a honeycomb-lattice Hamiltonian with four two-level systems per site, on-site blockade $V \to \infty$ (one excitation per site), uniform detuning $\Delta$, directional inter-site blockade $U$, and a homogeneous Rabi drive $\Omega$. It is shown to be a reduced version of the toric code blockade structure of Ref. [35], retaining its local $\mathbb{Z}_2$ plaquette symmetry while substantially shrinking the local Hilbert space (down to five states per site).

2. A full numerical phase diagram in the Rabi drive. Using infinite-size DMRG on cylinders of circumference $L_y = 2, 3, 4$ with an iMPS unit cell $L_x = 2$, the authors locate a single continuous phase transition at $\Omega_c/\Delta \approx 0.40$ for $U = K = 2\Delta$, separating the $\mathbb{Z}_2$ topologically ordered phase from a trivial paramagnet.

3. Identification of the transition's diagnostic signatures. The transition is characterized jointly by a closing of the excitation gap in the symmetric sector, a jump in the topological entanglement entropy from $\gamma_E \approx \ln 2$ to zero, and the vanishing of a Wilson loop operator expectation value that obeys a perimeter law in the topological phase.

4. A demonstration that the charge gap stays open. The gap between different symmetry sectors (the charge gap $\Delta_c$) remains finite through the transition, while charges become confined in the trivial phase — in contrast to the symmetric-sector gap, whose closing is consistent with flux-anyon condensation.

Main Findings

• Single continuous transition: A scan of the ground-state energy per unit cell across the transition in both directions shows smooth behavior with no hysteresis or metastability, consistent with a second-order transition at $\Omega_c/\Delta \approx 0.40$.

• Entanglement entropy signature: The entanglement entropy $S_0$ across the cylinder peaks sharply at the transition for all simulated circumferences. Deep in the topological phase ($\Omega \ll \Omega_c$) it approaches $S_0 = (L_y - 1)\ln 2$, the value expected for the toric code fixed point.

• Topological entanglement entropy: Extracted from the linear fit $S_0 = aL_y - \gamma_E$ over $L_y = 2, 3, 4$, the topological entanglement entropy is $\gamma_E \approx \ln 2$ for $\Omega < \Omega_c$ and drops to zero above $\Omega_c$; example fits are shown at $\Omega/\Delta = 0.3$ and $\Omega/\Delta = 0.6$.

• Wilson loop as an order parameter: $|\langle \Phi_y \rangle|$ is finite for $\Omega \leq \Omega_c$ and vanishes at the transition. Its finite-size scaling follows a perimeter law $|\langle \Phi_y \rangle| \sim e^{-L_y/\xi(\Omega)}$ with circumference, as expected for a topological phase.

• Symmetric-sector gap closes: The gap $\Delta_s$ for excitations with $U_p = +1$ on all plaquettes, computed for $L_y = 3$ at momenta $p_x = 0, \pi/2, 2\pi/3$, closes at the transition, consistent with condensation of flux anyons. It is a lower bound on the flux-anyon gap in the topological phase.

• Charge gap remains open: The charge gap $\Delta_c$ — the energy difference to the lowest state with two $U_p = -1$ plaquettes separated by $d = 2$ along $y$ — stays finite across the transition, and the inset compares this behavior for $U = {0.2, 2, 20}\Delta$ as a function of the rescaled distance from the transition. Charges become confined in the trivial phase.

• Topological phase is robust: The ordered phase survives up to Rabi drives comparable to the characteristic energy scale of the system, extending well beyond the weak-driving regime covered by the existing rigorous proofs.

• Note on truncation: The introduction states that the authors also derive the critical Rabi drive as a function of interaction strength $U$ using a reduced basis method; the detailed results of that scan (and the remainder of the Wilson-loop scaling discussion) lie beyond the truncated text supplied here and are not reported in this summary.

Methodology in Plain English

The authors treat each of the four two-level systems on a honeycomb site as an atom that can be in its ground or excited state, with strong repulsive "blockade" interactions preventing more than one excitation per site. Which atoms on neighboring sites block each other is decided by a rule that depends on the bond direction ($x$, $y$, or $z$), and this rule is what produces the local $\mathbb{Z}_2$ symmetry: a set of permutation operations that, multiplied around a plaquette, leave the Hamiltonian unchanged.

Because the on-site blockade is taken to infinity, each site's Hilbert space collapses to five states (empty, or one excitation in one of four positions). The authors then switch to a basis labeled by "Klein charges" in which the symmetry is diagonal, which makes it practical to enforce the symmetry in the numerics. An auxiliary term $K\sum_p(\mathbb{1}-U_p)$ is added; it does not change the eigenstates but separates the different symmetry sectors so the numerics can target the ground state and the charge excitations cleanly.

The simulations use infinite-system DMRG in the open-source TeNPy package, representing the ground state as an infinite matrix product state on a cylinder that is periodic in $y$. Bond dimensions up to $\chi = 1536$ were needed at $L_y = 4$ (with $\chi = 256$, or $\chi = 512$ near the transition, sufficient for $L_y = 2$), converging truncation errors below $10^{-8}$. Only the global charge $k_y = 0$ is conserved, deliberately not $k_x$, because for odd $L_y$ fixing $k_x$ would force the state into a superposition of topological sectors with $\langle\Phi_y\rangle = 0$ and spoil the extraction of the topological entanglement entropy. The residual $\langle H_K \rangle$ serves as a self-consistency check and as an estimate of the energy uncertainty.

Why This Matters

The paper bridges a rigorous but weak-driving mathematical result and the practical question of whether such topological order can be observed in a laboratory. It shows the topological phase occupies a large, physically accessible parameter regime and identifies concrete, measurable signatures — entanglement entropy scaling, Wilson loops, and gap behavior — that an experiment could target.

Applications (as motivated in the paper):

• Quantum memories: Topologically ordered phases are of interest as robust storage for quantum information because of their long-range entanglement.
• Topological quantum computation: Anyonic excitations in such phases are a candidate route to fault-tolerant quantum logic.
• Rydberg atom arrays and quantum simulators: Blockade interactions between Rydberg atoms are the physical origin of the model, making neutral-atom platforms a natural venue for realizing these phases.
• Alternative synthetic platforms: The authors explicitly state the door is opened for realizations on platforms beyond Rydberg arrays, since the construction relies only on two-body blockade interactions.

Industry relevance: The work is directly relevant to neutral-atom quantum computing and simulation efforts that engineer Rydberg blockade on programmable lattices, and to teams building symmetry-protected or topologically protected hardware-level encodings. It also provides benchmark-quality data for tensor-network methods applied to lattice gauge theories with local symmetries.

Future Directions

• Finite-size and bond-dimension extrapolation: Whether $\gamma_E \approx \ln 2$ and the Wilson loop perimeter law survive controlled extrapolation in $L_y$ and $\chi$ beyond the values used here is left to future, larger simulations.
• Critical exponents and universality class: The paper establishes that the transition is continuous but does not report a numerically extracted set of critical exponents or a universality-class assignment in the available text.
• Interaction-strength dependence: The full dependence of $\Omega_c$ on $U$ — promised via the reduced basis method — determines how much freedom experimentalists have in choosing blockade strengths without losing the topological phase.
• Experimental realization: Translating the four-two-level-systems-per-site, direction-dependent blockade graph into an actual Rydberg array layout, and identifying which of the predicted signatures (entanglement entropy, Wilson loop, charge confinement) is most accessible in a finite, noisy device.

Target Audience

Specialists in quantum many-body theory and tensor-network methods; researchers working on Rydberg atom arrays and neutral-atom quantum simulation; and anyone interested in routes to topologically ordered matter from two-body interactions rather than engineered multi-body terms. Physicists who already know the toric code and DMRG will get the most out of it, but the phase-diagram narrative and the diagnostic observables are readable by graduate students with a background in condensed matter or quantum information.

Authors’ abstract

We present a numerical study of a blockade Hamiltonian with a local $\mathbb{Z}_2$ symmetry. For such models, the existence of a topologically ordered ground state has been proven rigorously for weak Rabi driving [T. F. Maier et al., PRX Quantum 6, 030340 (2025)]. However, the phase diagram beyond weak driving is still an open question. We perform large-scale infinite-size density matrix renormalization group (iDMRG) simulations on a cylinder and map out the full phase diagram of the model. We find a single continuous phase transition between the topologically ordered phase and a trivial paramagnetic phase, and confirm numerically the stability of the topologically ordered phase up to Rabi drives comparable to the characteristic energy scale of the system. Furthermore, we compute the excitation gaps and show that the charge gap remains open across the phase transition and charges become confined in the trivial phase. On the other hand, the gap in the symmetric sector closes at the transition, consistent with the condensation of flux excitations.