Grand Canonical-like Thermalization of Quantum Many-body Scars

Here is an explanation of the paper "Grand Canonical-like Thermalization of Quantum Many-body Scars" using simple language and creative analogies.

The Big Picture: When Quantum Systems Break the Rules

Imagine you have a crowded dance floor (a quantum system). Usually, when you drop a beat (add energy), everyone starts dancing wildly and eventually, the whole floor settles into a uniform, chaotic groove. This is called thermalization. In physics, we have a rulebook for this called the Eigenstate Thermalization Hypothesis (ETH). It says: "If you look at any single dancer, they will eventually look just like the average of the whole crowd."

But, in some special quantum systems, there are a few dancers who refuse to blend in. They keep doing the exact same synchronized routine over and over again, ignoring the chaos around them. These are called Quantum Many-Body Scars (QMBS). They are the "rebellious teenagers" of the quantum world who break the rulebook.

For a long time, physicists were confused: Why do these rebels exist? How do we predict what they will do?

This paper by Wang and colleagues provides a new rulebook to explain them.

The New Perspective: The "Two-Variable" Thermostat

1. The Old Way vs. The New Way

The Old Way (Standard ETH):
Imagine trying to predict the temperature of a room. The old rulebook said you only need to know one thing: the total energy (how hot the room is). If you know the energy, you know everything about the room's behavior.

The New Way (This Paper):
The authors realized that in these "constrained" systems (where dancers can't move next to each other), energy isn't enough. You need a second piece of information: The Quasiparticle Number.

Think of it like a Grand Canonical Ensemble (a fancy physics term).

Energy ( $E$ ): How much "heat" or activity is in the system.
Quasiparticle Number ( $N$ ): Think of this as the "number of active dancers" or "information flow." In these systems, there are rules preventing certain moves (like two people standing too close). The "quasiparticle number" counts how many of these "forbidden moves" are being attempted or avoided.

The Analogy:
Imagine a busy coffee shop.

Standard ETH: You predict the noise level based only on the number of people (Energy).
This Paper: You realize the noise level also depends on how many people are talking loudly (Quasiparticle Number). Even if the shop is full (high energy), if everyone is whispering (low quasiparticle number), it's quiet. If a few people are shouting (high quasiparticle number), it's loud. You need both numbers to predict the atmosphere.

2. The "Open System" Trick

The authors used a clever trick to understand this. Instead of looking at the system as a closed box, they imagined it as an open system connected to a "ghost environment."

The Constraint: Imagine a bouncer at the door who stops people from entering a specific VIP area.
The Trick: Instead of just saying "No entry," the authors imagined the bouncer is actually swapping information with a ghost outside. Sometimes the ghost gives info back (negative dissipation).
The Result: This "ghost exchange" is mathematically equivalent to the bouncer's rule. It turns out that the "Quasiparticle Number" is actually a measure of how much information is being exchanged with this ghost.

Solving the Mystery of the "Scars"

Why do the "Scar" states (the rebellious dancers) behave so strangely?

The "Density of States" Map

In the old rulebook, physicists looked at a map of Energy to see how crowded the system was. They thought the rebels lived in a crowded area.

The authors drew a new map: The Energy vs. Quasiparticle Number Plane.

The Thermal States: These are like a dense forest. There are millions of trees (states) packed tightly together. If you pick a random spot, you are surrounded.
The Scar States: These are like a desert. They live in a region where there are almost no other states nearby.

The Analogy:
Imagine you are walking in a city.

Thermal states are in the downtown district. It's so crowded that if you drop a coin, it hits someone immediately. The system "thermalizes" (settles down) quickly because there are so many places to go.
Scar states are in a remote, empty valley. There are very few other people (states) around. Because the valley is so empty (Low Density of States), if you drop a coin, it takes a long time to hit anyone. The system gets "stuck" in a loop, creating those long, rhythmic oscillations we see in scars.

The "Cross Coherence Purity" (CCP)

To prove this, the authors invented a new tool called Cross Coherence Purity (CCP).

Think of CCP as a "fuzziness meter." It measures how much two different quantum states "blur" into each other.
The Discovery: They found that in the empty desert (low density of states), the "fuzziness" is high. The states are distinct and don't mix well. This lack of mixing is exactly why the scars don't thermalize—they stay in their own little loop.

The "Spectrum-Generating Algebra" (The Magic Ladder)

Physicists had noticed that scar states form a perfect ladder: State A, State B, State C, all spaced by the exact same amount of energy. It looked like magic.

The authors showed this isn't magic; it's geometry.
Because the scar states live in that empty desert (low density of states), there are very few "neighbors" to mess up their spacing.

Thermal states: In the crowded city, if you try to build a ladder, people bump into you, and the rungs get uneven.
Scar states: In the empty desert, you can build a perfect, straight ladder because no one is there to knock it over.

The "algebra" (the math rule that makes the ladder work) emerges naturally because the system is so sparse in that specific region.

Summary: What Does This Mean for Us?

Thermalization is more complex: We can't just look at energy. We have to look at "information flow" (quasiparticles) too.
Scars aren't broken; they are just lonely: The strange, non-thermal behavior of scars happens because they live in a "desert" of quantum states where there is very little competition.
A Unified Theory: This paper unifies the "normal" thermal behavior and the "weird" scar behavior into one single framework. It says: "If you are in a crowded place, you thermalize. If you are in an empty place, you keep oscillating."

The Takeaway:
The universe isn't just a chaotic mess. Even in the most complex quantum systems, there are hidden patterns. By looking at the system through a new lens (Energy + Quasiparticles), the authors showed that the "rebels" (scars) are actually following a very logical, albeit sparse, set of rules. They aren't breaking the laws of physics; they are just living in a very quiet neighborhood.

Here is a detailed technical summary of the paper "Grand Canonical-like Thermalization of Quantum Many-body Scars" by Wang et al.

1. Problem Statement

The Eigenstate Thermalization Hypothesis (ETH) is the standard framework explaining how isolated quantum many-body systems reach thermal equilibrium. It posits that:

Diagonal elements: Expectation values of local observables in energy eigenstates depend smoothly on energy ( $E$ ) and match the microcanonical average.
Off-diagonal elements: Matrix elements between different eigenstates are exponentially small, scaling as $\Omega(E)^{-1/2}$ (where $\Omega(E)$ is the density of states), leading to vanishing temporal fluctuations in the thermodynamic limit.

Quantum Many-Body Scars (QMBS) violate this hypothesis. In kinetically constrained systems (e.g., the PXP model), specific "scar" eigenstates exhibit:

Large deviations from thermal expectation values.
Anomalous, persistent temporal oscillations (quasi-periodic revivals) rather than thermalization.
A failure of the standard ETH to predict long-time averages or fluctuation magnitudes.

The core problem addressed is: How can we formulate a unified thermalization theory that accounts for both generic thermal states and the anomalous behavior of QMBS in kinetically constrained systems?

2. Methodology

The authors employ a multi-faceted approach combining open-system theory, statistical mechanics, and numerical simulation:

Open-System Mapping: They reinterpret kinetically constrained dynamics (which restrict the Hilbert space) as an effective open system. By coupling the system to an auxiliary environment via specific dissipation channels, they show that the constrained dynamics are equivalent to unitary evolution under a non-Hermitian Hamiltonian ( $H_N$ ) within the constrained subspace.
Quasiparticle Number Definition: Based on the open-system picture, they define a quasiparticle number operator ( $\hat{N}$ ) that quantifies the rate of information exchange between the system and the environment. This operator counts specific local patterns (e.g., $|y\rangle_{k,k+1}$ ) forbidden by the kinetic constraints.
Revised ETH Formulation:
- Diagonal: They propose that eigenstates are equilibrium states of a Grand Canonical-like ensemble, where the state is characterized by both Energy ( $E$ ) and Quasiparticle Number ( $N$ ). Thus, local observables depend on $O(E, N)$ rather than just $O(E)$ .
- Off-Diagonal: They generalize the Density of States (DOS) to a two-variable function $\Omega(E, N)$ defined on the energy-quasiparticle-number plane.
Cross Coherence Purity (CCP): To quantify off-diagonal matrix elements independent of specific observables, they introduce the CCP ( $T_{i,i'}$ ), defined as the purity of the reduced cross-density matrix between two eigenstates. They hypothesize $T_{i,i'} \propto \Omega(E, N)^{-1}$ .
Spectrum-Generating Algebra (SGA) Analysis: They investigate the algebraic structure (approximate ladder operators) responsible for the quasi-periodic dynamics of scars, analyzing its emergence in the context of the revised ETH.

3. Key Contributions

A. Grand Canonical-like Thermalization Framework

The paper establishes that kinetically constrained systems thermalize not to a canonical ensemble (fixed $E$ ) but to a grand canonical-like ensemble (fixed $E$ and $N$ ).

Result: The long-time average of any local observable $\bar{O}$ for an initial state is given by the ensemble average $\langle \hat{O} \rangle_{\beta, \mu}$ , where $\beta$ (inverse temperature) and $\mu$ (chemical potential) are fixed by the initial energy and quasiparticle number.
Significance: This explains why scar states (which have low $N$ ) deviate from canonical predictions but align perfectly with the grand canonical prediction.

B. Generalized Density of States and Fluctuations

The authors demonstrate that the suppression of off-diagonal matrix elements is governed by the DOS on the $(E, N)$ plane, $\Omega(E, N)$ , rather than just $\Omega(E)$ .

Finding: Scar states reside in regions of anomalously low DOS on the $(E, N)$ plane.
Mechanism: Because $\Omega(E, N)$ is small in these regions, the off-diagonal matrix elements (and thus temporal fluctuations) are enhanced (scaling as $\Omega(E, N)^{-1/2}$ ), preventing thermalization.
Verification: Numerical simulations show an inverse relationship between the Cross Coherence Purity (CCP) and the generalized DOS $\Omega(E, N)$ .

C. Emergent Spectrum-Generating Algebra (SGA)

The paper resolves the mystery of why scar states form an approximate algebraic structure (SGA) leading to quasi-periodic dynamics.

Insight: The SGA is not an independent anomaly but a natural consequence of the low-DOS regime.
Derivation: The condition for an approximate SGA (small deviation from exact ladder operators) depends on the denominator of a ratio involving off-diagonal matrix elements. Since these matrix elements scale with $\Omega(E, N)^{-1/2}$ , the low DOS in the scar region naturally suppresses the "noise" terms, allowing the algebraic structure to emerge.

4. Results

Unified Description of EEVs: Numerical data (Fig. 2) confirms that Eigenstate Expectation Values (EEVs) of local observables cannot be collapsed onto a single curve of Energy ( $E$ ) alone. However, they form a smooth surface when plotted against both $E$ and $N$ .
Accuracy of Predictions: The grand canonical-like ensemble accurately predicts the long-time averages and temporal fluctuations for both thermal states and scar states. In contrast, the standard canonical ensemble fails significantly for scar states.
CCP-DOS Relation: The inverse of the Cross Coherence Purity ($1/T_{grid} $) is linearly proportional to the generalized DOS$ \Omega(E, N)$, validating the revised off-diagonal ETH ansatz.
Origin of Oscillations: The large temporal fluctuations observed in scar states are directly attributed to the sparse distribution of eigenstates (low DOS) in the $(E, N)$ plane.
SGA Validation: The study shows that the approximate SGA observed in scars is a direct result of the low DOS, providing a thermodynamic explanation for the algebraic structure previously thought to be purely kinematic.

5. Significance

Theoretical Unification: This work bridges the gap between standard thermalization (ETH) and non-thermal phenomena (QMBS). It suggests that QMBS are not "exceptions" to thermalization but rather represent a specific phase of matter where thermalization is governed by an additional conserved quantity (quasiparticle number) and a sparse spectrum.
New Perspective on Constraints: By mapping kinetic constraints to dissipative open-system dynamics, the authors provide a powerful thermodynamic lens to analyze constrained quantum systems.
Predictive Power: The revised ETH allows for accurate prediction of long-time dynamics and fluctuations in systems where the standard ETH fails, offering a tool for understanding non-equilibrium dynamics in Rydberg atom arrays and other constrained platforms.
Fundamental Insight: It reframes the "violation" of ETH in scars as a consequence of the geometry of the spectrum in a generalized phase space ( $E, N$ ), rather than a breakdown of statistical mechanics principles.