Lee-Yang Zeros and Pseudocritical Drift in J-Q Néel-VBS Transitions

This study utilizes Lee-Yang zero scaling in stochastic series expansion quantum Monte Carlo simulations to demonstrate that the Néel-to-valence-bond-solid transition in square-lattice J-Q models is weakly first-order, characterized by a systematic pseudocritical drift and an effective order-parameter scaling dimension that violates unitary conformal field theory bounds in the thermodynamic limit.

The Gist

Imagine you are trying to figure out what happens when you slowly turn a dial on a complex machine. Does the machine change its behavior smoothly and continuously, like a dimmer switch fading a light from bright to dark? Or does it suddenly snap, like a light switch flipping from off to on?

In the world of quantum physics, scientists have been arguing about a specific "machine" called the J-Q model. This model describes how tiny magnets (spins) on a grid interact. The big question is: When these magnets switch from one organized pattern (called Néel) to another (called Valence-Bond Solid or VBS), do they do it smoothly (a "continuous" transition) or do they snap (a "first-order" transition)?

For years, the evidence was confusing. The machine looked like it was dimming smoothly, but the data was fuzzy, and the "snap" might just be hiding because the machine is too small to show the full picture.

This paper introduces a new, super-sensitive way to look at the machine: Lee–Yang Zeros.

The Analogy: The "Ghost" Switch

Think of the machine's behavior as a landscape. In normal physics, we look at the "real" ground (the real numbers) to see if there's a cliff (a phase transition). But sometimes, the cliff is hidden or the ground is too bumpy to see clearly.

The Lee–Yang method is like taking a magical flight into the "complex plane"—a hidden dimension where numbers can be imaginary. In this hidden dimension, the machine has "ghost switches" (the zeros).

• If the transition is smooth (Continuous): These ghost switches float around and slowly drift toward the real ground as the machine gets bigger. They act like a gentle slope.
• If the transition is a snap (First-Order): These ghost switches rush straight toward the real ground and pinch it tightly, like a pair of tongs closing. This pinch happens very fast and indicates a sudden, violent change.

What the Scientists Did

The researchers used a powerful computer simulation technique (called Quantum Monte Carlo) to track these "ghost switches" for different sizes of the machine.

1. The Test Run (The Benchmarks):
   First, they tested their method on machines where they already knew the answer.

   • The Smooth Machine: They looked at a model known to change smoothly. The ghost switches behaved exactly as expected for a smooth transition, drifting gently.
   • The Snapping Machine: They looked at a model known to snap. The ghost switches rushed in and pinched the ground, confirming the method works.

2. The Mystery Machine (The J-Q Models):
   Then, they applied this to the famous J-Q models (J-Q2 and J-Q3) that have been debated for years.

   • The Old View: Previous studies looked at the "real ground" and saw the machine acting like it was changing smoothly. It looked like a dimmer switch.
   • The New View: When the researchers tracked the "ghost switches," they saw something different. As they made the machine larger and larger, the switches didn't just drift gently. They started to drift downward and accelerate toward the "pinch" behavior.

The Big Discovery

The paper concludes that the J-Q models are not changing smoothly. Instead, they are undergoing a "weakly first-order" transition.

Here is the best way to visualize this:
Imagine a crowd of people trying to decide whether to stand up or sit down.

• Smooth Transition: Everyone slowly stands up one by one over a long time.
• Strong Snap: Everyone jumps up instantly at the same time.
• Weakly First-Order (The Discovery): It looks like everyone is slowly standing up for a long time (the "pseudocritical" regime). But if you wait long enough and watch a huge crowd, you realize they are actually all preparing to jump up together. The "ghost switches" in this paper are the early warning signs that the crowd is really about to jump, even though it looks like they are just stretching their legs.

Why This Matters

This is a big deal because it solves a long-standing mystery in physics.

• It proves that the "smooth" behavior seen in earlier experiments was an illusion caused by the machines being too small to show the final "snap."
• It establishes that the Lee–Yang Zero method is a super-powerful tool. It can see the "snap" coming even when everything else looks calm.
• It tells us that nature, in these specific quantum magnets, prefers a sudden change, even if that change is very subtle and takes a long time to reveal itself.

In short: The paper used a "ghost detector" to prove that a quantum magnet transition, which looked like a gentle fade, is actually a slow-motion explosion waiting to happen.

Technical Summary

1. Problem Statement

The paper addresses the long-standing debate regarding the nature of the quantum phase transition between the Néel antiferromagnetic (AFM) phase and the columnar valence-bond solid (VBS) phase in square-lattice $J$-$Q$ models.

• The Controversy: Conventional Landau-Ginzburg-Wilson theory suggests that a direct continuous transition between these two phases (breaking distinct symmetries) is unlikely without fine-tuning. However, the Deconfined Quantum Critical Point (DQCP) scenario posits a continuous transition governed by fractionalized spinons and emergent gauge fields.
• The Challenge: Numerical studies using Quantum Monte Carlo (QMC) have yielded ambiguous results. While Binder cumulants and correlation-length ratios often appear to cross at a single point (suggesting continuity), these crossings drift as system size increases. Furthermore, effective critical exponents depend heavily on fitting windows, and joint order-parameter distributions show emergent symmetry that could mimic a true critical point but actually arise from a weakly first-order transition with an extended pseudocritical regime.
• The Goal: The authors aim to definitively distinguish between a stable continuous critical point and a weakly first-order transition by analyzing the finite-size scaling of Lee–Yang zeros, which provide a direct probe of the singularity scale in the complex plane.

2. Methodology

The authors employ Stochastic Series Expansion (SSE) Quantum Monte Carlo combined with a Lee–Yang zero analysis in complex source fields.

• Lee–Yang Formalism:

  • Instead of observing real-axis observables, the authors introduce an imaginary source field $ig_Q$ coupled to an order parameter operator $O_Q$ (e.g., staggered magnetization for Néel, bond energy for VBS).
  • They compute the normalized generating function $G_{L,\beta}(ig_Q) = Z(ig_Q)/Z(0)$.
  • Lee–Yang Zeros: The values of $ig_Q$ where $G_{L,\beta}$ vanishes.
  • Scaling Distinction:
    • Continuous Transition: The leading zero $g_{LY}$ scales as $L^{-y_h}$, where $y_h$ is the renormalization group eigenvalue. This corresponds to a non-trivial anomalous dimension $\eta$.
    • First-Order Transition: The leading zero scales with the inverse spacetime volume, $g_{LY} \sim (\beta L^d)^{-1} \sim L^{-(d+z)}$. This corresponds to an effective scaling dimension $\Delta_\phi \to 0$ (or $\eta \to -1$).

• Implementation:

  • Reweighting: They sample configurations from a real, symmetry-preserving anchor Hamiltonian ($H_0$) and reweight them to evaluate the complex generating function.
  • Two-Size Estimator: To eliminate non-universal amplitudes, they construct an effective exponent estimator:
    $$\eta_{\text{eff}}(L) = d + z + 2 + \frac{2}{\ln 2} \ln \left[ \frac{g_{LY}(2L)}{g_{LY}(L)} \right]$$
    • For a continuous transition, $\eta_{\text{eff}} \to \eta$ (a constant).
    • For a first-order transition, $\eta_{\text{eff}} \to -1$.

• Models Studied:

  1. Benchmark 1 (Continuous): Columnar dimerized Heisenberg model (known $O(3)$ critical point).
  2. Benchmark 2 (First-Order): Checkerboard $J$-$Q$ (CBJQ) model (known first-order transition).
  3. Target Models: Uniform $J$-$Q_2$ and $J$-$Q_3$ models (candidates for DQCP).

3. Key Contributions

1. Development of a Symmetry-Resolved Lee–Yang Framework: The authors adapted SSE QMC to compute Lee–Yang zeros for specific symmetry channels (Néel and VBS) in quantum magnets, allowing for a direct diagnostic of transition order.
2. Validation of the Method: They successfully recovered the expected $O(3)$ scaling for the dimerized Heisenberg model and the clear $L^{-(d+z)}$ volume-law scaling for the CBJQ model, proving the method's reliability.
3. Resolution of the J-Q Debate: By applying this method to $J$-$Q_2$ and $J$-$Q_3$, they provided direct evidence that these transitions are weakly first-order, characterized by a slow drift toward the first-order scaling limit rather than a stable fixed point.

4. Results

• Benchmark Validation:

  • Dimerized Heisenberg: The effective exponent $\eta_{\text{eff}}$ converged to $\approx 0.0375$, matching the known 3D $O(3)$ value. The leading zeros followed the critical scaling law $L^{-y_h}$.
  • CBJQ Model: The effective exponent $\eta_{\text{eff}}$ drifted systematically toward $-1$ as system size increased, and the zeros followed the volume-law scaling $L^{-(d+z)}$, confirming the first-order nature.

• J-Q2 and J-Q3 Models:

  • Systematic Drift: In both Néel and VBS channels, the effective exponent $\eta_{\text{eff}}(L)$ showed a pronounced downward drift with increasing system size $L$.
  • Scaling Behavior: The leading zeros did not stabilize at a critical value but instead moved toward the first-order limit ($\eta \to -1$).
  • Effective Scaling Dimension: The data implies an effective scaling dimension $\Delta_\phi$ of the order-parameter field that decreases with size and vanishes in the thermodynamic limit.
  • Unitarity Bound Violation: The observed behavior lies below the scalar unitarity bound for any unitary relativistic conformal field theory in $2+1$ dimensions, which is only possible for a first-order transition.
  • Pseudocritical Regime: The drift is slow, indicating an extended pseudocritical regime where the system appears quasi-critical over accessible intermediate size scales before eventually bending toward the first-order coexistence behavior.

5. Significance

• Definitive Diagnosis: The study establishes Lee–Yang zeros as a sensitive, symmetry-resolved diagnostic tool capable of distinguishing true criticality from weakly first-order transitions, even when conventional real-axis observables (like Binder cumulants) appear stable or ambiguous.
• Resolution of DQCP Question: The results strongly support the interpretation that the Néel-VBS transition in uniform $J$-$Q$ models is weakly first-order, rather than a stable deconfined quantum critical point.
• Theoretical Implication: The vanishing of the effective scaling dimension in the thermodynamic limit enforces inverse spacetime-volume scaling, a hallmark of first-order transitions. This suggests that the "emergent symmetry" observed in these models is a finite-size artifact of the slow crossover rather than a sign of a true continuous phase transition.
• Methodological Impact: The work demonstrates the power of complexified partition function analysis (Lee–Yang theory) in resolving long-standing controversies in quantum many-body physics, offering a robust alternative to standard finite-size scaling analysis.