A universal approach to Renyi entropy of multiple disjoint intervals

The paper proposes a general theory for computing Renyi entropy of multiple disjoint intervals by leveraging the connection between the replica trick and swapping operations, which is validated through applications to the transverse-field Ising model that yield results consistent with conformal field theory at criticality and extendable beyond it.

The Gist

The Big Picture: Measuring the "Weirdness" of Quantum Connections

Imagine you have a giant, complex puzzle made of quantum particles. In the quantum world, pieces of this puzzle can be "entangled," meaning they are mysteriously connected even when they are far apart. If you change one piece, the other changes instantly, no matter the distance.

Scientists want to measure how strong this connection is. They use a tool called Entropy.

• Von Neumann Entropy is the standard ruler for measuring this connection.
• Rényi Entropy is a more flexible, multi-purpose ruler. It gives you a richer picture of the connection, like looking at a 3D hologram instead of a flat photo.

The Problem:
Usually, calculating this entropy is easy if the connected pieces are right next to each other. But what if the pieces are disjoint? Imagine you have a long string of beads, and you want to measure the connection between Bead #1 and Bead #10, while ignoring everything in between.

• If you have just two separate groups, it's hard but doable.
• If you have three, four, or ten separate groups scattered across the string, the math becomes a nightmare. It's like trying to solve a Rubik's cube that keeps changing its rules every time you touch it.

The Solution:
The authors of this paper (Han-Qing Shi and Hai-Qing Zhang) found a clever shortcut. They realized that instead of doing impossible math, you can use a "magic trick" called Swapping.

---

The Core Idea: The "Copy-Paste-Swap" Trick

To understand their method, let's use an analogy involving Photocopying and Swapping Cards.

1. The Setup: The Original and the Copy

Imagine you have a deck of cards representing your quantum system. You want to measure the connection between specific cards (let's call them the "Blue Cards" and the "Red Cards") that are separated by other cards.

In the old way (the "Replica Trick" used in physics), you would have to build a massive, multi-layered mathematical structure to calculate this. It's like trying to build a house of cards while standing on a trampoline.

The authors say: "Let's try something simpler."

• Take your original deck of cards.
• Make an exact photocopy of it. Now you have two identical decks: Deck A and Deck B.

2. The Swap: The "Magic Hand"

Now, imagine a magical hand that can reach into both decks at the same time.

• It grabs the Blue Cards from Deck A.
• It grabs the Blue Cards from Deck B.
• It swaps them. Now Deck A has the Blue Cards from Deck B, and Deck B has the Blue Cards from Deck A.
• The Red Cards (the parts you don't care about) stay exactly where they are.

3. The Measurement: The "Echo"

After the swap, you ask the decks: "How much did you change?"

• If the Blue Cards were totally independent of the Red Cards, the swap wouldn't matter much. The decks would look almost the same.
• If the Blue Cards were deeply entangled with the Red Cards, the swap would cause a huge "glitch" or "echo" in the system.

The authors discovered that the Rényi Entropy (the measure of connection) is directly related to how much the system "echoes" after this swap.

• High Echo = High Entropy (Strong connection).
• Low Echo = Low Entropy (Weak connection).

Why This is a Game-Changer

1. It Works for Any Number of Groups

Previously, if you wanted to measure the connection between two separate groups, you needed complex math. If you wanted to measure three or four groups, the math was so hard that no one had a formula for it.

With the Swapping Trick, it doesn't matter if you have two groups or twenty groups. You just make more copies of the deck (3 copies for 3rd-order entropy, 4 copies for 4th-order, etc.) and perform a specific "circular swap" (Deck 1 gives to Deck 2, Deck 2 gives to Deck 3, Deck 3 gives back to Deck 1).

It's like a game of musical chairs. No matter how many people are playing, the rule is the same: everyone passes their chair to the person on their right.

2. It Works Everywhere (Critical and Non-Critical)

The paper tested this on a famous model called the Ising Model (think of it as a row of tiny magnets that can point up or down).

• At the "Critical Point": This is the exact moment the magnets are deciding whether to all point up or all point down. The physics here is perfectly smooth and predictable (like a calm lake). The authors' method matched the perfect mathematical predictions of the experts.
• Away from the Critical Point: This is when the magnets are messy and chaotic (like a stormy sea). The old math formulas break down here. But the Swapping Trick still works perfectly! It can measure the entropy even when the system is messy and there is no simple formula to describe it.

The Analogy of the "Universal Translator"

Think of the old methods as trying to translate a book written in a language that only exists in a few ancient libraries. You need a specific dictionary for every single sentence.

The authors' method is like a Universal Translator.

• You don't need to know the specific rules of the language (the specific physics of the system).
• You just need to know the "Swap" rule.
• You apply the rule, and the machine tells you the answer, whether the text is from a calm, poetic section (Critical Point) or a chaotic, noisy section (Non-Critical).

Summary of Results

The authors tested this on a computer simulation of a 1D chain of magnets:

1. Two Disjoint Intervals: They checked their math against known theories. Result: Perfect match.
2. Three and Four Disjoint Intervals: They tried to measure connections between 3 and 4 separate groups of magnets. Result: They got clear, accurate numbers where no one had formulas before.
3. Different Magnetic Fields: They changed the strength of the magnetic field. The method worked smoothly as the system went from "calm" to "chaotic."

The Takeaway

This paper gives physicists a universal toolkit. Instead of struggling to invent new, complex math for every new type of quantum system or every new number of disconnected parts, they can now just use the Swapping Operation.

It's like realizing that to measure the distance between any two cities, you don't need a different map for every pair of cities. You just need a ruler and a straight line. The authors found the "ruler" for quantum entanglement that works for any number of disconnected islands in the quantum ocean.

Technical Summary

1. Problem Statement

The calculation of Rényi entropy ($S_m$) and entanglement entropy (von Neumann entropy) for quantum many-body systems is a fundamental challenge in quantum physics. While these quantities are well-understood for single contiguous subsystems or specific cases in Conformal Field Theory (CFT), significant difficulties arise when dealing with multiple disjoint intervals.

• Limitations of Current Methods:
  • CFT: For two disjoint intervals, CFT requires complex four-point correlation functions and model-dependent universal functions. For three or more disjoint intervals, the required higher-order correlation functions make analytical calculations intractable.
  • General Systems: Away from critical points (non-conformal regimes), there are no closed analytical formulas for the entropy of multiple disjoint subsystems.
  • Numerical Complexity: Direct computation via the density matrix is often impossible due to the exponential growth of the Hilbert space dimension.

The authors aim to develop a general, systematic, and analytical method to compute the Rényi entropy $S_m$ for an arbitrary number of disjoint intervals in general quantum systems, not limited to critical points.

2. Methodology: The Swapping Operator Approach

The core innovation of the paper is the proposal of a swapping operator method, which draws a theoretical parallel between the replica trick in Quantum Field Theory (QFT) and discrete state swapping in quantum mechanics.

A. Theoretical Framework

• Replica Trick Analogy: In QFT, the replica trick computes $\text{Tr}(\rho_A^m)$ by cutting the path integral into $m$ sheets and identifying boundaries cyclically ($\phi_k(A^+) = \phi_{k+1}(A^-)$).
• Swapping Operator Definition: The authors propose that this cyclic identification can be mimicked by a unitary swapping operator ($S_{\text{wap}}^{(m)}$) acting on $m$ copies of the quantum state $|\Psi\rangle$.
  • Consider $m$ copies of the system: $\bigotimes_{j=1}^m |\Psi_j\rangle$.
  • The subsystem $A$ consists of $n$ disjoint intervals $\{a_1, a_2, \dots, a_n\}$.
  • The operator $S_{\text{wap}}^{(m)}$ acts on the $m$ copies by cyclically exchanging the states of subsystem $A$ between adjacent copies:
    • States in $A$ of copy $j$ are swapped with states in $A$ of copy $j+1$.
    • The states in $A$ of the final copy ($m$) are swapped with the first copy ($1$).
  • States in the complement subsystem $B$ remain untouched.

B. Mathematical Derivation

The authors prove that the expectation value of this swapping operator is directly related to the trace of the $m$-th power of the reduced density matrix:
$$\langle S_{\text{wap}}^{(m)} \rangle = \text{Tr}(\rho_A^m)$$
Consequently, the $m$-th order Rényi entropy is given by:
$$S_m = \frac{1}{1-m} \ln \langle S_{\text{wap}}^{(m)} \rangle$$

• Proof Strategy: The proof involves expanding the state $|\Psi\rangle$ in a basis, applying the swapping operator to the composite state $\bigotimes |\Psi\rangle$, and showing that the resulting expectation value matches the definition of $\text{Tr}(\rho_A^m)$ derived from the reduced density matrix. This is demonstrated rigorously for:
  1. Two disjoint intervals ($n=2$).
  2. $n$ multiple disjoint intervals.
  3. General $m$-th order entropy.

3. Key Contributions

1. Universal Formula: Established a general formula linking Rényi entropy to the expectation value of a swapping operator, applicable to any number of disjoint intervals and any order $m$.
2. Bridging QFT and Lattice Models: Provided a conceptual bridge between the continuum replica trick in QFT and discrete lattice calculations, allowing the use of swapping operators in numerical simulations.
3. Beyond Criticality: Unlike CFT methods which are restricted to critical points, this method is valid for non-critical regimes (gapped phases) where no analytical closed-form solutions exist.
4. Computational Feasibility: Demonstrated that for finite spin systems, the expectation value $\langle S_{\text{wap}} \rangle$ can be computed exactly using standard linear algebra techniques (e.g., Singular Value Decomposition) without needing to construct the full reduced density matrix explicitly.

4. Results and Applications

The authors applied their theory to the one-dimensional Transverse-Field Quantum Ising Model (TFQIM) with periodic boundary conditions.

• System Setup:

  • Hamiltonian: $H = -J \sum (\sigma_i^z \sigma_{i+1}^z + h \sigma_i^x)$.
  • System sizes: Up to $N=24$ sites (Hilbert space dimension $\sim 10^8$).
  • States: Ground states ($|\Psi\rangle$) calculated via exact diagonalization.

• Case Studies:

  1. Two Disjoint Intervals:
     • Calculated $S_2$ varying the separation distance and interval sizes.
     • Validation: At the critical point ($h=1$), the results matched perfectly with analytical CFT predictions (using four-point correlation functions).
     • Non-Critical: For $h \neq 1$ (paramagnetic/ferromagnetic phases), the method provided results where CFT formulas are unavailable.
  2. Three and Four Disjoint Intervals:
     • Calculated $S_2$ for systems with 3 and 4 disjoint intervals.
     • Significance: Since no analytical CFT formulas exist for $n \geq 3$ intervals, this serves as a primary demonstration of the method's utility in previously intractable scenarios.
     • Observed that entropy peaks at the critical point ($h=1$) and decreases as the system moves away from criticality.

• Symmetry Observations: The results exhibited expected parity symmetries (e.g., $x \leftrightarrow L-x$) due to periodic boundary conditions, confirming the consistency of the numerical implementation.

5. Significance and Future Outlook

• Theoretical Impact: The work provides a robust, "universal" tool for studying multipartite entanglement, a crucial resource in quantum information and many-body physics.
• Practical Application: The method is not limited to ground states; it can be applied to dynamical cases (time-evolved states) and higher-dimensional systems.
• Experimental Relevance: Since the expectation value of a swapping operator can be measured in quantum Monte Carlo simulations or potentially in quantum simulators (as hinted by reference [15]), this approach offers a pathway to experimentally probe complex entanglement structures in systems with multiple disjoint regions.
• Future Work: The authors suggest extending this to higher-order Rényi entropies ($m > 2$) and applying it to higher-dimensional lattice models.

In summary, Shi and Zhang have successfully generalized the computation of Rényi entropy for multiple disjoint intervals by introducing a swapping operator formalism that unifies analytical insights from CFT with numerical capabilities for general quantum systems.