Specific-heat anomaly in frustrated magnets with vacancy defects

This paper demonstrates that vacancy defects in a two-dimensional geometrically frustrated magnet, such as the triangular-lattice antiferromagnetic Ising model, induce a specific-heat peak at a characteristic temperature determined by the defect concentration, arising from the thermal relaxation of entropy-constrained bulk spins.

The Gist

Imagine a crowded dance floor where everyone is trying to dance in perfect harmony, but the music is tricky. This is the world of frustrated magnets.

In these materials, the atoms (spins) want to align in a specific pattern, but the geometry of the crystal lattice makes it impossible for everyone to be happy at the same time. It's like a game of musical chairs where there are more chairs than people, but the rules say everyone must sit next to someone they dislike. Because of this "frustration," the atoms can't settle into a single, rigid order even when it gets very cold. Instead, they keep jiggling and shifting, creating a state of constant, chaotic potential energy.

Now, imagine you take a few dancers off the floor. You create vacancies (empty spots).

The Main Discovery: The "Vacancy Peak"

The paper by Sedik, Zhu, and Syzranov reveals a surprising thing happens when you remove dancers from this chaotic floor.

In a perfect, empty-free dance floor, as you cool the room down, the dancers slowly slow their movements. The "heat capacity" (a measure of how much energy the system absorbs or releases as it changes temperature) drops smoothly to zero.

But when you have empty spots (vacancies), something weird happens. As you cool the system down, the heat capacity doesn't just drop; it suddenly spikes at a specific low temperature, creating a mountain-like peak on a graph.

The Analogy: The "Frozen Crowd" and the "Thaw"

Here is a simple way to visualize why this peak happens:

1. The Freeze (Low Temperature):
   Imagine the dance floor is very cold. The dancers are moving very slowly. The empty spots (vacancies) act like traffic cops or anchors. Because the dancers are so close to each other, the empty spots force the remaining dancers into very specific, rigid patterns to avoid breaking the rules of the dance.

   • Result: The empty spots "freeze" the freedom of the surrounding dancers. They are stuck in place. The system has very few options for how to move, so its "entropy" (disorder/freedom) is low.

2. The Thaw (Rising Temperature):
   Now, imagine you slowly turn up the heat. The dancers start to wiggle more.

   • At a certain critical temperature (determined by how many empty spots there are), the "traffic cops" lose their grip. The constraints imposed by the empty spots suddenly relax.
   • The dancers, who were previously frozen in rigid patterns, suddenly realize, "Hey, I can move again!" They gain a massive burst of freedom.
   • Result: This sudden explosion of freedom creates a huge spike in entropy (disorder).

3. The Peak:
   In physics, when a system suddenly gains a lot of freedom (entropy) as it warms up, it has to absorb a lot of heat to do so. This absorption shows up as a peak in heat capacity.

   • Think of it like a dam breaking. Before the break, the water is held back (frozen degrees of freedom). When the dam breaks (at temperature $T_{imp}$), a massive rush of water (entropy) is released. The "peak" is the moment that rush happens.

Why is this important?

• It's a Fingerprint: This peak acts like a fingerprint for the amount of "dirt" (defects) in the material. By measuring the temperature where the peak happens, scientists can calculate exactly how many missing atoms are in the crystal, even if they can't see them directly.
• Quantum Spin Liquids: Scientists are hunting for a special state of matter called a "Quantum Spin Liquid," where atoms never stop moving, even at absolute zero. However, real-world materials always have defects. This paper explains that these defects might be hiding a secret signal (the peak) that tells us how the material is behaving, helping researchers distinguish between a true quantum liquid and a messy, defective one.

The "Magic Number"

The authors did the math (using a model called the Ising model on a triangular lattice) and found a simple formula for when this peak happens:
$$T_{imp} \approx \frac{-4J}{\ln(\text{amount of missing spots})}$$

• $J$: How strongly the dancers want to follow the rules.
• $\ln(\text{missing spots})$: The logarithm of how many empty seats there are.

In plain English: The more empty seats you have, the higher the temperature needs to be to "thaw" the frozen dancers. If you have very few empty seats, the peak happens at a very low temperature. If you have many, it happens at a higher temperature.

Summary

This paper explains that holes in a material don't just make it weaker; they change how it heats up. By freezing the movement of nearby atoms, vacancies create a "bottleneck" of order. When the material warms up enough to break that bottleneck, it releases a burst of chaos, creating a distinct, measurable spike in heat capacity. It's a beautiful example of how a small defect can create a large, predictable signal in the quantum world.

Technical Summary

1. Problem Statement

Geometrically frustrated magnets (GFMs), particularly those considered candidates for quantum spin liquids (QSLs), are highly sensitive to quenched disorder such as randomly located impurities and vacancy defects. While clean frustrated systems often lack long-range magnetic order at low temperatures, the introduction of vacancies imposes severe constraints on the bulk spin degrees of freedom.
The central problem addressed is understanding how these vacancy defects alter the thermodynamics of 2D frustrated systems. Specifically, the authors investigate whether vacancies induce new characteristic energy scales and anomalies in the heat capacity $C(T)$ that are distinct from the intrinsic energy scales of the clean material (such as the Curie-Weiss temperature $\theta_{CW}$ or the activation gap of the clean system).

2. Methodology

The authors employ an analytical approach based on the Kac-Ward formalism (loop representation) applied to the antiferromagnetic (AFM) Ising model on a triangular lattice.

• Loop Representation: The partition function $Z$ is expressed as a sum over closed loops on the lattice. The contribution of each loop depends on its length $r$ and a parity factor determined by the number of self-intersections and the rotation angle of the loop tangent.
• Defect-Free System: The authors first reproduce the known thermodynamics of the clean triangular lattice Ising model. They utilize a Fourier transform to diagonalize the transfer matrix (a $6 \times 6$ matrix in direction space) to calculate the entropy and heat capacity, confirming the well-known residual entropy $S_0 \approx 0.323$ per spin and the exponential suppression of heat capacity at low temperatures with an activation gap of $4J$.
• Single Vacancy Calculation: To model a vacancy at site $\rho_0$, the authors modify the partition function by subtracting the contributions of all loops that pass through $\rho_0$. This is achieved by averaging the loop sums over all possible vacancy locations, effectively replacing the sum over loops containing $\rho_0$ with a fraction proportional to the loop length.
• Dilute Limit: For a system with a low concentration of vacancies ($n_{imp} = N_{imp}/N \ll 1$), the authors assume the contributions of individual vacancies are additive, provided the temperature is high enough that the spin correlation length $\xi(T)$ is smaller than the average inter-vacancy distance.

3. Key Contributions and Results

A. Emergence of a Vacancy-Induced Heat Capacity Peak

The primary result is the prediction of a distinct peak in the heat capacity $C(T)$ at a specific temperature $T_{imp}$, driven solely by the presence of vacancy defects.

• Peak Temperature: The temperature of this anomaly is determined by the vacancy concentration $n_{imp}$ and the coupling $J$:
  $$T_{imp} \approx -\frac{4J}{\ln n_{imp}}$$
  This scale is distinct from the Curie-Weiss temperature ($\theta_{CW} \sim zJ$) and the intrinsic low-temperature peaks of clean Heisenberg models.
• Heat Capacity Behavior:
  • Low $T$ ($T \ll T_{imp}$): Vacancies constrain the bulk spins, effectively "freezing" degrees of freedom. The heat capacity contribution from vacancies vanishes as $T \to 0$.
  • Intermediate $T$ ($T \sim T_{imp}$): As temperature increases, the constraints imposed by vacancies relax. The correlation length $\xi(T) \sim \exp(2J/T)$ shrinks, eventually becoming smaller than the inter-vacancy distance. This relaxation releases entropy, causing a sharp rise in $C(T)$.
  • High $T$ ($T \gg T_{imp}$): The vacancy contribution decays as $C_{vac} \propto \frac{N_{imp}}{T^2} e^{2J/T}$.

B. Entropy Analysis

The authors analyze the entropy associated with this peak:

• The entropy released at the peak, $S_{vac} = \int_{peak} \frac{C_{vac}(T)}{T} dT$, scales as $\sqrt{N N_{imp}}$.
• This entropy originates from the lowest-energy degrees of freedom (ground states in the Ising limit).
• While vacancies reduce the total ground-state entropy of the system by $N_{imp} \ln 2$ (by lifting degeneracy), the redistribution of entropy creates a large, localized peak at $T_{imp}$. The peak entropy is significantly larger than the total entropy reduction caused by the missing sites.

C. Analytical Derivations

The paper provides explicit analytical expressions for the triangular lattice:

• Clean System Heat Capacity: $C_0(T) \propto \frac{J^2}{T^2} e^{-4J/T}$ (exponentially suppressed).
• Vacancy Contribution: For $T \gtrsim T_{imp}$, the vacancy contribution is:
  $$C_{vac}(T) \approx \frac{2J^2 N_{imp}}{T^2} e^{2J/T}$$
  The factor of $2J$ in the exponent (compared to $4J$ for the clean system) arises because the vacancy creates a local environment where the energy cost to flip spins is lower, effectively creating "fractionalized" excitations that are bound at low temperatures but become active at $T_{imp}$.

4. Physical Interpretation

The mechanism is interpreted as follows:

1. Constraint at $T=0$: A vacancy forces the surrounding spins into a specific staggered pattern to satisfy the AFM coupling, reducing the number of accessible ground states locally.
2. Relaxation at $T_{imp}$: As temperature rises, thermal fluctuations allow spins to overcome these local constraints. The "frozen" degrees of freedom associated with the vacancy become active.
3. Correlation Length Crossover: The anomaly occurs when the thermal correlation length $\xi(T)$ drops below the average distance between vacancies. Below this temperature, vacancies interact and constrain the bulk; above it, they act as isolated defects, releasing the constrained entropy.

5. Significance and Implications

• Experimental Signature: The paper predicts a specific "impurity peak" in the heat capacity of frustrated magnets (e.g., triangular lattice materials like $Ba_3NiSb_2O_9$ or organic salts) that is distinct from intrinsic QSL features. This provides a diagnostic tool to distinguish between intrinsic low-energy physics and disorder-induced effects.
• Universality: While derived for the Ising model, the authors argue this mechanism is generic for frustrated magnets, including Heisenberg models. In Heisenberg systems, which may already exhibit two peaks (one from transverse coupling and one from Curie-Weiss physics), vacancies would introduce a third peak at $T_{imp}$.
• Impact on QSL Searches: The results caution that low-temperature anomalies often attributed to exotic quantum states (like spin liquids) could instead be artifacts of vacancy defects. The presence of even a small concentration of vacancies can mask or mimic intrinsic low-temperature behavior.
• Theoretical Framework: The work extends the Kac-Ward loop method to include quenched disorder analytically, offering a rigorous framework for studying the thermodynamics of disordered frustrated systems without relying solely on numerical simulations.

In conclusion, the paper establishes that vacancy defects in geometrically frustrated magnets are not merely passive perturbations but active agents that reorganize the entropy landscape, creating a robust, concentration-dependent thermodynamic anomaly at low temperatures.