Ground-state Entropy of the Ising model on a Frustrated lattice

This paper reports the ground-state entropy of the 2D Ising model on the Shastry-Sutherland lattice and investigates a generalized version of the model where the constraint on zero-temperature configurations is continuously removed.

The Gist

Imagine a giant, endless floor made of tiles. But this isn't a normal floor; it's a special pattern of squares and triangles stuck together, known in physics as the Shastry-Sutherland lattice.

On this floor, we place tiny magnets (called "spins") at every corner. Each magnet can point either Up or Down. The rule of the game is simple: neighbors hate being the same. If two magnets are next to each other, they want to point in opposite directions to be "happy" (low energy). This is called an antiferromagnetic setup.

The Problem: The Frustrated Floor

Here's the catch: The floor is shaped in a way that makes it impossible for everyone to be happy at once. This is called frustration.

Imagine a triangle of three magnets. If Magnet A points Up and Magnet B points Down to satisfy their bond, Magnet C is stuck. It can't be opposite to both A and B at the same time. One bond will always be unhappy.

In this specific lattice, there are two types of connections:

1. The Sides: The edges of the squares and triangles.
2. The Diagonals: The lines cutting across the squares.

The paper studies what happens when the "diagonal" connections are very strong (stronger than the sides).

The Two Scenarios

Scenario A: The "Strict" Rule (High Strength)
When the diagonal connections are super strong, the magnets have a very easy time. They just pair up on every diagonal line: one Up, one Down. It's like a dance where every partner is strictly assigned.

• Result: There are many ways to arrange these pairs, but the rules are rigid. The "disorder" (or entropy) is easy to calculate.

Scenario B: The "Relaxed" Rule (The Sweet Spot)
The paper focuses on a specific moment where the diagonal strength is just right (a value called $\alpha = 1$). Suddenly, the rules loosen up. Now, magnets on the diagonal lines are allowed to point in the same direction (both Up or both Down), which was forbidden in the strict scenario.

• The Chaos: This tiny permission creates a massive explosion of possibilities. The magnets can now arrange themselves in countless different ways while still keeping the total energy at the lowest possible level.
• The Question: How many ways can they do this? In physics, we call this number the Ground-State Entropy. It's a measure of how "confused" or "disordered" the system is even when it's as cold as possible (absolute zero).

How the Authors Solved It

Calculating this number is like trying to count every possible way to arrange a deck of cards in a room the size of a galaxy. It's too big for a normal computer.

The authors used two clever tricks:

1. The "Row-by-Row" Method (Transfer Matrix): Imagine building the floor one row of magnets at a time. They created a mathematical machine that calculates how many ways you can add the next row based on the one before it. They ran this on small sections and used math to guess what happens on an infinite floor.
2. The "Corner" Method (CTMRG): This is like looking at a single spot on the floor and asking, "If I zoom out infinitely, what does the average neighborhood look like?" They used a modern, high-powered algorithm (Tensor Networks) to simulate this infinite zoom.

The Big Discovery

After running these complex calculations, the authors found the exact number for how "disordered" this system is at the sweet spot ($\alpha = 1$).

• The Number: The entropy is approximately 0.4588 (per magnet).
• Why it matters: Before this paper, scientists only knew a "lower bound" (a minimum guess). They knew it was at least this much, but didn't know the exact ceiling. This paper pins down the exact value.

The "Magic Dial"

To make sure their math was right, the authors introduced a "dial" (a parameter called $r$).

• Turn the dial to 0: You force the magnets to follow the strict rules (no parallel spins on diagonals). The system is simple, and the math is easy.
• Turn the dial to 1: You allow the relaxed rules. The system becomes complex and "frustrated."

They watched the entropy grow smoothly as they turned the dial from 0 to 1. This confirmed that their calculations were consistent and that the transition from the "strict" world to the "frustrated" world is continuous, not a sudden jump.

Summary

In simple terms, the authors solved a long-standing puzzle about a specific pattern of magnets. They figured out exactly how many different ways these magnets can arrange themselves when they are at their lowest energy state, but still stuck in a "frustrated" pattern where they can't all be happy. They found the answer is roughly 0.4588, a precise number that had been hiding in the math for years.

Technical Summary

Technical Summary: Ground-State Entropy of the Ising Model on a Frustrated Lattice

Problem Statement
This work addresses the calculation of the exact ground-state entropy for the antiferromagnetic Ising model defined on the Shastry-Sutherland (SS) lattice. While previous studies [1] established that the system exhibits a highly degenerate ground state for the parameter regime $\alpha \geq 1$ (where $\alpha$ represents the ratio of diagonal to square-bond exchange interactions), the precise value of the extensive entropy remained unknown. The specific challenge lies in the regime $\alpha = 1$, where the ground state manifold includes not only the spin-pair configurations found at $\alpha > 1$ but also new configurations involving parallel spins on diagonal bonds. The authors aim to compute this entropy exactly and investigate the continuous transition of the system's properties as the constraint on these parallel configurations is relaxed.

Methodology
The authors employ two primary computational approaches to solve the zero-temperature ($T=0$) statistical mechanics problem:

1. Transfer Matrix Formulation:
   The problem is mapped to a row-to-row transfer matrix $T_L$ acting on a chain of length $L$. The lattice is decomposed into alternating rows of type (A) and (B), corresponding to the diagonal orientations of the bonds. The partition function is expressed as $Z_{LM} = \text{Tr}(T_L^{M/2})$. To study the transition between the $\alpha > 1$ limit and the $\alpha = 1$ case, the authors introduce a continuous parameter $r$ (where $r=0$ corresponds to $\alpha > 1$ and $r=1$ to $\alpha = 1$). This parameter weights the "forbidden" parallel-spin configurations on diagonal bonds. The transfer matrix elements are constructed using local operators $V^{(A)}$ and $V^{(B)}$, which are expressed in terms of Pauli matrices and permutation operators.

2. Numerical and Variational Techniques:

   • Corner Transfer Matrix Renormalization Group (CTMRG): To access the thermodynamic limit ($L, M \to \infty$), the authors utilize the CTMRG algorithm. This tensor network method approximates the environment of a local tensor using corner and edge transfer matrices truncated to a dimension $\chi$. The partition function is computed as a contraction of a grid of tensors representing the Boltzmann weights of local clusters.
   • Exact Diagonalization and Variational Bounds: For finite systems, the largest eigenvalue of the transfer matrix is calculated using exact diagonalization (via the DiracQ package). Additionally, a variational lower bound is derived using a Rayleigh-Ritz ansatz. The authors construct a specific trial state $|\Psi_A\rangle$ based on the leading eigenvector of the operator $T_A$ and calculate the expectation value $\langle \Psi_A | T_B T_A | \Psi_A \rangle$ to establish a rigorous lower bound for the entropy.

Key Results

• Ground-State Entropy at $\alpha = 1$: Using CTMRG with a truncation dimension $\chi = 16$, the authors determine the entropy per site to be $\Sigma \approx 0.45877772$. This result is fully converged within machine precision.
• Variational Lower Bound: The variational approach yields an infinite-size lower bound of $\Sigma_{\text{bound}} \approx 0.45668258$. The finite-size transfer matrix calculations for $L=4, 6, 8, 10$ extrapolate to a value of $\Sigma \sim 0.4588 \pm 0.0002$, which is consistent with the CTMRG result.
• Generalized Model ($r$ dependence): By varying the parameter $r$ from 0 to 1, the authors map the continuous evolution of the system:
  • At $r=0$ (corresponding to $\alpha > 1$), the entropy is exactly $\Sigma = \frac{1}{2} \ln 2 \approx 0.3466$. The system exhibits "superstability," where the transfer matrix possesses a flat eigenfunction with equal weights for all configurations.
  • As $r$ increases to 1 ($\alpha = 1$), the entropy increases continuously to the calculated value of $\approx 0.4588$.
  • The fraction of ferromagnetic diagonal bonds ($n_{fmd}$) and the spin-spin correlation length ($\xi$) are also computed. The correlation length exhibits a weak logarithmic singularity as $r \to 0$.
• Asymptotic Behavior: Near $r=0$, the authors derive asymptotic formulas showing that the entropy scales as $\Sigma(r) \sim \frac{1}{2}\ln 2 + \frac{r}{8}$ and the fraction of parallel bonds scales as $n_{fmd}(r) \sim \frac{r}{4}$.

Significance and Claims
The paper claims to provide the first exact determination of the ground-state entropy for the antiferromagnetic Ising model on the Shastry-Sutherland lattice at the critical point $\alpha = 1$. The work resolves a long-standing open problem where only lower bounds were previously available.

The authors highlight that their variational lower bound corrects a previous numerical error found in earlier literature [1], where a value of 0.4812 was reported (approximately 5% higher than the exact value). By combining the CTMRG method with rigorous variational bounds and finite-size scaling, the authors establish a precise value for the entropy, characterizing the extensive degeneracy of the frustrated ground state. The study also elucidates the nature of the transition between the $\alpha > 1$ and $\alpha = 1$ regimes, demonstrating that while the transition in terms of the parameter $\alpha$ is abrupt (due to the sudden appearance of new configurations), the relaxation of the constraint via the parameter $r$ reveals a smooth, continuous evolution of thermodynamic quantities.