Complete finite-size scaling theory of Renyi thermal entropy for second, first and weak first order quantum phase transitions

This paper establishes a unified finite-size scaling framework based on Renyi thermal entropy and its derivative to accurately distinguish between second-order, first-order, and weak first-order quantum phase transitions, offering a robust numerical tool that resolves long-standing ambiguities in identifying weak first-order transitions within finite-size simulations.

The Gist

Imagine you are a detective trying to solve a mystery in the quantum world. The mystery? What kind of "phase transition" is happening when a material changes its state?

In the quantum realm, materials can suddenly switch from one behavior to another (like a magnet turning on or off). Physicists call these switches "phase transitions." There are three main types of suspects:

1. The Smooth Operator (Second-Order): A gentle, continuous change.
2. The Hard Switch (First-Order): A sudden, violent jump (like water boiling instantly into steam).
3. The Imposter (Weak First-Order): This is the tricky one. It looks like a smooth change, but deep down, it's actually a hard switch. It's a "wolf in sheep's clothing."

For decades, scientists have struggled to tell the difference between the Smooth Operator and the Imposter when they are running computer simulations. The problem is that the "Imposter" is so good at hiding that even with super-computers, it looks exactly like a smooth transition.

The New Detective Tool: The "Thermal Entropy" Mirror

The authors of this paper, a team from Westlake University, have invented a new magnifying glass to solve this case. They call it the Rényi Thermal Entropy (RTE) and its derivative, the DRTE.

Here is how it works, using a simple analogy:

The Analogy: The Noisy Crowd vs. The Whisper

Imagine you are trying to hear a specific whisper (the "singular" part of the physics that tells you the truth) in a massive, noisy stadium.

• Old Methods: If you just listen to the whole stadium, the noise (analytic contributions) is so loud that you can't hear the whisper. You might think the crowd is just murmuring smoothly, missing the fact that a riot is about to break out.
• The New Method (DRTE): The authors realized that if you take two different recordings of the crowd—one at a normal volume and one at double volume—and subtract them in a specific way, the loud noise cancels out perfectly.
• The Result: Suddenly, the background noise disappears, and you are left with only the whisper. This whisper tells you exactly what is happening.

How the New Tool Solves the Mystery

The paper shows that this new tool (DRTE) behaves differently for each type of transition, acting like a "smoking gun" for the detective:

1. The Smooth Operator (Second-Order)

• What it looks like: The DRTE curve goes up to a single peak and then comes down.
• The Clue: All the curves from different-sized simulations cross at the exact same point. This confirms it's a genuine, smooth transition.

2. The Hard Switch (First-Order)

• What it looks like: The DRTE curve does something wild. It shoots up to a high positive peak, crashes down to a deep negative valley, and crosses zero right in the middle.
• The Clue: This "Double-Peak" structure is the signature of a sudden jump. It's like seeing a car suddenly jump from 0 to 60 mph and then back to 0. It's impossible to miss.

3. The Imposter (Weak First-Order)

• The Mystery: For years, scientists looked at famous models (like the J-Q models) and thought they were seeing a Smooth Operator because the system was too small to see the jump.
• The Breakthrough: When the authors used their new DRTE tool, they saw the Double-Peak and the Zero Crossing even in these "imposter" cases.
• The Verdict: The tool revealed that these systems aren't smooth at all; they are actually Weak First-Order transitions. They were just hiding because the "correlation length" (the distance over which particles talk to each other) was so huge that the computer simulations were too small to see the whole picture.

Why This Matters

Think of this like a new type of X-ray machine. Before, doctors could only see bones (obvious transitions). Now, they can see soft tissue tumors (weak transitions) that were previously invisible.

• It's Universal: It works for all kinds of quantum materials.
• It's Efficient: You don't need a supercomputer the size of a city to find the answer; you just need the right math.
• It Settles Debates: There has been a huge argument in physics about whether certain exotic quantum states are "continuous" or "discontinuous." This paper says, "Stop arguing, here is the proof: they are discontinuous."

In a Nutshell

The authors found a mathematical trick to filter out the "static" in quantum simulations. By doing so, they created a clear, unmistakable signal that tells us exactly how a quantum material is changing its state. They used this to prove that some of the most famous "smooth" transitions in physics are actually sneaky, sudden jumps, finally solving a mystery that has puzzled scientists for years.

Technical Summary

1. Problem Statement

Determining the nature of a quantum phase transition (QPT)—specifically distinguishing between continuous (second-order), strong first-order, and weak first-order transitions—is a fundamental challenge in quantum many-body physics.

• The Core Difficulty: Weak first-order transitions often occur near a critical fixed point with an extremely large correlation length ($\xi$). In finite-size simulations (where system size $L \ll \xi$), these transitions mimic continuous critical behavior, making them indistinguishable from true second-order transitions using standard observables.
• Limitations of Current Methods: Conventional finite-size scaling theories often fail to provide a "smoking-gun" signature for weak first-order transitions. Standard order parameters or free energy derivatives often contain large analytic contributions that obscure the singular parts governing critical behavior.
• The Gap: There is currently no unified, effective finite-size scaling theory capable of robustly differentiating these three types of transitions within computationally accessible system sizes.

2. Methodology

The authors propose a unified framework based on the Rényi Thermal Entropy (RTE) and its derivative (DRTE) as primary observables.

• Rényi Thermal Entropy ($S^{(n)}$):
  Defined as $S^{(n)} = \frac{1}{1-n} \ln \left[ \frac{Z(n\beta)}{Z(\beta)^n} \right]$.
  The authors focus on the second-order case ($n=2$), which corresponds to the purity.

  • Theoretical Insight: By constructing the difference between partition functions at different temperatures ($\ln Z(2\beta) - 2\ln Z(\beta)$), the dominant analytic contributions to the free energy (which follow volume laws) cancel out. This isolates the singular part of the free energy, which contains the critical information.

• Derivative of RTE (DRTE):
  The derivative with respect to a tuning parameter $g$ (e.g., coupling constant) is calculated as:
  $$\frac{dS^{(n)}}{dg} = \frac{1}{1-n} \left[ -n\beta \left\langle \frac{dH}{dg} \right\rangle_{Z(n\beta)} + n\beta \left\langle \frac{dH}{dg} \right\rangle_{Z(\beta)} \right]$$
  This quantity naturally isolates the singular behavior and exhibits distinct scaling behaviors for different transition types.

• Numerical Implementation:

  • Algorithm: The authors utilize Quantum Monte Carlo (QMC) simulations, specifically the Stochastic Series Expansion (SSE) method combined with a bipartite reweight-annealing algorithm.
  • Efficiency: The RTE is computed via reweighting to handle the ratio of partition functions, while the DRTE is computed directly from expectation values of the Hamiltonian derivative in two separate ensembles ($Z(\beta)$ and $Z(2\beta)$).
  • Models Studied:
    1. Second-order: Bilayer Ising-Heisenberg model (Ising universality), Dimerized Anisotropic Heisenberg model (O(2) and O(3) universality).
    2. First-order: 2D Anisotropic Heisenberg model (symmetry-enhanced first-order), Checkerboard J-Q model.
    3. Weak First-order: $J-Q_2$ and $J-Q_3$ models (candidates for Deconfined Quantum Criticality).

3. Key Contributions: Unified Scaling Theory

The paper derives complete finite-size scaling forms for the DRTE across all three transition types:

A. Second-Order (Continuous) Transitions

• Scaling Form: The DRTE scales as $L^{1/\nu} \tilde{S}'(\Delta g L^{1/\nu})$.
• Signature:
  • The RTE curves for different system sizes intersect at a single critical point ($g_c$).
  • The DRTE exhibits a single peak that diverges with system size.
  • Result: Accurate extraction of critical exponents ($\nu$) and $g_c$ via data collapse.

B. Strong First-Order Transitions

• Scaling Form: The DRTE scales as $L^{d+z} f(\Delta g L^{d+z})$, where $d$ is spatial dimension and $z$ is the dynamical exponent.
• Signature:
  • Double-Peak Structure: A robust double-peak structure with one positive and one negative peak located on opposite sides of the transition.
  • Zero Crossing: The DRTE crosses exactly zero at the transition point ($g_c$).
  • Discontinuity: In the thermodynamic limit, the DRTE shows a discontinuity (jump) rather than a divergence, reflecting the latent heat and phase coexistence.

C. Weak First-Order Transitions

• Scaling Form: A hybrid scaling ansatz is proposed:
  $$\frac{\partial S^{(2)}}{\partial g} = \left[ 2L^{d+z}\Delta f'(g) + \tilde{S}'(\Delta g L^{1/\nu})L^{1/\nu} \right] [w(2\beta) - w(\beta)]$$
  Here, the first term represents the underlying first-order mechanism (suppressed by a small coefficient $\Delta f'$), and the second term represents the "singular" correction from the nearby critical fixed point.
• Signature:
  • Hybrid Behavior: For accessible system sizes ($L \ll \xi$), the data collapses using continuous exponents ($\nu$) due to the dominance of the second term.
  • Smoking-Gun: Despite the continuous-like scaling, the double-peak structure and the zero crossing at $g_c$ persist. This combination (continuous scaling + first-order topology) uniquely identifies a weak first-order transition.

4. Key Results

The authors validated their theory using high-precision QMC simulations:

1. Second-Order Verification:

   • Successfully demonstrated data collapse for (2+1)D Ising, O(2), and O(3) critical points.
   • Extracted critical exponents ($\nu$) and critical points ($g_c$) that match known theoretical values with high precision.

2. First-Order Verification:

   • Observed the characteristic double-peak structure and zero crossing in the 2D Anisotropic Heisenberg model and Checkerboard J-Q model.
   • Confirmed the $L^{d+z}$ volume-law scaling of the peak amplitude.

3. Weak First-Order Resolution (The Breakthrough):

   • Applied the theory to the debated $J-Q_2$ and $J-Q_3$ models, which are candidates for Deconfined Quantum Criticality (DQC).
   • Finding: While the DRTE data collapsed using continuous exponents (mimicking DQC), the double-peak structure and zero crossing were unambiguously observed.
   • Conclusion: The transitions in these models are weak first-order, not continuous. This resolves a long-standing ambiguity in the field, identifying the weak first-order nature at relatively modest system sizes ($L \le 56$), where previous methods failed.

5. Significance

• Resolution of Long-Standing Debates: The paper provides a definitive tool to distinguish between continuous and weak first-order transitions, settling controversies regarding the nature of Deconfined Quantum Criticality in $J-Q$ models.
• General and Unbiased Tool: The RTE/DRTE approach does not rely on specific order parameters (which may be unknown or difficult to define for topological transitions). It is applicable to a wide range of quantum systems.
• Computational Efficiency: The method allows for the identification of transition types in finite sizes that are much smaller than the correlation length, significantly reducing the computational resources required compared to traditional methods that need to reach the thermodynamic limit.
• Theoretical Framework: It establishes a complete finite-size scaling theory that unifies the description of second-order, strong first-order, and weak first-order transitions under a single observable framework.

In summary, this work introduces the Derivative of Rényi Thermal Entropy (DRTE) as a powerful, universal diagnostic tool that isolates singular free energy contributions, enabling the unambiguous detection of weak first-order transitions through a unique combination of continuous scaling behavior and first-order topological signatures (double-peak and zero-crossing).