Pseudocriticality in antiferromagnetic spin chains

By combining advanced quantum Monte Carlo simulations with a novel loop estimator for Rényi entanglement entropy, this study demonstrates that an SU($N$) generalization of the Heisenberg antiferromagnet in 1+1 dimensions exhibits weak first-order pseudocriticality driven by proximity to a complex conformal field theory, a finding that accurately recovers the real part of the complex central charge for $N>2$ and reinterprets the dimerized phase of the spin-1 chain as pseudocritical.

The Gist

Imagine you are a physicist trying to understand how a chain of tiny magnets (spins) behaves when you cool them down to absolute zero. Usually, these chains do one of two things: they either line up perfectly in an orderly pattern (like soldiers), or they stay in a chaotic, fluctuating state where everything is connected over long distances (like a crowd at a concert).

But sometimes, nature plays a trick. It creates a "fake" critical point. It looks like the system is on the verge of a massive, smooth transition, but it's actually just a very slow, weak snap into a different state. The authors of this paper call this "Pseudocriticality."

Here is a simple breakdown of what they did, using some everyday analogies.

1. The Setup: A Chain of Colorful Beads

The researchers studied a specific type of magnetic chain called the SU(N) Heisenberg model.

• The Analogy: Imagine a necklace made of beads. In a normal magnet, the beads are just "Up" or "Down." In this model, the beads have "colors" or "flavors" (like Red, Blue, Green, etc.). The number of colors is represented by the letter $N$.
• The Goal: They wanted to see what happens as they change the number of colors ($N$) from 2 to 3, 4, and so on.

2. The Mystery: The "Ghost" Phase

In physics, there is a famous concept called a Critical Point. This is the exact moment a material changes state (like ice melting to water). At this point, the system becomes "scale-invariant," meaning it looks the same whether you zoom in or out.

However, for certain values of $N$ (specifically when $N > 2$), the math suggests that the "perfect" critical point doesn't actually exist in our real world. Instead, it exists in a parallel universe of complex numbers (mathematical ghosts).

• The Analogy: Imagine you are driving toward a destination. Usually, you drive straight there. But in this case, the destination is a "ghost town" located in a parallel dimension. You can't actually get there, but because it's so close, your car (the physical system) starts driving very slowly, swerving, and acting like it might be arriving. This slow, swerving behavior is the "walking" or "pseudocritical" behavior.

3. The Problem: How Do You Measure a Ghost?

To prove this "ghost" critical point exists, the scientists needed to measure something called the Central Charge.

• The Analogy: Think of the Central Charge as the "volume knob" of the system's complexity. In a normal critical point, the volume is loud and clear. In a "ghost" critical point, the volume is a complex number (a mix of real and imaginary sound).
• The Challenge: You can't measure an imaginary sound directly. However, if you look at how the volume changes as you look at the chain from different distances (sizes), you can hear the "echo" of that ghost. The volume starts to drift or slide as the chain gets longer.

4. The Solution: The "Non-Equilibrium Work" Trick

The authors used a super-advanced computer simulation called Quantum Monte Carlo. But standard simulations were too slow and inaccurate for this tricky problem.

• The Innovation: They invented a new way to measure the "entanglement entropy" (a measure of how connected the beads are).
• The Analogy: Imagine you want to measure the weight of a feather, but the scale is too sensitive. Instead of just putting the feather on, you imagine a process where you slowly "inflate" a balloon attached to the feather and measure the work it takes to blow it up. By measuring the "work" done during this imaginary process, they could calculate the weight with incredible precision.
• The "Replica Shift": They also had to deal with a problem where the simulation got stuck in one pattern (like a car stuck in mud). They developed a "replica shift" technique, which is like giving the car a gentle nudge to the left or right to help it get unstuck and explore the whole road.

5. The Results: Catching the Ghost

What did they find?

1. Perfect Agreement: When they measured the "volume knob" (Central Charge) for different numbers of colors ($N$), the results matched the predictions of the "ghost" theory perfectly.
2. The Drift: For $N > 2$, they saw the volume knob slowly drift downward as the chain got longer. This drift is the signature of the system being near a complex, ghostly critical point.
3. The Spin-1 Chain: They discovered that a specific, well-known magnetic chain (the Spin-1 chain) is actually living right next to this ghost. This means the "dimerized" phase (where the beads pair up) isn't just a boring, stable state; it's actually a "pseudocritical" state, teetering on the edge of a complex transition.

Why Does This Matter?

This is a big deal for two reasons:

• It Solves a Puzzle: It explains why some magnetic materials behave strangely, showing "almost critical" behavior without actually being critical.
• It Connects Worlds: It proves that "ghost" theories (Complex Conformal Field Theories) aren't just math tricks; they have real, measurable effects on physical materials. It's like proving that a ghost town exists by measuring the fog rolling off it.

In a nutshell: The authors built a super-precise microscope to look at a chain of magnets. They found that for certain settings, the magnets aren't just behaving normally; they are "walking" toward a critical point that exists in a mathematical parallel universe. By measuring how the magnets "drift" as they get larger, they successfully mapped out the properties of this invisible, ghostly world.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the phenomenon of weak first-order pseudocriticality, where systems exhibit approximate scale invariance and slow renormalization group (RG) flows despite undergoing a first-order phase transition. This behavior is often attributed to the system's proximity to a complex Conformal Field Theory (CFT) fixed point located outside the physical parameter space (i.e., with complex-valued couplings).

While this mechanism has been proposed to explain "walking" behavior in 2+1D deconfined quantum critical points and the $Q > 4$ Potts model, rigorous verification in simpler 1+1D quantum models is lacking. The authors aim to:

• Establish the existence of complex-CFT-induced pseudocriticality in a 1+1D quantum spin model.
• Investigate the SU(N) generalization of the Heisenberg antiferromagnet, a model that interpolates between the spin-1/2 Heisenberg chain ($N=2$) and the spin-1 biquadratic chain ($N=3$).
• Determine if the model's proximity to a complex CFT can be tuned continuously via the parameter $N$, and whether the real part of the complex central charge can be recovered from finite-size drifts in the pseudocritical regime ($N > 2$).

2. Methodology

The authors employ state-of-the-art Quantum Monte Carlo (QMC) simulations combined with advanced entanglement entropy estimation techniques.

• Model: The study focuses on the SU(N) Heisenberg antiferromagnet on a 1D chain with fundamental representation on one sublattice and conjugate representation on the other. The Hamiltonian is defined as the negative sum of nearest-neighbor singlet projectors:
  $$H = -J \sum_{i=1}^{L} \sum_{\alpha, \beta=1}^{N} |\alpha_i \alpha_{i+1}\rangle \langle \beta_i \beta_{i+1}|$$
  This model maps to the quantum Potts chain with $Q=N^2$ and is equivalent to classical O(n) loop models.

• Continuous N: Unlike standard lattice models restricted to integer $N$, the authors utilize a Resummed Stochastic Series Expansion (RSSE) loop representation. This allows the symmetry group parameter $N$ to be treated as a continuous real variable ($N > 0$), enabling a fine-grained scan of the phase diagram.

• Rényi Entanglement Entropy (EE) Estimation:

  • The second Rényi entropy ($S_2$) is calculated as a free energy difference between two generalized partition functions, $Z_A$ (with a trace over subsystem $A$) and $Z_\emptyset$ (empty set).
  • Nonequilibrium Work Protocol: To compute $S_2$ with high precision, the authors introduce an interpolating ensemble $Z(\lambda)$ where $\lambda$ varies from 0 to 1. They employ a Jarzynski equality based protocol:
    $$e^{-S_2} = \langle e^{-W} \rangle$$
    where $W$ is the work done during a nonequilibrium process varying $\lambda$.
  • Loop Estimator: A new loop-based estimator is developed for the work calculation, leveraging the loop decomposition of the partition function to reduce stochastic errors significantly.

• Replica Shift Update: For $N > 2$, the system enters a dimerized phase with two degenerate ground states. Standard QMC updates suffer from long tunneling times (ergodicity breaking) between these patterns. The authors introduce a "replica shift" update, shifting the operator string by one lattice site, to ensure efficient sampling and accurate EE calculation in the gapped regime.

3. Key Contributions

1. Continuous Tuning of Pseudocriticality: The paper demonstrates that the SU(N) chain in 1+1D provides a tunable platform to study the transition from critical ($N \le 2$) to pseudocritical ($N > 2$) behavior by continuously varying $N$.
2. Novel Entanglement Algorithm: The development of a highly efficient, nonequilibrium loop protocol for calculating Rényi EE in continuous-$N$ QMC simulations, overcoming ergodicity issues in dimerized phases via the replica shift method.
3. Recovery of Complex CFT Data: The authors successfully extract the real part of the complex central charge ($\text{Re}(c)$) for $N > 2$ by analyzing the finite-size drift of the central charge, confirming the "walking" behavior predicted by RG flows near complex fixed points.
4. Spin-1 Chain Insight: By setting $N=3$, the work provides new evidence that the dimerized phase of the spin-1 biquadratic chain is not merely gapped but is pseudocritical, residing in close proximity to a complex CFT.

4. Results

• Central Charge ($c$) Behavior:

  • $N < 1$: The central charge is negative, consistent with CFT predictions for non-unitary theories.
  • $1 \le N \le 2$: The system is gapless. The extracted central charge $c$ converges to the CFT prediction $c = 1 - \frac{6(1-g)^2}{g}$ (where $N = -2\cos(\pi g)$). The results show excellent agreement with theory.
  • $N > 2$: The CFT central charge becomes complex ($c = c_R + i c_I$). In finite-size systems, the measured central charge $c(L)$ initially converges to the real part $c_R$ but eventually drifts toward zero (the thermodynamic limit of a gapped phase).
  • Drift Analysis: For $N \ge 2.6$, the authors fit the drift of $c(L)$ using the RG flow equation derived for complex fixed points:
    $$c(L) = c_R - \alpha \tan(\alpha^{1/3} \ln(L/L_0))$$
    This fitting allows them to recover $c_R$ with remarkable accuracy, matching the real part of the theoretical complex CFT central charge even up to $N=3.6$ ($Q \approx 13$ Potts model).

• Subsystem Oscillations:

  • The Rényi EE exhibits oscillatory terms decaying as $D(l_A, L)^{-x}$.
  • The decay exponent $x$ is found to match the Potts thermal exponent relation $x = 2 - y_T^{\text{Potts}}$, rather than the standard Luttinger parameter, providing a unique signature of the underlying CFT structure.
  • For $N > 2$, the oscillatory form breaks down, signaling the onset of complex scaling dimensions.

• Spin Wave Velocity: The authors also computed the spin wave velocity $v$, finding it decreases exponentially with $N$. For $N=2$, $v \approx \pi/2$, matching theoretical predictions.

5. Significance

• Validation of Complex CFT Theory: This work provides one of the most direct and precise numerical validations of the theory that weak first-order transitions and "walking" behavior are driven by proximity to complex CFT fixed points.
• Bridging Dimensions: It connects 1+1D exactly solvable models (via CFT) with the physics of 2+1D deconfined criticality, suggesting that the mechanisms observed in higher dimensions may have simpler, rigorous counterparts in 1D.
• Spin-1 Physics: The findings redefine the understanding of the spin-1 biquadratic chain's dimerized phase, suggesting it is a pseudocritical state rather than a simple gapped state, opening new avenues for studying non-Hermitian extensions and complex CFTs in field theories with topological terms.
• Methodological Advancement: The combination of continuous-$N$ QMC, nonequilibrium work protocols, and replica shifts sets a new standard for calculating entanglement properties in quantum many-body systems, particularly in regimes with strong dimerization or complex fixed points.