Prethermal gauge structure and surface growth in $\mathbb{Z}_2$ lattice gauge theories

This paper numerically demonstrates that experimentally feasible two-body Ising interactions in $(2+1)$D Rydberg atom arrays can stabilize a prethermal $\mathbb{Z}_2$ lattice gauge structure, which eventually breaks down via defect proliferation to reveal thermalization dynamics governed by the Kardar-Parisi-Zhang universality class.

The Gist

Imagine you are trying to keep a massive, chaotic dance floor perfectly organized. You have thousands of dancers (atoms), and you want them to follow a very strict set of rules: "If you are in this spot, your neighbors must be in that spot." In physics, we call these rules gauge symmetries. They are the "laws of the universe" for this little system.

The problem is, in the real world, things get messy. If you nudge the dancers, they might break the rules. Usually, once they break a rule, the whole dance floor collapses into chaos (thermalization) very quickly.

This paper is about a team of physicists who discovered a way to keep the dance floor organized for a surprisingly long time, even when the music is playing loudly. They also discovered that when the rules do eventually break, they don't just collapse randomly; they break in a very specific, beautiful pattern, like a bubble expanding in a liquid.

Here is the story of their discovery, broken down into simple parts:

1. The "Force Field" Shield (Prethermalization)

The researchers built a model using a honeycomb grid (like a beehive) where every intersection has a dancer. They wanted to enforce a rule called Gauss's Law: "The number of dancers in a specific pattern must always be even."

To keep this rule, they created a "force field" (a high energy cost). Imagine that if a dancer breaks the rule, they have to pay a huge fine (energy penalty).

• The Setup: They start with everyone following the rules perfectly.
• The Disturbance: They turn on a "nudge" (a laser or magnetic field) that tries to make the dancers spin and break the rules.
• The Result: Because the "fine" for breaking the rule is so high, the dancers resist. They enter a Prethermal Phase. This is like a "time-out" period where the system acts as if the rules are still perfect, even though the nudge is trying to break them. It's like a heavy door that is hard to push open; it stays shut for a long time even if you push it.

2. The "Bubble" Breakdown

Eventually, the nudge wins. The rules break. But here is the cool part: it doesn't happen everywhere at once.

• Imagine a drop of water on a hot pan. It doesn't instantly turn into steam; a little bubble forms, then grows.
• In this system, a single spot breaks the rule (a "defect"). This defect acts like a seed.
• Neighboring spots see this and start breaking the rules too. The "bad behavior" spreads out like a growing bubble or a stain spreading on a shirt.
• The researchers found that this spreading isn't random. It follows a specific mathematical pattern known as KPZ universality. Think of it like sand piling up on a beach. The pile gets rougher and rougher in a very predictable way. This pattern had never been seen in this specific type of quantum system before.

3. The "Crystal Ball" Problem (Why Computers Failed)

The team tried to predict this behavior using different types of computer simulations:

• The "Perfect" Computer (Exact Diagonalization): This is like simulating every single dancer's exact move. It works, but it's so slow that you can only simulate a tiny dance floor (like 4x4 dancers).
• The "Gambler" Computer (DTWA): This method tries to guess the future by adding random noise to the simulation. Usually, this works great for big systems. But it failed here. It couldn't see the "time-out" period. It thought the rules would break immediately.
• The "Captain" Computer (Mean-Field): This method treats the dancers as a single, smooth fluid rather than individuals. Surprisingly, this simple method worked better than the "Gambler" computer!

Why did the "Gambler" fail? The paper suggests that the "Gambler" method accidentally introduced tiny errors that acted like seeds for the rule-breaking. Because the system is so sensitive to these seeds (due to the hidden rules), the simulation collapsed too fast. It's like trying to balance a house of cards while someone is sneezing nearby; the "Gambler" method sneezed, and the house fell.

4. Why This Matters

This isn't just a theoretical game. The researchers showed that this system can be built using Rydberg atoms (super-excited atoms) in a lab.

• For Quantum Computers: Building a quantum computer is hard because these machines are very fragile and break rules easily. This paper shows a way to build a "shield" that protects the rules for a long time, giving us more time to do calculations before the system breaks down.
• For Physics: It reveals a hidden layer of order in how things break down. Even when a system is falling apart, it follows a universal "recipe" (the KPZ pattern) that connects it to other seemingly unrelated things, like how snowflakes grow or how traffic jams spread.

The Big Takeaway

The universe has a way of holding onto order for a long time, even when you try to shake it up. And when that order finally breaks, it doesn't just shatter; it spreads like a bubble, following a beautiful, predictable rhythm that we can now see and measure. This gives scientists a new tool to build better quantum simulators and understand the deep, hidden laws of how matter behaves.

Technical Summary

1. Problem Statement

The paper addresses two central challenges in non-equilibrium many-body physics:

1. Thermalization in Lattice Gauge Theories (LGTs): Understanding how isolated interacting systems reach thermal equilibrium is difficult, particularly in kinetically constrained models like LGTs where local constraints (Gauss' laws) restrict the Hilbert space.
2. Numerical Limitations: Studying the thermalization of (2+1)D LGTs beyond small system sizes is notoriously difficult due to the exponential scaling of the Hilbert space, limiting exact diagonalization (ED) to small systems and making long-time dynamics inaccessible.
3. Gauge Symmetry Protection: In experimental quantum simulators (e.g., Rydberg atoms), gauge symmetries are often approximate. The authors seek to understand how to stabilize a gauge-invariant subspace against symmetry-breaking perturbations and how the system eventually thermalizes when that protection fails.

2. Methodology

The authors employ a multi-faceted approach combining classical mean-field simulations, semi-classical approximations, and exact quantum benchmarks:

• Model System: A (2+1)D spin system on a honeycomb lattice representing a $Z_2$ LGT with dynamical matter.
  • Spins: Located on sites (matter) and links (electric fields).
  • Hamiltonian ($\hat{H} = \hat{H}_{Z2} + \hat{H}_V + \hat{H}_\Omega$):
    • $\hat{H}_{Z2}$: Gauge-invariant coupling between matter and fields.
    • $\hat{H}_V$: A local "pseudogenerator" term composed of two-body interactions that energetically penalizes violations of Gauss' law (creating an energy gap $\propto V$).
    • $\hat{H}_\Omega$: A transverse field drive that breaks the gauge symmetry.
• Numerical Techniques:
  • Mean-Field Dynamics: The primary tool for simulating large systems ($N \approx 2000$ spins, $20 \times 20$ plaquettes). Classical spins evolve via Poisson brackets ($\dot{\vec{\sigma}} = \{\vec{\sigma}, H\}$), allowing for efficient computation of long-time dynamics.
  • Discrete Time Wigner Approximation (DTWA): A semi-classical method used to benchmark small systems and test the role of quantum fluctuations.
  • Exact Diagonalization (ED): Used for small systems ($2 \times 2$ plaquettes) to provide the ground truth for quantum dynamics.
• Analysis:
  • Quantification of gauge symmetry breaking via the time-averaged error $\epsilon(t)$.
  • Analysis of the breakdown of the prethermal phase using surface growth models, specifically looking at the proliferation of Gauss' law defects.
  • Scaling analysis to identify universality classes (Family-Vicsek scaling).

3. Key Contributions

1. Stabilization of Prethermal $Z_2$ Gauge Structure: The authors demonstrate that for strong protection ($V \gg \Omega$), the system enters a long-lived prethermal plateau where the gauge symmetry is effectively restored, despite the presence of a symmetry-breaking drive.
2. Mechanism of Breakdown: They identify that the eventual breakdown of this metastable state is driven by the nucleation and proliferation of Gauss' law defects, analogous to bubble formation in false vacuum decay.
3. Discovery of KPZ Universality: A major finding is that the spatio-temporal correlations of the defect proliferation follow the Kardar-Parisi-Zhang (KPZ) universality class (specifically (1+1)D KPZ), characterized by non-linear surface growth exponents.
4. Benchmarking Approximations: The study reveals a critical failure of the DTWA method in this context. While mean-field dynamics qualitatively matches ED, DTWA fails to capture the prethermal plateau, highlighting that thermalization in LGTs is governed by emergent local constraints rather than simple scrambling of fluctuations.

4. Key Results

A. Prethermal Plateau and Timescales

• Dynamics: Starting from a gauge-invariant state, the system exhibits three regimes:
  1. Early time: Universal quadratic growth of error ($\epsilon(t) \propto t^2$).
  2. Prethermal Plateau: The error stabilizes at a low value $\epsilon_{pre} \propto (\Omega/V)^2$. The duration of this plateau scales linearly with the protection strength ($t_{crit} \propto V/\Omega$).
  3. Late Time: The gauge symmetry breaks down completely, and the system thermalizes to an infinite-temperature state ($\epsilon \to 1$).
• Critical Protection: Analytic derivation of a translationally invariant model without matter reveals a phase transition at a critical protection strength $(V/\Omega)_c \approx 3.33$. Below this, the system oscillates wildly; above it, a stable potential well forms, constraining the dynamics near $G=+1$.

B. Surface Growth and KPZ Scaling

• Defect Proliferation: Upon exiting the prethermal plateau, random fluctuations nucleate regions of gauge violation. These regions grow and merge, creating a "surface" separating defect-free and defective regions.
• Scaling Exponents: The authors analyze the height function $h(X,t)$ of this surface.
  • Ballistic Spreading: The mean height grows as $h(t) \propto t^{1.01}$.
  • Roughening: The surface width grows as $\delta h(t) \propto t^\beta$ with $\beta \approx 1/3$.
  • Family-Vicsek Scaling: A dynamical scaling analysis confirms the exponents $\alpha = 1/2$ (roughness) and $z = 3/2$ (dynamic), which are the defining characteristics of the (1+1)D KPZ universality class. This reveals a hidden feature in the thermalization of multi-point correlators.

C. Benchmarking and Quantum Relevance

• Mean-Field vs. ED: In small systems, the classical mean-field dynamics shows excellent qualitative agreement with Exact Diagonalization (ED), validating the use of mean-field theory for large-scale predictions in this regime.
• Failure of DTWA: DTWA fails to reproduce the prethermal plateau, predicting immediate thermalization. The authors attribute this to the stochastic noise in DTWA seeding gauge-symmetry-breaking sectors that are energetically suppressed in the true quantum evolution. This underscores that LGT thermalization is not merely a scrambling of fluctuations but is heavily constrained by emergent local symmetries.
• Experimental Feasibility: The model is directly implementable in Rydberg atom arrays. The authors propose a specific regime (facilitation-like regime with detuning $\Delta = V$) where the system naturally realizes the required $Z_2$ Gauss' law constraints without dynamical matter, stabilizing even/odd excitation sectors.

5. Significance

• Quantum Simulation: The work provides a concrete roadmap for simulating (2+1)D lattice gauge theories on current and near-future quantum simulators (e.g., Rydberg tweezers with thousands of qubits). It demonstrates that gauge constraints can be stabilized for experimentally relevant timescales.
• Universality in Non-Equilibrium Physics: By linking the breakdown of gauge symmetry to the KPZ universality class, the paper bridges the gap between high-energy physics concepts (gauge theories, false vacuum decay) and statistical mechanics (surface growth, non-equilibrium criticality).
• Methodological Insight: The failure of DTWA serves as a crucial warning for the community: semi-classical methods that rely on stochastic sampling of phase space may fail to capture the subtle, emergent constraints of gauge theories, necessitating the use of quantum simulators or more sophisticated numerical methods for accurate long-time predictions.
• Thermalization Theory: It offers a new perspective on thermalization in kinetically constrained systems, showing that the pathway to equilibrium is governed by the nucleation and growth of topological defects rather than simple ergodicity.