The toric code under antiferromagnetic isotropic Heisenberg interactions

This study employs neural-network quantum states and Schrieffer-Wolff transformations to demonstrate that isotropic antiferromagnetic Heisenberg perturbations drive the toric code through a quantum phase transition into a fourfold-degenerate Néel phase, characterized by a critical point identified via fidelity and entanglement diagnostics.

The Gist

Imagine you have a very special, magical tablecloth. This isn't just any tablecloth; it's a Toric Code. In the world of quantum physics, this tablecloth represents a state of matter that is incredibly robust. It holds its secrets (quantum information) not in the fabric itself, but in the knots and loops woven into it. Because these knots are global properties of the whole cloth, you can't untie them by just poking or prodding a single spot. This makes the Toric Code a perfect candidate for a quantum memory, a hard drive that never loses data because it's protected by the laws of topology.

However, in the real world, nothing is perfect. Your magical tablecloth is sitting in a room where the air is slightly turbulent. These air currents represent Heisenberg interactions—tiny, unavoidable magnetic forces between the atoms (spins) that make up the cloth. The researchers in this paper wanted to answer a crucial question: How much wind can this magical tablecloth take before the knots unravel and the memory is lost?

Here is a breakdown of their findings using simple analogies:

1. The Setup: The Perfect Knot vs. The Wind

• The Toric Code (The Tablecloth): Imagine a grid of tiny magnets. In the perfect "Toric Code" state, these magnets are arranged in a specific pattern where they don't care about their immediate neighbors; they care about the whole loop. It's like a dance where everyone is holding hands in a giant circle. As long as the circle is intact, the dance is stable.
• The Heisenberg Perturbation (The Wind): The researchers introduced a "wind" that pushes every magnet in all directions at once (up, down, left, right). This is the isotropic antiferromagnetic Heisenberg interaction. It's a messy, chaotic force that tries to scramble the neat dance.

2. The Investigation: Watching the Knots

The team used two main tools to study what happens when the wind blows:

• The "Schrieffer-Wolff" (SW) Map: Think of this as a theoretical map or a set of rules. It allows physicists to predict what happens when the wind is light. They found that for weak winds, the magic tablecloth is surprisingly tough. The wind just slightly "renormalizes" (tweaks) the local rules of the dance, but the big, global knots remain intact. The topological order is robust.
• Neural-Network Quantum States (NQS): This is their "super-computer eye." Since the math gets too hard to solve by hand once the wind gets strong, they used an AI (a neural network) trained to act like a quantum physicist. This AI learned to simulate the behavior of millions of magnets at once, allowing them to watch the system evolve as the wind got stronger.

3. The Turning Point: The Critical Wind Speed

As they cranked up the "wind" (the coupling strength $J$), they looked for the moment the tablecloth would snap.

• The Fidelity Susceptibility: Imagine measuring how much the dance changes when you add a tiny bit more wind. If the dance is stable, a little wind causes a little change. But right before the knots unravel, the dance becomes hypersensitive. The researchers found a "peak" in this sensitivity at a specific wind speed.
• The Result: They calculated that the topological order breaks down when the wind strength reaches about 0.164. It's like finding the exact speed of a hurricane that will finally tear the roof off a house.

4. What Happens After the Breakdown?

Once the wind passes that critical speed, the magical tablecloth doesn't just fall apart into chaos; it transforms into something else entirely.

• The N´eel Phase: The system settles into a new, ordered state called a N´eel phase.
• The Analogy: Imagine the dancers, who were previously holding hands in a giant, invisible circle (topological order), suddenly decide to stop dancing together and instead form a rigid, alternating pattern: "Up, Down, Up, Down." They are now locked in a local battle with their immediate neighbors.
• The "Cat State": Interestingly, because of the symmetry of the system, the new state is a "superposition" of four different patterns (like a cat that is simultaneously alive, dead, sleeping, and awake). It's a four-way tie between different magnetic arrangements. The researchers used a tool called staggered magnetization (checking if the "Up, Down" pattern is strong) to confirm this new phase had arrived.

5. Why This Matters

This paper is a victory for both theory and practice:

• For Theory: They proved that you can use a mathematical "map" (SW transformation) to understand the early stages of the breakdown, and then switch to an "AI eye" (NQS) to see the full picture. They showed that the topological protection is incredibly strong against local noise, only failing when the noise is strong enough to create a "loop" of errors that spans the entire system.
• For Technology: Real-world quantum computers (like those from Google or IBM) are trying to build these Toric Codes to store data. This study tells engineers: "Be careful! Even small, unavoidable magnetic interactions between your qubits can eventually destroy your memory if they get too strong. But as long as you keep them below this critical threshold, your topological memory is safe."

In a nutshell: The researchers took a perfect, knot-based quantum memory, blew a chaotic wind on it, and used AI and math to find exactly when the knots would untie. They discovered that the memory is very tough, but once the wind gets too strong, the system snaps into a rigid, alternating pattern, losing its magical topological protection.

Technical Summary

Here is a detailed technical summary of the paper "The toric code under antiferromagnetic isotropic Heisenberg interactions" by Won Jang, Robert Peters, and Thore Posske.

1. Problem Statement

The Toric Code is a paradigmatic model of topological order, characterized by long-range entanglement, ground state degeneracy dependent on manifold topology, and robustness against local perturbations. While its stability under uniform magnetic fields (Ising-type perturbations) is well-studied, the impact of isotropic antiferromagnetic (AFM) Heisenberg interactions remains less understood.

In realistic physical implementations (e.g., superconducting qubits or Mott insulators), residual two-body Heisenberg couplings inevitably arise due to crosstalk or imperfect fine-tuning. Unlike magnetic fields which selectively mobilize specific anyon species, the isotropic Heisenberg interaction couples all spin components ($\sigma^x, \sigma^y, \sigma^z$) simultaneously. This creates a complex competition between local stabilizer renormalization and non-local multi-spin loop processes. The central question is: How does this specific perturbation drive the breakdown of topological order, and what is the nature of the resulting phase?

2. Methodology

The authors employ a dual approach combining analytical perturbation theory and variational machine learning:

• Neural-Network Quantum States (NQS):

  • They use a variational ansatz based on Convolutional Neural Networks (CNNs) to approximate the ground state wavefunction.
  • Symmetry Adaptation: The network explicitly enforces translational symmetry (via weight sharing), $C_4$ rotational symmetry, and global parity constraints.
  • Marshall Gauge Transformation: To handle the sign problem inherent in the AFM Heisenberg model, they apply a Marshall sign transformation. This renders the Hamiltonian stoquastic (non-positive off-diagonal elements in the $\sigma^z$ basis), allowing the use of a real-valued, non-negative NQS ansatz, which significantly improves optimization stability.
  • Training: The network is optimized using Variational Monte Carlo (VMC) with Stochastic Reconfiguration (SR) as a preconditioner.

• Schrieffer-Wolff Transformation (SWT):

  • They derive an effective low-energy Hamiltonian by block-diagonalizing the full Hamiltonian into the toric-code ground state subspace.
  • This analytical framework allows them to track how the perturbation renormalizes local operators (stabilizers) and dresses non-local observables (Wilson loops) order-by-order in the coupling strength $J$.

• Diagnostics:

  • Fidelity Susceptibility ($\chi_F$): To detect the quantum phase transition and extract critical exponents via finite-size scaling.
  • Non-contractible Wilson Loops: Used to probe topological order. Specifically, they analyze the expectation values of single and double-winding loops.
  • Topological Entanglement Entropy (TEE): Calculated via the Kitaev-Preskill (KP) construction using the second Rényi entropy to quantify the loss of long-range entanglement.
  • Staggered Magnetization & Binder Cumulants: Used to characterize the symmetry-broken phase beyond the transition.

3. Key Contributions

A. Analytical Insights via SWT

• Renormalization vs. Mixing: The authors prove that at low orders, the Heisenberg perturbation only renormalizes local stabilizers ($A_v, B_p$) and non-contractible Wilson loops multiplicatively.
• Topological Sector Mixing: Crucially, they demonstrate that mixing between topological sectors (which breaks the degeneracy) only occurs at a perturbative order proportional to the system size ($L$). This confirms the exponential robustness of the topological phase against local perturbations.
• Parity Dependence: They derive distinct effective Hamiltonians for even and odd system sizes ($L$).
  • Odd $L$: Global parity constraints forbid single-loop terms, leading to a ground state dominated by correlated loop products (Bell-type states).
  • Even $L$: Single-loop terms are allowed, leading to a ground state polarized along the bisector of the logical $X$ and $Z$ axes.

B. Numerical Discovery of the Phase Transition

• Critical Point: Using finite-size scaling of the fidelity susceptibility and the logarithmic susceptibility of Wilson loops, they locate the critical coupling at $J_c \approx 0.164$.
• Critical Exponent: The scaling analysis yields a correlation-length exponent $\nu \approx 0.7$ (derived from $2/\nu \approx 2.85$), consistent with a continuous phase transition.
• Breakdown Mechanism: The transition is marked by the rapid decay of non-contractible Wilson loop expectation values and the deviation of the Topological Entanglement Entropy from its quantized value ($\gamma = -\ln 2$).

C. Characterization of the New Phase

• Emergent Phase: Beyond $J_c$, the system transitions into an antiferromagnetic Néel phase.
• Symmetry Breaking: This phase is characterized by a fourfold-degenerate symmetry-broken manifold.
• X/Z Néel Order: Unlike standard Heisenberg antiferromagnets which might exhibit SU(2) symmetry, the interplay with the toric code structure and duality leads to an easy-plane (X/Z) Néel order. The $Y$-component magnetization is strongly suppressed, while $X$ and $Z$ components develop order.
• Finite-Size Effects: In finite systems, the ground state remains a symmetric superposition (cat state) of the four degenerate sectors, resulting in vanishing linear magnetization but non-zero even moments (structure factors).

4. Key Results

1. Agreement in Weak Coupling: For $J \lesssim 0.1$, the NQS results for ground state energy and observables match the second-order SWT predictions perfectly, validating the perturbative breakdown of topological order.
2. Transition Signatures:
   • Fidelity Susceptibility: Shows a sharp peak at $J_c$ that grows with system size $L$.
   • Wilson Loops: The expectation values of non-contractible loops drop from $\sim 1$ (or $1/\sqrt{2}$depending on parity) to near zero as$J$ increases, signaling the loss of topological protection.
   • TEE: The topological entanglement entropy $\gamma_{KP}$ rises from $-\ln 2$ toward zero, confirming the transition to a trivial phase.
3. Nature of the Ordered Phase:
   • The staggered magnetization structure factors $S_x$ and $S_z$ grow significantly, while $S_y$ remains suppressed.
   • Binder cumulants converge to $1/3$for large$L$, identifying the order as a four-sector X/Z cat state in finite systems, which spontaneously breaks into a single X or Z Néel sector in the thermodynamic limit.

5. Significance

• Realistic Robustness: The study provides a rigorous assessment of the toric code's stability against the most physically relevant perturbation (isotropic Heisenberg exchange) found in real quantum hardware and condensed matter systems.
• Methodological Framework: It establishes a robust pipeline combining Schrieffer-Wolff transformations (for analytical control) and symmetry-adapted NQS (for non-perturbative numerical accuracy). This framework is applicable to a broad class of topological models and perturbations.
• Phase Transition Physics: It elucidates how local two-spin interactions can drive a topological phase transition, revealing a specific "easy-plane" Néel order that respects the underlying duality of the toric code.
• Finite-Size vs. Thermodynamic Limit: The work clarifies the subtle distinction between finite-size ground states (symmetric superpositions/cat states) and the thermodynamic limit (spontaneous symmetry breaking), a crucial consideration for interpreting numerical simulations of topological phases.

In summary, this paper demonstrates that while the toric code is robust against weak Heisenberg perturbations, it undergoes a continuous phase transition at a well-defined critical point into a magnetically ordered phase, providing a comprehensive map of the stability landscape for topological quantum memories.