🔢 mathematics

Breakdown of the thermodynamic limit in quantum spin and dimer models

This paper demonstrates that the thermodynamic limit can break down in quantum spin and dimer models by constructing Hamiltonians on square and square-octagon lattices where the ground state phases on diamond-shaped domains differ fundamentally from those on square domains, exhibiting geometry-dependent macroscopic regions with distinct ordering and correlation behaviors.

Original authors: Jeet Shah, Laura Shou, Jeremy Shuler, Victor Galitski

Published 2026-06-23

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jeet Shah, Laura Shou, Jeremy Shuler, Victor Galitski

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Idea: Shape Matters More Than You Think

Imagine you have a giant block of iron. If you cut it into a perfect cube or a perfect sphere, the inside of the metal feels exactly the same in both shapes. It's solid, metallic, and uniform. In physics, there is a golden rule called the Thermodynamic Limit. It says that if you make a system (like a collection of atoms) big enough, the shape of its edges shouldn't matter. The "bulk" (the middle part) should always behave the same way, regardless of whether the container is a square, a circle, or a diamond.

This paper says: "Not always."

The authors discovered specific quantum systems where the shape of the boundary completely changes what happens in the middle. If you take the exact same set of rules and put them in a square box, you get one type of behavior. If you put them in a diamond-shaped box, the middle splits into different zones with totally different behaviors. The "bulk" is no longer uniform; it depends entirely on the shape of the container.

The Characters: Quantum Spin and Dimers

To understand how they did this, we need to meet the "actors" in their play:

The Spins: Imagine tiny magnets on a grid. They can point up or down.
The Dimers: Think of these as dominoes. A dimer is a pair of neighboring spots on the grid that are "occupied" together. The rule is strict: No two dominoes can share a spot. Every spot must be covered by exactly one domino.
The Quantum Twist: In the real world, dominoes sit still. In this quantum world, the dominoes can "flip" and rearrange themselves instantly. The system exists in a superposition of all possible ways to cover the grid with dominoes at the same time.

Example 1: The Aztec Diamond (The "Frozen Corners")

The authors first looked at a grid shaped like a Diamond (specifically, an "Aztec Diamond").

The Square Box: If you play this game on a standard square grid, the dominoes are chaotic and jumbled everywhere. There is no order; it's like a liquid where everything is moving randomly.
The Diamond Box: When they played the exact same game on a diamond-shaped grid, something magical happened.
- The Corners: The four corners of the diamond became frozen. The dominoes there locked into a perfect, rigid pattern (like a crystal). They stopped moving entirely.
- The Center: The middle of the diamond remained liquid and chaotic, just like the square grid.
- The Arctic Circle: Separating the frozen corners from the liquid center is a perfect circle (called the "Arctic Circle").

The Analogy: Imagine a room full of people dancing wildly (the liquid center). But if the room is shaped like a diamond, the people in the four corners suddenly freeze in place, holding hands in a perfect line, while the people in the middle keep dancing. The shape of the room forced the corners to freeze, even though the music (the rules of the game) was the same everywhere.

Example 2: The Square-Octagon Fortress (The "Three-Zone City")

The authors then tried a more complex grid made of squares and octagons (like a stop sign surrounded by squares). They built this in a diamond shape, which they call a "Fortress."

Here, the breakdown of the "Thermodynamic Limit" is even more dramatic. The diamond fortress splits into three distinct zones:

The Frozen Corners: Just like the Aztec diamond, the corners are rigid and ordered.
The Critical Ring: A ring around the center where the dominoes are chaotic and "critical" (a state of high sensitivity and fluctuation).
The Gaseous Center: The very middle is a different kind of ordered state, where the dominoes are arranged in a specific, short-range pattern that is neither frozen nor chaotic.

The Analogy: Imagine a city built in a diamond shape.

The outskirts are a quiet, frozen suburb where everyone follows a strict routine.
The inner ring is a bustling, chaotic festival where anything can happen.
The downtown core is a highly organized, efficient business district with a different kind of order.
All these zones exist simultaneously in the same city, defined only by the diamond shape of the city limits.

How They Proved It (The "Mathematical Microscope")

Usually, to study these quantum systems, scientists use Monte Carlo simulations. This is like running a computer simulation where you randomly pick domino arrangements millions of times to see what usually happens. It's like trying to guess the weather by flipping a coin a million times. It's slow and often inaccurate for these specific shapes.

The authors invented a better way. They used a mathematical tool called the Kasteleyn Matrix.

The Analogy: Instead of flipping coins to guess the weather, they built a perfect, exact map of the atmosphere. They could calculate the exact probability of every single domino placement without guessing.
They used this to measure "correlators" (how much one domino's position affects another far away).
They found that in the center of the fortress, the dominoes were not in a "Quantum Spin Liquid" (a mysterious, highly entangled state physicists love). Instead, they were in a simple, ordered state. The "magic" of the quantum liquid was broken by the geometry.

The Conclusion

The paper proves that for certain quantum systems, geometry is destiny.

In the standard view of physics, if you make a system big enough, the edges don't matter. This paper shows that for quantum spin and dimer models, the edges dictate the entire story. A diamond shape creates a city with three different neighborhoods (frozen, chaotic, and ordered), while a square shape creates a uniform city.

This means that when we try to understand the "bulk" properties of matter, we can't just ignore the shape of the container. In these specific quantum worlds, the shape of the boundary writes the rules for the entire interior.

Technical Summary: Breakdown of the Thermodynamic Limit in Quantum Spin and Dimer Models

Problem Statement
The thermodynamic limit is a foundational assumption in statistical mechanics, positing that as a system size approaches infinity at fixed particle density, macroscopic bulk phases emerge that are independent of the system's boundary shape. While boundary effects are typically confined to a thin edge region in generic interacting many-body systems, this paper investigates whether this assumption holds for specific quantum spin and dimer models. The authors challenge the standard notion by constructing models where the shape of the boundary induces macroscopic regions with distinct bulk phases, causing local correlators deep in the interior to remain boundary-dependent even in the infinite-size limit.

Methodology
The authors employ a combination of exact analytical methods and high-precision numerical techniques to study quantum dimer models on two specific lattices: the square lattice with diamond-shaped boundaries (the "Aztec diamond") and the square-octagon lattice with diamond-shaped boundaries (the "square-octagon fortress").

Hamiltonian Construction: The authors define quantum spin Hamiltonians on these lattices. In the low-energy limit (strong coupling limit), these spin models map exactly to quantum dimer models. Specifically, they focus on the Rokhsar-Kivelson (RK) point, where the ground state is an equal-weight superposition of all valid dimer coverings (perfect matchings) of the graph.
Kasteleyn Matrix Formalism: To analyze the ground states, the authors utilize the Kasteleyn method. They construct signed adjacency matrices (Kasteleyn matrices) for the planar graphs. The inverse of these matrices allows for the exact calculation of dimer probabilities and correlators via Kenyon's formula.
Vison Correlator Calculation: A key methodological contribution is the development of an exact, efficient algorithm to calculate the vison correlator—a nonlocal string order parameter crucial for diagnosing quantum spin liquid phases.
- Instead of relying on standard Monte Carlo sampling, which suffers from slow convergence and finite sampling errors (typically $\sim 10^{-6}$ ), the authors use fraction-free sparse LU (FFLU) decomposition. This method relies on exact integer arithmetic to compute the required submatrices of the inverse Kasteleyn matrix, achieving precisions up to $10^{-21}$ .
- For the square-octagon lattice with rectangular boundaries, they further utilize the Toeplitz structure of the inverse Kasteleyn matrix to derive analytical asymptotic results for the vison correlator using the Strong Szegő limit theorem.
Correlator Analysis: The authors compute both the dimer-dimer connected correlator (to determine correlation length and criticality) and the vison correlator (to distinguish between deconfined spin liquid phases and confined ordered phases).

Key Contributions and Results

Breakdown of the Thermodynamic Limit (Aztec Diamond):
The authors demonstrate that for the quantum dimer model on the Aztec diamond (square lattice with 45° rotated boundaries), the ground state exhibits quantum phase separation.
- Frozen Corners: Regions near the four corners display a "frozen" phase with staggered dimer order (dimer occupation probability $\to 1$ or $0$).
- Critical Core: The interior region, bounded by the "Arctic circle," exhibits critical behavior with power-law decaying dimer-dimer correlators.
- Significance: This proves that the thermodynamic limit is not well-defined for this system without specifying the boundary shape. The phase separation is macroscopic (occupying $\approx 21\%$ of the area) and not merely a boundary effect.
Square-Octagon Lattice with Rectangular Boundaries:
The authors analyze the square-octagon lattice with standard square boundaries to establish a baseline.
- Gapped Ordered Phase: They find that the dimer-dimer correlator decays exponentially, indicating a gapped phase.
- Vison Correlator: Crucially, the vison correlator approaches a non-zero constant ( $\approx 0.7746$ ) along certain paths. This rules out a $Z_2$ quantum spin liquid phase (where visons would be gapped excitations and the correlator would decay to zero).
- Nature of Order: The system is identified as a short-range entangled ordered phase, characterized by a background of flippable square plaquettes interspersed with a finite density of flippable octagonal defects.
Robust Phase Separation (Square-Octagon Fortress):
The authors extend the study to the square-octagon lattice with diamond-shaped boundaries (the "fortress").
- Three Distinct Regions: The ground state exhibits three macroscopic regions separated by an "octic curve":
  1. Frozen Regions: Near the corners, exhibiting staggered dimer order.
  2. Critical Region: Between the frozen corners and the center, exhibiting power-law correlations.
  3. Gaseous (Central) Region: The innermost region, behaving exactly like the square-octagon lattice with rectangular boundaries (exponentially decaying dimer correlators and constant vison correlator).
- Stability: Unlike the Aztec diamond case where the critical region might be unstable to perturbations away from the RK point, the authors argue that the phase separation in the fortress is likely robust. This is because the frozen and central regions are both gapped (exponentially decaying correlations), suggesting that a critical region must persist between them to mediate the transition, even if the system is perturbed away from the exact RK point.

Significance and Claims
The paper claims to provide the first explicit examples of quantum spin models where the thermodynamic limit breaks down due to boundary shape dependence, resulting in macroscopic coexistence of distinct quantum phases.

Theoretical Implication: The work generalizes the concept of boundary dependence from classical dimer models (where it is well-known) to quantum many-body systems. It demonstrates that the shape of the boundary can dictate the bulk phase structure, challenging the standard assumption that bulk observables become independent of boundary conditions in the thermodynamic limit.
Methodological Advancement: The paper highlights the superiority of Kasteleyn-based exact methods over Monte Carlo simulations for calculating vison correlators in dimer models, offering a tool to diagnose quantum phases with arbitrary precision.
Phase Coexistence: The identification of a stable, three-region phase separation in the square-octagon fortress suggests that such boundary-induced phase coexistence is not restricted to fine-tuned critical points but may persist over extended regions of the phase diagram.

The authors conclude that these findings open new avenues for understanding how geometry influences quantum matter, particularly in the context of quantum spin liquids and topological orders.